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.

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.

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.

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.

HDAC8 mutations in Cornelia de Lange syndrome affect the cohesin acetylation cycle.

Cornelia de Lange syndrome (CdLS) is a dominantly inherited congenital malformation disorder, caused by mutations in the cohesin-loading protein NIPBL for nearly 60% of individuals with classical CdLS, and by mutations in the core cohesin components SMC1A (~5%) and SMC3 (<1%) for a smaller fraction of probands. In humans, the multisubunit complex cohesin is made up of SMC1, SMC3, RAD21 and a STAG protein. These form a ring structure that is proposed to encircle sister chromatids to mediate sister chromatid cohesion and also has key roles in gene regulation. SMC3 is acetylated during S-phase to establish cohesiveness of chromatin-loaded cohesin, and in yeast, the class I histone deacetylase Hos1 deacetylates SMC3 during anaphase. Here we identify HDAC8 as the vertebrate SMC3 deacetylase, as well as loss-of-function HDAC8 mutations in six CdLS probands. Loss of HDAC8 activity results in increased SMC3 acetylation and inefficient dissolution of the 'used' cohesin complex released from chromatin in both prophase and anaphase. SMC3 with retained acetylation is loaded onto chromatin, and chromatin immunoprecipitation sequencing analysis demonstrates decreased occupancy of cohesin localization sites that results in a consistent pattern of altered transcription seen in CdLS cell lines with either NIPBL or HDAC8 mutations.

Role for LSM genes in the regulation of circadian rhythms.

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.

Loss-of-function HDAC8 mutations cause a phenotypic spectrum of Cornelia de Lange syndrome-like features, ocular hypertelorism, large fontanelle and X-linked inheritance.

Cornelia de Lange syndrome (CdLS) is a multisystem genetic disorder with distinct facies, growth failure, intellectual disability, distal limb anomalies, gastrointestinal and neurological disease. Mutations in NIPBL, encoding a cohesin regulatory protein, account for >80% of cases with typical facies. Mutations in the core cohesin complex proteins, encoded by the SMC1A, SMC3 and RAD21 genes, together account for ∼5% of subjects, often with atypical CdLS features. Recently, we identified mutations in the X-linked gene HDAC8 as the cause of a small number of CdLS cases. Here, we report a cohort of 38 individuals with an emerging spectrum of features caused by HDAC8 mutations. For several individuals, the diagnosis of CdLS was not considered prior to genomic testing. Most mutations identified are missense and de novo. Many cases are heterozygous females, each with marked skewing of X-inactivation in peripheral blood DNA. We also identified eight hemizygous males who are more severely affected. The craniofacial appearance caused by HDAC8 mutations overlaps that of typical CdLS but often displays delayed anterior fontanelle closure, ocular hypertelorism, hooding of the eyelids, a broader nose and dental anomalies, which may be useful discriminating features. HDAC8 encodes the lysine deacetylase for the cohesin subunit SMC3 and analysis of the functional consequences of the missense mutations indicates that all cause a loss of enzymatic function. These data demonstrate that loss-of-function mutations in HDAC8 cause a range of overlapping human developmental phenotypes, including a phenotypically distinct subgroup of CdLS.

PECONPI: a novel software for uncovering pathogenic copy number variations in non-syndromic sensorineural hearing loss and other genetically heterogeneous disorders.

This report describes an algorithm developed to predict the pathogenicity of copy number variants (CNVs) in large sample cohorts. CNVs (genomic deletions and duplications) are found in healthy individuals and in individuals with genetic diagnoses, and differentiation of these two classes of CNVs can be challenging and usually requires extensive manual curation. We have developed PECONPI, an algorithm to assess the pathogenicity of CNVs based on gene content and CNV frequency. This software was applied to a large cohort of patients with genetically heterogeneous non-syndromic hearing loss to score and rank each CNV based on its relative pathogenicity. Of 636 individuals tested, we identified the likely underlying etiology of the hearing loss in 14 (2%) of the patients (1 with a homozygous deletion, 7 with a deletion of a known hearing loss gene and a point mutation on the trans allele and 6 with a deletion larger than 1 Mb). We also identified two probands with smaller deletions encompassing genes that may be functionally related to their hearing loss. The ability of PECONPI to determine the pathogenicity of CNVs was tested on a second genetically heterogeneous cohort with congenital heart defects (CHDs). It successfully identified a likely etiology in 6 of 355 individuals (2%). We believe this tool is useful for researchers with large genetically heterogeneous cohorts to help identify known pathogenic causes and novel disease genes.

PSMD11 modulates circadian clock function through PER and CRY nuclear translocation.

The molecular circadian clock is regulated by a transcriptional translational feedback loop. However, the post-translational control mechanisms are less understood. The NRON complex is a large ribonucleoprotein complex, consisting of a lncRNA and several proteins. Components of the complex play a distinct role in regulating protein phosphorylation, synthesis, stability, and translocation in cellular processes. This includes the NFAT and the circadian clock pathway. PSMD11 is a component of the NRON complex and a lid component of the 26S proteasome. Among the PSMD family members, PSMD11 has a more specific role in circadian clock function. Here, we used cell and biochemical approaches and characterized the role of PSMD11 in regulating the stability and nuclear translocation of circadian clock proteins. We used size exclusion chromatography to enrich the NRON complex in the cytosolic and nuclear fractions. More specifically, PSMD11 knockdown affected the abundance of PER2 and CRY2 proteins and the nuclear translocation of CRY1. This changed the relative abundance of CRY1 and CRY2 in the nucleus. Thus, this work defines the role of PSMD11 in the NRON complex regulating the nuclear translocation of circadian repressors, thereby enabling cellular circadian oscillations.

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.

