Michael E. Hughes

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

Significance We generated high-resolution multiorgan expression data showing that nearly half of all genes in the mouse genome oscillate with circadian rhythm somewhere in the body. Such widespread transcriptional oscillations have not been previously reported in mammals. Applying pathway analysis, we observed new clock-mediated spatiotemporal relationships. Moreover, we found a majority of best-selling drugs in the United States target circadian gene products. Many of these drugs have relatively short half-lives, and our data predict which may benefit from timed dosing. To characterize the role of the circadian clock in mouse physiology and behavior, we used RNA-seq and DNA arrays to quantify the transcriptomes of 12 mouse organs over time. We found 43% of all protein coding genes showed circadian rhythms in transcription somewhere in the body, largely in an organ-specific manner. In most organs, we noticed the expression of many oscillating genes peaked during transcriptional “rush hours” preceding dawn and dusk. Looking at the genomic landscape of rhythmic genes, we saw that they clustered together, were longer, and had more spliceforms than nonoscillating genes. Systems-level analysis revealed intricate rhythmic orchestration of gene pathways throughout the body. We also found oscillations in the expression of more than 1,000 known and novel noncoding RNAs (ncRNAs). Supporting their potential role in mediating clock function, ncRNAs conserved between mouse and human showed rhythmic expression in similar proportions as protein coding genes. Importantly, we also found that the majority of best-selling drugs and World Health Organization essential medicines directly target the products of rhythmic genes. Many of these drugs have short half-lives and may benefit from timed dosage. In sum, this study highlights critical, systemic, and surprising roles of the mammalian circadian clock and provides a blueprint for advancement in chronotherapy.

Harmonics of Circadian Gene Transcription in Mammals

The circadian clock is a molecular and cellular oscillator found in most mammalian tissues that regulates rhythmic physiology and behavior. Numerous investigations have addressed the contribution of circadian rhythmicity to cellular, organ, and organismal physiology. We recently developed a method to look at transcriptional oscillations with unprecedented precision and accuracy using high-density time sampling. Here, we report a comparison of oscillating transcription from mouse liver, NIH3T3, and U2OS cells. Several surprising observations resulted from this study, including a 100-fold difference in the number of cycling transcripts in autonomous cellular models of the oscillator versus tissues harvested from intact mice. Strikingly, we found two clusters of genes that cycle at the second and third harmonic of circadian rhythmicity in liver, but not cultured cells. Validation experiments show that 12-hour oscillatory transcripts occur in several other peripheral tissues as well including heart, kidney, and lungs. These harmonics are lost ex vivo, as well as under restricted feeding conditions. Taken in sum, these studies illustrate the importance of time sampling with respect to multiple testing, suggest caution in use of autonomous cellular models to study clock output, and demonstrate the existence of harmonics of circadian gene expression in the mouse.

JTK_CYCLE: An Efficient Nonparametric Algorithm for Detecting Rhythmic Components in Genome-Scale Data Sets

Circadian rhythms are oscillations of physiology, behavior, and metabolism that have period lengths near 24 hours. In several model organisms and humans, circadian clock genes have been characterized and found to be transcription factors. Because of this, researchers have used microarrays to characterize global regulation of gene expression and algorithmic approaches to detect cycling. This article presents a new algorithm, JTK_CYCLE, designed to efficiently identify and characterize cycling variables in large data sets. Compared with COSOPT and the Fisher’s G test, two commonly used methods for detecting cycling transcripts, JTK_CYCLE distinguishes between rhythmic and nonrhythmic transcripts more reliably and efficiently. JTK_CYCLE’s increased resistance to outliers results in considerably greater sensitivity and specificity. Moreover, JTK_CYCLE accurately measures the period, phase, and amplitude of cycling transcripts, facilitating downstream analyses. Finally, JTK_CYCLE is several orders of magnitude faster than COSOPT, making it ideal for large-scale data sets. JTK_CYCLE was used to analyze legacy data sets including NIH3T3 cells, which have comparatively low amplitude oscillations. JTK_CYCLE’s improved power led to the identification of a novel cluster of RNA-interacting genes whose abundance is under clear circadian regulation. These data suggest that JTK_CYCLE is an ideal tool for identifying and characterizing oscillations in genome-scale data sets.

