Mutations of the circadian clock genes Cry, Per, or Bmal1 have different effects on the transcribed and nontranscribed strands of cycling genes.

One important goal of circadian medicine is to apply time-of-day dosing to improve the efficacy of chemotherapy. However, limited knowledge of how the circadian clock regulates DNA repair presents a challenge to mechanism-based clinical application. We studied time-series genome-wide nucleotide excision repair in liver and kidney of wild type and three different clock mutant genotypes (Cry1-/-Cry2-/-, Per1-/-Per2-/-, and Bmal1-/-). Rhythmic repair on the nontranscribed strand was lost in all three clock mutants. Conversely, rhythmic repair of hundreds of genes on the transcribed strand (TSs) persisted in the livers of Cry1-/-Cry2-/- and Per1-/-Per2-/- mice. We identified a tissue-specific, promoter element-driven repair mode on TSs of collagen and angiogenesis genes in the absence of clock activators or repressors. Furthermore, repair on TSs of thousands of genes was altered when the circadian clock is disrupted. These data contribute to a better understanding of the regulatory role of the circadian clock on nucleotide excision repair in mammals and may be invaluable toward the design of time-aware platinum-based interventions in cancer.

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.

Ontogeny and function of the circadian clock in intestinal organoids.

Circadian rhythms regulate diverse aspects of gastrointestinal physiology ranging from the composition of microbiota to motility. However, development of the intestinal circadian clock and detailed mechanisms regulating circadian physiology of the intestine remain largely unknown. In this report, we show that both pluripotent stem cell-derived human intestinal organoids engrafted into mice and patient-derived human intestinal enteroids possess circadian rhythms and demonstrate circadian phase-dependent necrotic cell death responses to Clostridium difficile toxin B (TcdB). Intriguingly, mouse and human enteroids demonstrate anti-phasic necrotic cell death responses to TcdB. RNA-Seq analysis shows that ~3-10% of the detectable transcripts are rhythmically expressed in mouse and human enteroids. Remarkably, we observe anti-phasic gene expression of Rac1, a small GTPase directly inactivated by TcdB, between mouse and human enteroids, and disruption of Rac1 abolishes clock-dependent necrotic cell death responses. Our findings uncover robust functions of circadian rhythms regulating clock-controlled genes in both mouse and human enteroids governing organism-specific, circadian phase-dependent necrotic cell death responses, and lay a foundation for human organ- and disease-specific investigation of clock functions using human organoids for translational applications.

Obstructive sleep apnea in a mouse model is associated with tissue-specific transcriptomic changes in circadian rhythmicity and mean 24-hour gene expression.

Intermittent hypoxia (IH) is a major clinical feature of obstructive sleep apnea (OSA). The mechanisms that become dysregulated after periods of exposure to IH are unclear, particularly in the early stages of disease. The circadian clock governs a wide array of biological functions and is intimately associated with stabilization of hypoxia-inducible factors (HIFs) under hypoxic conditions. In patients, IH occurs during the sleep phase of the 24-hour sleep-wake cycle, potentially affecting their circadian rhythms. Alterations in the circadian clock have the potential to accelerate pathological processes, including other comorbid conditions that can be associated with chronic, untreated OSA. We hypothesized that changes in the circadian clock would manifest differently in those organs and systems known to be impacted by OSA. Using an IH model to represent OSA, we evaluated circadian rhythmicity and mean 24-hour expression of the transcriptome in 6 different mouse tissues, including the liver, lung, kidney, muscle, heart, and cerebellum, after a 7-day exposure to IH. We found that transcriptomic changes within cardiopulmonary tissues were more affected by IH than other tissues. Also, IH exposure resulted in an overall increase in core body temperature. Our findings demonstrate a relationship between early exposure to IH and changes in specific physiological outcomes. This study provides insight into the early pathophysiological mechanisms associated with IH.

The circadian clock is disrupted in pancreatic cancer.

