Lauren J. Francey

CYCLOPS reveals human transcriptional rhythms in health and disease

Significance Circadian rhythms influence most aspects of physiology and behavior. However, how do we apply this knowledge in medicine? Identifying molecular mechanisms in humans is challenging as existing large-scale datasets rarely include time of day. To address this problem, we combine understanding of periodic structure, evolutionary conservation, and unsupervised machine learning to order unordered human biopsy data along a periodic cycle. We show this works using ordered mouse and human data and that it gives consistent results when applied to populations on different continents. Then, we investigate molecular rhythms in normal human lung and liver and cancerous liver. Finally, we demonstrate proof of concept by finding the best time to administer a chemotherapeutic drug in an animal model. 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.

Role for LSM genes in the regulation of circadian rhythms.

Significance There is increasing evidence that previously considered core constituents of multi-subunit complexes involved in RNA processing play regulatory rather than passive roles in the control of gene expression, but specific signaling pathways in which they participate are not known. Here we show that SM-like (LSM) genes, which encode core components of the spliceosome, are regulated by the circadian clock and control clock function in plants and mammals, revealing convergent evolutionary mechanisms mediating posttranscriptional regulation of circadian networks across kingdoms. Growing evidence suggests that core spliceosomal components differentially affect RNA processing of specific genes; however, whether changes in the levels or activities of these factors control specific signaling pathways is largely unknown. Here we show that some SM-like (LSM) genes, which encode core components of the spliceosomal U6 small nuclear ribonucleoprotein complex, regulate circadian rhythms in plants and mammals. We found that the circadian clock regulates the expression of LSM5 in Arabidopsis plants and several LSM genes in mouse suprachiasmatic nucleus. Further, mutations in LSM5 or LSM4 in Arabidopsis, or down-regulation of LSM3, LSM5, or LSM7 expression in human cells, lengthens the circadian period. Although we identified changes in the expression and alternative splicing of some core clock genes in Arabidopsis lsm5 mutants, the precise molecular mechanism causing period lengthening remains to be identified. Genome-wide expression analysis of either a weak lsm5 or a strong lsm4 mutant allele in Arabidopsis revealed larger effects on alternative splicing than on constitutive splicing. Remarkably, large splicing defects were not observed in most of the introns evaluated using RNA-seq in the strong lsm4 mutant allele used in this study. These findings support the idea that some LSM genes play both regulatory and constitutive roles in RNA processing, contributing to the fine-tuning of specific signaling pathways.

Guidelines for Genome-Scale Analysis of Biological Rhythms

Genome biology approaches have made enormous contributions to our understanding of biological rhythms, particularly in identifying outputs of the clock, including RNAs, proteins, and metabolites, whose abundance oscillates throughout the day. These methods hold significant promise for future discovery, particularly when combined with computational modeling. However, genome-scale experiments are costly and laborious, yielding “big data” that are conceptually and statistically difficult to analyze. There is no obvious consensus regarding design or analysis. Here we discuss the relevant technical considerations to generate reproducible, statistically sound, and broadly useful genome-scale data. Rather than suggest a set of rigid rules, we aim to codify principles by which investigators, reviewers, and readers of the primary literature can evaluate the suitability of different experimental designs for measuring different aspects of biological rhythms. We introduce CircaInSilico, a web-based application for generating synthetic genome biology data to benchmark statistical methods for studying biological rhythms. Finally, we discuss several unmet analytical needs, including applications to clinical medicine, and suggest productive avenues to address them.

KPNB1 mediates PER/CRY nuclear translocation and circadian clock function

Regulated nuclear translocation of the PER/CRY repressor complex is critical for negative feedback regulation of the circadian clock of mammals. However, the precise molecular mechanism is not fully understood. Here, we report that KPNB1, an importin β component of the ncRNA repressor of nuclear factor of activated T cells (NRON) ribonucleoprotein complex, mediates nuclear translocation and repressor function of the PER/CRY complex. RNAi depletion of KPNB1 traps the PER/CRY complex in the cytoplasm by blocking nuclear entry of PER proteins in human cells. KPNB1 interacts mainly with PER proteins and directs PER/CRY nuclear transport in a circadian fashion. Interestingly, KPNB1 regulates the PER/CRY nuclear entry and repressor function, independently of importin α, its classical partner. Moreover, inducible inhibition of the conserved Drosophila importin β in lateral neurons abolishes behavioral rhythms in flies. Collectively, these data show that KPNB1 is required for timely nuclear import of PER/CRY in the negative feedback regulation of the circadian clock. DOI: http://dx.doi.org/10.7554/eLife.08647.001

The Local Edge Machine: inference of dynamic models of gene regulation

We present a novel approach, the Local Edge Machine, for the inference of regulatory interactions directly from time-series gene expression data. We demonstrate its performance, robustness, and scalability on in silico datasets with varying behaviors, sizes, and degrees of complexity. Moreover, we demonstrate its ability to incorporate biological prior information and make informative predictions on a well-characterized in vivo system using data from budding yeast that have been synchronized in the cell cycle. Finally, we use an atlas of transcription data in a mammalian circadian system to illustrate how the method can be used for discovery in the context of large complex networks.

Population level rhythms in human skin: implications for circadian medicine

