skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Detecting Rhythmic Gene Expression in Single Cell Transcriptomics
Abstract An autonomous, environmentally-synchronizable circadian rhythm is a ubiquitous feature of life on Earth. In multicellular organisms, this rhythm is generated by a transcription–translation feedback loop present in nearly every cell that drives daily expression of thousands of genes in a tissue–dependent manner. Identifying the genes that are under circadian control can elucidate the mechanisms by which physiological processes are coordinated in multicellular organisms. Today, transcriptomic profiling at the single-cell level provides an unprecedented opportunity to understand the function of cell-level clocks. However, while many cycling detection algorithms have been developed to identify genes under circadian control in bulk transcriptomic data, it is not known how best to adapt these algorithms to single-cell RNAseq data. Here, we benchmark commonly used circadian detection methods on their reliability and efficiency when applied to single cell RNAseq data. Our results provide guidance on adapting existing cycling detection methods to the single-cell domain, and elucidate opportunities for more robust and efficient rhythm detection in single-cell data. We also propose a subsampling procedure combined with harmonic regression as an efficient, reliable strategy to detect circadian genes in the single–cell setting.  more » « less
Award ID(s):
1764421
PAR ID:
10525843
Author(s) / Creator(s):
;
Publisher / Repository:
bioRxiv
Date Published:
Format(s):
Medium: X
Institution:
bioRxiv
Sponsoring Org:
National Science Foundation
More Like this
  1. The circadian rhythm drives the oscillatory expression of thousands of genes across all tissues, coordinating physiological processes. The effect of this rhythm on health has generated increasing interest in discovering genes under circadian control by searching for periodic patterns in transcriptomic time-series experiments. While algorithms for detecting cycling transcripts have advanced, there remains little guidance quantifying the effect of experimental design and analysis choices on cycling detection accuracy. We present TimeTrial, a user-friendly benchmarking framework using both real and synthetic data to investigate cycle detection algorithms’ performance and improve circadian experimental design. Results show that the optimal choice of analysis method depends on the sampling scheme, noise level, and shape of the waveform of interest and provides guidance on the impact of sampling frequency and duration on cycling detection accuracy. The TimeTrial software is freely available for download and may also be accessed through a web interface. By supplying a tool to vary and optimize experimental design considerations, TimeTrial will enhance circadian transcriptomics studies. 
    more » « less
  2. Abstract As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher’s performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research. 
    more » « less
  3. Motivation The circadian rhythm drives the oscillatory expression of thousands of genes across all tissues. The recent revolution in high-throughput transcriptomics, coupled with the significant implications of the circadian clock for human health, has sparked an interest in circadian profiling studies to discover genes under circadian control. Result We present TimeCycle: a topology-based rhythm detection method designed to identify cycling transcripts. For a given time-series, the method reconstructs the state space using time-delay embedding, a data transformation technique from dynamical systems theory. In the embedded space, Takens’ theorem proves that the dynamics of a rhythmic signal will exhibit circular patterns. The degree of circularity of the embedding is calculated as a persistence score using persistent homology, an algebraic method for discerning the topological features of data. By comparing the persistence scores to a bootstrapped null distribution, cycling genes are identified. Results in both synthetic and biological data highlight Time-Cycle’s ability to identify cycling genes across a range of sampling schemes, number of replicates, and missing data. Comparison to competing methods highlights their relative strengths, providing guidance as to the optimal choice of cycling detection method. Availability and Implementation A fully documented open-source R package implementing Time-Cycle is available at: https://nesscoder.github.io/TimeCycle/. 
    more » « less
  4. Abstract In mammals, T-cell migration is under circadian control, likely to anticipate daily rhythms in infection risk. Glucocorticoids are a major controller of circadian processes and malnutrition is associated with increased glucocorticoid secretion. Previous studies suggest malnutrition may impart a “super-quiescent” phenotype to T-cells, enabling a greater number of naïve T-cells to survive short-term malnutrition albeit with diminished function. Thus, we hypothesize that malnourished T-cells may conserve energy by disengaging from rhythmic migration under circadian control and/or foregoing migration to reside in the bone marrow instead. To test this hypothesis, the total number of nucleated cells and naïve CD4+ and CD8+ T-cells in the blood, spleen, bone marrow, and brachial and mesenteric lymph nodes were enumerated by flow cytometry every four hours over the course of one day from control and malnourished mice. Additionally, expression levels of CD127 and CXCR4 in both T-cell populations and the concentration of glucocorticoids in the blood were assessed. A better understanding of how malnutrition affects the circadian rhythm of T-cell migration will not only help identify the mechanisms of how circadian rhythms work, but also how organisms’ circadian rhythms change in response to malnutrition. This knowledge of how malnutrition disrupts the circadian rhythm of T-cells may help improve vaccination strategies in malnourished children. Supported by NSF-MRI [DBI- 1920116] NSF -RUI [IOS-1951881] 
    more » « less
  5. Circadian rhythms are internal processes repeating approximately every 24 hours in living organisms. The dominant circadian pacemaker is synchronized to the environmental light-dark cycle. Other circadian pacemakers, which can have noncanonical circadian mechanisms, are revealed by arousing stimuli, such as scheduled feeding, palatable meals and running wheel access, or methamphetamine administration. Organisms also have ultradian rhythms, which have periods shorter than circadian rhythms. However, the biological mechanism, origin, and functional significance of ultradian rhythms are not well-elucidated. The dominant circadian rhythm often masks ultradian rhythms; therefore, we disabled the canonical circadian clock of mice by knocking out Per1/2/3 genes, where Per1 and Per2 are essential components of the mammalian light-sensitive circadian mechanism. Furthermore, we recorded wheel-running activity every minute under constant darkness for 272 days. We then investigated rhythmic components in the absence of external influences, applying unique multiscale time-resolved methods to analyze the oscillatory dynamics with time-varying frequencies. We found four rhythmic components with periods of ∼17 h, ∼8 h, ∼4 h, and ∼20 min. When the ∼17-h rhythm was prominent, the ∼8-h rhythm was of low amplitude. This phenomenon occurred periodically approximately every 2-3 weeks. We found that the ∼4-h and ∼20-min rhythms were harmonics of the ∼8-h rhythm. Coupling analysis of the ridge-extracted instantaneous frequencies revealed strong and stable phase coupling from the slower oscillations (∼17, ∼8, and ∼4 h) to the faster oscillations (∼20 min), and weak and less stable phase coupling in the reverse direction and between the slower oscillations. Together, this study elucidated the relationship between the oscillators in the absence of the canonical circadian clock, which is critical for understanding their functional significance. These studies are essential as disruption of circadian rhythms contributes to diseases, such as cancer and obesity, as well as mood disorders. 
    more » « less