Bromodomain Proteins Epigenetically Regulate the Mitotically Associated lncRNA MANCR in Triple Negative Breast Cancer Cells.

Long non-coding RNA (lncRNA)-mediated control of gene expression contributes to regulation of biological processes that include proliferation and phenotype, as well as compromised expression of genes that are functionally linked to cancer initiation and tumor progression. lncRNAs have emerged as novel targets and biomarkers in breast cancer. We have shown that mitotically associated lncRNA MANCR is expressed in triple-negative breast cancer (TNBC) cells and that it serves a critical role in promoting genome stability and survival in aggressive breast cancer cells. Using an siRNA strategy, we selectively depleted BRD2, BRD3, and BRD4, singly and in combination, to establish which bromodomain proteins regulate MANCR expression in TNBC cells. Our findings were confirmed by using in situ hybridization combined with immunofluorescence analysis that revealed BRD4, either alone or with BRD2 and BRD3, can support MANCR regulation of TNBC cells. Here we provide evidence for MANCR-responsive epigenetic control of super enhancers by histone modifications that are required for gene transcription to support cell survival and expression of the epithelial tumor phenotype in triple negative breast cancer cells.

Thyroid hormone dependent transcriptional programming by TRß requires SWI/SNF chromatin remodelers.

Transcriptional regulation in response to thyroid hormone (3,5,3'-triiodo-l-thyronine, T3) is a dynamic and cell-type specific process that maintains cellular homeostasis and identity in all tissues. However, our understanding of the mechanisms of thyroid hormone receptor (TR) actions at the molecular level are actively being refined. We used an integrated genomics approach to profile and characterize the cistrome of TRß, map changes in chromatin accessibility, and capture the transcriptomic changes in response to T3 in normal human thyroid cells. There are significant shifts in TRß genomic occupancy in response to T3, which are associated with differential chromatin accessibility, and differential recruitment of SWI/SNF chromatin remodelers. We further demonstrate selective recruitment of BAF and PBAF SWI/SNF complexes to TRß binding sites, revealing novel differential functions in regulating chromatin accessibility and gene expression. Our findings highlight three distinct modes of TRß interaction with chromatin and coordination of coregulator activity.

peaksat: an R package for ChIP-seq peak saturation analysis.

BACKGROUND: Epigenomic profiling assays such as ChIP-seq have been widely used to map the genome-wide enrichment profiles of chromatin-associated proteins and posttranslational histone modifications. Sequencing depth is a key parameter in experimental design and quality control. However, due to variable sequencing depth requirements across experimental conditions, it can be challenging to determine optimal sequencing depth, particularly for projects involving multiple targets or cell types. RESULTS: We developed the peaksat R package to provide target read depth estimates for epigenomic experiments based on the analysis of peak saturation curves. We applied peaksat to establish the distinctive read depth requirements for ChIP-seq studies of histone modifications in different cell lines. Using peaksat, we were able to estimate the target read depth required per library to obtain high-quality peak calls for downstream analysis. In addition, peaksat was applied to other sequence-enrichment methods including CUT&RUN and ATAC-seq. CONCLUSION: peaksat addresses a need for researchers to make informed decisions about whether their sequencing data has been generated to an adequate depth and subsequently sufficient meaningful peaks, and failing that, how many more reads would be required per library. peaksat is applicable to other sequence-based methods that include calling peaks in their analysis.

Spatiotemporal Epigenetic Control of the Histone Gene Chromatin Landscape during the Cell Cycle.

Higher-order genomic organization supports the activation of histone genes in response to cell cycle regulatory cues that epigenetically mediates stringent control of transcription at the G1/S-phase transition. Histone locus bodies (HLBs) are dynamic, non-membranous, phase-separated nuclear domains where the regulatory machinery for histone gene expression is organized and assembled to support spatiotemporal epigenetic control of histone genes. HLBs provide molecular hubs that support synthesis and processing of DNA replication-dependent histone mRNAs. These regulatory microenvironments support long-range genomic interactions among non-contiguous histone genes within a single topologically associating domain (TAD). HLBs respond to activation of the cyclin E/CDK2/NPAT/HINFP pathway at the G1/S transition. HINFP and its coactivator NPAT form a complex within HLBs that controls histone mRNA transcription to support histone protein synthesis and packaging of newly replicated DNA. Loss of HINFP compromises H4 gene expression and chromatin formation, which may result in DNA damage and impede cell cycle progression. HLBs provide a paradigm for higher-order genomic organization of a subnuclear domain that executes an obligatory cell cycle-controlled function in response to cyclin E/CDK2 signaling. Understanding the coordinately and spatiotemporally organized regulatory programs in focally defined nuclear domains provides insight into molecular infrastructure for responsiveness to cell signaling pathways that mediate biological control of growth, differentiation phenotype, and are compromised in cancer.

Giving Time Purpose: The Synechococcus elongatus Clock in a Broader Network Context.

Early research on the cyanobacterial clock focused on characterizing the genes needed to keep, entrain, and convey time within the cell. As the scope of assays used in molecular genetics has expanded to capture systems-level properties (e.g., RNA-seq, ChIP-seq, metabolomics, high-throughput screening of genetic variants), so has our understanding of how the clock fits within and influences a broader cellular context. Here we review the work that has established a global perspective of the clock, with a focus on (a) an emerging network-centric view of clock architecture, (b) mechanistic insights into how temporal and environmental cues are transmitted and integrated within this network,

Domain-invariant features for mechanism of action prediction in a multi-cell-line drug screen.

MOTIVATION: High-content screening is an important tool in drug discovery and characterization. Often, high-content drug screens are performed on one single-cell line. Yet, a single-cell line cannot be thought of as a perfect disease model. Many diseases feature an important molecular heterogeneity. Consequently, a drug may be effective against one molecular subtype of a disease, but less so against another. To characterize drugs with respect to their effect not only on one cell line but on a panel of cell lines is therefore a promising strategy to streamline the drug discovery process. RESULTS: The contribution of this article is 2-fold. First, we investigate whether we can predict drug mechanism of action (MOA) at the molecular level without optimization of the MOA classes to the screen specificities. To this end, we benchmark a set of algorithms within a conventional pipeline, and evaluate their MOA prediction performance according to a statistically rigorous framework. Second, we extend this conventional pipeline to the simultaneous analysis of multiple cell lines, each manifesting potentially different morphological baselines. For this, we propose multi-task autoencoders, including a domain-adaptive model used to construct domain-invariant feature representations across cell lines. We apply these methods to a pilot screen of two triple negative breast cancer cell lines as models for two different molecular subtypes of the disease. AVAILABILITY AND IMPLEMENTATION: https://github.com/jcboyd/multi-cell-line or https://zenodo.org/record/2677923. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RUNX1 and RUNX2 transcription factors function in opposing roles to regulate breast cancer stem cells.

Breast cancer stem cells (BCSCs) are competent to initiate tumor formation and growth and refractory to conventional therapies. Consequently BCSCs are implicated in tumor recurrence. Many signaling cascades associated with BCSCs are critical for epithelial-to-mesenchymal transition (EMT). We developed a model system to mechanistically examine BCSCs in basal-like breast cancer using MCF10AT1 FACS sorted for CD24 (negative/low in BCSCs) and CD44 (positive/high in BCSCs). Ingenuity Pathway Analysis comparing RNA-seq on the CD24-/low versus CD24+/high MCF10AT1 indicates that the top activated upstream regulators include TWIST1, TGFß1, OCT4, and other factors known to be increased in BCSCs and during EMT. The top inhibited upstream regulators include ESR1, TP63, and FAS. Consistent with our results, many genes previously demonstrated to be regulated by RUNX factors are altered in BCSCs. The RUNX2 interaction network is the top significant pathway altered between CD24-/low and CD24+/high MCF10AT1. RUNX1 is higher in expression at the RNA level than RUNX2. RUNX3 is not expressed. While, human-specific quantitative polymerase chain reaction primers demonstrate that RUNX1 and CDH1 decrease in human MCF10CA1a cells that have grown tumors within the murine mammary fat pad microenvironment, RUNX2 and VIM increase. Treatment with an inhibitor of RUNX binding to CBFß for 5 days followed by a 7-day recovery period results in EMT suggesting that loss of RUNX1, rather than increase in RUNX2, is a driver of EMT in early stage breast cancer. Increased understanding of RUNX regulation on BCSCs and EMT will provide novel insight into therapeutic strategies to prevent recurrence.

RUNX1-dependent mechanisms in biological control and dysregulation in cancer.

The RUNX1 transcription factor has recently been shown to be obligatory for normal development. RUNX1 controls the expression of genes essential for proper development in many cell lineages and tissues including blood, bone, cartilage, hair follicles, and mammary glands. Compromised RUNX1 regulation is associated with many cancers. In this review, we highlight evidence for RUNX1 control in both invertebrate and mammalian development and recent novel findings of perturbed RUNX1 control in breast cancer that has implications for other solid tumors. As RUNX1 is essential for definitive hematopoiesis, RUNX1 mutations in hematopoietic lineage cells have been implicated in the etiology of several leukemias. Studies of solid tumors have revealed a context-dependent function for RUNX1 either as an oncogene or a tumor suppressor. These RUNX1 functions have been reported for breast, prostate, lung, and skin cancers that are related to cancer subtypes and different stages of tumor development. Growing evidence suggests that RUNX1 suppresses aggressiveness in most breast cancer subtypes particularly in the early stage of tumorigenesis. Several studies have identified RUNX1 suppression of the breast cancer epithelial-to-mesenchymal transition. Most recently, RUNX1 repression of cancer stem cells and tumorsphere formation was reported for breast cancer. It is anticipated that these new discoveries of the context-dependent diversity of RUNX1 functions will lead to innovative therapeutic strategies for the intervention of cancer and other abnormalities of normal tissues.

Higher order genomic organization and regulatory compartmentalization for cell cycle control at the G1/S-phase transition.

Fidelity of histone gene regulation, and ultimately of histone protein biosynthesis, is obligatory for packaging of newly replicated DNA into chromatin. Control of histone gene expression within the 3-dimensional context of nuclear organization is reflected by two well documented observations. DNA replication-dependent histone mRNAs are synthesized at specialized subnuclear domains designated histone locus bodies (HLBs), in response to activation of the growth factor dependent Cyclin E/CDK2/HINFP/NPAT pathway at the G1/S transition in mammalian cells. Complete loss of the histone gene regulatory factors HINFP or NPAT disrupts HLB integrity that is necessary for coordinate control of DNA replication and histone gene transcription. Here we review the molecular histone-related requirements for G1/S-phase progression during the cell cycle. Recently developed experimental strategies, now enable us to explore mechanisms involved in dynamic control of histone gene expression in the context of the temporal (cell cycle) and spatial (HLBs) remodeling of the histone gene loci.

Intranuclear and higher-order chromatin organization of the major histone gene cluster in breast cancer.

Alterations in nuclear morphology are common in cancer progression. However, the degree to which gross morphological abnormalities translate into compromised higher-order chromatin organization is poorly understood. To explore the functional links between gene expression and chromatin structure in breast cancer, we performed RNA-seq gene expression analysis on the basal breast cancer progression model based on human MCF10A cells. Positional gene enrichment identified the major histone gene cluster at chromosome 6p22 as one of the most significantly upregulated (and not amplified) clusters of genes from the normal-like MCF10A to premalignant MCF10AT1 and metastatic MCF10CA1a cells. This cluster is subdivided into three sub-clusters of histone genes that are organized into hierarchical topologically associating domains (TADs). Interestingly, the sub-clusters of histone genes are located at TAD boundaries and interact more frequently with each other than the regions in-between them, suggesting that the histone sub-clusters form an active chromatin hub. The anchor sites of loops within this hub are occupied by CTCF, a known chromatin organizer. These histone genes are transcribed and processed at a specific sub-nuclear microenvironment termed the major histone locus body (HLB). While the overall chromatin structure of the major HLB is maintained across breast cancer progression, we detected alterations in its structure that may relate to gene expression. Importantly, breast tumor specimens also exhibit a coordinate pattern of upregulation across the major histone gene cluster. Our results provide a novel insight into the connection between the higher-order chromatin organization of the major HLB and its regulation during breast cancer progression.

Regulation of osteogenesis by long noncoding RNAs: An epigenetic mechanism contributing to bone formation.

Long noncoding RNAs (lncRNAs) have recently emerged as novel regulators of lineage commitment, differentiation, development, viability, and disease progression. Few studies have examined their role in osteogenesis; however, given their critical and wide-ranging roles in other tissues, lncRNAs are most likely vital regulators of osteogenesis. In this study, we extensively characterized lncRNA expression in mesenchymal cells during commitment and differentiation to the osteoblast lineage using a whole transcriptome sequencing approach (RNA-Seq). Using mouse primary mesenchymal stromal cells (mMSC), we identified 1438 annotated lncRNAs expressed during MSC differentiation, 462 of which are differentially expressed. We performed guilt-by-association analysis using lncRNA and mRNA expression profiles to identify lncRNAs influencing MSC commitment and differentiation. These findings open novel dimensions for exploring lncRNAs in regulating normal bone formation and in skeletal disorders.

Cross-talk and regulatory interactions between the essential response regulator RpaB and cyanobacterial circadian clock output.

The response regulator RpaB (regulator of phycobilisome associated B), part of an essential two-component system conserved in cyanobacteria that responds to multiple environmental signals, has recently been implicated in the control of cell dimensions and of circadian rhythms of gene expression in the model cyanobacterium Synechococcus elongatus PCC 7942. However, little is known of the molecular mechanisms that underlie RpaB functions. In this study we show that the regulation of phenotypes by RpaB is intimately connected with the activity of RpaA (regulator of phycobilisome associated A), the master regulator of circadian transcription patterns. RpaB affects RpaA activity both through control of gene expression, a function requiring an intact effector domain, and via altering RpaA phosphorylation, a function mediated through the N-terminal receiver domain of RpaB. Thus, both phosphorylation cross-talk and coregulation of target genes play a role in the genetic interactions between the RpaA and RpaB pathways. In addition, RpaB∼P levels appear critical for survival under light:dark cycles, conditions in which RpaB phosphorylation is environmentally driven independent of the circadian clock. We propose that the complex regulatory interactions between the essential and environmentally sensitive NblS-RpaB system and the SasA-RpaA clock output system integrate relevant extra- and intracellular signals to the circadian clock.

Enhanced biomass and oxidative stress tolerance of Synechococcus elongatus PCC 7942 overexpressing the DHAR gene from Brassica juncea.

OBJECTIVES: To improve the oxidative stress tolerance, biomass yield, and ascorbate/dehydroascorbate (AsA/DHA) ratio of Synechococcus elongatus PCC 7942 in the presence of H2O2, by heterologous expression of the dehydroascorbate reductase (DHAR) gene from Brassica juncea (BrDHAR). RESULTS: Under H2O2 stress, overexpression of BrDHAR in the transgenic strain (BrD) of S. elongatus greatly increased the AsA/DHA ratio. As part of the AsA recycling system, the oxidative stress response induced by reactive oxygen species was enhanced, and intracellular H2O2 level decreased. In addition, under H2O2 stress conditions, the BrD strain displayed increased growth rate and biomass, as well as higher chlorophyll content and deeper pigmentation than did wild-type and control strains. CONCLUSION: By maintaining the AsA pool and redox homeostasis, the heterologous expression of BrDHAR increased S. elongatus tolerance to H2O2 stress, improving the biomass yield under these conditions. The results suggest that the BrD strain of S. elongatus, with its ability to attenuate the deleterious effects of ROS caused by environmental stressors, could be a promising platform for the generation of biofuels and other valuable bioproducts.

Single mutations in sasA enable a simpler ΔcikA gene network architecture with equivalent circadian properties.

The circadian input kinase of the cyanobacterium Synechococcus elongatus PCC 7942 (CikA) is important both for synchronizing circadian rhythms with external environmental cycles and for transferring temporal information between the oscillator and the global transcriptional regulator RpaA (regulator of phycobilisome-associated A). KOs of cikA result in one of the most severely altered but still rhythmic circadian phenotypes observed. We chemically mutagenized a cikA-null S. elongatus strain and screened for second-site suppressor mutations that could restore normal circadian rhythms. We identified two independent mutations in the Synechococcus adaptive sensor A (sasA) gene that produce nearly WT rhythms of gene expression, likely because they compensate for the loss of CikA on the temporal phosphorylation of RpaA. Additionally, these mutations restore the ability to reset the clock after a short dark pulse through an output-independent pathway, suggesting that SasA can influence entrainment through direct interactions with KaiC, a property previously unattributed to it. These experiments question the evolutionary advantage of integrating CikA into the cyanobacterial clock, challenge the conventional construct of separable input and output pathways, and show how easily the cell can adapt to restore phenotype in a severely compromised genetic network.

A Combined Computational and Genetic Approach Uncovers Network Interactions of the Cyanobacterial Circadian Clock.

UNLABELLED: Two-component systems (TCS) that employ histidine kinases (HK) and response regulators (RR) are critical mediators of cellular signaling in bacteria. In the model cyanobacterium Synechococcus elongatus PCC 7942, TCSs control global rhythms of transcription that reflect an integration of time information from the circadian clock with a variety of cellular and environmental inputs. The HK CikA and the SasA/RpaA TCS transduce time information from the circadian oscillator to modulate downstream cellular processes. Despite immense progress in understanding of the circadian clock itself, many of the connections between the clock and other cellular signaling systems have remained enigmatic. To narrow the search for additional TCS components that connect to the clock, we utilized direct-coupling analysis (DCA), a statistical analysis of covariant residues among related amino acid sequences, to infer coevolution of new and known clock TCS components. DCA revealed a high degree of interaction specificity between SasA and CikA with RpaA, as expected, but also with the phosphate-responsive response regulator SphR. Coevolutionary analysis also predicted strong specificity between RpaA and a previously undescribed kinase, HK0480 (herein CikB). A knockout of the gene for CikB (cikB) in a sasA cikA null background eliminated the RpaA phosphorylation and RpaA-controlled transcription that is otherwise present in that background and suppressed cell elongation, supporting the notion that CikB is an interactor with RpaA and the clock network. This study demonstrates the power of DCA to identify subnetworks and key interactions in signaling pathways and of combinatorial mutagenesis to explore the phenotypic consequences. Such a combined strategy is broadly applicable to other prokaryotic systems. IMPORTANCE: Signaling networks are complex and extensive, comprising multiple integrated pathways that respond to cellular and environmental cues. A TCS interaction model, based on DCA, independently confirmed known interactions and revealed a core set of subnetworks within the larger HK-RR set. We validated high-scoring candidate proteins via combinatorial genetics, demonstrating that DCA can be utilized to reduce the search space of complex protein networks and to infer undiscovered specific interactions for signaling proteins in vivo Significantly, new interactions that link circadian response to cell division and fitness in a light/dark cycle were uncovered. The combined analysis also uncovered a more basic core clock, illustrating the synergy and applicability of a combined computational and genetic approach for investigating prokaryotic signaling networks.

Oncofetal Epigenetic Bivalency in Breast Cancer Cells: H3K4 and H3K27 Tri-Methylation as a Biomarker for Phenotypic Plasticity.

Alterations in the epigenetic landscape are fundamental drivers of aberrant gene expression that contribute to cancer progression and pathology. Understanding specific modes of epigenetic regulation can be used to identify novel biomarkers or targets for therapeutic intervention to clinically treat solid tumors and leukemias. The bivalent marking of gene promoters by H3K4me3 and H3K27me3 is a primary mechanism to poise genes for expression in pluripotent embryonic stem cells (ESC). In this study we interrogated three well-established mammary cell lines to model epigenetic programming observed among breast cancer subtypes. Evidence is provided for a distinct bivalent signature, activating and repressive histone marks co-residing at the same gene promoter, in the MCF7 (ESR/PGR+) luminal breast cancer cell line. We identified a subset of genes, enriched for developmental pathways that regulate cellular phenotype and signaling, and partially recapitulate the bivalent character observed in ESC. We validated the biological relevance of this "oncofetal epigenetic" signature using data from ESR/PGR+ tumor samples from breast cancer patients. This signature of oncofetal epigenetic control is an informative biomarker and may provide novel therapeutic targets, selective for both recurring and treatment-resistant cancers. J. Cell. Physiol. 231: 2474-2481, 2016. © 2016 Wiley Periodicals, Inc.

Expression of Ribosomal RNA and Protein Genes in Human Embryonic Stem Cells Is Associated With the Activating H3K4me3 Histone Mark.

Embryonic stem cells (ESCs) exhibit unrestricted and indefinite, but stringently controlled, proliferation, and can differentiate into any lineage in the body. In the current study, we test the hypothesis that expression of ribosomal RNA (rRNA) and ribosomal protein genes (RPGs) contribute to the ability of hESCs to proliferate indefinitely. Consistent with the accelerated growth rate of hESCs, we find that hESC lines H1 and H9 both exhibit significantly higher levels of rRNA when compared to a panel of normal and cancer human cell lines. Although many RPGs are expressed at levels that comparable to other human cell lines, a few RPGs also exhibit higher expression levels. In situ nuclear run-on assays reveal that both nucleoli in hESCs actively transcribe nascent rRNA. Employing genome-wide chromatin immunoprecipitation-deep sequencing and bioinformatics approaches, we discovered that, RPGs are dominantly marked by the activating H3K4me3 histone mark in the G1, M, and G2 phases of the cell cycle. Interestingly, the rDNA repeats are marked by the activating H3K4me3 only in the M phase, and repressive H3K27me3 histone mark in all three cell cycle phases. Bioinformatics analyses also reveal that Myc, a known regulator of cell growth and proliferation, occupies both the rRNA genes and RPGs. Functionally, down-regulation of Myc expression by siRNA results in a concomitant decrease in rRNA levels. Together, our results show that expression of rRNA, which is regulated by the Myc pluripotency transcription factor, and of RPGs in hESCs is associated with the activating H3K4me3 modification. J. Cell. Physiol. 231: 2007-2013, 2016. © 2016 Wiley Periodicals, Inc.

MicroRNA-378-mediated suppression of Runx1 alleviates the aggressive phenotype of triple-negative MDA-MB-231 human breast cancer cells.

The Runx1 transcription factor, known for its essential role in normal hematopoiesis, was reported in limited studies to be mutated or associated with human breast tumor tissues. Runx1 increases concomitantly with disease progression in the MMTV-PyMT transgenic mouse model of breast cancer. Compelling questions relate to mechanisms that regulate Runx1 expression in breast cancer. Here, we tested the hypothesis that dysregulation of Runx1-targeting microRNAs (miRNAs) allows for pathologic increase of Runx1 during breast cancer progression. Microarray profiling of the MMTV-PyMT model revealed significant downregulation of numerous miRNAs predicted to target Runx1. One of these, miR-378, was inversely correlated with Runx1 expression during breast cancer progression in mice and in human breast cancer cell lines MCF7 and triple-negative MDA-MB-231 that represent early- and late-stage diseases, respectively. MiR-378 is nearly absent in MDA-MB-231 cells. Luciferase reporter assays revealed that miR-378 binds the Runx1 3' untranslated region (3'UTR) and inhibits Runx1 expression. Functionally, we demonstrated that ectopic expression of miR-378 in MDA-MB-231 cells inhibited Runx1 and suppressed migration and invasion, while inhibition of miR-378 in MCF7 cells increased Runx1 levels and cell migration. Depletion of Runx1 in late-stage breast cancer cells resulted in increased expression of both the miR-378 host gene PPARGC1B and pre-miR-378, suggesting a feedback loop. Taken together, our study identifies a novel and clinically relevant mechanism for regulation of Runx1 in breast cancer that is mediated by a PPARGC1B-miR-378-Runx1 regulatory pathway. Our results highlight the translational potential of miRNA replacement therapy for inhibiting Runx1 in breast cancer.