Skin is the largest organ in the body and serves important barrier, regulatory, and sensory functions. Like other tissues, skin is subject to temporal fluctuations in physiological responses under both homeostatic and stressed states. To gain insight into these fluctuations, we investigated the role of the circadian clock in the transcriptional regulation of epidermis using a hybrid experimental design, where a limited set of human subjects (n=20) were sampled throughout the 24 h cycle and a larger population (n=219) were sampled once. By looking at pairwise correlations of core clock genes in 298 skin samples, we found a robust circadian oscillator in skin at the population level. Encouraged by this, we used CYCLOPS to reconstruct the temporal order of all samples and identified hundreds of rhythmically-expressed genes at the population level in human skin. We compared these results with published time-series skin data from mouse and show strong concordance in circadian phase across species for both transcripts and pathways. Further, like blood, skin is readily accessible and a potential source of biomarkers. Using ZeitZeiger, we identified a biomarker set for human skin that is capable of reporting circadian phase to within 3 h from a single sample. In summary, we show rhythms in human skin that persist at the population scale and a path to develop robust single-sample circadian biomarkers. One Sentence Summary Human epidermis shows strong circadian rhythms at the population scale and provides a better source for developing robust, single-sample circadian phase biomarkers than human blood.

A population-based human enCYCLOPedia for circadian medicine

The discovery that half of the mammalian protein-coding genome is clock-regulated has clear implications for medicine. Indeed, recent studies demonstrate time-of-day impact on therapeutic outcomes in human heart disease and cancer. Yet biological time is rarely given clinical consideration. A key barrier is the absence of information on the what and where of molecular rhythms in the human body. Here, we have applied CYCLOPS, an algorithm designed to reconstruct sample order in the absence of time-of-day information, to the GTEx collection of 632 human donors contributing 4,292 RNA-seq samples from 13 distinct human tissue types. We identify rhythms in expression across the body that persist at the population-level. This includes a set of ‘ubiquitous cyclers’ comprised of well-established circadian clock factors but also many genes without prior circadian context. Among thousands of tissue-divergent rhythms, we discover a set of genes robustly oscillating in cardiovascular tissue, including key drug targets relevant to heart disease. These results also have implications for genetic studies where circadian variability may have masked genetic influence. It is our hope that the human enCYCLOPedia helps drive the translation of circadian biology into prospective clinical trials in cardiology and many other therapeutic areas. One Sentence Summary Bioinformatic analyses on thousands of human tissue samples reveals an enCYCLOPedia of rhythmic gene expression across the body and identifies key translational opportunities for circadian medicine in cardiovascular disease.

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

Bioinformatic analyses of human tissue samples were used to build a database of rhythmically expressed genes across the body. Body timing Although the existence of circadian clock–dependent modulation of gene expression in humans has been known for more than a decade, the relevance of the circadian clock in drug response and therapeutic outcome has been only recently appreciated. Now, Ruben et al. used an algorithm called cyclic ordering by periodic structure (CYCLOPS) to create a database of cycling genes in 13 human tissues. The authors show that several rhythmically expressed genes code for known drug targets or for proteins involved in drug transport and metabolism. The data represent a useful resource for circadian medicine and strengthen the notion that circadian rhythms should be considered when determining therapeutic interventions. 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.

Computational and experimental insights into the circadian effects of SIRT1

Significance Circadian rhythms are oscillations with a period of 24 h inherent in numerous biological processes. The prevailing model describing the molecular machinery of the mammalian circadian rhythms is governed by a core set of genes including BMAL1, CLOCK, and PER2. However, emerging evidence highlighted the importance of a new but contradictory role of SIRT1 on the core circadian machinery. Briefly, the two major contradictions are that (i) BMAL1 and PER2 are directly deacetylated by SIRT1 and (ii) loss of SIRT1 leads to dampening and robustness of the circadian amplitude. These contradictions remain without explanation or resolution. Our findings provide support for PER2 as a direct target of SIRT1 and identify a potential role for PGC1α in the circadian network. 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α.

Circadian Dysregulation: The Next Frontier in Obstructive Sleep Apnea Research

Objective To review the effects of the circadian clock on homeostasis, the functional interaction between the circadian clock and hypoxia-inducible factors, and the role of circadian dysregulation in the progression of cardiopulmonary disease in obstructive sleep apnea (OSA). Data Sources The MEDLINE database was accessed through PubMed. Review Methods A general review is presented on molecular pathways disrupted in OSA, circadian rhythms and the role of the circadian clock, hypoxia signaling, crosstalk between the circadian and hypoxia systems, the role of the circadian clock in cardiovascular disease, and implications for practice. Studies included in this State of the Art Review demonstrate the potential contribution of the circadian clock and hypoxia in animal models or human disease. Conclusions Molecular crosstalk between the circadian clock and hypoxia-inducible factors has not been evaluated in disease models of OSA. Implications for Practice Pediatric OSA is highly prevalent and, if left untreated, may lead to cardiopulmonary sequelae. Changes in inflammatory markers that normally demonstrate circadian rhythmicity are also seen among patients with OSA. Hypoxia-inducible transcription factors interact with core circadian clock transcription factors; however, the interplay between these pathways has not been elucidated in the cardiopulmonary system. This gap in knowledge hinders our ability to identify potential biomarkers of OSA and develop alternative therapeutic strategies. A deeper understanding of the mechanisms by which OSA impinges on clock function and the impact of clock dysregulation on the cardiopulmonary system may lead to future advancements for the care of patients with OSA. The aim of this review is to shed light on this important clinical topic.