Circadian distribution of epileptiform discharges in epilepsy: Candidate mechanisms of variability.

Epilepsy is a serious neurological disorder characterised by a tendency to have recurrent, spontaneous, seizures. Classically, seizures are assumed to occur at random. However, recent research has uncovered underlying rhythms both in seizures and in key signatures of epilepsy-so-called interictal epileptiform activity-with timescales that vary from hours and days through to months. Understanding the physiological mechanisms that determine these rhythmic patterns of epileptiform discharges remains an open question. Many people with epilepsy identify precipitants of their seizures, the most common of which include stress, sleep deprivation and fatigue. To quantify the impact of these physiological factors, we analysed 24-hour EEG recordings from a cohort of 107 people with idiopathic generalized epilepsy. We found two subgroups with distinct distributions of epileptiform discharges: one with highest incidence during sleep and the other during day-time. We interrogated these data using a mathematical model that describes the transitions between background and epileptiform activity in large-scale brain networks. This model was extended to include a time-dependent forcing term, where the excitability of nodes within the network could be modulated by other factors. We calibrated this forcing term using independently-collected human cortisol (the primary stress-responsive hormone characterised by circadian and ultradian patterns of secretion) data and sleep-staged EEG from healthy human participants. We found that either the dynamics of cortisol or sleep stage transition, or a combination of both, could explain most of the observed distributions of epileptiform discharges. Our findings provide conceptual evidence for the existence of underlying physiological drivers of rhythms of epileptiform discharges. These findings should motivate future research to explore these mechanisms in carefully designed experiments using animal models or people with epilepsy.

Modelling the dynamic interaction of systemic inflammation and the hypothalamic-pituitary-adrenal (HPA) axis during and after cardiac surgery.

Major surgery and critical illness produce a potentially life-threatening systemic inflammatory response. The hypothalamic-pituitary-adrenal (HPA) axis is one of the key physiological systems that counterbalances this systemic inflammation through changes in adrenocorticotrophic hormone (ACTH) and cortisol. These hormones normally exhibit highly correlated ultradian pulsatility with an amplitude modulated by circadian processes. However, these dynamics are disrupted by major surgery and critical illness. In this work, we characterize the inflammatory, ACTH and cortisol responses of patients undergoing cardiac surgery and show that the HPA axis response can be classified into one of three phenotypes: single-pulse, two-pulse and multiple-pulse dynamics. We develop a mathematical model of cortisol secretion and metabolism that predicts the physiological mechanisms responsible for these different phenotypes. We show that the effects of inflammatory mediators are important only in the single-pulse pattern in which normal pulsatility is lost-suggesting that this phenotype could be indicative of the greatest inflammatory response. Investigating whether and how these phenotypes are correlated with clinical outcomes will be critical to patient prognosis and designing interventions to improve recovery.

Fast dynamics in the HPA axis: Insight from mathematical and experimental studies.

The activity of the hypothalamic-pituitary-adrenal (HPA) axis is characterised by complex dynamics spanning several timescales. This ranges from slow circadian rhythms in blood hormone concentration to faster ultradian pulses of hormone secretion and even more rapid oscillations in electrical and calcium activity in neuroendocrine cells of the hypothalamus and pituitary gland. Here, we focus on the system's oscillations on the short timescale. We highlight some of the mathematical modelling and experimental work that has been carried out to characterise the mechanisms regulating this highly dynamic mode of neuroendocrine signalling and discuss some future directions that may be explored to enhance understanding of HPA function.

Dynamic Hormone Control of Stress and Fertility.

