Behavioral, Neurobiological, and Genetic Analysis of the Circadian Mutant Duper

The recently discovered circadian mutant hamster duper has a short period of ~23 hours and exhibits exaggerated phase shifts in response to a 15-min light pulse. To increase the understanding of the duper mutation, I performed behavioral, neurobiological, and genetic experiments. Behavioral studies using photic and non-photic stimuli found that large phase shifts exhibited by duper hamsters are specific to photic cues, but not to phase. Additionally, 2/3 of duper hamsters, but no WTs, displayed transient ultradian wheel-running patterns when transferred from light to dark at CT 18. This suggests that the mutation may weaken coupling among components of the circadian pacemaker. Anatomical and immunocytochemical analysis of the SCN was used to examine the neurobiological mechanisms of large light-induced phase shifts in dupers. Brains were collected from duper and WT hamsters at CT 12 and 15 as well as 1, 2, 3, 6 and 9 hours following a light pulse, or control handling, at CT 15. Surprisingly, the only difference in PER1 (a core clock protein) expression in the SCN between dupers and WTs was seen 2-hours after a light pulse; duper hamsters displayed a significantly greater percentage of retinorecipient VIP cells co-labeled with PER1 compared to WTs. Additional differences between genotypes occurred 9 hours after CT15 (controls). In the SCN, the number of PER1-ir cells was significantly greater in WT than duper hamsters, however this finding was reversed in the PVN. This anatomical mismatch suggests the mutation may affect signaling between the SCN and extra-SCN oscillators. Finally, to identify the genetic basis of the duper phenotype, I crossed dupers with a novel ecotype in order to perform fast homozygosity mapping. Duper transmitted onto the novel ecotype with the predicted Mendelian inheritance of phenotype. I collected DNA from F2 duper hamsters, and expect fast homozygosity mapping will identify candidate genetic regions of the duper mutation. Additional behavioral experiments in F2 dupers demonstrated that duper hamsters are resistant to jet lag. As duper is a unique circadian mutation, understanding of the behavioral phenotype, neurobiological mechanism, and genetic basis of the duper mutation will greatly increase our knowledge of the circadian system

Modeling gene regulatory networks through data integration

Modeling gene regulatory networks has become a problem of great interest in biology and medical research. Most common methods for learning regulatory dependencies rely on observations in the form of gene expression data. In this dissertation, computational models for gene regulation have been developed based on constrained regression by integrating comprehensive gene expression data for M. tuberculosis with genome-scale ChIP-Seq interaction data. The resulting models confirmed predictive power for expression in independent stress conditions and identified mechanisms driving hypoxic adaptation and lipid metabolism in M. tuberculosis. I then used the regulatory network model for M. tuberculosis to identify factors responding to stress conditions and drug treatments, revealing drug synergies and conditions that potentiate drug treatments. These results can guide and optimize design of drug treatments for this pathogen. I took the next step in this direction, by proposing a new probabilistic framework for learning modular structures in gene regulatory networks from gene expression and protein-DNA interaction data, combining the ideas of module networks and stochastic blockmodels. These models also capture combinatorial interactions between regulators. Comparisons with other network modeling methods that rely solely on expression data, showed the essentiality of integrating ChIP-Seq data in identifying direct regulatory links in M. tuberculosis. Moreover, this work demonstrates the theoretical advantages of integrating ChIP-Seq data for the class of widely-used module network models. The systems approach and statistical modeling presented in this dissertation can also be applied to problems in other organisms. A similar approach was taken to model the regulatory network controlling genes with circadian gene expression in Neurospora crassa, through integrating time-course expression data with ChIP-Seq data. The models explained combinatorial regulations leading to different phase differences in circadian rhythms. The Neurospora crassa network model also works as a tool to manipulate the phases of target genes

Quantifying the Life Stages of a Biomolecule: Implications for the Circadian Transcriptome