Genome-wide SNP genotyping identifies the Stereocilin (STRC) gene as a major contributor to pediatric bilateral sensorineural hearing impairment.

Hearing loss is the most prevalent sensory perception deficit in humans, affecting 1/500 newborns, can be syndromic or nonsyndromic and is genetically heterogeneous. Nearly 80% of inherited nonsyndromic bilateral sensorineural hearing loss (NBSNHI) is autosomal recessive. Although many causal genes have been identified, most are minor contributors, except for GJB2, which accounts for nearly 50% of all recessive cases of severe to profound congenital NBSNHI in some populations. More than 60% of children with a NBSNHI do not have an identifiable genetic cause. To identify genetic contributors, we genotyped 659 GJB2 mutation negative pediatric probands with NBSNHI and assayed for copy number variants (CNVs). After identifying 8 mild-moderate NBSNHI probands with a Chr15q15.3 deletion encompassing the Stereocilin (STRC) gene amongst this cohort, sequencing of STRC was undertaken in these probands as well as 50 probands and 14 siblings with mild-moderate NBSNHI and 40 probands with moderately severe-profound NBSNHI who were GJB2 mutation negative. The existence of a STRC pseudogene that is 99.6% homologous to the STRC coding region has made the sequencing interpretation complicated. We identified 7/50 probands in the mild-moderate cohort to have biallelic alterations in STRC, not including the 8 previously identified deletions. We also identified 2/40 probands to have biallelic alterations in the moderately severe-profound NBSNHI cohort, notably no large deletions in combination with another variant were found in this cohort. The data suggest that STRC may be a common contributor to NBSNHI among GJB2 mutation negative probands, especially in those with mild to moderate hearing impairment.

Mechanistic insights into the events that lead to synergistic induction of interleukin 6 transcription upon activation of the aryl hydrocarbon receptor and inflammatory signaling.

The aryl hydrocarbon receptor (AHR) is the ligand-activated transcription factor responsible for mediating the toxicological effects of dioxin and xenobiotic metabolism. However, recent evidence has implicated the AHR in additional, nonmetabolic physiological processes, including immune regulation. Certain tumor cells are largely nonresponsive to cytokine-mediated induction of the pro-survival cytokine interleukin (IL) 6. We have demonstrated that multiple nonresponsive tumor lines are able to undergo synergistic induction of IL6 following combinatorial treatment with IL1beta and the AHR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin. Such data implicate the AHR in tumor expansion, although the mechanistic basis for the AHR-dependent synergistic induction of IL6 has not been determined. Here, we demonstrate that ligand-activated AHR is involved in priming the IL6 promoter through binding to nonconsensus dioxin response elements located upstream of the IL6 start site. Such binding appears to render the promoter more permissive to IL1beta-induced binding of NF-kappaB components. The nature of the AHR-dependent increases in IL6 promoter transcriptional potential has been shown to involve a reorganization of repressive complexes as exemplified by the presence of HDAC1 and HDAC3. Dismissal of these HDACs correlates with post-translational modifications of promoter-bound NF-kappaB components in a time-dependent manner. Thus the AHR plays a role in derepressing the IL6 promoter, leading to synergistic IL6 expression in the presence of inflammatory signals. These observations may explain the association between enhanced expression of AHR and tumor aggressiveness. It is likely that AHR-mediated priming is not restricted to the IL6 promoter and may contribute to the expression of a variety of genes, which do not have consensus dioxin response elements.

Intermittent Hypoxia Alters the Circadian Expression of Clock Genes in Mouse Brain and Liver.

At least one-third of adults in the United States experience intermittent hypoxia (IH) due to health or living conditions. The majority of these adults suffer with sleep breathing conditions and associated circadian rhythm disorders. The impact of IH on the circadian clock is not well characterized. In the current study, we used an IH mouse model to understand the impact of IH on the circadian gene expression of the canonical clock genes in the central (the brain) and peripheral (the liver) tissues. Gene expression was measured using a Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR). CircaCompare was used to evaluate the differential rhythmicity between normoxia and IH. Our observations suggested that the circadian clock in the liver was less sensitive to IH compared to the circadian clock in the brain.

Short-term exposure to intermittent hypoxia leads to changes in gene expression seen in chronic pulmonary disease.

Obstructive sleep apnea (OSA) results from episodes of airway collapse and intermittent hypoxia (IH) and is associated with a host of health complications. Although the lung is the first organ to sense changes in oxygen levels, little is known about the consequences of IH to the lung hypoxia-inducible factor-responsive pathways. We hypothesized that exposure to IH would lead to cell-specific up- and downregulation of diverse expression pathways. We identified changes in circadian and immune pathways in lungs from mice exposed to IH. Among all cell types, endothelial cells showed the most prominent transcriptional changes. Upregulated genes in myofibroblast cells were enriched for genes associated with pulmonary hypertension and included targets of several drugs currently used to treat chronic pulmonary diseases. A better understanding of the pathophysiologic mechanisms underlying diseases associated with OSA could improve our therapeutic approaches, directing therapies to the most relevant cells and molecular pathways.