Tumor circadian clock strength influences metastatic potential and predicts patient prognosis in luminal A breast cancer.

Studies in shift workers and model organisms link circadian disruption to breast cancer. However, molecular circadian rhythms in noncancerous and cancerous human breast tissues and their clinical relevance are largely unknown. We reconstructed rhythms informatically, integrating locally collected, time-stamped biopsies with public datasets. For noncancerous breast tissue, inflammatory, epithelial-mesenchymal transition (EMT), and estrogen responsiveness pathways show circadian modulation. Among tumors, clock correlation analysis demonstrates subtype-specific changes in circadian organization. Luminal A organoids and informatic ordering of luminal A samples exhibit continued, albeit dampened and reprogrammed rhythms. However, CYCLOPS magnitude, a measure of global rhythm strength, varied widely among luminal A samples. Cycling of EMT pathway genes was markedly increased in high-magnitude luminal A tumors. Surprisingly, patients with high-magnitude tumors had reduced 5-y survival. Correspondingly, 3D luminal A cultures show reduced invasion following molecular clock disruption. This study links subtype-specific circadian disruption in breast cancer to EMT, metastatic potential, and prognosis.

Exploring phylogeny to find the function of sleep.

During sleep, animals do not eat, reproduce or forage. Sleeping animals are vulnerable to predation. Yet, the persistence of sleep despite evolutionary pressures, and the deleterious effects of sleep deprivation, indicate that sleep serves a function or functions that cannot easily be bypassed. Recent research demonstrates sleep to be phylogenetically far more pervasive than previously appreciated; it is possible that the very first animals slept. Here, we give an overview of sleep across various species, with the aim of determining its original purpose. Sleep exists in animals without cephalized nervous systems and can be influenced by non-neuronal signals, including those associated with metabolic rhythms. Together, these observations support the notion that sleep serves metabolic functions in neural and non-neural tissues.

An in silico genome-wide screen for circadian clock strength in human samples.