Viele biologische Prozesse im Verhalten von ganzen Organismen, aber auch in den Prozessen und der biochemischen Zusammensetzung von Zellen zeigen einen zirkadianen Rhythmus, also einen Rhythmus mit einer Periode von etwa 24 Stunden. Diese 24-Stunden-Rhythmen sind in der Genexpression auf allen Ebenen zu finden: von der Tran- skriptionsinitiation bis zur Proteindegradation. Auf Transkriptebene, zirkadiane mRNA-Produktion und mRNA-Abundanz ist umfassend gemessen. Auf der anderen Seite, zirkadiane posttranskriptionelle Regulation ist weit weniger verstanden. In dieser Arbeit untersuche ich, wie bisher ungemessene, posttranskriptionelle Prozesse die rhythmischen Eigenschaften der Genexpression beeinflussen. Dazu beschreibe ich die Lebensstadien eines Bio-Moleküls mit einem Modell-Motiv, einer einfachen Differentialgleichung mit zeitabhängigen, rhythmischen Raten. Als erstes diskutiere ich die Einschränkungen von Phase und Amplitude zirkadianer Transkripte, die nur von konstanter PTR beeinflusst werden. Bei vielen gemessenen Transkripten sind diese Einschränkungen verletzt. In diesen Fällen muss es eine rhythmische PTR geben. Ich untersuche, welche rhythmische PTR diese Fälle erklären können und führe einen statistischen Test ein, der auf unbeobachtete, rhythmische PTR testet. Durch die Analyse zweier Datensätze von Mausleber und -niere finde ich, dass 18% aller zirkadianen Gene in Niere und 34% in Leber rhythmisch posttranskriptionell reguliert sind. Im zweiten Teil analysiere ich weitere Aspekte von PTR in einem Hypothesen-getriebenen Ansatz. Ich zeige, dass Spleißen mit einem Rhythmus von 24 Stunden 12 Stunden-Rhythmen in der Abundanz von mRNA erzeugen kann. Als nächstes schlage ich ein Modell vor, das rhythmische Degradation von Mitgliedern der zirkadianen Uhr beschreibt. Schließlich erweitere ich das Modell-Grundmotiv zu einer partiellen Differentialgleichung (PDG), die das “Altern” von Molekülen beschreibt.In almost all organisms on Earth, many behavioral, physiological, and biochemical activities oscillate with a circadian rhythm, a rhythm with a period of about 24 hours. In gene expression, the 24-hour-rhythm can be found on all stages: from transcription initiation to protein degradation. On the transcript level, circadian mRNA production and mRNA abundance are comprehensively charted through numerous genome-wide high throughput studies. Circadian post-transcriptional regulation, however, is less well understood. In this thesis, I will investigate how unobserved post-transcriptional processes influence rhythmic properties of gene expression. To this end, I quantify the life-stages of biomolecules using one modeling motif, a simple ordinary differential equation describing production and degradation with time-dependent rhythmic rates. This basic modeling motif is systematically varied to examine and discuss various influences of post-transcriptional regulation (PTR) on circadian mRNA expression. I first discuss the restrictions of rhythmic phase and amplitude of circadian transcripts influenced by non-rhythmic PTR. For many genes these restrictions are violated and we have to assume the existence of a rhythmic PTR. I discuss which rhythmic PTR can explain these findings and further introduce a statistical test to quantify the extent of unobserved rhythmic PTR. Analyzing two data sets on mouse liver and kidney, I find that 18% of circadian genes in kidney and 34% in liver are under rhythmic post-transcriptional control. In a second part, I analyze more specific aspects of PTR in a hypothesis-driven approach. Firstly, I find that splicing with a rhythm of 24 hours is able to generate 12-hour rhythms in abundance of mature mRNA. Secondly, I propose and analyze a model to investigate rhythmic degradation of core clock genes. And finally, I extend the core modeling motif to a partial differential equation (PDE) model that accounts for the “aging” process of molecules

Modeling complex cell regulation in the zebrafish circadian clock

The interdisciplinary "systems biology" approach of combining traditional biological investigations with tools from the mathematical and computer sciences has enabled novel insights into many highly complex and dynamic biological systems. The use of models has, for instance, revealed much about the intricate feedback mechanisms and acute importance of gene regulatory networks, and one such network of special note is our internal time keeper, or circadian clock. The circadian clock plays a pivotal role in modulating critical physiological processes, and has also been implicated, either directly or indirectly, in a whole range of pathological states. This research project investigates how the underlying dynamics of the circadian clock in the zebrafish model organism may be captured by a mathematical model, considering in particular the entrainment effect due to external cues such as light. Simulated data is contrasted with experimental results from different light regime experiments to validate the model and guide its refinement. Furthermore, various statistical methods are implemented to process the raw data and support its analysis. Extending the initial deterministic approach to take into account stochastic effects and additive population level effects emerges as a powerful means of representing the circadian signal decay in prolonged darkness, as well as light initiated re-synchronization as a strong component of entrainment. Consequently, it emerges that stochastic effects may be considered an essential feature of the circadian clock in zebrafish. A further cornerstone of the project is the implementation of an integrated simulation environment, including a Sequential Monte Carlo parameter estimation function, which succeeds in predicting a range of previously determined and also novel suitable parameter values. However, considerable difficulties in obtaining parameter values that satisfy the entire range of important target values simultaneously highlights the inherent complexity of accurately simulating the circadian clock