Abstract:Background Recent circadian clock studies using gene expression microarray in two different tissues of mouse have revealed not all circadian-related genes are synchronized in phase or peak expression times across tissues in vivo . Instead, some circadian-related genes may be delayed by 4–8 hrs in peak expression in one tissue relative to the other. These interesting biological observations prompt a statistical question regarding how to distinguish the synchronized genes from genes that are systematically lagged in phase/peak expression time across two tissues. Results We propose a set of techniques from circular statistics to analyze phase angles of circadian-related genes in two tissues. We first estimate the phases of a cycling gene separately in each tissue, which are then used to estimate the paired angular difference of the phase angles of the gene in the two tissues. These differences are modeled as a mixture of two von Mises distributions which enables us to cluster genes into two groups; one group having synchronized transcripts with the same phase in the two tissues, the other containing transcripts with a discrepancy in phase between the two tissues. For each cluster of genes we assess the association of phases across the tissue types using circular-circular regression. We also develop a bootstrap methodology based on a circular-circular regression model to evaluate the improvement in fit provided by allowing two components versus a one-component von-Mises model. Conclusion We applied our proposed methodologies to the circadian-related genes common to heart and liver tissues in Storch et al . [ 2 ], and found that an estimated 80% of circadian-related transcripts common to heart and liver tissues were synchronized in phase, and the other 20% of transcripts were lagged about 8 hours in liver relative to heart. The bootstrap p -value for being one cluster is 0.063, which suggests the possibility of two clusters. Our methodologies can be extended to analyze peak expression times of circadian-related genes across more than two tissues, for example, kidney, heart, liver, and the suprachiasmatic nuclei (SCN) of the hypothalamus.

CircaKB: a comprehensive knowledgebase of circadian genes across multiple species

CirGRDB: a Database for the Genome-Wide Deciphering Circadian Genes and Regulators.

CircFunBase: a database for functional circular RNAs.

CGDB: a database of circadian genes in eukaryotes.

A circadian gene expression atlas in mammals: Implications for biology and medicine

Identification of Human Circadian Genes Based on Time Course Gene Expression Profiles by Using a Deep Learning Method.

Circ2Disease: a Manually Curated Database of Experimentally Validated Circrnas in Human Disease

Genome-wide circadian rhythm detection methods: systematic evaluations and practical guidelines

Genetic and epigeneticregulations of mammalian circadian rhythms.

What time is it? Deep learning approaches for circadian rhythms

A Computational Model for Functional Mapping of Genes That Regulate Intra-Cellular Circadian Rhythms

Inferring circadian gene regulatory relationships from gene expression data with a hybrid framework

CIRCpedia V2: an Updated Database for Comprehensive Circular RNA Annotation and Expression Comparison.

Dual-Approach Co-expression Analysis Framework (D-CAF) Enables Identification of Novel Circadian Regulation From Multi-Omic Timeseries Data

Assembly of a Comprehensive Regulatory Network for the Mammalian Circadian Clock: A Bioinformatics Approach

High-throughput Discovery of Genetic Determinants of Circadian Misalignment.

Phase analysis of circadian-related genes in two tissues

Clocks not winding down: unravelling circadian networks

Guidelines for Genome-Scale Analysis of Biological Rhythms

A novel bioinformatics approach to reveal the role of circadian oscillations in AD development

Statistical mechanics of clock gene networks underlying circadian rhythms