Disruption of the circadian clock is linked to cancer development and progression. Establishing this connection has proven beneficial for understanding cancer pathogenesis, determining prognosis, and uncovering novel therapeutic targets. However, barriers to characterizing the circadian clock in human pancreas and human pancreatic cancer-one of the deadliest malignancies-have hindered an appreciation of its role in this cancer. Here, we employed normalized coefficient of variation (nCV) and clock correlation analysis in human population-level data to determine the functioning of the circadian clock in pancreas cancer and adjacent normal tissue. We found a substantially attenuated clock in the pancreatic cancer tissue. Then we exploited our existing mouse pancreatic transcriptome data to perform an analysis of the human normal and pancreas cancer samples using a machine learning method, cyclic ordering by periodic structure (CYCLOPS). Through CYCLOPS ordering, we confirmed the nCV and clock correlation findings of an intact circadian clock in normal pancreas with robust cycling of several core clock genes. However, in pancreas cancer, there was a loss of rhythmicity of many core clock genes with an inability to effectively order the cancer samples, providing substantive evidence of a dysregulated clock. The implications of clock disruption were further assessed with a Bmal1 knockout pancreas cancer model, which revealed that an arrhythmic clock caused accelerated cancer growth and worse survival, accompanied by chemoresistance and enrichment of key cancer-related pathways. These findings provide strong evidence for clock disruption in human pancreas cancer and demonstrate a link between circadian disruption and pancreas cancer progression.

duper is a null mutation of Cryptochrome 1 in Syrian hamsters.

The duper mutation is a recessive mutation that shortens the period length of the circadian rhythm in Syrian hamsters. These animals show a large phase shift when responding to light pulses. Limited genetic resources for the Syrian hamster (Mesocricetus auratus) presented a major obstacle to cloning duper. This caused the duper mutation to remain unknown for over a decade. In this study, we did a de novo genome assembly of Syrian hamsters with long-read sequencing data from two different platforms, Pacific Biosciences and Oxford Nanopore Technologies. Using two distinct ecotypes and a fast homozygosity mapping strategy, we identified duper as an early nonsense allele of Cryptochrome 1 (Cry1) leading to a short, unstable protein. CRY1 is known as a highly conserved component of the repressive limb of the core circadian clock. The genome assembly and other genomic datasets generated in this study will facilitate the use of the Syrian hamster in biomedical research.

A genome-wide RNAi screen for modifiers of the circadian clock in human cells.

Two decades of research identified more than a dozen clock genes and defined a biochemical feedback mechanism of circadian oscillator function. To identify additional clock genes and modifiers, we conducted a genome-wide small interfering RNA screen in a human cellular clock model. Knockdown of nearly 1000 genes reduced rhythm amplitude. Potent effects on period length or increased amplitude were less frequent; we found hundreds of these and confirmed them in secondary screens. Characterization of a subset of these genes demonstrated a dosage-dependent effect on oscillator function. Protein interaction network analysis showed that dozens of gene products directly or indirectly associate with known clock components. Pathway analysis revealed these genes are overrepresented for components of insulin and hedgehog signaling, the cell cycle, and the folate metabolism. Coupled with data showing many of these pathways are clock regulated, we conclude the clock is interconnected with many aspects of cellular function.

NF-κB modifies the mammalian circadian clock through interaction with the core clock protein BMAL1.