Neuroendocrine axes display a remarkable diversity of dynamic signaling processes relaying information between the brain, endocrine glands, and peripheral target tissues. These dynamic processes include oscillations, elastic responses to perturbations, and plastic long term changes observed from the cellular to the systems level. While small transient dynamic changes can be considered physiological, larger and longer disruptions are common in pathological scenarios involving more than one neuroendocrine axes, suggesting that a robust control of hormone dynamics would require the coordination of multiple neuroendocrine clocks. The idea of apparently different axes being in fact exquisitely intertwined through neuroendocrine signals can be investigated in the regulation of stress and fertility. The stress response and the reproductive cycle are controlled by the Hypothalamic-Pituitary-Adrenal (HPA) axis and the Hypothalamic-Pituitary-Gonadal (HPG) axis, respectively. Despite the evidence surrounding the effects of stress on fertility, as well as of the reproductive cycle on stress hormone dynamics, there is a limited understanding on how perturbations in one neuroendocrine axis propagate to the other. We hypothesize that the links between stress and fertility can be better understood by considering the HPA and HPG axes as coupled systems. In this manuscript, we investigate neuroendocrine rhythms associated to the stress response and reproduction by mathematically modeling the HPA and HPG axes as a network of interlocked oscillators. We postulate a network architecture based on physiological data and use the model to predict responses to stress perturbations under different hormonal contexts: normal physiological, gonadectomy, hormone replacement with estradiol or corticosterone (CORT), and high excess CORT (hiCORT) similar to hypercortisolism in humans. We validate our model predictions against experiments in rodents, and show how the dynamic responses of these endocrine axes are consistent with our postulated network architecture. Importantly, our model also predicts the conditions that ensure robustness of fertility to stress perturbations, and how chronodisruptions in glucocorticoid hormones can affect the reproductive axis' ability to withstand stress. This insight is key to understand how chronodisruption leads to disease, and to design interventions to restore normal rhythmicity and health.

Diving into the brain: deep-brain imaging techniques in conscious animals.

In most species, survival relies on the hypothalamic control of endocrine axes that regulate critical functions such as reproduction, growth, and metabolism. For decades, the complexity and inaccessibility of the hypothalamic-pituitary axis has prevented researchers from elucidating the relationship between the activity of endocrine hypothalamic neurons and pituitary hormone secretion. Indeed, the study of central control of endocrine function has been largely dominated by 'traditional' techniques that consist of studying in vitro or ex vivo isolated cell types without taking into account the complexity of regulatory mechanisms at the level of the brain, pituitary and periphery. Nowadays, by exploiting modern neuronal transfection and imaging techniques, it is possible to study hypothalamic neuron activity in situ, in real time, and in conscious animals. Deep-brain imaging of calcium activity can be performed through gradient-index lenses that are chronically implanted and offer a 'window into the brain' to image multiple neurons at single-cell resolution. With this review, we aim to highlight deep-brain imaging techniques that enable the study of neuroendocrine neurons in awake animals whilst maintaining the integrity of regulatory loops between the brain, pituitary and peripheral glands. Furthermore, to assist researchers in setting up these techniques, we discuss the equipment required and include a practical step-by-step guide to performing these deep-brain imaging studies.

Ion channel noise shapes the electrical activity of endocrine cells.

Endocrine cells in the pituitary gland typically display either spiking or bursting electrical activity, which is related to the level of hormone secretion. Recent work, which combines mathematical modelling with dynamic clamp experiments, suggests the difference is due to the presence or absence of a few large-conductance potassium channels. Since endocrine cells only contain a handful of these channels, it is likely that stochastic effects play an important role in the pattern of electrical activity. Here, for the first time, we explicitly determine the effect of such noise by studying a mathematical model that includes the realistic noisy opening and closing of ion channels. This allows us to investigate how noise affects the electrical activity, examine the origin of spiking and bursting, and determine which channel types are responsible for the greatest noise. Further, for the first time, we address the role of cell size in endocrine cell electrical activity, finding that larger cells typically display more bursting, while the smallest cells almost always only exhibit spiking behaviour.

A Novel Automated System Yields Reproducible Temporal Feeding Patterns in Laboratory Rodents.

