Systems biology of the clock in Neurospora crassa

The clock in growing hyphae and their synchronization in Neurospora crassa

Utilizing a microfluidic chip with serpentine channels, we inoculated the chip with an agar plug with Neurospora crassa mycelium and successfully captured individual hyphae in channels. For the first time, we report the presence of an autonomous clock in hyphae. Fluorescence of a mCherry reporter gene driven by a clock-controlled gene-2 promoter (ccg-2p) was measured simultaneously along hyphae every half an hour for at least 6 days. We entrained single hyphae to light over a wide range of day lengths, including 6,12, 24, and 36 h days. Hyphae tracked in individual serpentine channels were highly synchronized (K = 0.60-0.78). Furthermore, hyphae also displayed temperature compensation properties, where the oscillation period was stable over a physiological range of temperatures from 24 °C to 30 °C (Q10 = 1.00-1.10). A Clock Tube Model developed could mimic hyphal growth observed in the serpentine chip and provides a mechanism for the stable banding patterns seen in race tubes at the macroscopic scale and synchronization through molecules riding the growth wave in the device.

Designs for the simultaneous inference of concentration-response curves

BACKGROUND

An important problem in toxicology in the context of gene expression data is the simultaneous inference of a large number of concentration-response relationships. The quality of the inference substantially depends on the choice of design of the experiments, in particular, on the set of different concentrations, at which observations are taken for the different genes under consideration. As this set has to be the same for all genes, the efficient planning of such experiments is very challenging. We address this problem by determining efficient designs for the simultaneous inference of a large number of concentration-response models. For that purpose, we both construct a D-optimality criterion for simultaneous inference and a K-means procedure which clusters the support points of the locally D-optimal designs of the individual models.

RESULTS

We show that a planning of experiments that addresses the simultaneous inference of a large number of concentration-response relationships yields a substantially more accurate statistical analysis. In particular, we compare the performance of the constructed designs to the ones of other commonly used designs in terms of D-efficiencies and in terms of the quality of the resulting model fits using a real data example dealing with valproic acid. For the quality comparison we perform an extensive simulation study.

CONCLUSIONS

The design maximizing the D-optimality criterion for simultaneous inference improves the inference of the different concentration-response relationships substantially. The design based on the K-means procedure also performs well, whereas a log-equidistant design, which was also included in the analysis, performs poorly in terms of the quality of the simultaneous inference. Based on our findings, the D-optimal design for simultaneous inference should be used for upcoming analyses dealing with high-dimensional gene expression data.

Light Stress in Yeasts: Signaling and Responses in Creatures of the Night

Living organisms on the surface biosphere are periodically yet consistently exposed to light. The adaptive or protective evolution caused by this source of energy has led to the biological systems present in a large variety of organisms, including fungi. Among fungi, yeasts have developed essential protective responses against the deleterious effects of light. Stress generated by light exposure is propagated through the synthesis of hydrogen peroxide and mediated by regulatory factors that are also involved in the response to other stressors. These have included Msn2/4, Crz1, Yap1, and Mga2, thus suggesting that light stress is a common factor in the yeast environmental response.

The macroscopic limit to synchronization of cellular clocks in single cells of Neurospora crassa

We determined the macroscopic limit for phase synchronization of cellular clocks in an artificial tissue created by a “big chamber” microfluidic device to be about 150,000 cells or less. The dimensions of the microfluidic chamber allowed us to calculate an upper limit on the radius of a hypothesized quorum sensing signal molecule of 13.05 nm using a diffusion approximation for signal travel within the device. The use of a second microwell microfluidic device allowed the refinement of the macroscopic limit to a cell density of 2166 cells per fixed area of the device for phase synchronization. The measurement of averages over single cell trajectories in the microwell device supported a deterministic quorum sensing model identified by ensemble methods for clock phase synchronization. A strong inference framework was used to test the communication mechanism in phase synchronization of quorum sensing versus cell-to-cell contact, suggesting support for quorum sensing. Further evidence came from showing phase synchronization was density-dependent.

Characterizing the gene-environment interaction underlying natural morphological variation in Neurospora crassa conidiophores using high-throughput phenomics and transcriptomics

Neurospora crassa propagates through dissemination of conidia, which develop through specialized structures called conidiophores. Recent work has identified striking variation in conidiophore morphology, using a wild population collection from Louisiana, United States of America to classify 3 distinct phenotypes: Wild-Type, Wrap, and Bulky. Little is known about the impact of these phenotypes on sporulation or germination later in the N. crassa life cycle, or about the genetic variation that underlies them. In this study, we show that conidiophore morphology likely affects colonization capacity of wild N. crassa isolates through both sporulation distance and germination on different carbon sources. We generated and crossed homokaryotic strains belonging to each phenotypic group to more robustly fit a model for and estimate heritability of the complex trait, conidiophore architecture. Our fitted model suggests at least 3 genes and 2 epistatic interactions contribute to conidiophore phenotype, which has an estimated heritability of 0.47. To uncover genes contributing to these phenotypes, we performed RNA-sequencing on mycelia and conidiophores of strains representing each of the 3 phenotypes. Our results show that the Bulky strain had a distinct transcriptional profile from that of Wild-Type and Wrap, exhibiting differential expression patterns in clock-controlled genes (ccgs), the conidiation-specific gene con-6, and genes implicated in metabolism and communication. Combined, these results present novel ecological impacts of and differential gene expression underlying natural conidiophore morphological variation, a complex trait that has not yet been thoroughly explored.

The Resonance and Adaptation of Neurospora crassa Circadian and Conidiation Rhyth ms to Short Light-Dark Cycles

Circadian clocks control the physiological and behavioral rhythms to adapt to the environment with a period of ~24 h. However, the influences and mechanisms of the extreme light/dark cycles on the circadian clock remain unclear. We showed that, in Neurospora crassa, both the growth and the microconidia production contribute to adaptation in LD12:12 (12 h light/12 h dark, periodically). Mathematical modeling and experiments demonstrate that in short LD cycles, the expression of the core clock protein FREQUENCY was entrained to the LD cycles when LD > 3:3 while it free ran when T ≤ LD3:3. The conidial rhythmicity can resonate with a series of different LD conditions. Moreover, we demonstrate that the existence of unknown blue light photoreceptor(s) and the circadian clock might promote the conidiation rhythms that resonate with the environment. The ubiquitin E3 ligase FWD-1 and the previously described CRY-dependent oscillator system were implicated in regulating conidiation under short LD conditions. These findings shed new light on the resonance of Neurospora circadian clock and conidiation rhythms to short LD cycles, which may benefit the understandings of both the basic regulatory aspects of circadian clock and the adaptation of physiological rhythms to the extreme conditions.