An allele of the crm gene blocks cyanobacterial circadian rhythms.

The SasA-RpaA two-component system constitutes a key output pathway of the cyanobacterial Kai circadian oscillator. To date, rhythm of phycobilisome associated (rpaA) is the only gene other than kaiA, kaiB, and kaiC, which encode the oscillator itself, whose mutation causes completely arrhythmic gene expression. Here we report a unique transposon insertion allele in a small ORF located immediately upstream of rpaA in Synechococcus elongatus PCC 7942 termed crm (for circadian rhythmicity modulator), which results in arrhythmic promoter activity but does not affect steady-state levels of RpaA. The crm ORF complements the defect when expressed in trans, but only if it can be translated, suggesting that crm encodes a small protein. The crm1 insertion allele phenotypes are distinct from those of an rpaA null; crm1 mutants are able to grow in a light:dark cycle and have no detectable oscillations of KaiC phosphorylation, whereas low-amplitude KaiC phosphorylation rhythms persist in the absence of RpaA. Levels of phosphorylated RpaA in vivo measured over time are significantly altered compared with WT in the crm1 mutant as well as in the absence of KaiC. Taken together, these results are consistent with the hypothesis that the Crm polypeptide modulates a circadian-specific activity of RpaA.

Active output state of the Synechococcus Kai circadian oscillator.