BACKGROUND: The impact of temporal feeding patterns remains a major unanswered question in nutritional science. Progress has been hampered by the absence of a reliable method to impose temporal feeding in laboratory rodents, without the confounding influence of food-hoarding behavior. OBJECTIVE: The aim of this study was to develop and validate a reliable method for supplying crushed diets to laboratory rodents in consistent, relevant feeding patterns for prolonged periods. METHODS: We programmed our experimental feeding station to deliver a standard diet [StD; Atwater Fuel Energy (AFE) 13.9% fat] or high-fat diet (HFD; AFE 45% fat) during nocturnal grazing [providing 1/24th of the total daily food intake (tdF/I) of ad libitum-fed controls every 30 min] and meal-fed (3 × 1-h periods of ad libitum feeding) patterns in male rats (Sprague-Dawley: 4 wk old, 72-119 g) and mice [C57/Bl6J wild-type (WT): 6 mo old, 29-37 g], and ghrelin-null littermates (Ghr-/-; 27-34 g). RESULTS: Grazing yielded accurate, consistent feeding events in rats, with an approximately linear rise in nocturnal cumulative food intake [tdF/I (StD): 97.4 ± 1.5% accurate compared with manual measurement; R2 = 0.86; tdF/I (HFD): 99.0 ± 1.4% accurate; R2 = 0.86]. Meal-feeding produced 3 nocturnal meals of equal size and duration in StD-fed rats (tdF/I: 97.4 ± 0.9% accurate; R2 = 0.90), whereas the second meal size increased progressively in HFD-fed rats (44% higher on day 35 than on day 14; P < 0.01). Importantly, cumulative food intake in grazing and meal-fed rats was identical. Similar results were obtained in WT mice except that less restricted grazing induced hyperphagia (compared with meal-fed WT mice; P < 0.05 from day 1). This difference was abolished in Ghr-/- mice, with meal initiation delayed and meal duration enhanced. Neither pattern elevated corticosterone secretion in rats, but meal-feeding aligned ultradian pulses. CONCLUSIONS: We have established a consistent, measurable, researcher-defined, stress-free method for imposing temporal feeding patterns in rats and mice. This approach will facilitate progress in understanding the physiologic impact of feeding patterns.

Pulsatility of glucocorticoid hormones in pregnancy: Changes with gestation and obesity.

OBJECTIVE: Hypothalamic-pituitary-adrenal axis (HPA) activity is decreased in obese pregnancy and associates with increased foetal size. Pulsatile release of glucocorticoid hormones regulates their action in target tissues. Glucocorticoids are essential for normal foetal growth, but little is known about glucocorticoid pulsatility in pregnancy. We aimed to investigate the ultradian rhythm of glucocorticoid secretion during obese and lean pregnancy and nonpregnancy. DESIGN: Serum cortisol, cortisone, corticosterone and 11-dehydrocorticosterone were measured by LC-MS/MS from samples obtained at 10-minute intervals between 08.00-11.00 hours and 16.00-19.00 hours, from 8 lean (BMI <25 kg/m2 ) and 7 obese (BMI > 35 kg/m2 ) pregnant women between 16-24 weeks gestation and again at 30-36 weeks), and nonpregnant controls (lean n = 3, obese n = 4) during the luteal phase of their menstrual cycle. Interstitial fluid cortisol was measured by ELISA, from samples obtained using a portable microdialysis and automated collection device at 20-minute intervals over 24 hours. RESULTS: Serum cortisol AUC, highest peak and lowest trough increased significantly with gestation in lean and obese pregnant compared with nonpregnant subjects. Pulsatility of cortisol was detected in interstitial fluid. In pregnant subjects, interstitial fluid pulse frequency was significantly lower with advancing gestation in obese, but not in lean. CONCLUSIONS: We demonstrate cortisol pulsatility in interstitial fluid. Pulse frequency is altered with increased gestation and BMI. This may be a novel mechanism to explain decreased HPA activity in obese pregnancy.

Dynamic responses of the adrenal steroidogenic regulatory network.

The hypothalamic-pituitary-adrenal axis is a dynamic system regulating glucocorticoid hormone synthesis in the adrenal glands. Many key factors within the adrenal steroidogenic pathway have been identified and studied, but little is known about how these factors function collectively as a dynamic network of interacting components. To investigate this, we developed a mathematical model of the adrenal steroidogenic regulatory network that accounts for key regulatory processes occurring at different timescales. We used our model to predict the time evolution of steroidogenesis in response to physiological adrenocorticotropic hormone (ACTH) perturbations, ranging from basal pulses to larger stress-like stimulations (e.g., inflammatory stress). Testing these predictions experimentally in the rat, our results show that the steroidogenic regulatory network architecture is sufficient to respond to both small and large ACTH perturbations, but coupling this regulatory network with the immune pathway is necessary to explain the dissociated dynamics between ACTH and glucocorticoids observed under conditions of inflammatory stress.

Heterogeneity of Calcium Responses to Secretagogues in Corticotrophs From Male Rats.

