Differential Phasing between Circadian Clocks in the Brain and Peripheral Organs in Humans

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/.

Experimental design and power calculation in omics circadian rhythmicity detection using the cosinor model

Circadian clocks are 24-h endogenous oscillators in physiological and behavioral processes. Though recent transcriptomic studies have been successful in revealing the circadian rhythmicity in gene expression, the power calculation for omics circadian analysis have not been fully explored. In this paper, we develop a statistical method, namely CircaPower, to perform power calculation for circadian pattern detection. Our theoretical framework is determined by three key factors in circadian gene detection: sample size, intrinsic effect size and sampling design. Via simulations, we systematically investigate the impact of these key factors on circadian power calculation. We not only demonstrate that CircaPower is fast and accurate, but also show its underlying cosinor model is robust against variety of violations of model assumptions. In real applications, we demonstrate the performance of CircaPower using mouse pan-tissue data and human post-mortem brain data, and illustrate how to perform circadian power calculation using mouse skeleton muscle RNA-Seq pilot as case study. Our method CircaPower has been implemented in an R package, which is made publicly available on GitHub ( https://github.com/circaPower/circaPower).

Disruption of day-to-night changes in circadian gene expression with chronic tendinopathy

Overuse injury in tendon tissue (tendinopathy) is a frequent and costly musculoskeletal disorder and represents a major clinical problem with unsolved pathogenesis. Studies in mice have demonstrated that circadian clock-controlled genes are vital for protein homeostasis and important in the development of tendinopathy. We performed RNA sequencing, collagen content and ultrastructural analyses on human tendon biopsies obtained 12 h apart in healthy individuals to establish whether human tendon is a peripheral clock tissue and we performed RNA sequencing on patients with chronic tendinopathy to examine the expression of circadian clock genes in tendinopathic tissues. We found time-dependent expression of 280 RNAs including 11 conserved circadian clock genes in healthy tendons and markedly fewer (23) differential RNAs with chronic tendinopathy. Further, the expression of COL1A1 and COL1A2 was reduced at night but was not circadian rhythmic in synchronised human tenocyte cultures. In conclusion, day-to-night changes in gene expression in healthy human patellar tendons indicate a conserved circadian clock as well as the existence of a night reduction in collagen I expression. KEY POINTS: Tendinopathy is a major clinical problem with unsolved pathogenesis. Previous work in mice has shown that a robust circadian rhythm is required for collagen homeostasis in tendons. The use of circadian medicine in the diagnosis and treatment of tendinopathy has been stifled by the lack of studies on human tissue. Here, we establish that the expression of circadian clock genes in human tendons is time dependent, and now we have data to corroborate that circadian output is reduced in diseased tendon tissues. We consider our findings to be of significance in advancing the use of the tendon circadian clock as a therapeutic target or preclinical biomarker for tendinopathy.

TimeTeller for timing health: The potential of circadian medicine to improve performance, prevent disease and optimize treatment

Circadian medicine, the study of the effects of time on health and disease has seen an uprising in recent years as a means to enhance health and performance, and optimize treatment timing. Our endogenous time generating system -the circadian clock- regulates behavioural, physiological and cellular processes. Disruptions of the clock, via external factors like shift work or jet lag, or internal perturbations such as genetic alterations, are linked to an increased risk of various diseases like obesity, diabetes, cardiovascular diseases and cancer. By aligning an individual's circadian clock with optimal times for performing daily routines, physical and mental performance, and also the effectiveness of certain therapies can be improved. Despite the benefits of circadian medicine, the lack of non-invasive tools for characterizing the clock limits the potential of the field. TimeTeller is a non-invasive molecular/digital tool for the characterization of circadian rhythms and prediction of daily routines, including treatment timing, to unlock the potential of circadian medicine and implementing it in various settings. Given the multiple known and potentially yet unknown dependent health factors of individual circadian rhythms, the utility of this emerging biomarker is best exploited in data driven, personalized medicine use cases, using health information across lifestyle, care, and research settings.

Diurnal and circadian rhythmicity of the human blood transcriptome overlaps with organ- and tissue-specific expression of a non-human primate

Background

Twenty-four-hour rhythmicity in mammalian tissues and organs is driven by local circadian oscillators, systemic factors, the central circadian pacemaker and light-dark cycles. At the physiological level, the neural and endocrine systems synchronise gene expression in peripheral tissues and organs to the 24-h-day cycle, and disruption of such regulation has been shown to lead to pathological conditions. Thus, monitoring rhythmicity in tissues/organs holds promise for circadian medicine; however, most tissues and organs are not easily accessible in humans and alternative approaches to quantify circadian rhythmicity are needed. We investigated the overlap between rhythmic transcripts in human blood and transcripts shown to be rhythmic in 64 tissues/organs of the baboon, how these rhythms are aligned with light-dark cycles and each other, and whether timing of tissue-specific rhythmicity can be predicted from a blood sample.

Results

We compared rhythmicity in transcriptomic time series collected from humans and baboons using set logic, circular cross-correlation, circular clustering, functional enrichment analyses, and least squares regression. Of the 759 orthologous genes that were rhythmic in human blood, 652 (86%) were also rhythmic in at least one baboon tissue and most of these genes were associated with basic processes such as transcription and protein homeostasis. In total, 109 (17%) of the 652 overlapping rhythmic genes were reported as rhythmic in only one baboon tissue or organ and several of these genes have tissue/organ-specific functions. The timing of human and baboon rhythmic transcripts displayed prominent ‘night’ and ‘day’ clusters, with genes in the dark cluster associated with translation. Alignment between baboon rhythmic transcriptomes and the overlapping human blood transcriptome was significantly closer when light onset, rather than midpoint of light, or end of light period, was used as phase reference point. The timing of overlapping human and baboon rhythmic transcriptomes was significantly correlated in 25 tissue/organs with an average earlier timing of 3.21 h (SD 2.47 h) in human blood.

Conclusions

The human blood transcriptome contains sets of rhythmic genes that overlap with rhythmic genes of tissues/organs in baboon. The rhythmic sets vary across tissues/organs, but the timing of most rhythmic genes is similar in human blood and baboon tissues/organs. These results have implications for development of blood transcriptome-based biomarkers for circadian rhythmicity in tissues and organs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01258-7.

Expression patterns of clock genes in the kidney of two Lasiopodomys species

Previous studies showed that the kidney has its own molecular circadian clock expression regulation that maintains the homeostasis of physiological processes. However, limited information is available on the molecular mechanisms of the kidney circadian rhythm in subterranean rodents. Here, we report circadian gene expression in the kidney of subterranean Mandarin voles and the related aboveground Brandt’s voles, reared under 12L:12D (LD) or dark (DD) conditions, respectively. The results showed that the rhythmic genes were represented in Brandt’s voles in higher numbers under LD than DD conditions, but the number of rhythmic genes in Mandarin voles was similar between the two treatment conditions. The gene expression levels at different timepoints all showed reduced results under DD conditions compared with those in the LD cycle in Brandt’s voles, whereas the expression levels of the tested genes at certain Zeitgeber timepoints showed higher results than in the LD cycle in Mandarin voles. The gene expression peak showed chaotic resetting under DD conditions in both voles. We thus suggest that Mandarin and Brandt’s voles have different molecular circadian clock expression adjustment patterns in the kidney as an adaptation to different living environments. Mandarin voles seem to be more adapted to the dark environment, while Brandt’s voles are more dependent on external light conditions.

Diurnal shift of mouse activity by the deficiency of an aging-related gene Lmna

Nuclear lamina is a fundamental structure of the cell nucleus and regulates a wide range of molecular pathways. Defects of components of the nuclear lamina cause aging-like physiological disorders, called laminopathy. Generally, aging and diseases are often associated with perturbation of various time-of-day-dependent regulations, but it remains still elusive whether laminopathy induces any changes of the circadian clock and physiological rhythms. Here we demonstrated that deficiency of Lmna gene in mice caused an obvious shift of locomotor activities to the daytime. The abnormal activity profile was accompanied by a remarkable change in phase-angle between the central clock in the suprachiasmatic nucleus (SCN) and lung peripheral clocks, leaving the phase of the SCN clock unaffected by the mutation. These observations suggest that Lmna deficiency causes a change of the habitat from nocturnal to diurnal behaviors. On the other hand, molecular oscillation and its phase resetting mechanism were intact in both the Lmna-deficient cells and progeria-mimicking cells. Intriguingly, high-fat diet feeding extended the short lifespan and ameliorated the abnormalities of the behaviors and the phase of the peripheral clock in the Lmna-deficient mice. The present study supports the important contribution of the energy conditions to a shift between the diurnal and nocturnal activities.

Likelihood-based tests for detecting circadian rhythmicity and differential circadian patterns in transcriptomic applications

Circadian rhythmicity in transcriptomic profiles has been shown in many physiological processes, and the disruption of circadian patterns has been found to associate with several diseases. In this paper, we developed a series of likelihood-based methods to detect (i) circadian rhythmicity (denoted as LR_rhythmicity) and (ii) differential circadian patterns comparing two experimental conditions (denoted as LR_diff). In terms of circadian rhythmicity detection, we demonstrated that our proposed LR_rhythmicity could better control the type I error rate compared to existing methods under a wide variety of simulation settings. In terms of differential circadian patterns, we developed methods in detecting differential amplitude, differential phase, differential basal level and differential fit, which also successfully controlled the type I error rate. In addition, we demonstrated that the proposed LR_diff could achieve higher statistical power in detecting differential fit, compared to existing methods. The superior performance of LR_rhythmicity and LR_diff was demonstrated in four real data applications, including a brain aging data (gene expression microarray data of human postmortem brain), a time-restricted feeding data (RNA sequencing data of human skeletal muscles) and a scRNAseq data (single cell RNA sequencing data of mouse suprachiasmatic nucleus). An R package for our methods is publicly available on GitHub https://github.com/diffCircadian/diffCircadian.

CircadiOmic medicine and aging

The earth displays daily, seasonal and annual environmental cycles that have led to evolutionarily adapted ultradian, circadian and infradian rhythmicities in the entire biosphere. All biological organisms must adapt to these cycles that synchronize the function of their circadiome. The objective of this review is to discuss the latest knowledge regarding the role of circadiomics in health and aging. The biological timekeepers are responsive to the environmental cues at microsecond to seasonal time-scales and act with precision of a clock machinery. The robustness of these rhythms is essential to normal daily function of cells, tissues and organs. Mis-alignment of circadian rhythms makes the individual prone to aging, sleep disorders, cancer, diabetes, and neuro-degenerative diseases. Circadian and CircadiOmic medicine are emerging fields that leverage our in-depth understanding of health issues, that arise as a result of disturbances in circadian rhythms, towards establishing better therapeutic approaches in personalized medicine and for geroprotection.

Genes Relevant to Tissue Response to Cancer Therapy Display Diurnal Variation in mRNA Expression in Human Oral Mucosa

Background:

To address a critical gap for application of cancer chronotherapy of when would be the best time(s) for treating an individual cancer patient, we conducted a pilot study to characterize diurnal variations of gene expression in oral mucosal tissue, which is vulnerable to damage from cancer therapies.

Methods:

We conducted RNA-seq assay on individual oral mucosal samples collected from 11 healthy volunteers every 4 hours (6 time points). Using a cosine-based method, we estimated the individual and average values of peak-time and amplitude for each gene. Correlations between gene expression peak-times and age was examined, adjusting for individual’s sleep timing.

Results:

Among candidate gene pathways that are relevant to treatment response, 7 of 16 genes (PER3, CIART, TEF, PER1, PER2, CRY2, ARNTL) involved in circadian regulation and 1 of 118 genes (WEE1) involved in cell cycle regulation achieved p-value ≤ 0.1 and relative amplitude>0.1. The average peak times were approximately 10:15 for PER3, CIART and TEF, 10:45 for PER1, 13:00 for WEE1, PER2 and CRY2, and 19:30 for ARNTL. Ranges in peak times across individuals differed by gene (e.g., 8 hours for PER1; 16.7 hours for WEE1). Older people had later peak times for PER1 (r = 0.77, p = 0.03) and PER3 (r = 0.69, p-value = 0.06).

Conclusion:

In oral mucosa, expression of some genes relevant to treatment response displayed diurnal variation. These genes may be candidates for development of biomarkers for optimizing individual timing of cancer therapy using non-invasively collected oral mucosa.

Clocks, Rhythms, Sex, and Hearts: How Disrupted Circadian Rhythms, Time-of-Day, and Sex Influence Cardiovascular Health

Cardiovascular diseases are the top cause of mortality in the United States, and ischemic heart disease accounts for 16% of all deaths around the world. Modifiable risk factors such as diet and exercise have often been primary targets in addressing these conditions. However, mounting evidence suggests that environmental factors that disrupt physiological rhythms might contribute to the development of these diseases, as well as contribute to increasing other risk factors that are typically associated with cardiovascular disease. Exposure to light at night, transmeridian travel, and social jetlag disrupt endogenous circadian rhythms, which, in turn, alter carefully orchestrated bodily functioning, and elevate the risk of disease and injury. Research into how disrupted circadian rhythms affect physiology and behavior has begun to reveal the intricacies of how seemingly innocuous environmental and social factors have dramatic consequences on mammalian physiology and behavior. Despite the new focus on the importance of circadian rhythms, and how disrupted circadian rhythms contribute to cardiovascular diseases, many questions in this field remain unanswered. Further, neither time-of-day nor sex as a biological variable have been consistently and thoroughly taken into account in previous studies of circadian rhythm disruption and cardiovascular disease. In this review, we will first discuss biological rhythms and the master temporal regulator that controls these rhythms, focusing on the cardiovascular system, its rhythms, and the pathology associated with its disruption, while emphasizing the importance of the time-of-day as a variable that directly affects outcomes in controlled studies, and how temporal data will inform clinical practice and influence personalized medicine. Finally, we will discuss evidence supporting the existence of sex differences in cardiovascular function and outcomes following an injury, and highlight the need for consistent inclusion of both sexes in studies that aim to understand cardiovascular function and improve cardiovascular health.

Extensive tissue-specific expression variation and novel regulators underlying circadian behavior

Natural genetic variation affects circadian rhythms across the evolutionary tree, but the underlying molecular mechanisms are poorly understood. We investigated population-level, molecular circadian clock variation by generating >700 tissue-specific transcriptomes of Drosophila melanogaster (w1118 ) and 141 Drosophila Genetic Reference Panel (DGRP) lines. This comprehensive circadian gene expression atlas contains >1700 cycling genes including previously unknown central circadian clock components and tissue-specific regulators. Furthermore, >30% of DGRP lines exhibited aberrant circadian gene expression, revealing abundant genetic variation-mediated, intertissue circadian expression desynchrony. Genetic analysis of one line with the strongest deviating circadian expression uncovered a novel cry mutation that, as shown by protein structural modeling and brain immunohistochemistry, disrupts the light-driven flavin adenine dinucleotide cofactor photoreduction, providing in vivo support for the importance of this conserved photoentrainment mechanism. Together, our study revealed pervasive tissue-specific circadian expression variation with genetic variants acting upon tissue-specific regulatory networks to generate local gene expression oscillations.

Minimally Invasive Ways of Determining Circadian Rhythms in Humans

Circadian rhythm exerts a critical role in mammalian health and disease. A malfunctioning circadian clock can be a consequence, as well as the cause of several pathophysiologies. Clinical therapies and research may also be influenced by the clock. Since the most suitable manner of revealing this rhythm in humans is not yet established, we discuss existing methods and seek to determine the most feasible ones.

A Characteristic Expression Pattern of Core Circadian Genes in the Diurnal Tree Shrew

The day-active tree shrew may serve as an animal model of human-like diurnal rhythms. However, the molecular basis for circadian rhythms in this species has remained unclear. In the present study, we investigated the expression patterns of core circadian genes involved in transcriptional/translational feedback loops (TTFLs) in both central and peripheral tissues of the tree shrew. The expression of 12 core circadian genes exhibited similar rhythmic patterns in the olfactory bulb, prefrontal cortex, hippocampus, and cerebellum, while the hypothalamus exhibited the weakest oscillations. The rhythms in peripheral tissues, especially the liver, were much more robust than those in brain tissues. ARNTL and NPAS2 were weakly rhythmic in brain tissues but exhibited almost the strongest rhythmicity in peripheral tissues. CLOCK and CRY2 exhibited the weakest rhythms in both central and peripheral tissues, while NR1D1 and CIART exhibited robust rhythms in both tissues. Most of these circadian genes were highly expressed at light/dark transitions in both brain and peripheral tissues, such as ARNTL and NPAS2 peaking at dusk while PERs peaking at dawn. Additionally, the peripheral clock was phase-advanced relative to the brain clock, as there was a significant advance (2-4 h) for PER3, DBP, NR1D1 and NR1D2. Furthermore, these genes exhibited an anti-phasic relationship between the diurnal tree shrew and the nocturnal mouse (i.e., 12-h phasing differential). Collectively, our findings demonstrate a characteristic expression pattern of core circadian genes in the tree shrew, which may provide a means for elucidating molecular mechanisms of diurnal rhythms.

Effects of photoperiod on circadian clock genes in skin contribute to the regulation of hair follicle cycling of Inner Mongolia white cashmere goats

This study investigated how the effects of photoperiod on circadian clock genes in the skin contribute to the regulation of hair follicle cycling of Inner Mongolia white cashmere goats. Twenty-four female (non-pregnant) Inner Mongolia white cashmere goats aged 1 - 1.5 years old with similar live weights (mean, 20.36 ± 2.63 kg) were randomly allocated into two groups: a natural daily photoperiod group (NDPP group:10 - 16 hr Light, n = 12) and a short daily photoperiod group (SDPP group: 7 hr Light:17 hr Dark, n = 12). All goats were housed in individual pens from May 15 to October 15, 2015 and were fed the same diets. We detected the mRNA expression of brain and muscle arnt-like protein-1 (Bmal1), circadian locomotor output control kaput (Clock), cryptochrome-1 (Cry1), period homolog-1 (Per1) and Rev-erbα genes in the goat skin. ANOVA revealed a significant 24 hr (10:00 hr, 14:00 hr, 18:00 hr, 22:00 hr, 02:00 hr, 06:00 hr, 10:00 hr) variation between the SDPP and NDPP groups for three months (July, September, and October). In summary, the current results confirm that an intrinsic oscillating molecular clock exists in goat skin, and that the clock is important for potential timing mechanisms at the anagen phase of hair follicles, which would contribute to the regulation mechanisms of hair follicle cell proliferation.

Ten-Second Electrophysiology: Evaluation of the 3DEP Platform for high-speed, high-accuracy cell analysis

Electrical correlates of the physiological state of a cell, such as membrane conductance and capacitance, as well as cytoplasm conductivity, contain vital information about cellular function, ion transport across the membrane, and propagation of electrical signals. They are, however, difficult to measure; gold-standard techniques are typically unable to measure more than a few cells per day, making widespread adoption difficult and limiting statistical reproducibility. We have developed a dielectrophoretic platform using a disposable 3D electrode geometry that accurately (r² > 0.99) measures mean electrical properties of populations of ~20,000 cells, by taking parallel ensemble measurements of cells at 20 frequencies up to 45 MHz, in (typically) ten seconds. This allows acquisition of ultra-high-resolution (100-point) DEP spectra in under two minutes. Data acquired from a wide range of cells – from platelets to large cardiac cells - benchmark well with patch-clamp-data. These advantages are collectively demonstrated in a longitudinal (same-animal) study of rapidly-changing phenomena such as ultradian (2–3 hour) rhythmicity in whole blood samples of the common vole (Microtus arvalis), taken from 10 µl tail-nick blood samples and avoiding sacrifice of the animal that is typically required in these studies.

The importance of determining circadian parameters in pharmacological studies

In mammals, most molecular and cellular processes show circadian changes, leading to daily variations in physiology and ultimately in behaviour. Such daily variations induce a temporal coordination of processes that is essential to ensure homeostasis and health. Thus, it is of no surprise that pharmacokinetics (PK) and pharmacodynamics (PD) of many drugs are also subject to circadian variations, profoundly affecting their efficacy and tolerability. Understanding how circadian rhythms influence drug PK, PD, and toxicity might significantly improve treatment efficacy and decrease related side effects. Therefore, it is essential to take circadian variations into account and to determine circadian parameters in pharmacological studies, especially when drugs have a short half-life or target rhythmic processes. This review provides an overview of the current knowledge on circadian rhythms and their relevance to the field of pharmacology. Methodologies to evaluate circadian rhythms in vitro, in rodent models and in humans, from experimental to computational approaches, are described and discussed. Lastly, we aim at alerting the scientific, medical, and regulatory communities to the relevance of the physiological time, as a key parameter to be considered when designing pharmacological studies. This will eventually lead to more successful preclinical and clinical trials and pave the way to a more personalized treatment to the benefit of the patients.

Coupling the Circadian Clock to Homeostasis: The Role of Period in Timing Physiology

A plethora of physiological processes show stable and synchronized daily oscillations that are either driven or modulated by biological clocks. A circadian pacemaker located in the suprachiasmatic nucleus of the ventral hypothalamus coordinates 24-hour oscillations of central and peripheral physiology with the environment. The circadian clockwork involved in driving rhythmic physiology is composed of various clock genes that are interlocked via a complex feedback loop to generate precise yet plastic oscillations of ∼24 hours. This review focuses on the specific role of the core clockwork gene Period1 and its paralogs on intra-oscillator and extra-oscillator functions, including, but not limited to, hippocampus-dependent processes, cardiovascular function, appetite control, as well as glucose and lipid homeostasis. Alterations in Period gene function have been implicated in a wide range of physical and mental disorders. At the same time, a variety of conditions including metabolic disorders also impact clock gene expression, resulting in circadian disruptions, which in turn often exacerbates the disease state.

LimoRhyde: A Flexible Approach for Differential Analysis of Rhythmic Transcriptome Data

Unraveling the effects of genetic or environmental perturbations on biological rhythms requires detecting changes in rhythmicity across multiple conditions. Although methods to detect rhythmicity in genome-scale data are well established, methods to detect changes in rhythmicity or changes in average expression between experimental conditions are often ad hoc and statistically unreliable. Here we present LimoRhyde (linear models for rhythmicity, design), a flexible approach for analyzing transcriptome data from circadian systems. Borrowing from cosinor regression, LimoRhyde decomposes circadian or zeitgeber time into multiple components to fit a linear model to the expression of each gene. The linear model can accommodate any number of additional experimental variables, whether discrete or continuous, making it straightforward to detect differential rhythmicity and differential expression using state-of-the-art methods for analyzing microarray and RNA-seq data. In this approach, differential rhythmicity corresponds to a statistical interaction between an experimental variable and circadian time, whereas differential expression corresponds to the main effect of an experimental variable while accounting for circadian time. To validate LimoRhyde’s performance, we applied it to simulated data. To demonstrate LimoRhyde’s versatility, we applied it to murine and human circadian transcriptome datasets acquired under various experimental designs. Our results show how LimoRhyde systematizes the analysis of such data, and suggest that LimoRhyde could prove valuable for assessing how circadian systems respond to perturbations.

Rhythms of the Genome: Circadian Dynamics from Chromatin Topology, Tissue-Specific Gene Expression, to Behavior

Circadian rhythms in physiology and behavior evolved to resonate with daily cycles in the external environment. In mammals, organs orchestrate temporal physiology over the 24-h day, which requires extensive gene expression rhythms targeted to the right tissue. Although a core set of gene products oscillates across virtually all cell types, gene expression profiling across tissues over the 24-h day showed that rhythmic gene expression programs are tissue specific. We highlight recent progress in uncovering how the circadian clock interweaves with tissue-specific gene regulatory networks involving functions such as xenobiotic metabolism, glucose homeostasis, and sleep. This progress hinges on not only comprehensive experimental approaches but also computational methods for multivariate analysis of periodic functional genomics data. We emphasize dynamic chromatin interactions as a novel regulatory layer underlying circadian gene transcription, core clock functions, and ultimately behavior. Finally, we discuss perspectives on extending the knowledge of the circadian clock in animals to human chronobiology.

A database of tissue-specific rhythmically expressed human genes has potential applications in circadian medicine

The discovery that half of the mammalian protein-coding genome is regulated by the circadian clock has clear implications for medicine. Recent studies demonstrated that the circadian clock influences therapeutic outcomes in human heart disease and cancer. However, biological time is rarely given clinical consideration. A key barrier is the absence of information on tissue-specific molecular rhythms in the human body. We have applied the cyclic ordering by periodic structure (CYCLOPS) algorithm, designed to reconstruct sample temporal order in the absence of time-of-day information, to the gene expression collection of 13 tissues from 632 human donors. We identified rhythms in gene expression across the body; nearly half of protein-coding genes were shown to be cycling in at least 1 of the 13 tissues analyzed. One thousand of these cycling genes encode proteins that either transport or metabolize drugs or are themselves drug targets. These results provide a useful resource for studying the role of circadian rhythms in medicine and support the idea that biological time might play a role in determining drug response.

Timing of eating in adults across the weight spectrum: Metabolic factors and potential circadian mechanisms

Timing of eating is recognized as a significant contributor to body weight regulation. Disruption of sleep-wake cycles from a predominantly diurnal (daytime) to a delayed (evening) lifestyle leads to altered circadian rhythms and metabolic dysfunction. This article reviews current evidence for timed and delayed eating in individuals of normal weight and those with overweight or obesity: although some findings indicate a benefit of eating earlier in the daytime on weight and/or metabolic outcomes, results have not been uniformly consistent, and more rigorous and longer-duration studies are needed. We also review potential circadian mechanisms underlying the metabolic- and weight-related changes resulting from timed and delayed eating. Further identification of such mechanisms using deep phenotyping is required to determine targets for medical interventions for obesity and for prevention of metabolic syndrome and diabetes, and to inform clinical guidelines regarding eating schedules for management of weight and metabolic disease.

The circadian dynamics of small nucleolar RNA in the mouse liver

The circadian regulation of gene expression allows plants and animals to anticipate predictable environmental changes. While the influence of the circadian clock has recently been shown to extend to ribosome biogenesis, the dynamics and regulation of the many small nucleolar RNA that are required in pre-ribosomal RNA folding and modification are unknown. Using a novel computational method, we show that 18S and 28S pre-rRNA are subject to circadian regulation in a nuclear RNA sequencing time course. A population of snoRNA with circadian expression is identified that is functionally associated with rRNA modification. More generally, we find the abundance of snoRNA known to modify 18S and 28S to be inversely correlated with the abundance of their target. Cyclic patterns in the expression of a number of snoRNA indicate a coordination with rRNA maturation, potentially through an upregulation in their biogenesis, or their release from mature rRNA at the end of the previous cycle of rRNA maturation, in antiphase with the diurnal peak in pre-rRNA. Few cyclic snoRNA have cyclic host genes, indicating the action of regulatory mechanisms in addition to transcriptional activation of the host gene. For highly expressed independently transcribed snoRNA, we find a characteristic RNA polymerase II and H3K4me3 signature that correlates with mean snoRNA expression over the day.

Evidence for widespread dysregulation of circadian clock progression in human cancer

The ubiquitous daily rhythms in mammalian physiology are guided by progression of the circadian clock. In mice, systemic disruption of the clock can promote tumor growth. In vitro, multiple oncogenes can disrupt the clock. However, due to the difficulties of studying circadian rhythms in solid tissues in humans, whether the clock is disrupted within human tumors has remained unknown. We sought to determine the state of the circadian clock in human cancer using publicly available transcriptome data. We developed a method, called the clock correlation distance (CCD), to infer circadian clock progression in a group of samples based on the co-expression of 12 clock genes. Our method can be applied to modestly sized datasets in which samples are not labeled with time of day and coverage of the circadian cycle is incomplete. We used the method to define a signature of clock gene co-expression in healthy mouse organs, then validated the signature in healthy human tissues. By then comparing human tumor and non-tumor samples from twenty datasets of a range of cancer types, we discovered that clock gene co-expression in tumors is consistently perturbed. Subsequent analysis of data from clock gene knockouts in mice suggested that perturbed clock gene co-expression in human cancer is not caused solely by the inactivation of clock genes. Furthermore, focusing on lung cancer, we found that human lung tumors showed systematic changes in expression in a large set of genes previously inferred to be rhythmic in healthy lung. Our findings suggest that clock progression is dysregulated in many solid human cancers and that this dysregulation could have broad effects on circadian physiology within tumors. In addition, our approach opens the door to using publicly available data to infer circadian clock progression in a multitude of human phenotypes.

The retinol-binding protein receptor STRA6 regulates diurnal insulin responses

It has long been appreciated that insulin action is closely tied to circadian rhythms. However, the mechanisms that dictate diurnal insulin sensitivity in metabolic tissues are not well understood. Retinol-binding protein 4 (RBP4) has been implicated as a driver of insulin resistance in rodents and humans, and it has become an attractive drug target in type II diabetes. RBP4 is synthesized primarily in the liver where it binds retinol and transports it to tissues throughout the body. The retinol-RBP4 complex (holo-RBP) can be recognized by a cell-surface receptor known as stimulated by retinoic acid 6 (STRA6), which transports retinol into cells. Coupled to retinol transport, holo-RBP can activate STRA6-driven Janus kinase (JAK) signaling and downstream induction of signal transducer and activator of transcription (STAT) target genes. STRA6 signaling in white adipose tissue has been shown to inhibit insulin receptor responses. Here, we examined diurnal rhythmicity of the RBP4/STRA6 signaling axis and investigated whether STRA6 is necessary for diurnal variations in insulin sensitivity. We show that adipose tissue STRA6 undergoes circadian patterning driven in part by the nuclear transcription factor REV-ERBα. Furthermore, STRA6 is necessary for diurnal rhythmicity of insulin action and JAK/STAT signaling in adipose tissue. These findings establish that holo-RBP and its receptor STRA6 are potent regulators of diurnal insulin responses and suggest that the holo-RBP/STRA6 signaling axis may represent a novel therapeutic target in type II diabetes.

Machine learning identifies a compact gene set for monitoring the circadian clock in human blood

BACKGROUND

The circadian clock and the daily rhythms it produces are crucial for human health, but are often disrupted by the modern environment. At the same time, circadian rhythms may influence the efficacy and toxicity of therapeutics and the metabolic response to food intake. Developing treatments for circadian dysfunction, as well as optimizing the daily timing of treatments for other health conditions, will require a simple and accurate method to monitor the molecular state of the circadian clock.

METHODS

Here we used a recently developed method called ZeitZeiger to predict circadian time (CT, time of day according to the circadian clock) from genome-wide gene expression in human blood.

RESULTS

In cross-validation on 498 samples from 60 individuals across three publicly available datasets, ZeitZeiger predicted CT in single samples with a median absolute error of 2.1 h. The predictor trained on all 498 samples used 15 genes, only two of which are part of the core circadian clock. By then applying ZeitZeiger to 475 additional samples from the same three datasets, we quantified how the circadian clock in the blood was affected by various perturbations to the sleep-wake and light-dark cycles. Finally, we extended ZeitZeiger (1) to handle intra-individual variation by making predictions based on multiple samples taken a known time apart, and (2) to handle inter-individual variation by personalizing predictions based on samples from the respective individual. Each of these strategies improved prediction of CT by ~20%.

CONCLUSIONS

Our results are an important step towards precision circadian medicine. In addition, our generalizable extensions to ZeitZeiger may be applicable to the growing number of biological datasets that contain multiple observations per individual.

Comparative Analysis of Vertebrate Diurnal/Circadian Transcriptomes

From photosynthetic bacteria to mammals, the circadian clock evolved to track diurnal rhythms and enable organisms to anticipate daily recurring changes such as temperature and light. It orchestrates a broad spectrum of physiology such as the sleep/wake and eating/fasting cycles. While we have made tremendous advances in our understanding of the molecular details of the circadian clock mechanism and how it is synchronized with the environment, we still have rudimentary knowledge regarding its connection to help regulate diurnal physiology. One potential reason is the sheer size of the output network. Diurnal/circadian transcriptomic studies are reporting that around 10% of the expressed genome is rhythmically controlled. Zebrafish is an important model system for the study of the core circadian mechanism in vertebrate. As Zebrafish share more than 70% of its genes with human, it could also be an additional model in addition to rodent for exploring the diurnal/circadian output with potential for translational relevance. Here we performed comparative diurnal/circadian transcriptome analysis with established mouse liver and other tissue datasets. First, by combining liver tissue sampling in a 48h time series, transcription profiling using oligonucleotide arrays and bioinformatics analysis, we profiled rhythmic transcripts and identified 2609 rhythmic genes. The comparative analysis revealed interesting features of the output network regarding number of rhythmic genes, proportion of tissue specific genes and the extent of transcription factor family expression. Undoubtedly, the Zebrafish model system will help identify new vertebrate outputs and their regulators and provides leads for further characterization of the diurnal cis-regulatory network.