Deep sequencing the circadian and diurnal transcriptome of Drosophila brain

Eukaryotic circadian clocks include transcriptional/translational feedback loops that drive 24-h rhythms of transcription. These transcriptional rhythms underlie oscillations of protein abundance, thereby mediating circadian rhythms of behavior, physiology, and metabolism. Numerous studies over the last decade have used microarrays to profile circadian transcriptional rhythms in various organisms and tissues. Here we use RNA sequencing (RNA-seq) to profile the circadian transcriptome of Drosophila melanogaster brain from wild-type and period-null clock-defective animals. We identify several hundred transcripts whose abundance oscillates with 24-h periods in either constant darkness or 12 h light/dark diurnal cycles, including several noncoding RNAs (ncRNAs) that were not identified in previous microarray studies. Of particular interest are U snoRNA host genes (Uhgs), a family of diurnal cycling noncoding RNAs that encode the precursors of more than 50 box-C/D small nucleolar RNAs, key regulators of ribosomal biogenesis. Transcriptional profiling at the level of individual exons reveals alternative splice isoforms for many genes whose relative abundances are regulated by either period or circadian time, although the effect of circadian time is muted in comparison to that of period. Interestingly, period loss of function significantly alters the frequency of RNA editing at several editing sites, suggesting an unexpected link between a key circadian gene and RNA editing. We also identify tens of thousands of novel splicing events beyond those previously annotated by the modENCODE Consortium, including several that affect key circadian genes. These studies demonstrate extensive circadian control of ncRNA expression, reveal the extent of clock control of alternative splicing and RNA editing, and provide a novel, genome-wide map of splicing in Drosophila brain.

Circadian expression of clock genes in mouse macrophages, dendritic cells, and B cells

In mammals, circadian and daily rhythms influence nearly all aspects of physiology, ranging from behavior to gene expression. Functional molecular clocks have been described in the murine spleen and splenic NK cells. The aim of our study was to investigate the existence of molecular clock mechanisms in other immune cells. Therefore, we measured the circadian changes in gene expression of clock genes (Per1, Per2, Bmal1, and Clock) and clock-controlled transcription factors (Rev-erbα and Dbp) in splenic enriched macrophages, dendritic cells, and B cells in both mice entrained to a light-dark cycle and under constant environmental conditions. Our study reveals the existence of functional molecular clock mechanisms in splenic macrophages, dendritic cells, and B cells.

High-resolution Time Course Analysis of Gene Expression from Pituitary