Heterogeneity in homotypic cellular responses is an important feature of many biological systems, and it has been shown to be prominent in most anterior pituitary hormonal cell types. In this study, we analyze heterogeneity in the responses to hypothalamic secretagogues in the corticotroph cell population of adult male rats. Using the genetically encoded calcium indicator GCaMP6s, we determined the intracellular calcium responses of these cells to corticotropin-releasing hormone and arginine-vasopressin. Our experiments revealed marked population heterogeneity in the response to these peptides, in terms of amplitude and dynamics of the responses, as well as the sensitivity to different concentrations and duration of stimuli. However, repeated stimuli to the same cell produced remarkably consistent responses, indicating that these are deterministic on a cell-by-cell level. We also describe similar heterogeneity in the sensitivity of cells to inhibition by corticosterone. In summary, our results highlight a large degree of heterogeneity in the cellular mechanisms that govern corticotroph responses to their physiological stimuli; this could provide a mechanism to extend the dynamic range of the responses at the population level to allow adaptation to different physiological challenges.

60 YEARS OF NEUROENDOCRINOLOGY: Glucocorticoid dynamics: insights from mathematical, experimental and clinical studies.

A pulsatile pattern of secretion is a characteristic of many hormonal systems, including the glucocorticoid-producing hypothalamic-pituitary-adrenal (HPA) axis. Despite recent evidence supporting its importance for behavioral, neuroendocrine and transcriptional effects of glucocorticoids, there has been a paucity of information regarding the origin of glucocorticoid pulsatility. In this review we discuss the mechanisms regulating pulsatile dynamics of the HPA axis, and how these dynamics become disrupted in disease. Our recent mathematical, experimental and clinical studies show that glucocorticoid pulsatility can be generated and maintained by dynamic processes at the level of the pituitary-adrenal axis, and that an intra-adrenal negative feedback may contribute to these dynamics.

Dynamic pituitary-adrenal interactions in response to cardiac surgery.

OBJECTIVES: To characterize the dynamics of the pituitary-adrenal interaction during the course of coronary artery bypass grafting both on and off pump. Since our data pointed to a major change in adrenal responsiveness to adrenocorticotropic hormone, we used a reverse translation approach to investigate the molecular mechanisms underlying this change in a rat model of critical illness. CLINICAL STUDIES: Prospective observational study. ANIMAL STUDIES: Controlled experimental study. CLINICAL STUDIES: Cardiac surgery operating rooms and critical care units. ANIMAL STUDIES: University research laboratory. CLINICAL STUDIES: Twenty, male patients. ANIMAL STUDIES: Adult, male Sprague-Dawley rats. CLINICAL STUDIES: Coronary artery bypass graft-both on and off pump. ANIMAL STUDIES: Injection of either lipopolysaccharide or saline (controls) via a jugular vein cannula. CLINICAL STUDIES: Blood samples were taken for 24 hours from placement of the first venous access. Cortisol and adrenocorticotropic hormone were measured every 10 and 60 minutes, respectively, and corticosteroid-binding globulin was measured at the beginning and end of the 24-hour period and at the end of operation. There was an initial rise in both levels of adrenocorticotropic hormone and cortisol to supranormal values at around the end of surgery. Adrenocorticotropic hormone levels then returned toward preoperative values. Ultradian pulsatility of both adrenocorticotropic hormone and cortisol was maintained throughout the perioperative period in all individuals. The sensitivity of the adrenal gland to adrenocorticotropic hormone increased markedly at around 8 hours after surgery maintaining very high levels of cortisol in the face of "basal" levels of adrenocorticotropic hormone. This sensitivity began to return toward preoperative values at the end of the 24-hour sampling period. ANIMAL STUDIES: Adult, male Sprague-Dawley rats were given either lipopolysaccharide or sterile saline via a jugular vein cannula. Hourly blood samples were subsequently collected for adrenocorticotropic hormone and corticosterone measurement. Rats were killed 6 hours after the injection, and the adrenal glands were collected for measurement of steroidogenic acute regulatory protein, steroidogenic factor 1, and dosage-sensitive sex reversal, adrenal hypoplasia critical region, on chromosome X, gene 1 messenger RNAs and protein using real-time quantitative polymerase chain reaction and Western immunoblotting, respectively. Adrenal levels of the adrenocorticotropic hormone receptor (melanocortin type 2 receptor) messenger RNA and its accessory protein (melanocortin type 2 receptor accessory protein) were also measured by real-time quantitative polymerase chain reaction. In response to lipopolysaccharide, rats showed a pattern of adrenocorticotropic hormone and corticosterone that was similar to patients undergoing coronary artery bypass grafting. We were also able to demonstrate increased intra-adrenal corticosterone levels and an increase in steroidogenic acute regulatory protein, steroidogenic factor 1, and melanocortin type 2 receptor accessory protein messenger RNAs and steroidogenic acute regulatory protein, and a reduction in dosage-sensitive sex reversal, adrenal hypoplasia critical region, on chromosome X, gene 1 and melanocortin type 2 receptor messenger RNAs, 6 hours after lipopolysaccharide injection. CONCLUSIONS: Severe inflammatory stimuli activate the hypothalamic-pituitary-adrenal axis resulting in increased steroidogenic activity in the adrenal cortex and an elevation of cortisol levels in the blood. Following coronary artery bypass grafting, there is a massive increase in both adrenocorticotropic hormone and cortisol secretion. Despite a subsequent fall of adrenocorticotropic hormone to basal levels, cortisol remains elevated and coordinated adrenocorticotropic hormone-cortisol pulsatility is maintained. This suggested that there is an increase in adrenal sensitivity to adrenocorticotropic hormone, which we confirmed in our animal model of immune activation of the hypothalamic-pituitary-adrenal axis. Using this model, we were able to show that this increased adrenal sensitivity results from changes in the regulation of both stimulatory and inhibitory intra-adrenal signaling pathways. Increased understanding of the dynamics of normal hypothalamic-pituitary-adrenal responses to major surgery will provide us with a more rational approach to glucocorticoid therapy in critically ill patients.