A global search for novel transcription factors impacting the Neurospora crassa circadian clock

Eukaryotic circadian oscillators share a common circuit architecture, a negative feedback loop in which a positive element activates the transcription of a negative one that then represses the action of the former, inhibiting its own expression. While studies in mammals and insects have revealed additional transcriptional inputs modulating the expression of core clock components, this has been less characterized in the model Neurospora crassa, where the participation of other transcriptional components impacting circadian clock dynamics remains rather unexplored. Thus, we sought to identify additional transcriptional regulators modulating the N. crassa clock, following a reverse genetic screen based on luminescent circadian reporters and a collection of transcription factors (TFs) knockouts, successfully covering close to 60% of them. Besides the canonical core clock components WC-1 and -2, none of the tested transcriptional regulators proved to be essential for rhythmicity. Nevertheless, we identified a set of 23 TFs that when absent lead to discrete, but significant, changes in circadian period. While the current level of analysis does not provide mechanistic information about how these new players modulate circadian parameters, the results of this screen reveal that an important number of light and clock-regulated TFs, involved in a plethora of processes, are capable of modulating the clockworks. This partial reverse genetic clock screen also exemplifies how the N. crassa knockout collection continues to serve as an expedite platform to address broad biological questions.

Circadian Interactomics: How Research Into Protein-Protein Interactions Beyond the Core Clock Has Influenced the Model of Circadian Timekeeping

The circadian clock is the broadly conserved, protein-based, timekeeping mechanism that synchronizes biology to the Earth's 24-h light-dark cycle. Studies of the mechanisms of circadian timekeeping have placed great focus on the role that individual protein-protein interactions play in the creation of the timekeeping loop. However, research has shown that clock proteins most commonly act as part of large macromolecular protein complexes to facilitate circadian control over physiology. The formation of these complexes has led to the large-scale study of the proteins that comprise these complexes, termed here "circadian interactomics." Circadian interactomic studies of the macromolecular protein complexes that comprise the circadian clock have uncovered many basic principles of circadian timekeeping as well as mechanisms of circadian control over cellular physiology. In this review, we examine the wealth of knowledge accumulated using circadian interactomics approaches to investigate the macromolecular complexes of the core circadian clock, including insights into the core mechanisms that impart circadian timing and the clock's regulation of many physiological processes. We examine data acquired from the investigation of the macromolecular complexes centered on both the activating and repressing arm of the circadian clock and from many circadian model organisms.

Identifying a stochastic clock network with light entrainment for single cells of Neurospora crassa

Stochastic networks for the clock were identified by ensemble methods using genetic algorithms that captured the amplitude and period variation in single cell oscillators of Neurospora crassa. The genetic algorithms were at least an order of magnitude faster than ensemble methods using parallel tempering and appeared to provide a globally optimum solution from a random start in the initial guess of model parameters (i.e., rate constants and initial counts of molecules in a cell). The resulting goodness of fit [Formula: see text] was roughly halved versus solutions produced by ensemble methods using parallel tempering, and the resulting [Formula: see text] per data point was only [Formula: see text] = 2,708.05/953 = 2.84. The fitted model ensemble was robust to variation in proxies for "cell size". The fitted neutral models without cellular communication between single cells isolated by microfluidics provided evidence for only one Stochastic Resonance at one common level of stochastic intracellular noise across days from 6 to 36 h of light/dark (L/D) or in a D/D experiment. When the light-driven phase synchronization was strong as measured by the Kuramoto (K), there was degradation in the single cell oscillations away from the stochastic resonance. The rate constants for the stochastic clock network are consistent with those determined on a macroscopic scale of 10⁷ cells.

Circadian misalignment on submarines and other non-24-h environments - from research to application

Circadian clocks have important physiological and behavioral functions in humans and other organisms, which enable organisms to anticipate and respond to periodic environmental changes. Disturbances in circadian rhythms impair sleep, metabolism, and behavior. People with jet lag, night workers and shift workers are vulnerable to circadian misalignment. In addition, non-24-h cycles influence circadian rhythms and cause misalignment and disorders in different species, since these periods are beyond the entrainment ranges. In certain special conditions, e.g., on submarines and commercial ships, non-24-h watch schedules are often employed, which have also been demonstrated to be deleterious to circadian rhythms. Personnel working under such conditions suffer from circadian misalignment with their on-watch hours, leading to increased health risks and decreased cognitive performance. In this review, we summarize the research progress and knowledge concerning circadian rhythms on submarines and other environments in which non-24-h watch schedules are employed.

What is Phase in Cellular Clocks?

