Browse
Submit Data
Databases
API
Help

Dataset Information

14 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

ICON: An Adaptation of Infinite HMMs for Time Traces with Drift.

ABSTRACT: Bayesian nonparametric methods have recently transformed emerging areas within data science. One such promising method, the infinite hidden Markov model (iHMM), generalizes the HMM that itself has become a workhorse in single molecule data analysis. The iHMM goes beyond the HMM by self-consistently learning all parameters learned by the HMM in addition to learning the number of states without recourse to any model selection steps. Despite its generality, simple features (such as drift), common to single molecule time traces, result in an overinterpretation of drift and the introduction of artifact states. Here we present an adaptation of the iHMM that can treat data with drift originating from one or many traces (e.g., Förster resonance energy transfer). Our fully Bayesian method couples the iHMM to a continuous control process (drift) self-consistently learned while learning all other quantities determined by the iHMM (including state numbers). A key advantage of this method is that all traces-regardless of drift or states visited across traces-may now be treated on an equal footing, thereby eliminating user-dependent trace selection (based on drift levels), preprocessing to remove drift, and postprocessing model selection based on state number.

SUBMITTER: Sgouralis I

PROVIDER: S-EPMC5443975 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

An Introduction to Infinite HMMs for Single-Molecule Data Analysis.

Project description:The hidden Markov model (HMM) has been a workhorse of single-molecule data analysis and is now commonly used as a stand-alone tool in time series analysis or in conjunction with other analysis methods such as tracking. Here, we provide a conceptual introduction to an important generalization of the HMM, which is poised to have a deep impact across the field of biophysics: the infinite HMM (iHMM). As a modeling tool, iHMMs can analyze sequential data without a priori setting a specific number of states as required for the traditional (finite) HMM. Although the current literature on the iHMM is primarily intended for audiences in statistics, the idea is powerful and the iHMM's breadth in applicability outside machine learning and data science warrants a careful exposition. Here, we explain the key ideas underlying the iHMM, with a special emphasis on implementation, and provide a description of a code we are making freely available. In a companion article, we provide an important extension of the iHMM to accommodate complications such as drift.

| S-EPMC5448313 | biostudies-literature

Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift.

Project description:Resurrection studies are a useful tool to measure how phenotypic traits have changed in populations through time. If these trait modifications correlate with the environmental changes that occurred during the time period, it suggests that the phenotypic changes could be a response to selection. Selfing, through its reduction of effective size, could challenge the ability of a population to adapt to environmental changes. Here, we used a resurrection study to test for adaptation in a selfing population of Medicago truncatula, by comparing the genetic composition and flowering times across 22 generations. We found evidence for evolution toward earlier flowering times by about two days and a peculiar genetic structure, typical of highly selfing populations, where some multilocus genotypes (MLGs) are persistent through time. We used the change in frequency of the MLGs through time as a multilocus fitness measure and built a selection gradient that suggests evolution toward earlier flowering times. Yet, a simulation model revealed that the observed change in flowering time could be explained by drift alone, provided the effective size of the population is small enough (<150). These analyses suffer from the difficulty to estimate the effective size in a highly selfing population, where effective recombination is severely reduced.

| S-EPMC8794724 | biostudies-literature

Continuous-state HMMs for modeling time-series single-cell RNA-Seq data.

Project description:MOTIVATION:Methods for reconstructing developmental trajectories from time-series single-cell RNA-Seq (scRNA-Seq) data can be largely divided into two categories. The first, often referred to as pseudotime ordering methods are deterministic and rely on dimensionality reduction followed by an ordering step. The second learns a probabilistic branching model to represent the developmental process. While both types have been successful, each suffers from shortcomings that can impact their accuracy. RESULTS:We developed a new method based on continuous-state HMMs (CSHMMs) for representing and modeling time-series scRNA-Seq data. We define the CSHMM model and provide efficient learning and inference algorithms which allow the method to determine both the structure of the branching process and the assignment of cells to these branches. Analyzing several developmental single-cell datasets, we show that the CSHMM method accurately infers branching topology and correctly and continuously assign cells to paths, improving upon prior methods proposed for this task. Analysis of genes based on the continuous cell assignment identifies known and novel markers for different cell types. AVAILABILITY AND IMPLEMENTATION:Software and Supporting website: www.andrew.cmu.edu/user/chiehl1/CSHMM/. SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.