HPA axis-rhythms.

The hypothalamic-pituitary-adrenal (HPA) axis regulates circulating levels of glucocorticoid hormones, and is the major neuroendocrine system in mammals that provides a rapid response and defense against stress. Under basal (i.e., unstressed) conditions, glucocorticoids are released with a pronounced circadian rhythm, characterized by peak levels of glucocorticoids during the active phase, that is daytime in humans and nighttime in nocturnal animals such as mice and rats. When studied in more detail, it becomes clear that the circadian rhythm of the HPA axis is characterized by a pulsatile release of glucocorticoids from the adrenal gland that results in rapid ultradian oscillations of hormone levels both in the blood and within target tissues, including the brain. In this review, we discuss the regulation of these circadian and ultradian HPA rhythms, how these rhythms change in health and disease, and how they affect the physiology and behavior of the organism.

Constant light disrupts the circadian rhythm of steroidogenic proteins in the rat adrenal gland.

The circadian rhythm of corticosterone (CORT) secretion from the adrenal cortex is regulated by the suprachiasmatic nucleus (SCN), which is entrained to the light-dark cycle. Since the circadian CORT rhythm is associated with circadian expression of the steroidogenic acute regulatory (StAR) protein, we investigated the 24h pattern of hormonal secretion (ACTH and CORT), steroidogenic gene expression (StAR, SF-1, DAX1 and Nurr77) and the expression of genes involved in ACTH signalling (MC2R and MRAP) in rats entrained to a normal light-dark cycle. We found that circadian changes in ACTH and CORT were associated with the circadian expression of all gene targets; with SF-1, Nurr77 and MRAP peaking in the evening, and DAX1 and MC2R peaking in the morning. Since disruption of normal SCN activity by exposure to constant light abolishes the circadian rhythm of CORT in the rat, we also investigated whether the AM-PM variation of our target genes was also disrupted in rats exposed to constant light conditions for 5weeks. We found that the disruption of the AM-PM variation of ACTH and CORT secretion in rats exposed to constant light was accompanied by a loss of AM-PM variation in StAR, SF-1 and DAX1, and a reversed AM-PM variation in Nurr77, MC2R and MRAP. Our data suggest that circadian expression of StAR is regulated by the circadian expression of nuclear receptors and proteins involved in both ACTH signalling and StAR transcription. We propose that ACTH regulates the secretion of CORT via the circadian control of steroidogenic gene pathways that become dysregulated under the influence of constant light.

Ultradian corticosterone secretion is maintained in the absence of circadian cues.