In mammals, the circadian clock coordinates cell physiological processes including inflammation. Recent studies suggested a crosstalk between these two pathways. However, the mechanism of how inflammation affects the clock is not well understood. Here, we investigated the role of the proinflammatory transcription factor NF-κB in regulating clock function. Using a combination of genetic and pharmacological approaches, we show that perturbation of the canonical NF-κB subunit RELA in the human U2OS cellular model altered core clock gene expression. While RELA activation shortened period length and dampened amplitude, its inhibition lengthened period length and caused amplitude phenotypes. NF-κB perturbation also altered circadian rhythms in the master suprachiasmatic nucleus (SCN) clock and locomotor activity behavior under different light/dark conditions. We show that RELA, like the clock repressor CRY1, repressed the transcriptional activity of BMAL1/CLOCK at the circadian E-box cis-element. Biochemical and biophysical analysis showed that RELA binds to the transactivation domain of BMAL1. These data support a model in which NF-kB competes with CRY1 and coactivator CBP/p300 for BMAL1 binding to affect circadian transcription. This is further supported by chromatin immunoprecipitation analysis showing that binding of RELA, BMAL1 and CLOCK converges on the E-boxes of clock genes. Taken together, these data support a significant role for NF-κB in directly regulating the circadian clock and highlight mutual regulation between the circadian and inflammatory pathways.

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.

Normalized coefficient of variation (nCV): a method to evaluate circadian clock robustness in population scale data.

