Order restricted inference in chronobiology

CIRCUST: A novel methodology for temporal order reconstruction of molecular rhythms; validation and application towards a daily rhythm gene expression atlas in humans

The circadian system drives near-24-h oscillations in behaviors and biological processes. The underlying core molecular clock regulates the expression of other genes, and it has been shown that the expression of more than 50 percent of genes in mammals displays 24-h rhythmic patterns, with the specific genes that cycle varying from one tissue to another. Determining rhythmic gene expression patterns in human tissues sampled as single timepoints has several challenges, including the reconstruction of temporal order of highly noisy data. Previous methodologies have attempted to address these challenges in one or a small number of tissues for which rhythmic gene evolutionary conservation is assumed to be preserved. Here we introduce CIRCUST, a novel CIRCular-robUST methodology for analyzing molecular rhythms, that relies on circular statistics, is robust against noise, and requires fewer assumptions than existing methodologies. Next, we validated the method against four controlled experiments in which sampling times were known, and finally, CIRCUST was applied to 34 tissues from the Genotype-Tissue Expression (GTEx) dataset with the aim towards building a comprehensive daily rhythm gene expression atlas in humans. The validation and application shown here indicate that CIRCUST provides a flexible framework to formulate and solve the issues related to the analysis of molecular rhythms in human tissues. CIRCUST methodology is publicly available at https://github.com/yolandalago/CIRCUST/.

Frequency Modulated Möbius Model Accurately Predicts Rhythmic Signals in Biological and Physical Sciences

Motivated by applications in physical and biological sciences, we developed a Frequency Modulated Möbius (FMM) model to describe rhythmic patterns in oscillatory systems. Unlike standard symmetric sinusoidal models, FMM is a flexible parametric model that allows deformations to sinusoidal shape to accommodate commonly seen asymmetries in applications. FMM model parameters are easy to estimate and the model is easy to interpret complex rhythmic data. We illustrate FMM model in three disparate applications, namely, circadian clock gene expression, corticoptropin levels in depressed patients and the temporal light intensity patterns of distant stars. In each case, FMM model is demonstrated to be flexible, scientifically plausible and easy to interpret. Analysis of synthetic data derived from patterns of real data, suggest that FMM model fits the data very well both visually as well as in terms of the goodness of fit measure total mean squared error. An R language based software for implementing FMM model is available.