Plasma levels of corticosterone exhibit both circadian and ultradian rhythms. The circadian component of these rhythms is regulated by the suprachiasmatic nucleus (SCN). Our studies investigate the importance of the SCN in regulating ultradian rhythmicity. Two approaches were used to dissociate the hypothalamic-pituitary-adrenal (HPA) axis from normal circadian input in rats: (i) exposure to a constant light (LL) environment and (ii) electrolytic lesioning of the SCN. Blood was sampled using an automated sampling system. As expected, both treatments resulted in a loss of the circadian pattern of corticosterone secretion. Ultradian pulsatile secretion of corticosterone however, was maintained across the 24 h in all animals. Furthermore, the loss of SCN input revealed an underlying relationship between locomotor and HPA activity. In control (LD) rats there was no clear correlation between ultradian locomotor activity and hormone secretion, whereas, in LL rats, episodes of ultradian activity were consistently followed by periods of increased pulsatile hormone secretion. These data clearly demonstrate that the ultradian rhythm of corticosterone secretion is generated through a mechanism independent of the SCN input, supporting recent evidence for a sub-hypothalamic pulse generator.

The origin of glucocorticoid hormone oscillations.

Oscillating levels of adrenal glucocorticoid hormones are essential for optimal gene expression, and for maintaining physiological and behavioural responsiveness to stress. The biological basis for these oscillations is not known, but a neuronal "pulse generator" within the hypothalamus has remained a popular hypothesis. We demonstrate that pulsatile hypothalamic activity is not required for generating ultradian glucocorticoid oscillations. We show that a constant level of corticotrophin-releasing hormone (CRH) can activate a dynamic pituitary-adrenal peripheral network to produce ultradian adrenocorticotrophic hormone and glucocorticoid oscillations with a physiological frequency. This oscillatory response to CRH is dose dependent and becomes disrupted for higher levels of CRH. These data suggest that glucocorticoid oscillations result from a sub-hypothalamic pituitary-adrenal system, which functions as a deterministic peripheral hormone oscillator with a characteristic ultradian frequency. This constitutes a novel mechanism by which the level, rather than the pattern, of CRH determines the dynamics of glucocorticoid hormone secretion.

Characterizing dynamic interactions between ultradian glucocorticoid rhythmicity and acute stress using the phase response curve.

The hypothalamic-pituitary-adrenal (HPA) axis is a dynamic oscillatory hormone signalling system that regulates the pulsatile secretion of glucocorticoids from the adrenal glands. In addition to regulation of basal levels of glucocorticoids, the HPA axis provides a rapid hormonal response to stress that is vitally important for homeostasis. Recently it has become clear that glucocorticoid pulses encode an important biological signal that regulates receptor signalling both in the central nervous system and in peripheral tissues. It is therefore important to understand how stressful stimuli disrupt the pulsatile dynamics of this system. Using a computational model that incorporates the crucial feed-forward and feedback components of the axis, we provide novel insight into experimental observations that the size of the stress-induced hormonal response is critically dependent on the timing of the stress. Further, we employ the theory of Phase Response Curves to show that an acute stressor acts as a phase-resetting mechanism for the ultradian rhythm of glucocorticoid secretion. Using our model, we demonstrate that the magnitude of an acute stress is a critical factor in determining whether the system resets via a Type 1 or Type 0 mechanism. By fitting our model to our in vivo stress-response data, we show that the glucocorticoid response to an acute noise stress in rats is governed by a Type 0 phase-resetting curve. Our results provide additional evidence for the concept of a deterministic sub-hypothalamic oscillator regulating the ultradian glucocorticoid rhythm, which constitutes a highly responsive peripheral hormone system that interacts dynamically with hypothalamic inputs to regulate the overall hormonal response to stress.

Origin of ultradian pulsatility in the hypothalamic-pituitary-adrenal axis.

The hypothalamic-pituitary-adrenal (HPA) axis is a neuroendocrine system that regulates the circulating levels of vital glucocorticoid hormones. The activity of the HPA axis is characterized not only by a classic circadian rhythm, but also by an ultradian pattern of discrete pulsatile release of glucocorticoids. A number of psychiatric and metabolic diseases are associated with changes in glucocorticoid pulsatility, and it is now clear that glucocorticoid responsive genes respond to these rapid fluctuations in a biologically meaningful way. Theoretical modelling has enabled us to identify and explore potential mechanisms underlying the ultradian activity in this axis, which to date have not been identified successfully. We demonstrate that the combination of delay with feed-forward and feedback loops in the pituitary-adrenal system is sufficient to give rise to ultradian pulsatility in the absence of an ultradian source from a supra-pituitary site. Moreover, our model enables us to predict the different patterns of glucocorticoid release mediated by changes in hypophysial-portal corticotrophin-releasing hormone levels, with results that parallel our experimental in vivo data.