| S-EPMC6853676 | biostudies-literature

Generalizing HMMs to Continuous Time for Fast Kinetics: Hidden Markov Jump Processes.

Project description:The hidden Markov model (HMM) is a framework for time series analysis widely applied to single-molecule experiments. Although initially developed for applications outside the natural sciences, the HMM has traditionally been used to interpret signals generated by physical systems, such as single molecules, evolving in a discrete state space observed at discrete time levels dictated by the data acquisition rate. Within the HMM framework, transitions between states are modeled as occurring at the end of each data acquisition period and are described using transition probabilities. Yet, whereas measurements are often performed at discrete time levels in the natural sciences, physical systems evolve in continuous time according to transition rates. It then follows that the modeling assumptions underlying the HMM are justified if the transition rates of a physical process from state to state are small as compared to the data acquisition rate. In other words, HMMs apply to slow kinetics. The problem is, because the transition rates are unknown in principle, it is unclear, a priori, whether the HMM applies to a particular system. For this reason, we must generalize HMMs for physical systems, such as single molecules, because these switch between discrete states in "continuous time". We do so by exploiting recent mathematical tools developed in the context of inferring Markov jump processes and propose the hidden Markov jump process. We explicitly show in what limit the hidden Markov jump process reduces to the HMM. Resolving the discrete time discrepancy of the HMM has clear implications: we no longer need to assume that processes, such as molecular events, must occur on timescales slower than data acquisition and can learn transition rates even if these are on the same timescale or otherwise exceed data acquisition rates.

| S-EPMC7896036 | biostudies-literature

Time makes histone H3 modifications drift

Project description:To disclose the epigenetic drift of time passing, we determined the genome-wide distributions of mono- and tri-methylated lysine 4 and acetylated and tri-methylated lysine 27 of histone H3 in the livers of healthy 3, 6 and 12 months old C57BL/6 mice. The comparison of different age profiles of histone H3 marks revealed global redistribution of histone H3 modifications with time, in particular in intergenic regions and near transcription start sites, as well as altered correlation between the profiles of different histone modifications. Moreover, feeding mice with caloric restriction diet, a treatment known to retard aging, preserved younger state of histone H3 in these genomic regions.

2022-07-13 | GSE179090 | GEO

Evaluation of drift time specific collision energies

Project description:Comparison among the following data independent acquisition modes provided on the Waters Synapt G2-S instrument: MSE, MSE with Ion Mobility (HDMSE), and MSE with Ion Mobility and drift time specific collision energies (UDMSE). The following on-column loads and LC gradient lengths have been used for the DDA data (1-20130710-63647): * 200 ng, 90 minutes gradient * 300 ng, 180 minutes gradient * 1000 ng, 180 minutes gradient. DDA data has been processed by using MaxQuant (v.1.3.0.5) software, searching against a organism specific (Homo Sapiens) Uniprot/Swissprot database. FDR has been controlled by using a pseudo-reverse (KR positions kept) database at both protein and peptide levels to 1%. Proteins with less than two peptides (minimum of 6 amino acids length, up to 2 missed cleavages and 3 variable modifications allowed. Variable identifications were : methionine oxidation, fixed modifications were : cysteine carbamidomethylation ) identified have been discarded. A match between runs of a 2 minutes windows among the three technical replicates has been performed.

2013-12-05 | PXD000303 | Pride

Decoding Single Molecule Time Traces with Dynamic Disorder.