In both the suprachiasmatic nucleus (SCN) and peripheral tissues, the circadian oscillator drives rhythmic transcription of downstream target genes. Recently, a number of studies have used DNA microarrays to systematically identify oscillating transcripts in plants, fruit flies, rats, and mice. These studies have identified several dozen to many hundred rhythmically expressed genes by sampling tissues every 4 hours for 1, 2, or more days. To extend this work, we have performed DNA microarray analysis on RNA derived from the mouse pituitary sampled every hour for 2 days. COSOPT and Fishers G-test were used at a false-discovery rate of less than 5% to identify more than 250 genes in the pituitary that oscillate with a 24-hour period length. We found that increasing the frequency of sampling across the circadian day dramatically increased the statistical power of both COSOPT and Fishers G-test, resulting in considerably more high-confidence identifications of rhythmic transcripts than previously described. Finally, to extend the utility of these data sets, a Web-based resource has been constructed (at http://wasabi.itmat.upenn.edu/circa/mouse ) that is freely available to the research community.

Blood Contrast Enhancement with a Novel, Non-Gaseous Nanoparticle Contrast Agent

Modern ultrasound contrast agents primarily comprise microbubble formulations that circulate in the intravascular compartment and are designed to enhance acoustic signals reflected from the blood pool. A variety of shell materials have been utilized to stabilize gas bubbles of the order of 1–10 microns in diameter. Reflectivity from microbubbles is enhanced by resonance and non-linear physical effects. However, the overall efficacy of bubbles as contrast agents must be considered in light of their marked instability to insonification pressures, marked attenuation artifacts, “blooming” effects, and their short circulatory half-life. Low molecular weight gaseous perfluorocarbon formulations have been utilized in vivo because they may offer advantages in formulation and reflectivity. In contrast, higher molecular weight perfluorocarbon emulsions that are liquid at body temperature have been formulated as nongaseous nanoparticle preparations (diameters 100– 300 nanometers), originally for use as blood substitutes. Unfortunately they exhibit low inherent echogenicity and are poor blood pool contrast agents under conditions of conventional 2-D echocardiography or harmonic imaging, or when imaged with color flow or spectral Doppler. Nevertheless, these nanoparticle formulations are chemically inert, manifest long circulatory half-lives, are not destroyed by ultrasonic imaging, and they possess low acoustic attenuation. Such features might still render them of interest as blood pool contrast agents if properly formulated and imaged. Recently, a new ultrasonic imaging modality, Power Doppler Harmonic Imaging (PDHI), has been introduced (4). This technique color-encodes changes in acoustic signal amplitude and motion of ultrasonic scatterers between insonifying pulses. PDHI has been used in a number of clinical studies to assess coronary artery bypass graft patency, tumor blood flow, and myocardial perfusion. In view of the exquisite sensitivity of Doppler for detecting the presence of small scatterers with limited scattering cross-sections as compared to microbubbles (e.g., red blood cells), and the enhanced ability of PDHI to register backscatter power, we hypothesized that certain liquid perfluorocarbon nanoparticle emulsions (5) might be more efficiently detected with this new imaging modality. Furthermore, although we have demonstrated previously that the liquid nanoparticle emulsions do not manifest any appreciable resonance behavior at clinically relevant imaging frequencies, they have performed well as targeted imaging agents in vitro and in vivo over a very broad range of frequencies (5–50 MHz)(6–8). Thus we anticipated that the PDHI method might permit imaging of these nanoparticles in the blood pool without reliance on any intrinsic resonance behavior.

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

Two independent studies, one of them using a computational approach, identified CHRONO, a gene shown to modulate the activity of circadian transcription factors and alter circadian behavior in mice.

Age-Associated Disruption of Molecular Clock Expression in Skeletal Muscle of the Spontaneously Hypertensive Rat

It is well known that spontaneously hypertensive rats (SHR) develop muscle pathologies with hypertension and heart failure, though the mechanism remains poorly understood. Woon et al. (2007) linked the circadian clock gene Bmal1 to hypertension and metabolic dysfunction in the SHR. Building on these findings, we compared the expression pattern of several core-clock genes in the gastrocnemius muscle of aged SHR (80 weeks; overt heart failure) compared to aged-matched control WKY strain. Heart failure was associated with marked effects on the expression of Bmal1, Clock and Rora in addition to several non-circadian genes important in regulating skeletal muscle phenotype including Mck, Ttn and Mef2c. We next performed circadian time-course collections at a young age (8 weeks; pre-hypertensive) and adult age (22 weeks; hypertensive) to determine if clock gene expression was disrupted in gastrocnemius, heart and liver tissues prior to or after the rats became hypertensive. We found that hypertensive/hypertrophic SHR showed a dampening of peak Bmal1 and Rev-erb expression in the liver, and the clock-controlled gene Pgc1α in the gastrocnemius. In addition, the core-clock gene Clock and the muscle-specific, clock-controlled gene Myod1, no longer maintained a circadian pattern of expression in gastrocnemius from the hypertensive SHR. These findings provide a framework to suggest a mechanism whereby chronic heart failure leads to skeletal muscle pathologies; prolonged dysregulation of the molecular clock in skeletal muscle results in altered Clock, Pgc1α and Myod1 expression which in turn leads to the mis-regulation of target genes important for mechanical and metabolic function of skeletal muscle.