SUMMARY: Robust oscillation of clock genes is a core feature of the circadian system. Relative amplitude (rAMP) measures the robustness of clock gene oscillations but only works for longitudinal samples. We lack a method for estimating robust oscillations from human samples without labeled time. We show that the normalized coefficient of variation (nCV) of 10 clock genes is linearly correlated with their normalized rAMP, independent of time labels. We found that the mean nCV of clock genes are consistently decreased in tumors compared to nontumors, suggesting a new therapeutic target in cancer treatment by enhancing clock robustness. nCV can provide a simple measure of the clock robustness in population-level datasets. AVAILABILITY AND IMPLEMENTATION: The nCV package (https://github.com/gangwug/nCV) and web application (https://github.com/gangwug/nCVapp) are available on the GitHub repository. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A large-scale study reveals 24-h operational rhythms in hospital treatment.

Hospitals operate 24 h a day, and it is assumed that important clinical decisions occur continuously around the clock. However, many aspects of hospital operation occur at specific times of day, including medical team rounding and shift changes. It is unclear whether this impacts patient care, as no studies have addressed this. We analyzed the daily distribution of ∼500,000 doses of 12 separate drugs in 1,546 inpatients at a major children's hospital in the United States from 2010 to 2017. We tracked both order time (when a care provider places an electronic request for a drug) and dosing time (when the patient receives the drug). Order times were time-of-day-dependent, marked by distinct morning-time surges and overnight lulls. Nearly one-third of all 103,847 orders for treatment were placed between 8:00 AM and 12:00 PM. First doses from each order were also rhythmic but shifted by 2 h. These 24-h rhythms in orders and first doses were remarkably consistent across drugs, diagnosis, and hospital units. This rhythm in hospital medicine coincided with medical team rounding time, not necessarily immediate medical need. Lastly, we show that the clinical response to hydralazine, an acute antihypertensive, is dosing time-dependent and greatest at night, when the fewest doses were administered. The prevailing dogma is that hospital treatment is administered as needed regardless of time of day. Our findings challenge this notion and reveal a potential operational barrier to best clinical care.

CAMPAREE: a robust and configurable RNA expression simulator.

BACKGROUND: The accurate interpretation of RNA-Seq data presents a moving target as scientists continue to introduce new experimental techniques and analysis algorithms. Simulated datasets are an invaluable tool to accurately assess the performance of RNA-Seq analysis methods. However, existing RNA-Seq simulators focus on modeling the technical biases and artifacts of sequencing, rather than on simulating the original RNA samples. A first step in simulating RNA-Seq is to simulate RNA. RESULTS: To fill this need, we developed the Configurable And Modular Program Allowing RNA Expression Emulation (CAMPAREE), a simulator using empirical data to simulate diploid RNA samples at the level of individual molecules. We demonstrated CAMPAREE's use for generating idealized coverage plots from real data, and for adding the ability to generate allele-specific data to existing RNA-Seq simulators that do not natively support this feature. CONCLUSIONS: Separating input sample modeling from library preparation/sequencing offers added flexibility for both users and developers to mix-and-match different sample and sequencing simulators to suit their specific needs. Furthermore, the ability to maintain sample and sequencing simulators independently provides greater agility to incorporate new biological findings about transcriptomics and new developments in sequencing technologies. Additionally, by simulating at the level of individual molecules, CAMPAREE has the potential to model molecules transcribed from the same genes as a heterogeneous population of transcripts with different states of degradation and processing (splicing, editing, etc.). CAMPAREE was developed in Python, is open source, and freely available at https://github.com/itmat/CAMPAREE .

Polycystin-1 regulates bone development through an interaction with the transcriptional coactivator TAZ.

Polycystin-1 (PC1), encoded by the PKD1 gene that is mutated in the autosomal dominant polycystic kidney disease, regulates a number of processes including bone development. Activity of the transcription factor RunX2, which controls osteoblast differentiation, is reduced in Pkd1 mutant mice but the mechanism governing PC1 activation of RunX2 is unclear. PC1 undergoes regulated cleavage that releases its C-terminal tail (CTT), which translocates to the nucleus to modulate transcriptional pathways involved in proliferation and apoptosis. We find that the cleaved CTT of PC1 (PC1-CTT) stimulates the transcriptional coactivator TAZ (Wwtr1), an essential coactivator of RunX2. PC1-CTT physically interacts with TAZ, stimulating RunX2 transcriptional activity in pre-osteoblast cells in a TAZ-dependent manner. The PC1-CTT increases the interaction between TAZ and RunX2 and enhances the recruitment of the p300 transcriptional co-regulatory protein to the TAZ/RunX2/PC1-CTT complex. Zebrafish injected with morpholinos directed against pkd1 manifest severe bone calcification defects and a curly tail phenotype. Injection of messenger RNA (mRNA) encoding the PC1-CTT into pkd1-morphant fish restores bone mineralization and reduces the severity of the curly tail phenotype. These effects are abolished by co-injection of morpholinos directed against TAZ. Injection of mRNA encoding a dominant-active TAZ construct is sufficient to rescue both the curly tail phenotype and the skeletal defects observed in pkd1-morpholino treated fish. Thus, TAZ constitutes a key mechanistic link through which PC1 mediates its physiological functions.

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

Cisplatin-DNA adduct repair of transcribed genes is controlled by two circadian programs in mouse tissues.

Cisplatin is a major cancer chemotherapeutic drug. It kills cancer cells by damaging their DNA, mainly in the form of Pt-d(GpG) diadducts. However, it also has serious side effects, including nephrotoxicity and hepatotoxicity that limit its usefulness. Chronotherapy is taking circadian time into account during therapy to improve the therapeutic index, by improving efficacy and/or limiting toxicity. To this end, we tested the impact of clock time on excision repair of cisplatin-induced DNA damage at single-nucleotide resolution across the genome in mouse kidney and liver. We found that genome repair is controlled by two circadian programs. Repair of the transcribed strand (TS) of active, circadian-controlled genes is dictated by each gene's phase of transcription, which falls across the circadian cycle with prominent peaks at dawn and dusk. In contrast, repair of the nontranscribed strand (NTS) of all genes, repair of intergenic DNA, and global repair overall peaks at Zeitgeber time ZT08, as basal repair capacity, which is controlled by the circadian clock, peaks at this circadian time. Consequently, the TS and NTS of many genes are repaired out of phase. As most cancers are thought to have defective circadian rhythms, these results suggest that future research on timed dosage of cisplatin could potentially reduce damage to healthy tissue and improve its therapeutic index.

mTOR signaling regulates central and peripheral circadian clock function.

The circadian clock coordinates physiology and metabolism. mTOR (mammalian/mechanistic target of rapamycin) is a major intracellular sensor that integrates nutrient and energy status to regulate protein synthesis, metabolism, and cell growth. Previous studies have identified a key role for mTOR in regulating photic entrainment and synchrony of the central circadian clock in the suprachiasmatic nucleus (SCN). Given that mTOR activities exhibit robust circadian oscillations in a variety of tissues and cells including the SCN, here we continued to investigate the role of mTOR in orchestrating autonomous clock functions in central and peripheral circadian oscillators. Using a combination of genetic and pharmacological approaches we show that mTOR regulates intrinsic clock properties including period and amplitude. In peripheral clock models of hepatocytes and adipocytes, mTOR inhibition lengthens period and dampens amplitude, whereas mTOR activation shortens period and augments amplitude. Constitutive activation of mTOR in Tsc2-/-fibroblasts elevates levels of core clock proteins, including CRY1, BMAL1 and CLOCK. Serum stimulation induces CRY1 upregulation in fibroblasts in an mTOR-dependent but Bmal1- and Period-independent manner. Consistent with results from cellular clock models, mTOR perturbation also regulates period and amplitude in the ex vivo SCN and liver clocks. Further, mTOR heterozygous mice show lengthened circadian period of locomotor activity in both constant darkness and constant light. Together, these results support a significant role for mTOR in circadian timekeeping and in linking metabolic states to circadian clock functions.

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.

Genome-wide effect of pulmonary airway epithelial cell-specific Bmal1 deletion.

Pulmonary airway epithelial cells (AECs) form a critical interface between host and environment. We investigated the role of the circadian clock using mice bearing targeted deletion of the circadian gene brain and muscle ARNT-like 1 (Bmal1) in AECs. Pulmonary neutrophil infiltration, biomechanical function, and responses to influenza infection were all disrupted. A circadian time-series RNA sequencing study of laser-captured AECs revealed widespread disruption in genes of the core circadian clock and output pathways regulating cell metabolism (lipids and xenobiotics), extracellular matrix, and chemokine signaling, but strikingly also the gain of a novel rhythmic transcriptome in Bmal1-targeted cells. Many of the rhythmic components were replicated in primary AECs cultured in air-liquid interface, indicating significant cell autonomy for control of pulmonary circadian physiology. Finally, we found that metabolic cues dictate phasing of the pulmonary clock and circadian responses to immunologic challenges. Thus, the local circadian clock in AECs is vital in lung health by coordinating major cell processes such as metabolism and immunity.-Zhang, Z. Hunter, L., Wu, G., Maidstone, R., Mizoro, Y., Vonslow, R., Fife, M., Hopwood, T., Begley, N., Saer, B., Wang, P., Cunningham, P., Baxter, M., Durrington, H., Blaikley, J. F., Hussell, T., Rattray, M., Hogenesch, J. B., Gibbs, J., Ray, D. W., Loudon, A. S. I. Genome-wide effect of pulmonary airway epithelial cell-specific Bmal1 deletion.

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.

Experimental and statistical reevaluation provides no evidence for Drosophila courtship song rhythms.

From 1980 to 1992, a series of influential papers reported on the discovery, genetics, and evolution of a periodic cycling of the interval between Drosophila male courtship song pulses. The molecular mechanisms underlying this periodicity were never described. To reinitiate investigation of this phenomenon, we previously performed automated segmentation of songs but failed to detect the proposed rhythm [Arthur BJ, et al. (2013) BMC Biol 11:11; Stern DL (2014) BMC Biol 12:38]. Kyriacou et al. [Kyriacou CP, et al. (2017) Proc Natl Acad Sci USA 114:1970-1975] report that we failed to detect song rhythms because (i) our flies did not sing enough and (ii) our segmenter did not identify many of the song pulses. Kyriacou et al. manually annotated a subset of our recordings and reported that two strains displayed rhythms with genotype-specific periodicity, in agreement with their original reports. We cannot replicate this finding and show that the manually annotated data, the original automatically segmented data, and a new dataset provide no evidence for either the existence of song rhythms or song periodicity differences between genotypes. Furthermore, we have reexamined our methods and analysis and find that our automated segmentation method was not biased to prevent detection of putative song periodicity. We conclude that there is no evidence for the existence of Drosophila courtship song rhythms.