The mechanisms by which cellular oscillators keep time and transmit temporal information are poorly understood. In cyanobacteria, the timekeeping aspect of the circadian oscillator, composed of the KaiA, KaiB, and KaiC proteins, involves a cyclic progression of phosphorylation states at Ser431 and Thr432 of KaiC. Elucidating the mechanism that uses this temporal information to modulate gene expression is complicated by unknowns regarding the number, structure, and regulatory effects of output components. To identify oscillator signaling states without a complete description of the output machinery, we defined a simple metric, Kai-complex output activity (KOA), that represents the difference in expression of reporter genes between strains that carry specific variants of KaiC and baseline strains that lack KaiC. In the absence of the oscillator, expression of the class 1 paradigm promoter P(kaiBC) was locked at its usual peak level; conversely, that of the class 2 paradigm promoter P(purF) was locked at its trough level. However, for both classes of promoters, peak KOA in wild-type strains coincided late in the circadian cycle near subjective dawn, when KaiC-pST becomes most prevalent (Ser431 is phosphorylated and Thr432 is not). Analogously, peak KOA was detected specifically for the phosphomimetic of KaiC-pST (KaiC-ET). Notably, peak KOA required KaiB, indicating that a KaiBC complex is involved in the output activity. We also found evidence that phosphorylated RpaA (regulator of phycobilisome associated) represses an RpaA-independent output of KOA. A simple mathematical expression successfully simulated two key features of the oscillator-the time of peak KOA and the peak-to-trough amplitude changes.