MOTIVATION: Years of time-series gene expression studies have built a strong understanding of clock-controlled pathways across species. However, comparatively little is known about how 'non-clock' pathways influence clock function. We need a strong understanding of clock-coupled pathways in human tissues to better appreciate the links between disease and clock function. RESULTS: We developed a new computational approach to explore candidate pathways coupled to the clock in human tissues. This method, termed LTM, is an in silico screen to infer genetic influences on circadian clock function. LTM uses natural variation in gene expression in human data and directly links gene expression variation to clock strength independent of longitudinal data. We applied LTM to three human skin and one melanoma datasets and found that the cell cycle is the top candidate clock-coupled pathway in healthy skin. In addition, we applied LTM to thousands of tumor samples from 11 cancer types in the TCGA database and found that extracellular matrix organization-related pathways are tightly associated with the clock strength in humans. Further analysis shows that clock strength in tumor samples is correlated with the proportion of cancer-associated fibroblasts and endothelial cells. Therefore, we show both the power of LTM in predicting clock-coupled pathways and classify factors associated with clock strength in human tissues. AVAILABILITY AND IMPLEMENTATION: LTM is available on GitHub (https://github.com/gangwug/LTMR) and figshare (https://figshare.com/articles/software/LTMR/21217604) to facilitate its use. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational and experimental insights into the circadian effects of SIRT1.

The circadian clock orchestrates 24-h rhythms in physiology in most living organisms. At the molecular level, the dogma is that circadian oscillations are based on a negative transcriptional feedback loop. Recent studies found the NAD+-dependent histone deacetylase, SIRT1, directly regulates acetylation status of clock components and influences circadian amplitude in cells. While Nakahata et al. [Nakahata Y, Kaluzova M (2008) Cell 134:329-340] reported that loss of SIRT1 increases amplitude through BMAL1 acetylation, Asher et al. [Asher G, Gatfield D (2008) Cell 134:317-328] reported that loss of SIRT1 decreases amplitude through an increase in acetylated PER2. To address this SIRT1 paradox, we developed a circadian enzymatic model. Predictions from this model and experimental validation strongly align with the findings of Asher et al., with PER2 as the primary target of SIRT1. Further, the model suggested SIRT1 influences BMAL1 expression through actions on PGC1α. We validated this finding experimentally. Thus, our computational and experimental approaches suggest SIRT1 positively regulates clock function through actions on PER2 and PGC1α.

Population-level rhythms in human skin with implications for circadian medicine.

Skin is the largest organ in the body and serves important barrier, regulatory, and sensory functions. The epidermal layer shows rhythmic physiological responses to daily environmental variation (e.g., DNA repair). We investigated the role of the circadian clock in the transcriptional regulation of epidermis using a hybrid experimental design, in which a limited set of human subjects (n = 20) were sampled throughout the 24-h cycle and a larger population (n = 219) were sampled once. We found a robust circadian oscillator in human epidermis at the population level using pairwise correlations of clock and clock-associated genes in 298 epidermis samples. We then used CYCLOPS to reconstruct the temporal order of all samples, and identified hundreds of rhythmically expressed genes at the population level in human epidermis. We compared these results with published time-series skin data from mice and found a strong concordance in circadian phase across species for both transcripts and pathways. Furthermore, like blood, epidermis is readily accessible and a potential source of biomarkers. Using ZeitZeiger, we identified a biomarker set for human epidermis that is capable of reporting circadian phase to within 3 hours from a single sample. In summary, we show rhythms in human epidermis that persist at the population scale and describe a path to develop robust single-sample circadian biomarkers.

CYCLOPS reveals human transcriptional rhythms in health and disease.

Circadian rhythms modulate many aspects of physiology. Knowledge of the molecular basis of these rhythms has exploded in the last 20 years. However, most of these data are from model organisms, and translation to clinical practice has been limited. Here, we present an approach to identify molecular rhythms in humans from thousands of unordered expression measurements. Our algorithm, cyclic ordering by periodic structure (CYCLOPS), uses evolutionary conservation and machine learning to identify elliptical structure in high-dimensional data. From this structure, CYCLOPS estimates the phase of each sample. We validated CYCLOPS using temporally ordered mouse and human data and demonstrated its consistency on human data from two independent research sites. We used this approach to identify rhythmic transcripts in human liver and lung, including hundreds of drug targets and disease genes. Importantly, for many genes, the circadian variation in expression exceeded variation from genetic and other environmental factors. We also analyzed hepatocellular carcinoma samples and show these solid tumors maintain circadian function but with aberrant output. Finally, to show how this method can catalyze medical translation, we show that dosage time can temporally segregate efficacy from dose-limiting toxicity of streptozocin, a chemotherapeutic drug. In sum, these data show the power of CYCLOPS and temporal reconstruction in bridging basic circadian research and clinical medicine.

MetaCycle: an integrated R package to evaluate periodicity in large scale data.

Detecting periodicity in large scale data remains a challenge. While efforts have been made to identify best of breed algorithms, relatively little research has gone into integrating these methods in a generalizable method. Here, we present MetaCycle, an R package that incorporates ARSER, JTK_CYCLE and Lomb-Scargle to conveniently evaluate periodicity in time-series data. MetaCycle has two functions, meta2d and meta3d, designed to analyze two-dimensional and three-dimensional time-series datasets, respectively. Meta2d implements N-version programming concepts using a suite of algorithms and integrating their results. AVAILABILITY AND IMPLEMENTATION: MetaCycle package is available on the CRAN repository (https://cran.r-project.org/web/packages/MetaCycle/index.html) and GitHub (https://github.com/gangwug/MetaCycle). CONTACT: hogenesch@gmail.comSupplementary information: Supplementary data are available at Bioinformatics online.

Machine learning helps identify CHRONO as a circadian clock component.

Over the last decades, researchers have characterized a set of "clock genes" that drive daily rhythms in physiology and behavior. This arduous work has yielded results with far-reaching consequences in metabolic, psychiatric, and neoplastic disorders. Recent attempts to expand our understanding of circadian regulation have moved beyond the mutagenesis screens that identified the first clock components, employing higher throughput genomic and proteomic techniques. In order to further accelerate clock gene discovery, we utilized a computer-assisted approach to identify and prioritize candidate clock components. We used a simple form of probabilistic machine learning to integrate biologically relevant, genome-scale data and ranked genes on their similarity to known clock components. We then used a secondary experimental screen to characterize the top candidates. We found that several physically interact with known clock components in a mammalian two-hybrid screen and modulate in vitro cellular rhythms in an immortalized mouse fibroblast line (NIH 3T3). One candidate, Gene Model 129, interacts with BMAL1 and functionally represses the key driver of molecular rhythms, the BMAL1/CLOCK transcriptional complex. Given these results, we have renamed the gene CHRONO (computationally highlighted repressor of the network oscillator). Bi-molecular fluorescence complementation and co-immunoprecipitation demonstrate that CHRONO represses by abrogating the binding of BMAL1 to its transcriptional co-activator CBP. Most importantly, CHRONO knockout mice display a prolonged free-running circadian period similar to, or more drastic than, six other clock components. We conclude that CHRONO is a functional clock component providing a new layer of control on circadian molecular dynamics.

Design and analysis of large-scale biological rhythm studies: a comparison of algorithms for detecting periodic signals in biological data.

MOTIVATION: To discover and study periodic processes in biological systems, we sought to identify periodic patterns in their gene expression data. We surveyed a large number of available methods for identifying periodicity in time series data and chose representatives of different mathematical perspectives that performed well on both synthetic data and biological data. Synthetic data were used to evaluate how each algorithm responds to different curve shapes, periods, phase shifts, noise levels and sampling rates. The biological datasets we tested represent a variety of periodic processes from different organisms, including the cell cycle and metabolic cycle in Saccharomyces cerevisiae, circadian rhythms in Mus musculus and the root clock in Arabidopsis thaliana. RESULTS: From these results, we discovered that each algorithm had different strengths. Based on our findings, we make recommendations for selecting and applying these methods depending on the nature of the data and the periodic patterns of interest. Additionally, these results can also be used to inform the design of large-scale biological rhythm experiments so that the resulting data can be used with these algorithms to detect periodic signals more effectively.

Circadian transcriptome of pancreatic adenocarcinoma unravels chronotherapeutic targets.

Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer characterized by a poor outcome and an increasing incidence. A significant majority (>80%) of newly diagnosed cases are deemed unresectable, leaving chemotherapy as the sole viable option, though with only moderate success. This necessitates the identification of improved therapeutic options for PDA. We hypothesized that there are temporal variations in cancer-relevant processes within PDA tumors, offering insights into the optimal timing of drug administration - a concept termed chronotherapy. In this study, we explored the presence of the circadian transcriptome in PDA using patient-derived organoids and validated these findings by comparing PDA data from The Cancer Genome Atlas with noncancerous healthy pancreas data from GTEx. Several PDA-associated pathways (cell cycle, stress response, Rho GTPase signaling) and cancer driver hub genes (EGFR and JUN) exhibited a cancer-specific rhythmic pattern intricately linked to the circadian clock. Through the integration of multiple functional measurements for rhythmic cancer driver genes, we identified top chronotherapy targets and validated key findings in molecularly divergent pancreatic cancer cell lines. Testing the chemotherapeutic efficacy of clinically relevant drugs further revealed temporal variations that correlated with drug-target cycling. Collectively, our study unravels the PDA circadian transcriptome and highlights a potential approach for optimizing chrono-chemotherapeutic efficacy.

A glial circadian gene expression atlas reveals cell type and disease-specific reprogramming in response to amyloid pathology or aging.

While circadian rhythm disruption may promote neurodegenerative disease, how aging and neurodegenerative pathology impact circadian gene expression patterns in different brain cell types is unknown. Here, we used translating ribosome affinity purification methods to define the circadian translatomes of astrocytes, microglia, and bulk cerebral cortex, in healthy mouse brain and in the settings of amyloid-beta plaque pathology or aging. Our data reveal that glial circadian translatomes are highly cell type-specific and exhibit profound, context-dependent reprogramming of rhythmic transcripts in response to amyloid pathology or aging. Transcripts involved in glial activation, immunometabolism, and proteostasis, as well as nearly half of all Alzheimer Disease (AD)-associated risk genes, displayed circadian oscillations, many of which were altered by pathology. Amyloid-related differential gene expression was also dependent on time of day. Thus, circadian rhythms in gene expression are cell- and context dependent and provide important insights into glial gene regulation in health, AD, and aging.

Sleep is not just for the brain: transcriptional responses to sleep in peripheral tissues.

BACKGROUND: Many have assumed that the primary function of sleep is for the brain. We evaluated the molecular consequences of sleep and sleep deprivation outside the brain, in heart and lung. Using microarrays we compared gene expression in tissue from sleeping and sleep deprived mice euthanized at the same diurnal times. RESULTS: In each tissue, nearly two thousand genes demonstrated statistically significant differential expression as a function of sleep/wake behavioral state. To mitigate the influence of an artificial deprivation protocol, we identified a subset of these transcripts as specifically sleep-enhanced or sleep-repressed by requiring that their expression also change over the course of unperturbed sleep. 3% and 6% of the assayed transcripts showed "sleep specific" changes in the lung and heart respectively. Sleep specific transcripts in these tissues demonstrated highly significant overlap and shared temporal dynamics. Markers of cellular stress and the unfolded protein response were reduced during sleep in both tissues. These results mirror previous findings in brain. Sleep-enhanced pathways reflected the unique metabolic functions of each tissue. Transcripts related to carbohydrate and sulfur metabolic processes were enhanced by sleep in the lung, and collectively favor buffering from oxidative stress. DNA repair and protein metabolism annotations were significantly enriched among the sleep-enhanced transcripts in the heart. Our results also suggest that sleep may provide a Zeitgeber, or synchronizing cue, in the lung as a large cluster of transcripts demonstrated systematic changes in inter-animal variability as a function of both sleep duration and circadian time. CONCLUSION: Our data support the notion that the molecular consequences of sleep/wake behavioral state extend beyond the brain to include peripheral tissues. Sleep state induces a highly overlapping response in both heart and lung. We conclude that sleep enhances organ specific molecular functions and that it has a ubiquitous role in reducing cellular metabolic stress in both brain and peripheral tissues. Finally, our data suggest a novel role for sleep in synchronizing transcription in peripheral tissues.

Developmental emergence of sleep rhythms enables long-term memory in Drosophila.

In adulthood, sleep-wake rhythms are one of the most prominent behaviors under circadian control. However, during early life, sleep is spread across the 24-hour day. The mechanism through which sleep rhythms emerge, and consequent advantage conferred to a juvenile animal, is unknown. In the second-instar Drosophila larvae (L2), like in human infants, sleep is not under circadian control. We identify the precise developmental time point when the clock begins to regulate sleep in Drosophila, leading to emergence of sleep rhythms in early third-instars (L3). At this stage, a cellular connection forms between DN1a clock neurons and arousal-promoting Dh44 neurons, bringing arousal under clock control to drive emergence of circadian sleep. Last, we demonstrate that L3 but not L2 larvae exhibit long-term memory (LTM) of aversive cues and that this LTM depends upon deep sleep generated once sleep rhythms begin. We propose that the developmental emergence of circadian sleep enables more complex cognitive processes, including the onset of enduring memories.

Integration of genome-scale data identifies candidate sleep regulators.

STUDY OBJECTIVES: Genetics impacts sleep, yet, the molecular mechanisms underlying sleep regulation remain elusive. In this study, we built machine learning models to predict sleep genes based on their similarity to genes that are known to regulate sleep. METHODS: We trained a prediction model on thousands of published datasets, representing circadian, immune, sleep deprivation, and many other processes, using a manually curated list of 109 sleep genes. RESULTS: Our predictions fit with prior knowledge of sleep regulation and identified key genes and pathways to pursue in follow-up studies. As an example, we focused on the NF-κB pathway and showed that chronic activation of NF-κB in a genetic mouse model impacted the sleep-wake patterns. CONCLUSION: Our study highlights the power of machine learning in integrating prior knowledge and genome-wide data to study genetic regulation of complex behaviors such as sleep.

Subtype-specific circadian clock dysregulation modulates breast cancer biology, invasiveness, and prognosis.

Studies in shift workers and model organisms link circadian disruption to breast cancer. However, molecular rhythms in non-cancerous and cancerous human breast tissues are largely unknown. We reconstructed rhythms informatically, integrating locally collected, time-stamped biopsies with public datasets. For non-cancerous tissue, the inferred order of core-circadian genes matches established physiology. Inflammatory, epithelial-mesenchymal transition (EMT), and estrogen responsiveness pathways show circadian modulation. Among tumors, clock correlation analysis demonstrates subtype-specific changes in circadian organization. Luminal A organoids and informatic ordering of Luminal A samples exhibit continued, albeit disrupted rhythms. However, CYCLOPS magnitude, a measure of global rhythm strength, varied widely among Luminal A samples. Cycling of EMT pathway genes was markedly increased in high-magnitude Luminal A tumors. Patients with high-magnitude tumors had reduced 5-year survival. Correspondingly, 3D Luminal A cultures show reduced invasion following molecular clock disruption. This study links subtype-specific circadian disruption in breast cancer to EMT, metastatic potential, and prognosis.

Multi-tissue transcriptional changes and core circadian clock disruption following intensive care.

Objective: Both critical illness and current care have been hypothesized to upset daily rhythms and impair molecular circadian function. However, the influence of critical illness on clock function in different tissues and on circadian output genes are unknown. Here we evaluate the effect of critical care and illness on transcription, focusing on the functional organization of the core circadian oscillator. Methods: We downloaded RNAseq count data from the Genotype-Tissue Expression (GTEx) project. Treating mechanical ventilation as a marker for intensive care, we stratified samples into acute death (AD) and intensive care (IC) groups based on the documented Hardy Death Scale. We restricted our analysis to the 25 tissues with >50 samples in each group. Using the edgeR package and controlling for collection center, gender, and age, we identified transcripts differentially expressed between the AD and IC groups. Overrepresentation and enrichment methods were used to identify gene sets modulated by intensive care across tissues. For each tissue, we then calculated the delta clock correlation distance (ΔCCD), a comparative measure of the functional organization of the core circadian oscillator, in the both the AD and IC groups. The statistical significance of the ΔCCD was assessed by permutation, modifying a pre-existing R package to control for confounding variables. Results: Intensive care, as marked by ventilation, significantly modulated the expression of thousands of genes. Transcripts that were modulated in ≥75% of tissues were enriched for genes involved in mitochondrial energetics, cellular stress, metabolism, and notably circadian regulation. Transcripts that were more markedly affected, in ≥10 tissues, were enriched for inflammation, complement and immune pathways. Oscillator organization, as assessed by ΔCCD, was significantly reduced in the intensive care group in 11/25 tissues. Conclusion: Our findings support the hypothesis that patients in intensive care have impaired molecular circadian rhythms. Tissues involved in metabolism and energetics demonstrated the most marked changes in oscillator organization. In adipose tissue, there was a significant overlap between transcripts previously established to be modulated by sleep deprivation and fasting with those modulated by critical care. This work suggests that intensive care protocols that restore sleep/wake and nutritional rhythms may be of benefit.

Treatment of a patient with a circadian sleep-wake disorder using a combination of melatonin and metoprolol.

CITATION: Circadian rhythm sleep-wake disorders result from the lack of synchronization between endogenous circadian rhythms and daily environmental or behavioral cycles. Current treatment of circadian rhythm sleep-wake disorders relies on strengthening normal zeitgebers, or temporal cues, through the combination of strict behavioral modification, controlled light exposure, and supplemental melatonin or melatonin receptor agonists. These therapies can be difficult to maintain and are supported with only limited clinical outcome data. The effectiveness of exogenous melatonin, in particular, may be reduced by the patient's continued production of endogenous melatonin with a temporal pattern that is not conducive to the desired sleep schedule. Here we describe the case of a single, sighted patient with a circadian rhythm sleep-wake disorder who benefited from the combined use of a beta blocker to suppress endogenous melatonin secretion along with the timed administration of exogenous melatonin. We suggest that the positive results obtained justify further study of this mechanism-guided approach. CITATION: Gehrman PR, Anafi RC. Treatment of a patient with a circadian sleep-wake disorder using a combination of melatonin and metoprolol. J Clin Sleep Med. 2021;17(10):2121-2124.

CRY1-CBS binding regulates circadian clock function and metabolism.

Circadian disruption influences metabolic health. Metabolism modulates circadian function. However, the mechanisms coupling circadian rhythms and metabolism remain poorly understood. Here, we report that cystathionine ß-synthase (CBS), a central enzyme in one-carbon metabolism, functionally interacts with the core circadian protein cryptochrome 1 (CRY1). In cells, CBS augments CRY1-mediated repression of the CLOCK/BMAL1 complex and shortens circadian period. Notably, we find that mutant CBS-I278T protein, the most common cause of homocystinuria, does not bind CRY1 or regulate its repressor activity. Transgenic CbsZn/Zn  mice, while maintaining circadian locomotor activity period, exhibit reduced circadian power and increased expression of E-BOX outputs. CBS function is reciprocally influenced by CRY1 binding. CRY1 modulates enzymatic activity of the CBS. Liver extracts from Cry1-/- mice show reduced CBS activity that normalizes after the addition of exogenous wild-type (WT) CRY1. Metabolomic analysis of WT, CbsZn/Zn , Cry1-/- , and Cry2-/- samples highlights the metabolic importance of endogenous CRY1. We observed temporal variation in one-carbon and transsulfuration pathways attributable to CRY1-induced CBS activation. CBS-CRY1 binding provides a post-translational switch to modulate cellular circadian physiology and metabolic control.