Project description:Single molecule time trajectories of biomolecules provide glimpses into complex folding landscapes that are difficult to visualize using conventional ensemble measurements. Recent experiments and theoretical analyses have highlighted dynamic disorder in certain classes of biomolecules, whose dynamic pattern of conformational transitions is affected by slower transition dynamics of internal state hidden in a low dimensional projection. A systematic means to analyze such data is, however, currently not well developed. Here we report a new algorithm-Variational Bayes-double chain Markov model (VB-DCMM)-to analyze single molecule time trajectories that display dynamic disorder. The proposed analysis employing VB-DCMM allows us to detect the presence of dynamic disorder, if any, in each trajectory, identify the number of internal states, and estimate transition rates between the internal states as well as the rates of conformational transition within each internal state. Applying VB-DCMM algorithm to single molecule FRET data of H-DNA in 100 mM-Na+ solution, followed by data clustering, we show that at least 6 kinetic paths linking 4 distinct internal states are required to correctly interpret the duplex-triplex transitions of H-DNA.

| S-EPMC5226833 | biostudies-literature

Important role of genetic drift in rapid polygenic adaptation.

Project description:We analyzed a model to determine the factors that facilitate or limit rapid polygenic adaptation. This model includes population genetic terms of mutation and both directional and stabilizing selection on a highly polygenic trait in a diploid population of finite size. First, we derived the equilibrium distribution of the allele frequencies of the multilocus model by diffusion approximation. This formula describing the equilibrium allele frequencies as a mutation-selection-drift balance was examined by computer simulation using parameter values inferred for human height, a well-studied polygenic trait. Second, assuming that a sudden environmental shift of the fitness optimum occurs while the population is in equilibrium, we analyzed the adaptation of the trait to the new optimum. The speed at which the trait mean approaches the new optimum increases with the equilibrium genetic variance. Thus, large population size and/or large mutation rate may facilitate rapid adaptation. Third, the contribution of an individual locus i to polygenic adaptation depends on the compound parameter ?ipi0qi0 , where ?i is the effect size, pi0 the equilibrium frequency of the trait-increasing allele of this locus, and qi0=1-pi0 . Thus, only loci with large values of this parameter contribute coherently to polygenic adaptation. Given that mutation rates are relatively small, this is more likely in large populations, in which the effects of drift are limited.

| S-EPMC7029068 | biostudies-literature

Clustering gene expression time series data using an infinite Gaussian process mixture model

Project description:In order to identify and characterize novel human gene expression responses to glucocorticoids, we exposed the human lung adenocarcinoma cell line, A549, to the synthetic glucocorticoid dexamethasone for 1, 3, 5, 7, 9, and 11 hrs in duration as well as to a paired vehicle control, ethanol. We assayed gene expression with RNA-seq and clustered gene expression profiles using an infinite Gaussian process mixture model.

2017-10-08 | GSE104714 | GEO

Mountain chickadees from different elevations sing different songs: acoustic adaptation, temporal drift or signal of local adaptation?

Project description:Song in songbirds is widely thought to function in mate choice and male-male competition. Song is also phenotypically plastic and typically learned from local adults; therefore, it varies across geographical space and can serve as a cue for an individual's location of origin, with females commonly preferring males from their respective location. Geographical variation in song dialect may reflect acoustic adaptation to different environments and/or serve as a signal of local adaptation. In montane environments, environmental differences can occur over an elevation gradient, favouring local adaptations across small spatial scales. We tested whether food caching mountain chickadees, known to exhibit elevation-related differences in food caching intensity, spatial memory and the hippocampus, also sing different dialects despite continuous distribution and close proximity. Male songs were collected from high and low elevations at two different mountains (separated by 35 km) to test whether song differs between elevations and/or between adjacent populations at each mountain. Song structure varied significantly between high and low elevation adjacent populations from the same mountain and between populations from different mountains at the same elevations, despite a continuous distribution across each mountain slope. These results suggest that elevation-related differences in song structure in chickadees might serve as a signal for local adaptation.

| S-EPMC4448878 | biostudies-other

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data