Four inter-related measures of phase are described to study the phase synchronization of cellular oscillators, and computation of these measures is described and illustrated on single cell fluorescence data from the model filamentous fungus, Neurospora crassa. One of these four measures is the phase shift ϕ in a sinusoid of the form x(t) = A(cos(ωt + ϕ), where t is time. The other measures arise by creating a replica of the periodic process x(t) called the Hilbert transform x̃(t), which is 90 degrees out of phase with the original process x(t). The second phase measure is the phase angle F^H(t) between the replica x̃(t) and x(t), taking values between -π and π. At extreme values the Hilbert Phase is discontinuous, and a continuous form F^C(t) of the Hilbert Phase is used, measuring time on the nonnegative real axis (t). The continuous Hilbert Phase F^C(t) is used to define the phase M^C(t₁,t₀) for an experiment beginning at time t₀ and ending at time t₁. In that phase differences at time t₀ are often of ancillary interest, the Hilbert Phase F^C(t₀) is subtracted from F^C(t₁). This difference is divided by 2π to obtain the phase M^C(t₁,t₀) in cycles. Both the Hilbert Phase F^C(t) and the phase M^C(t₁,t₀) are functions of time and useful in studying when oscillators phase-synchronize in time in signal processing and circadian rhythms in particular. The phase of cellular clocks is fundamentally different from circadian clocks at the macroscopic scale because there is an hourly cycle superimposed on the circadian cycle.

Ensemble methods for stochastic networks with special reference to the biological clock of Neurospora crassa

A major challenge in systems biology is to infer the parameters of regulatory networks that operate in a noisy environment, such as in a single cell. In a stochastic regime it is hard to distinguish noise from the real signal and to infer the noise contribution to the dynamical behavior. When the genetic network displays oscillatory dynamics, it is even harder to infer the parameters that produce the oscillations. To address this issue we introduce a new estimation method built on a combination of stochastic simulations, mass action kinetics and ensemble network simulations in which we match the average periodogram and phase of the model to that of the data. The method is relatively fast (compared to Metropolis-Hastings Monte Carlo Methods), easy to parallelize, applicable to large oscillatory networks and large (~2000 cells) single cell expression data sets, and it quantifies the noise impact on the observed dynamics. Standard errors of estimated rate coefficients are typically two orders of magnitude smaller than the mean from single cell experiments with on the order of ~1000 cells. We also provide a method to assess the goodness of fit of the stochastic network using the Hilbert phase of single cells. An analysis of phase departures from the null model with no communication between cells is consistent with a hypothesis of Stochastic Resonance describing single cell oscillators. Stochastic Resonance provides a physical mechanism whereby intracellular noise plays a positive role in establishing oscillatory behavior, but may require model parameters, such as rate coefficients, that differ substantially from those extracted at the macroscopic level from measurements on populations of millions of communicating, synchronized cells.

Maximally informative next experiments for nonlinear models

Mathematical modeling is a powerful tool in systems biology; we focus here on improving the reliability of model predictions by reducing the uncertainty in model dynamics through experimental design. Model-based experimental design is a process by which experiments can be systematically chosen to reduce dynamic uncertainty in a given model. We discuss the Maximally Informative Next Experiment (MINE) method for group-wise selection of points in an experimental design and present a convergence result for MINE with nonlinear models. As an application, we illustrate the method on polynomial regression and an ODE model for immune system dynamics. The MINE criterion sequentially determines experiments that can be conducted to best refine model dynamics.

Diurnal regulation of RNA polymerase III transcription is under the control of both the feeding-fasting response and the circadian clock

RNA polymerase III (Pol III) synthesizes short noncoding RNAs, many of which are essential for translation. Accordingly, Pol III activity is tightly regulated with cell growth and proliferation by factors such as MYC, RB1, TRP53, and MAF1. MAF1 is a repressor of Pol III transcription whose activity is controlled by phosphorylation; in particular, it is inactivated through phosphorylation by the TORC1 kinase complex, a sensor of nutrient availability. Pol III regulation is thus sensitive to environmental cues, yet a diurnal profile of Pol III transcription activity is so far lacking. Here, we first use gene expression arrays to measure mRNA accumulation during the diurnal cycle in the livers of (1) wild-type mice, (2) arrhythmic Arntl knockout mice, (3) mice fed at regular intervals during both night and day, and (4) mice lacking the Maf1 gene, and so provide a comprehensive view of the changes in cyclic mRNA accumulation occurring in these different systems. We then show that Pol III occupancy of its target genes rises before the onset of the night, stays high during the night, when mice normally ingest food and when translation is known to be increased, and decreases in daytime. Whereas higher Pol III occupancy during the night reflects a MAF1-dependent response to feeding, the rise of Pol III occupancy before the onset of the night reflects a circadian clock-dependent response. Thus, Pol III transcription during the diurnal cycle is regulated both in response to nutrients and by the circadian clock, which allows anticipatory Pol III transcription.

Predicting network modules of cell cycle regulators using relative protein abundance statistics

BACKGROUND

Parameter estimation in systems biology is typically done by enforcing experimental observations through an objective function as the parameter space of a model is explored by numerical simulations. Past studies have shown that one usually finds a set of "feasible" parameter vectors that fit the available experimental data equally well, and that these alternative vectors can make different predictions under novel experimental conditions. In this study, we characterize the feasible region of a complex model of the budding yeast cell cycle under a large set of discrete experimental constraints in order to test whether the statistical features of relative protein abundance predictions are influenced by the topology of the cell cycle regulatory network.

RESULTS

Using differential evolution, we generate an ensemble of feasible parameter vectors that reproduce the phenotypes (viable or inviable) of wild-type yeast cells and 110 mutant strains. We use this ensemble to predict the phenotypes of 129 mutant strains for which experimental data is not available. We identify 86 novel mutants that are predicted to be viable and then rank the cell cycle proteins in terms of their contributions to cumulative variability of relative protein abundance predictions. Proteins involved in "regulation of cell size" and "regulation of G1/S transition" contribute most to predictive variability, whereas proteins involved in "positive regulation of transcription involved in exit from mitosis," "mitotic spindle assembly checkpoint" and "negative regulation of cyclin-dependent protein kinase by cyclin degradation" contribute the least. These results suggest that the statistics of these predictions may be generating patterns specific to individual network modules (START, S/G2/M, and EXIT). To test this hypothesis, we develop random forest models for predicting the network modules of cell cycle regulators using relative abundance statistics as model inputs. Predictive performance is assessed by the areas under receiver operating characteristics curves (AUC). Our models generate an AUC range of 0.83-0.87 as opposed to randomized models with AUC values around 0.50.

CONCLUSIONS

By using differential evolution and random forest modeling, we show that the model prediction statistics generate distinct network module-specific patterns within the cell cycle network.

The Neurospora Transcription Factor ADV-1 Transduces Light Signals and Temporal Information to Control Rhythmic Expression of Genes Involved in Cell Fusion

Light and the circadian clock have a profound effect on the biology of organisms through the regulation of large sets of genes. Toward understanding how light and the circadian clock regulate gene expression, we used genome-wide approaches to identify the direct and indirect targets of the light-responsive and clock-controlled transcription factor ADV-1 in Neurospora crassa A large proportion of ADV-1 targets were found to be light- and/or clock-controlled, and enriched for genes involved in development, metabolism, cell growth, and cell fusion. We show that ADV-1 is necessary for transducing light and/or temporal information to its immediate downstream targets, including controlling rhythms in genes critical to somatic cell fusion. However, while ADV-1 targets are altered in predictable ways in Δadv-1 cells in response to light, this is not always the case for rhythmic target gene expression. These data suggest that a complex regulatory network downstream of ADV-1 functions to generate distinct temporal dynamics of target gene expression relative to the central clock mechanism.

lac-1 and lag-1 with ras-1 affect aging and the biological clock in Neurospora crassa

Using an automated cell counting technique developed previously (Case et al., Ecology and Evolution 2014; 4: 3494), we explore the lifespan effects of lac‐1, a ceramide synthase gene paralogous to lag‐1 in Neurospora crassa in conjunction with the band bd (ras‐1) gene. We find that the replicative lifespan of a lac‐1^KObd double mutants is short, about one race tube cycle, and this double mutant lacks a strong ~21‐hr clock cycle as shown by race tube and fluorometer analysis of fluorescent strains including lac‐1^KO. This short replicative lifespan phenotype is contrasted with a very long estimated chronological lifespan for lac‐1^KObd double mutants from 247 to 462 days based on our regression analyses on log viability, and for the single mutant lac‐1^KO, 161 days. Both of these estimated lifespans are much higher than that of previously studied WT and bd single mutant strains. In a lac‐1 rescue and induction experiment, the expression of lac‐1⁺ as driven by a quinic acid‐dependent promoter actually decreases the median chronological lifespan of cells down to only 7 days, much lower than the 34‐day median lifespan found in control bd conidia also grown on quinic acid media, which we interpret as an effect of balancing selection acting on ceramide levels based on previous findings from the literature. Prior work has shown phytoceramides can act as a signal for apoptosis in stressed N. crassa cells. To test this hypothesis of balancing selection on phytoceramide levels, we examine the viability of WT, lag‐1^KObd, and lac‐1^KObd strains following the dual stresses of heat and glycolysis inhibition, along with phytoceramide treatments of different dosages. We find that the phytoceramide dosage–response curve is altered in the lag‐1^KObd mutant, but not in the lac‐1^KObd mutant. We conclude that phytoceramide production is responsible for the previously reported longevity effects in the lag‐1^KObd mutant, but a different ceramide may be responsible for the longevity effect observed in the lac‐1^KObd mutant.

Synchronizing stochastic circadian oscillators in single cells of Neurospora crassa

The synchronization of stochastic coupled oscillators is a central problem in physics and an emerging problem in biology, particularly in the context of circadian rhythms. Most measurements on the biological clock are made at the macroscopic level of millions of cells. Here measurements are made on the oscillators in single cells of the model fungal system, Neurospora crassa, with droplet microfluidics and the use of a fluorescent recorder hooked up to a promoter on a clock controlled gene-2 (ccg-2). The oscillators of individual cells are stochastic with a period near 21 hours (h), and using a stochastic clock network ensemble fitted by Markov Chain Monte Carlo implemented on general-purpose graphical processing units (or GPGPUs) we estimated that >94% of the variation in ccg-2 expression was stochastic (as opposed to experimental error). To overcome this stochasticity at the macroscopic level, cells must synchronize their oscillators. Using a classic measure of similarity in cell trajectories within droplets, the intraclass correlation (ICC), the synchronization surface ICC is measured on >25,000 cells as a function of the number of neighboring cells within a droplet and of time. The synchronization surface provides evidence that cells communicate, and synchronization varies with genotype.

Modulation of Circadian Gene Expression and Metabolic Compensation by the RCO-1 Corepressor of Neurospora crassa

Neurospora crassa is a model organism for the study of circadian clocks, molecular machineries that confer ∼24-hr rhythms to different processes at the cellular and organismal levels. The FREQUENCY (FRQ) protein is a central component of the Neurospora core clock, a transcription/translation negative feedback loop that controls genome-wide rhythmic gene expression. A genetic screen aimed at determining new components involved in the latter process identified regulation of conidiation 1 (rco-1), the ortholog of the Saccharomyces cerevisiae Tup1 corepressor, as affecting period length. By employing bioluminescent transcriptional and translational fusion reporters, we evaluated frq and FRQ expression levels in the rco-1 mutant background observing that, in contrast to prior reports, frq and FRQ expression are robustly rhythmic in the absence of RCO-1, although both amplitude and period length of the core clock are affected. Moreover, we detected a defect in metabolic compensation, such that high-glucose concentrations in the medium result in a significant decrease in period when RCO-1 is absent. Proteins physically interacting with RCO-1 were identified through co-immunoprecipitation and mass spectrometry; these include several components involved in chromatin remodeling and transcription, some of which, when absent, lead to a slight change in period. In the aggregate, these results indicate a dual role for RCO-1: although it is not essential for core-clock function, it regulates proper period and amplitude of core-clock dynamics and is also required for the rhythmic regulation of several clock-controlled genes.

Assessing the relevance of light for fungi: Implications and insights into the network of signal transmission

Light represents an important environmental cue, which provides information enabling fungi to prepare and react to the different ambient conditions between day and night. This adaptation requires both anticipation of the changing conditions, which is accomplished by daily rhythmicity of gene expression brought about by the circadian clock, and reaction to sudden illumination. Besides perception of the light signal, also integration of this signal with other environmental cues, most importantly nutrient availability, necessitates light-dependent regulation of signal transduction pathways and metabolic pathways. An influence of light and/or the circadian clock is known for the cAMP pathway, heterotrimeric G-protein signaling, mitogen-activated protein kinases, two-component phosphorelays, and Ca(2+) signaling. Moreover, also the target of rapamycin signaling pathway and reactive oxygen species as signal transducing elements are assumed to be connected to the light-response pathway. The interplay of the light-response pathway with signaling cascades results in light-dependent regulation of primary and secondary metabolism, morphology, development, biocontrol activity, and virulence. The frequent use of fungi in biotechnology as well as analysis of fungi in the artificial environment of a laboratory therefore requires careful consideration of still operative evolutionary heritage of these organisms. This review summarizes the diverse effects of light on fungi and the mechanisms they apply to deal both with the information content and with the harmful properties of light. Additionally, the implications of the reaction of fungi to light in a laboratory environment for experimental work and industrial applications are discussed.

The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species

The genus Trichoderma contains fungi with high relevance for humans, with applications in enzyme production for plant cell wall degradation and use in biocontrol. Here, we provide a broad, comprehensive overview of the genomic content of these species for "hot topic" research aspects, including CAZymes, transport, transcription factors, and development, along with a detailed analysis and annotation of less-studied topics, such as signal transduction, genome integrity, chromatin, photobiology, or lipid, sulfur, and nitrogen metabolism in T. reesei, T. atroviride, and T. virens, and we open up new perspectives to those topics discussed previously. In total, we covered more than 2,000 of the predicted 9,000 to 11,000 genes of each Trichoderma species discussed, which is >20% of the respective gene content. Additionally, we considered available transcriptome data for the annotated genes. Highlights of our analyses include overall carbohydrate cleavage preferences due to the different genomic contents and regulation of the respective genes. We found light regulation of many sulfur metabolic genes. Additionally, a new Golgi 1,2-mannosidase likely involved in N-linked glycosylation was detected, as were indications for the ability of Trichoderma spp. to generate hybrid galactose-containing N-linked glycans. The genomic inventory of effector proteins revealed numerous compounds unique to Trichoderma, and these warrant further investigation. We found interesting expansions in the Trichoderma genus in several signaling pathways, such as G-protein-coupled receptors, RAS GTPases, and casein kinases. A particularly interesting feature absolutely unique to T. atroviride is the duplication of the alternative sulfur amino acid synthesis pathway.

Circadian Control of Global Transcription

Circadian rhythms exist in most if not all organisms on the Earth and manifest in various aspects of physiology and behavior. These rhythmic processes are believed to be driven by endogenous molecular clocks that regulate rhythmic expression of clock-controlled genes (CCGs). CCGs consist of a significant portion of the genome and are involved in diverse biological pathways. The transcription of CCGs is tuned by rhythmic actions of transcription factors and circadian alterations in chromatin. Here, we review the circadian control of CCG transcription in five model organisms that are widely used, including cyanobacterium, fungus, plant, fruit fly, and mouse. Comparing the similarity and differences in the five organisms could help us better understand the function of the circadian clock, as well as its output mechanisms adapted to meet the demands of diverse environmental conditions.

Seeing the world differently: variability in the photosensory mechanisms of two model fungi

Light plays an important role for most organisms on this planet, serving either as a source of energy or information for the adaptation of biological processes to specific times of day. The fungal kingdom is estimated to contain well over a million species, possibly 10-fold more, and it is estimated that a majority of the fungi respond to light, eliciting changes in several physiological characteristics including pathogenesis, development and secondary metabolism. Two model organisms for photobiological studies have taken centre-stage over the last few decades--Neurospora crassa and Aspergillus nidulans. In this review, we will first discuss our understanding of the light response in N. crassa, about which the most is known, and will then juxtapose N. crassa with A. nidulans, which, as will be described below, provides an excellent template for understanding photosensory cross-talk. Finally, we will end with a commentary on the variability of the light response among other relevant fungi, and how our molecular understanding in the aforementioned model organisms still provides a strong base for dissecting light responses in such species.

The circadian system as an organizer of metabolism

The regulation of metabolism by circadian systems is believed to be a key reason for the extensive representation of circadian rhythms within the tree of life. Despite this, surprisingly little work has focused on the link between metabolism and the clock in Neurospora, a key model system in circadian research. The analysis that has been performed has focused on the unidirectional control from the clock to metabolism and largely ignored the feedback from metabolism on the clock. Recent efforts to understand these links have broken new ground, revealing bidirectional control from the clock to metabolism and vise-versa, showing just how strongly interconnected these two cellular systems can be in fungi. This review describes both well understood and emerging links between the clock and metabolic output of fungi as well as the role that metabolism plays in influencing the rhythm set by the clock.

Efficient Optimization of Stimuli for Model-Based Design of Experiments to Resolve Dynamical Uncertainty

This model-based design of experiments (MBDOE) method determines the input magnitudes of an experimental stimuli to apply and the associated measurements that should be taken to optimally constrain the uncertain dynamics of a biological system under study. The ideal global solution for this experiment design problem is generally computationally intractable because of parametric uncertainties in the mathematical model of the biological system. Others have addressed this issue by limiting the solution to a local estimate of the model parameters. Here we present an approach that is independent of the local parameter constraint. This approach is made computationally efficient and tractable by the use of: (1) sparse grid interpolation that approximates the biological system dynamics, (2) representative parameters that uniformly represent the data-consistent dynamical space, and (3) probability weights of the represented experimentally distinguishable dynamics. Our approach identifies data-consistent representative parameters using sparse grid interpolants, constructs the optimal input sequence from a greedy search, and defines the associated optimal measurements using a scenario tree. We explore the optimality of this MBDOE algorithm using a 3-dimensional Hes1 model and a 19-dimensional T-cell receptor model. The 19-dimensional T-cell model also demonstrates the MBDOE algorithm’s scalability to higher dimensions. In both cases, the dynamical uncertainty region that bounds the trajectories of the target system states were reduced by as much as 86% and 99% respectively after completing the designed experiments in silico. Our results suggest that for resolving dynamical uncertainty, the ability to design an input sequence paired with its associated measurements is particularly important when limited by the number of measurements.

Many mathematical models that have been developed for biological systems are limited because the complex systems are not well understood, the parameters are not known, and available data is limited and noisy. On the other hand, experiments to support model development are limited in terms of costs and time, feasible inputs and feasible measurements. MBDOE combines the mathematical models with experiment design to strategically design optimal experiments to obtain data that will contribute to the understanding of the systems. Our approach extends current capabilities of existing MBDOE techniques to make them more useful for scientists to resolve the trajectories of the system under study. It identifies the optimal conditions for stimuli and measurements that yield the most information about the system given the practical limitations. Exploration of the input space is not a trivial extension to MBDOE methods used for determining optimal measurements due to the nonlinear nature of many biological system models. The exploration of the system dynamics elicited by different inputs requires a computationally efficient and tractable approach. Our approach plans optimal experiments to reduce dynamical uncertainty in the output of selected target states of the biological system.

Analysis of Circadian Rhythms in the Basal Filamentous Ascomycete Pyronema confluens

Many organisms use circadian clocks to adapt to daily changes in the environment. Major insights into the molecular mechanisms of circadian oscillators have been gained through studies of the model organism Neurospora crassa; however, little is known about molecular components of circadian clocks in other fungi. An important part of the N. crassa circadian clock is the frequency (frq) gene, homologs of which can be found in Sordariomycetes, Dothideomycetes, and Leotiomycetes, but not Eurotiomycetes. Recently, we identified a frq homolog in Pyronema confluens, a member of the early-diverging Pezizomycete lineage of filamentous ascomycetes. The P. confluens FRQ shares many conserved domains with the N. crassa FRQ. However, there is no known morphological phenotype showing overt circadian rhythmicity in P. confluens. To investigate whether a molecular clock is present, we analyzed frq transcription in constant darkness, and found circadian oscillation of frq with a peak in the subjective morning. This rhythm was temperature compensated. To identify additional clock-controlled genes, we performed RNA sequencing of two time points (subjective morning and evening). Circadian expression of two morning-specific genes was verified by reverse transcription quantitative polymerase chain reaction (RT-qPCR) over a full time course, whereas expression of two putative morning-specific and five putative evening-specific genes could not be verified as circadian. frq expression was synchronized, but not entrained by light. In summary, we have found evidence for two of the three main properties of circadian rhythms (free-running rhythm, temperature compensation) in P. confluens, suggesting that a circadian clock with rhythmically expressed frq is present in this basal filamentous ascomycete.

Quantitative Circadian Phosphoproteomic Analysis of Arabidopsis Reveals Extensive Clock Control of Key Components in Physiological, Metabolic, and Signaling Pathways

The circadian clock provides adaptive advantages to an organism, resulting in increased fitness and survival. The phosphorylation events that regulate circadian-dependent signaling and the processes which post-translationally respond to clock-gated signals are largely unknown. To better elucidate post-translational events tied to the circadian system we carried out a survey of circadian-regulated protein phosphorylation events in Arabidopsis seedlings. A large-scale mass spectrometry-based quantitative phosphoproteomics approach employing TiO2-based phosphopeptide enrichment techniques identified and quantified 1586 phosphopeptides on 1080 protein groups. A total of 102 phosphopeptides displayed significant changes in abundance, enabling the identification of specific patterns of response to circadian rhythms. Our approach was sensitive enough to quantitate oscillations in the phosphorylation of low abundance clock proteins (early flowering4; ELF4 and pseudoresponse regulator3; PRR3) as well as other transcription factors and kinases. During constant light, extensive cyclic changes in phosphorylation status occurred in critical regulators, implicating direct or indirect regulation by the circadian system. These included proteins influencing transcriptional regulation, translation, metabolism, stress and phytohormones-mediated responses. We validated our analysis using the elf4-211 allele, in which an S45L transition removes the phosphorylation herein identified. We show that removal of this phosphorylatable site diminishes interaction with early flowering3 (ELF3), a key partner in a tripartite evening complex required for circadian cycling. elf4-211 lengthens period, which increases with increasing temperature, relative to the wild type, resulting in a more stable temperature compensation of circadian period over a wider temperature range.

Ascomycete fungal actins differentially support plant spatial cell and organ development

Actin interacts with a wide variety of cytoplasmic and nuclear proteins to support spatial development in nearly all eukaryotes. Null mutations in plant vegetative actins produce dramatically altered cell, tissue, and organ morphologies. Animal cytoplasmic actins (e.g., human HsACTB, HsACTG1) and some ancestral protist actins fully suppress these mutant phenotypes suggesting that some animal, plant, and protist actins share functional competence for spatial development. Considering that fungi have a phylogenetic origin closer to animals than plants, we were interested to explore whether the fungal actins may have this same capacity to function in plants and support development. We ectopically expressed actins from four highly divergent ascomycete fungi in two different Arabidopsis double vegetative actin null mutants. We found that expression of actin from the earliest diverging ascomycete subphyla, the archiascomycete Schizosaccharomyces pombe, qualitatively and quantitatively suppressed the root cell polarity and root organ developmental defects of act8/act7 mutants and the root-hairless cell elongation phenotype of act2/act8 mutants. Interestingly, the actin from the pyrenomycete Neurospora crassa was modestly effective in the suppression of vegetative actin mutant phenotypes. In contrast, actins from the saccharomycetes Saccharomyces cerevisiae and Candida albicans were unable to support any aspect of plant development, and moreover induced severe dwarfism and sterility. These data imply that basal fungi inherited an actin with full competence for spatial development from their protist ancestor and maintained it via non-progressive sequence evolution, while the later more derived fungal species lost these activities.

Microarray hybridization analysis of light-dependent gene expression in Penicillium chrysogenum identifies bZIP transcription factor PcAtfA

The fungal velvet complex is a light-dependent master regulator of secondary metabolism and development in the major penicillin producer, Penicillium chrysogenum. However, the light-dependent mechanism is unclear. To identify velvet-dependent transcriptional regulators that show light-regulated expression, we performed microarray hybridizations with RNA isolated from P. chrysogenum ΔPcku70 cultures grown under 13 different long-term, light-dependent growth conditions. We compared these expression data to data from two velvet complex deletion mutants; one lacked a subunit of the velvet complex (ΔPcvelA), and the other lacked a velvet-associated protein (ΔPclaeA). We sought to identify genes that were up-regulated in light, but down-regulated in ΔPcvelA and ΔPclaeA. We identified 148 co-regulated genes that displayed this regulatory pattern. In silico analyses of the co-regulated genes identified six proteins with fungal-specific transcription factor domains. Among these, we selected the bZIP transcription factor, PcAtfA, for functional characterization in deletion and complementation strains. Our data clearly indicates that PcAtfA governs spore germination. This comparative analysis of different microarray hybridization data sets provided results that may be useful for identifying genes for future functional analyses.

Dissecting the mechanisms of the clock in Neurospora

The circadian clock exists to synchronize inner physiology with the external world, allowing life to anticipate and adapt to the continual changes that occur in an organism's environment. The clock architecture is highly conserved, present in almost all major branches of life. Within eukaryotes, the filamentous fungus Neurospora crassa has consistently been used as an excellent model organism to uncover the basic circadian physiology and molecular biology. The Neurospora model has elucidated our fundamental understanding of the clock as nested positive and negative feedback loop, regulated by transcriptional and posttranscriptional processes. This review will examine the basics of circadian rhythms in the model filamentous fungus N. crassa as well as highlight the output of the clock in Neurospora and the reasons that N. crassa has continued to be a strong model for the study of circadian rhythms. It will also synopsize classical and emerging methods in the study of the circadian clock.

Selection and evaluation of reference genes for expression studies with quantitative PCR in the model fungus Neurospora crassa under different environmental conditions in continuous culture

Neurospora crassa has served as a model organism for studying circadian pathways and more recently has gained attention in the biofuel industry due to its enhanced capacity for cellulase production. However, in order to optimize N. crassa for biotechnological applications, metabolic pathways during growth under different environmental conditions must be addressed. Reverse-transcription quantitative PCR (RT-qPCR) is a technique that provides a high-throughput platform from which to measure the expression of a large set of genes over time. The selection of a suitable reference gene is critical for gene expression studies using relative quantification, as this strategy is based on normalization of target gene expression to a reference gene whose expression is stable under the experimental conditions. This study evaluated twelve candidate reference genes for use with N. crassa when grown in continuous culture bioreactors under different light and temperature conditions. Based on combined stability values from NormFinder and Best Keeper software packages, the following are the most appropriate reference genes under conditions of: (1) light/dark cycling: btl, asl, and vma1; (2) all-dark growth: btl, tbp, vma1, and vma2; (3) temperature flux: btl, vma1, act, and asl; (4) all conditions combined: vma1, vma2, tbp, and btl. Since N. crassa exists as different cell types (uni- or multi-nucleated), expression changes in a subset of the candidate genes was further assessed using absolute quantification. A strong negative correlation was found to exist between ratio and threshold cycle (CT) values, demonstrating that CT changes serve as a reliable reflection of transcript, and not gene copy number, fluctuations. The results of this study identified genes that are appropriate for use as reference genes in RT-qPCR studies with N. crassa and demonstrated that even with the presence of different cell types, relative quantification is an acceptable method for measuring gene expression changes during growth in bioreactors.

Analysis of clock-regulated genes in Neurospora reveals widespread posttranscriptional control of metabolic potential

Neurospora crassa has been for decades a principal model for filamentous fungal genetics and physiology as well as for understanding the mechanism of circadian clocks. Eukaryotic fungal and animal clocks comprise transcription-translation-based feedback loops that control rhythmic transcription of a substantial fraction of these transcriptomes, yielding the changes in protein abundance that mediate circadian regulation of physiology and metabolism: Understanding circadian control of gene expression is key to understanding eukaryotic, including fungal, physiology. Indeed, the isolation of clock-controlled genes (ccgs) was pioneered in Neurospora where circadian output begins with binding of the core circadian transcription factor WCC to a subset of ccg promoters, including those of many transcription factors. High temporal resolution (2-h) sampling over 48 h using RNA sequencing (RNA-Seq) identified circadianly expressed genes in Neurospora, revealing that from ∼10% to as much 40% of the transcriptome can be expressed under circadian control. Functional classifications of these genes revealed strong enrichment in pathways involving metabolism, protein synthesis, and stress responses; in broad terms, daytime metabolic potential favors catabolism, energy production, and precursor assembly, whereas night activities favor biosynthesis of cellular components and growth. Discriminative regular expression motif elicitation (DREME) identified key promoter motifs highly correlated with the temporal regulation of ccgs. Correlations between ccg abundance from RNA-Seq, the degree of ccg-promoter activation as reported by ccg-promoter-luciferase fusions, and binding of WCC as measured by ChIP-Seq, are not strong. Therefore, although circadian activation is critical to ccg rhythmicity, posttranscriptional regulation plays a major role in determining rhythmicity at the mRNA level.

A MINE alternative to D-optimal designs for the linear model

Doing large-scale genomics experiments can be expensive, and so experimenters want to get the most information out of each experiment. To this end the Maximally Informative Next Experiment (MINE) criterion for experimental design was developed. Here we explore this idea in a simplified context, the linear model. Four variations of the MINE method for the linear model were created: MINE-like, MINE, MINE with random orthonormal basis, and MINE with random rotation. Each method varies in how it maximizes the MINE criterion. Theorem 1 establishes sufficient conditions for the maximization of the MINE criterion under the linear model. Theorem 2 establishes when the MINE criterion is equivalent to the classic design criterion of D-optimality. By simulation under the linear model, we establish that the MINE with random orthonormal basis and MINE with random rotation are faster to discover the true linear relation with regression coefficients and observations when . We also establish in simulations with , , and 1000 replicates that these two variations of MINE also display a lower false positive rate than the MINE-like method and additionally, for a majority of the experiments, for the MINE method.

Fungal photobiology: visible light as a signal for stress, space and time

Visible light is an important source of energy and information for much of life on this planet. Though fungi are neither photosynthetic nor capable of observing adjacent objects, it is estimated that the majority of fungal species display some form of light response, ranging from developmental decision-making to metabolic reprogramming to pathogenesis. As such, advances in our understanding of fungal photobiology will likely reach the broad fields impacted by these organisms, including agriculture, industry and medicine. In this review, we will first describe the mechanisms by which fungi sense light and then discuss the selective advantages likely imparted by their ability to do so.

The aging biological clock in Neurospora crassa

The biological clock affects aging through ras-1 (bd) and lag-1, and these two longevity genes together affect a clock phenotype and the clock oscillator in Neurospora crassa. Using an automated cell-counting technique for measuring conidial longevity, we show that the clock-associated genes lag-1 and ras-1 (bd) are true chronological longevity genes. For example, wild type (WT) has an estimated median life span of 24 days, while the double mutant lag-1, ras-1 (bd) has an estimated median life span of 120 days for macroconidia. We establish the biochemical function of lag-1 by complementing LAG1 and LAC1 in Saccharomyces cerevisiae with lag-1 in N. crassa. Longevity genes can affect the clock as well in that, the double mutant lag-1, ras-1 (bd) can stop the circadian rhythm in asexual reproduction (i.e., banding in race tubes) and lengthen the period of the frequency oscillator to 41 h. In contrast to the ras-1 (bd), lag-1 effects on chronological longevity, we find that this double mutant undergoes replicative senescence (i.e., the loss of replication function with time), unlike WT or the single mutants, lag-1 and ras-1 (bd). These results support the hypothesis that sphingolipid metabolism links aging and the biological clock through a common stress response

Genome-wide characterization of light-regulated genes in Neurospora crassa

The filamentous fungus Neurospora crassa responds to light in complex ways. To thoroughly study the transcriptional response of this organism to light, RNA-seq was used to analyze capped and polyadenylated mRNA prepared from mycelium grown for 24 hr in the dark and then exposed to light for 0 (control) 15, 60, 120, and 240 min. More than three-quarters of all defined protein coding genes (79%) were expressed in these cells. The increased sensitivity of RNA-seq compared with previous microarray studies revealed that the RNA levels for 31% of expressed genes were affected two-fold or more by exposure to light. Additionally, a large class of mRNAs, enriched for transcripts specifying products involved in rRNA metabolism, showed decreased expression in response to light, indicating a heretofore undocumented effect of light on this pathway. Based on measured changes in mRNA levels, light generally increases cellular metabolism and at the same time causes significant oxidative stress to the organism. To deal with this stress, protective photopigments are made, antioxidants are produced, and genes involved in ribosome biogenesis are transiently repressed.

Interconnections of reactive oxygen species homeostasis and circadian rhythm in Neurospora crassa

SIGNIFICANCE

Both circadian rhythm and the production of reactive oxygen species (ROS) are fundamental features of aerobic eukaryotic cells. The circadian clock enhances the fitness of organisms by enabling them to anticipate cycling changes in the surroundings. ROS generation in the cell is often altered in response to environmental changes, but oscillations in ROS levels may also reflect endogenous metabolic fluctuations governed by the circadian clock. On the other hand, an effective regulation and timing of antioxidant mechanisms may be crucial in the defense of cellular integrity. Thus, an interaction between the circadian timekeeping machinery and ROS homeostasis or signaling in both directions may be of advantage at all phylogenetic levels.

RECENT ADVANCES

The Frequency-White Collar-1 and White Collar-2 oscillator (FWO) of the filamentous fungus Neurospora crassa is well characterized at the molecular level. Several members of the ROS homeostasis were found to be controlled by the circadian clock, and ROS levels display circadian rhythm in Neurospora. On the other hand, multiple data indicate that ROS affect the molecular oscillator.

CRITICAL ISSUES

Increasing evidence suggests the interplay between ROS homeostasis and oscillators that may be partially or fully independent of the FWO. In addition, ROS may be part of a complex cellular network synchronizing non-transcriptional oscillators with timekeeping machineries based on the classical transcription-translation feedback mechanism.

FUTURE DIRECTIONS

Further investigations are needed to clarify how the different layers of the bidirectional interactions between ROS homeostasis and circadian regulation are interconnected.

Model-based analysis for qualitative data: an application in Drosophila germline stem cell regulation

Discovery in developmental biology is often driven by intuition that relies on the integration of multiple types of data such as fluorescent images, phenotypes, and the outcomes of biochemical assays. Mathematical modeling helps elucidate the biological mechanisms at play as the networks become increasingly large and complex. However, the available data is frequently under-utilized due to incompatibility with quantitative model tuning techniques. This is the case for stem cell regulation mechanisms explored in the Drosophila germarium through fluorescent immunohistochemistry. To enable better integration of biological data with modeling in this and similar situations, we have developed a general parameter estimation process to quantitatively optimize models with qualitative data. The process employs a modified version of the Optimal Scaling method from social and behavioral sciences, and multi-objective optimization to evaluate the trade-off between fitting different datasets (e.g. wild type vs. mutant). Using only published imaging data in the germarium, we first evaluated support for a published intracellular regulatory network by considering alternative connections of the same regulatory players. Simply screening networks against wild type data identified hundreds of feasible alternatives. Of these, five parsimonious variants were found and compared by multi-objective analysis including mutant data and dynamic constraints. With these data, the current model is supported over the alternatives, but support for a biochemically observed feedback element is weak (i.e. these data do not measure the feedback effect well). When also comparing new hypothetical models, the available data do not discriminate. To begin addressing the limitations in data, we performed a model-based experiment design and provide recommendations for experiments to refine model parameters and discriminate increasingly complex hypotheses.

We developed a process to quantitatively fit mathematical models using qualitative data, and applied it in the study of how stem cells are regulated in the fruit fly ovary. The available published data we collected are fluorescent images of protein and mRNA expression from genetic experiments. Despite lacking quantitative data, the new process makes available quantitative model analysis techniques to reliably compare different models and guide future experiments. We found that the current consensus regulatory model is supported, but that the data are indeed insufficient to address more complex hypotheses. With the quantitatively fit models, we evaluated hypothetical experiments and estimated which future measurements should best refine or test models. The model fitting process we have developed is applicable to many biological studies where qualitative data are common, and can accelerate progress through more efficient experimentation.

Circadian rhythms synchronize mitosis in Neurospora crassa

The cell cycle and the circadian clock communicate with each other, resulting in circadian-gated cell division cycles. Alterations in this network may lead to diseases such as cancer. Therefore, it is critical to identify molecular components that connect these two oscillators. However, molecular mechanisms between the clock and the cell cycle remain largely unknown. A model filamentous fungus, Neurospora crassa, is a multinucleate system used to elucidate molecular mechanisms of circadian rhythms, but not used to investigate the molecular coupling between these two oscillators. In this report, we show that a conserved coupling between the circadian clock and the cell cycle exists via serine/threonine protein kinase-29 (STK-29), the Neurospora homolog of mammalian WEE1 kinase. Based on this finding, we established a mathematical model that predicts circadian oscillations of cell cycle components and circadian clock-dependent synchronized nuclear divisions. We experimentally demonstrate that G1 and G2 cyclins, CLN-1 and CLB-1, respectively, oscillate in a circadian manner with bioluminescence reporters. The oscillations of clb-1 and stk-29 gene expression are abolished in a circadian arrhythmic frq(ko) mutant. Additionally, we show the light-induced phase shifts of a core circadian component, frq, as well as the gene expression of the cell cycle components clb-1 and stk-29, which may alter the timing of divisions. We then used a histone hH1-GFP reporter to observe nuclear divisions over time, and show that a large number of nuclear divisions occur in the evening. Our findings demonstrate the circadian clock-dependent molecular dynamics of cell cycle components that result in synchronized nuclear divisions in Neurospora.

Specifying informative experiment stimulation conditions for resolving dynamical uncertainty in biological systems

A computationally efficient model-based design of experiments (MBDOE) strategy is developed to plan an optimal experiment by specifying the experimental stimulation magnitudes and measurement points. The strategy is extended from previous work which optimized the experimental design over a space of measurable species and time points. We include system inputs (stimulation conditions) in the experiment design search to investigate if the addition of perturbations enhances the ability of the MBDOE method to resolve uncertainties in system dynamics. The MBDOE problem is made computationally tractable by using a sparse-grid approximation of the model output dynamics, pre-specifying the time points at which the input or experimental perturbations can be applied, and creating scenario trees to explore the endogenous uncertainty. Consecutive scenario trees are used to determine the best input magnitudes and select the optimal associated measurement species and time points. We demonstrate the effectiveness of this strategy on a T-Cell Receptor (TCR) signaling pathway model.

Circadian yin-yang regulation and its manipulation to globally reprogram gene expression

BACKGROUND

The cyanobacterial circadian program exerts genome-wide control of gene expression. KaiC undergoes rhythms of phosphorylation that are regulated by interactions with KaiA and KaiB. The phosphorylation status of KaiC is thought to mediate global transcription via output factors SasA, CikA, LabA, RpaA, and RpaB. Overexpression of kaiC has been reported to globally repress gene expression.

RESULTS

Here, we show that the positive circadian component KaiA upregulates "subjective dusk" genes and that its overexpression deactivates rhythmic gene expression without significantly affecting growth rates in constant light. We analyze the global patterns of expression that are regulated by KaiA versus KaiC and find in contrast to the previous report of KaiC repression that there is a "yin-yang" regulation of gene expression whereby kaiA overexpression activates "dusk genes" and represses "dawn genes," whereas kaiC overexpression complementarily activates dawn genes and represses dusk genes. Moreover, continuous induction of kaiA latched KaiABC-regulated gene expression to provide constitutively increased transcript levels of diverse endogenous and heterologous genes that are expressed in the predominant subjective dusk phase. In addition to analyzing KaiA regulation of endogenous gene expression, we apply these insights to the expression of heterologous proteins whose products are of potential value, namely human proinsulin, foreign luciferase, and exogenous hydrogenase.

CONCLUSIONS

Both KaiC and KaiA complementarily contribute to the regulation of circadian gene expression via yin-yang switching. Circadian patterns can be reprogrammed by overexpression of kaiA or kaiC to constitutively enhance gene expression, and this reprogramming can improve 24/7 production of heterologous proteins that are useful as pharmaceuticals or biofuels.