Heart rate variability biofeedback as a method for assessing baroreflex function: a preliminary study of resonance in the cardiovascular system

This study describes the use of a biofeedback method for the noninvasive study of baroreflex mechanisms. Five previously untrained healthy male participants learned to control oscillations in heart rate using biofeedback training to modify their heart rate variability at specific frequencies. They were instructed to match computer-generated sinusoidal oscillations with oscillations in heart rate at seven frequencies within the range of 0.01–0.14 Hz. All participants successfully produced high-amplitude target-frequency oscillations in both heart rate and blood pressure. Stable and predictable transfer functions between heart rate and blood pressure were obtained in all participants. The highest oscillation amplitudes were produced in the range of 0.055–0.11 Hz for heart rate and 0.02–0.055 Hz for blood pressure. Transfer functions were calculated among sinusoidal oscillations in the target stimuli, heart rate, blood pressure, and respiration for frequencies at which subjects received training. High and low target-frequency oscillation amplitudes at specific frequencies could be explained by resonance among various oscillatory processes in the cardiovascular system. The exact resonant frequencies differed among individuals. Changes in heart rate oscillations could not be completely explained by changes in breathing. The biofeedback method also allowed us to quantity characteristics of inertia, delay, and speed sensitivity in baroreflex system. We discuss the implications of these findings for using heart rate variability biofeedback as an aid in diagnosing various autonomic and cardiovascular system disorders and as a method for treating these disorders.     

I. Cytochrome oscillations in continuous culture: Action of chloramphenicol on Saccharomyces carlsbergensis and Candida utilis

Oscillations in the A, B and C group cytochromes of Saccharomyces carsbergensis and Candida utilis have been observed in continuous culture following the addition and removal of chloramphenicol. With S. carlsbergensis reproduction of these oscillations proved difficult due to the fermentative nature of the yeast. With C. utilis at D = 0.1 h⁻¹ the oscillations in teh A group cytochromes were statistically significant and 3–6 times the standard deviation observed in steady sate. The oscillations provide evidence of a feedback loop controlling cytochrome synthesis and have a period of the same order of that predicted by Goodwin for the protein (enzyme) component of his model (viz. 6–8 h). With C. utilis the period of these oscillations was a function of the growth rate varying from approx. 10 h at a growth rate of 0.1 h⁻¹ (mean generation time 6.9 h) to 3 h at a growth rate of 0.3 h⁻¹ (mean generation time 2.3 h).     

A G-Protein Subunit Translocation Embedded Network Motif Underlies GPCR Regulation of Calcium Oscillations

G-protein βγ subunits translocate reversibly from the plasma membrane to internal membranes on receptor activation. Translocation rates differ depending on the γ subunit type. There is limited understanding of the role of the differential rates of G_βγ translocation in modulating signaling dynamics in a cell. Bifurcation analysis of the calcium oscillatory network structure predicts that the translocation rate of a signaling protein can regulate the damping of system oscillation. Here, we examined whether the G_βγ translocation rate regulates calcium oscillations induced by G-protein-coupled receptor activation. Oscillations in HeLa cells expressing γ subunit types with different translocation rates were imaged and quantitated. The results show that differential G_βγ translocation rates can underlie the diversity in damping characteristics of calcium oscillations among cells. Mathematical modeling shows that a translocation embedded motif regulates damping of G-protein-mediated calcium oscillations consistent with experimental data. The current study indicates that such a motif may act as a tuning mechanism to design oscillations with varying damping patterns by using intracellular translocation of a signaling component.     

UNDAMPED OSCILLATIONS OF PYRIDINE NUCLEOTIDE AND OXYGEN TENSION IN CHEMOSTAT CULTURES OF KLEBSIELLA AEROGENES

Klebsiella aerogenes was grown in chemostat culture with the pH controlled to ±0.01 and temperature to ±0.1°C. The oxygen tension of the culture was regulated by changing the partial pressure of oxygen in the gas phase and recorded by means of an oxygen electrode. Reduced pyridine nucleotide was monitored continuously in the culture by means of direct fluorimetry. On applying an anaerobic shock to the culture, damped oscillations in pyridine nucleotide fluorescence were obtained. Further anaerobic shocks decreased the damping and eventually gave rise to undamped oscillations of a 2–3 min period which continued for several days. These oscillations were paralleled by oscillations of the same frequency in respiration rate. The amplitude of the oscillations in the respiration rate was equivalent to only 1% of the total steady-state respiration, whereas that of pyridine nucleotide oscillations was equivalent to 10% of the total aerobic/anaerobic fluorescence response. The oscillations ceased on interrupting the glucose feed but restarted on adding excess glucose to the culture. Addition of succinate also restarted the oscillations so that they appear not to be of glycolytic origin. The frequency of oscillations varied with growth rate and conditions. Oscillations of much lower frequency were obtained under limited-oxygen and anaerobic conditions than under fully aerobic conditions. Under glucose-limited conditions, fluctuations were found in adenosine triphosphate (ATP) content which were in phase with the pyridine nucleotide oscillations, but under nitrogen-limited growth conditions no such fluctuations in ATP were observed. The primary oscillating pathway could not be identified but the mechanism would appear to be quite different from that involved in oscillations observed in yeast cells. The synchronization of oscillations and observations of negative damping could be explained by a syntalysis effect.     

Estimating population trends using population viability analyses for the conservation ofCapra pyrenaica

Large herbivore populations can suffer important oscillations with considerable effects on ecosystem functions and services, yet our capacity to predict population fate is limited and conditional upon the availability of data. This study investigated the interannual variation in the growth rate of populations ofCapra pyrenaica Schinz, 1838, and its extinction risk by comparing the dynamics of populations that were stable for more than two decades (Gredos and Tortosa-Beceite), populations that had increased recently (Tejeda-Almijara), and populations that were in decline (Cazorla-Segura) or extinct (the Pyrenees population; hereafter, bucardo). To estimate quasi-extinction threshold assessments (50% of population extinct in this study), which have implications for the conservation of the species, we used empirical data and the predictions derived from several theoretical models. The results indicate that when variance of log population growth rate reaches a specific threshold, the probability of quasi-extinction increased drastically. ForC. pyrenaica, we recommend keeping population variance < 0.05, which will reduce the likelihood that the irruptive oscillations caused by environmental and demographic stochasticity will put the population at risk. Models to predict the dynamics ofC. pyrenaica populations should incorporate temporal stochasticity because, in this study, it strongly increased the likelihood that a population declined.     

Interplay of Intrinsic and Synaptic Conductances in the Generation of High-Frequency Oscillations in Interneuronal Networks with Irregular Spiking

High-frequency oscillations (above 30 Hz) have been observed in sensory and higher-order brain areas, and are believed to constitute a general hallmark of functional neuronal activation. Fast inhibition in interneuronal networks has been suggested as a general mechanism for the generation of high-frequency oscillations. Certain classes of interneurons exhibit subthreshold oscillations, but the effect of this intrinsic neuronal property on the population rhythm is not completely understood. We study the influence of intrinsic damped subthreshold oscillations in the emergence of collective high-frequency oscillations, and elucidate the dynamical mechanisms that underlie this phenomenon. We simulate neuronal networks composed of either Integrate-and-Fire (IF) or Generalized Integrate-and-Fire (GIF) neurons. The IF model displays purely passive subthreshold dynamics, while the GIF model exhibits subthreshold damped oscillations. Individual neurons receive inhibitory synaptic currents mediated by spiking activity in their neighbors as well as noisy synaptic bombardment, and fire irregularly at a lower rate than population frequency. We identify three factors that affect the influence of single-neuron properties on synchronization mediated by inhibition: i) the firing rate response to the noisy background input, ii) the membrane potential distribution, and iii) the shape of Inhibitory Post-Synaptic Potentials (IPSPs). For hyperpolarizing inhibition, the GIF IPSP profile (factor iii)) exhibits post-inhibitory rebound, which induces a coherent spike-mediated depolarization across cells that greatly facilitates synchronous oscillations. This effect dominates the network dynamics, hence GIF networks display stronger oscillations than IF networks. However, the restorative current in the GIF neuron lowers firing rates and narrows the membrane potential distribution (factors i) and ii), respectively), which tend to decrease synchrony. If inhibition is shunting instead of hyperpolarizing, post-inhibitory rebound is not elicited and factors i) and ii) dominate, yielding lower synchrony in GIF networks than in IF networks.     

A G-Protein Subunit Translocation Embedded Network Motif Underlies GPCR Regulation of Calcium Oscillations

G-protein βγ subunits translocate reversibly from the plasma membrane to internal membranes on receptor activation. Translocation rates differ depending on the γ subunit type. There is limited understanding of the role of the differential rates of G_βγ translocation in modulating signaling dynamics in a cell. Bifurcation analysis of the calcium oscillatory network structure predicts that the translocation rate of a signaling protein can regulate the damping of system oscillation. Here, we examined whether the G_βγ translocation rate regulates calcium oscillations induced by G-protein-coupled receptor activation. Oscillations in HeLa cells expressing γ subunit types with different translocation rates were imaged and quantitated. The results show that differential G_βγ translocation rates can underlie the diversity in damping characteristics of calcium oscillations among cells. Mathematical modeling shows that a translocation embedded motif regulates damping of G-protein-mediated calcium oscillations consistent with experimental data. The current study indicates that such a motif may act as a tuning mechanism to design oscillations with varying damping patterns by using intracellular translocation of a signaling component.     

Growth dynamics of Agrobacterium tumefaciens in chemostat cultures limited by carbon source and mineral nutrients

Agrobacterium tumefaciens was grown in a chemostat in a chemically-defined medium which had α-methyl d-glucoside, magnesium, manganese, phosphate or urea as the growth-limiting nutrient. Steady-state biomass concentrations were dependent on the specific growth rate of the organism when α-methyl d-glucoside, manganese or phosphate were growth-limiting nutrients. During magnesium-limited growth, large undamped oscillations in biomass concentration occurred. In all chemostat cultures a variant organism was selected which had lost the ability to grow in the medium supplied, but survived on products of carbon metabolism derived from the wild-type.     

Modeling light adaptation in circadian clock: prediction of the response that stabilizes entrainment

Periods of biological clocks are close to but often different from the rotation period of the earth. Thus, the clocks of organisms must be adjusted to synchronize with day-night cycles. The primary signal that adjusts the clocks is light. In Neurospora, light transiently up-regulates the expression of specific clock genes. This molecular response to light is called light adaptation. Does light adaptation occur in other organisms? Using published experimental data, we first estimated the time course of the up-regulation rate of gene expression by light. Intriguingly, the estimated up-regulation rate was transient during light period in mice as well as Neurospora. Next, we constructed a computational model to consider how light adaptation had an effect on the entrainment of circadian oscillation to 24-h light-dark cycles. We found that cellular oscillations are more likely to be destabilized without light adaption especially when light intensity is very high. From the present results, we predict that the instability of circadian oscillations under 24-h light-dark cycles can be experimentally observed if light adaptation is altered. We conclude that the functional consequence of light adaptation is to increase the adjustability to 24-h light-dark cycles and then adapt to fluctuating environments in nature.     

Oscillations in a size-structured prey-predator model

This article introduces a predator–prey model with the prey structured by body size, based on reports in the literature that predation rates are prey-size specific. The model is built on the foundation of the one-species physiologically structured models studied earlier. Three types of equilibria are found: extinction, multiple prey-only equilibria and possibly multiple predator–prey coexistence equilibria. The stabilities of the equilibria are investigated. Comparison is made with the underlying ODE Lotka–Volterra model. It turns out that the ODE model can exhibit sustain oscillations if there is an Allee effect in the net reproduction rate, that is the net reproduction rate grows for some range of the prey’s population size. In contrast, it is shown that the structured PDE model can exhibit sustain oscillations even if the net reproductive rate is strictly declining with prey population size. We find that predation, even size-non-specific linear predation can destabilize a stable prey-only equilibrium, if reproduction is size specific and limited to individuals of large enough size. Furthermore, we show that size-specific predation can also destabilize the predator–prey equilibrium in the PDE model. We surmise that size-specific predation allows for temporary prey escape which is responsible for destabilization in the predator–prey dynamics.     

Physiology of spontaneous [Ca2+]i oscillations in the isolated vasopressin and oxytocin neurones of the rat supraoptic nucleus

The magnocellular vasopressin (AVP) and oxytocin (OT) neurones exhibit specific electrophysiological behaviour, synthesise AVP and OT peptides and secrete them into the neurohypophysial system in response to various physiological stimulations. The activity of these neurones is regulated by the very same peptides released either somato-dendritically or when applied to supraoptic nucleus (SON) preparations in vitro. The AVP and OT, secreted somato-dendritically (i.e. in the SON proper) act through specific autoreceptors, induce distinct Ca²⁺ signals and regulate cellular events. Here, we demonstrate that about 70% of freshly isolated individual SON neurones from the adult non-transgenic or transgenic rats bearing AVP (AVP-eGFP) or OT (OT-mRFP1) markers, produce distinct spontaneous [Ca²⁺]_i oscillations. In the neurones identified (through specific fluorescence), about 80% of AVP neurones and about 60% of OT neurones exhibited these oscillations. Exposure to AVP triggered [Ca²⁺]_i oscillations in silent AVP neurones, or modified the oscillatory pattern in spontaneously active cells. Hyper- and hypo-osmotic stimuli (325 or 275 mOsmol/l) respectively intensified or inhibited spontaneous [Ca²⁺]_i dynamics. In rats dehydrated for 3 or 5 days almost 90% of neurones displayed spontaneous [Ca²⁺]_i oscillations. More than 80% of OT-mRFP1 neurones from 3 to 6-day-lactating rats were oscillatory vs. about 44% (OT-mRFP1 neurones) in virgins. Together, these results unveil for the first time that both AVP and OT neurones maintain, via Ca²⁺ signals, their remarkable intrinsic in vivo physiological properties in an isolated condition.     

Role of Spontaneous Current Oscillations during High-Efficiency Electrotransformation of Thermophilic Anaerobes

Current oscillations at about 24 MHz were observed during electrotransformation (ET) of the thermophilic anaerobes Clostridium thermocellum ATCC 27405, C. thermocellum DSM 1313, and Thermoanaerobacterium saccharolyticum YS 485, using a pulse gated by a square signal generated by a custom generator. In experiments in which only the field strength was varied, all three of these strains resulted in a one-to-one correspondence between the appearance of current oscillations and successful ET. Oscillations accompanied ET of both C. thermocellum strains only at field strengths of ≥12 kV/cm, and ET was only observed above the same threshold. Similarly, for T. saccharolyticum, oscillations were only observed at field strengths of ≥10 kV/cm, and ET was only observed above the same threshold. When a passive electrical filter consisting of an inductor and resistor in parallel was added to the system to prevent the development of oscillations, ET efficiencies were reduced dramatically for all three strains at all field strengths tested. The maximum tested field strength, 25 kV/cm, resulted in the maximum measured transformation efficiency for all three strains. At this field strength, the efficiency of ET in the absence of oscillations was decreased compared to that observed in the presence of oscillations by 500-fold for C. thermocellum ATCC 27405, 2,500-fold for C. thermocellum DSM 1313, and 280-fold for T. saccharolyticum. Controls using the same apparatus with Escherichia coli cells or a resistor with a value representative of the direct current resistance of typical cell samples did not develop oscillations, and ET efficiencies obtained with E. coli were the same with or without the electrical filter included in the pulse generator circuit. The results are interpreted to indicate that spontaneously arising oscillations have a large beneficial effect on transformation efficiency in the system employed here and that the development of oscillations in this system is affected by the cell species present.     

Dynamic oscillatory circuit of regulation of capillary hemodynamics

We used laser Doppler flowmetry with wavelet analysis of blood flow oscillations, computer capillaroscopy, and thermometry of the nail bed in 30 subjects to show an important role of the oscillatory circuit in the regulation of capillary hemodynamics, number of functioning capillaries, and linear and volumetric velocity of blood flow. The number of functioning capillaries is regulated by oscillations of myogenic and sensory peptidergic origin. The appearance of sensory oscillations, especially high-amplitude oscillations, is an adaptive neurotrophic mechanism that significantly increases the number of functioning capillaries and intensity of blood flow from arterioles to capillaries. The linear velocity of blood flow depends on both the tone of microvessels and changes in the dynamic component of blood pressure. Under conditions of skin hypoperfusion, the mean linear velocity of capillary blood flow may be inversely related to the extracapillary perfusion, including the amplitude of heart rate (A _h) and oscillations of the tone of precapillary sphincters, whereas under conditions of vasodilation and increased skin perfusion, it may be inversely related to the amplitude of arteriolar oscillations of endothelial or neurogenic sympathetic origin (A _maxe + n) and the shunting index. The A _h affects the linear velocity of blood flow in the arterial part of capillaries, whereas the A _maxe + n influences the same factor in the venous part. The contribution of oscillations to the regulation of the linear velocity varies depending on the perfusion and skin temperature. The resultant tone of distributing microvessels is determined by the competition between the stationary and oscillatory components. In addition to changes in the amplitude, the frequency of vasomotions may also be important. The regulatory importance of the oscillatory circuit is increased with a decrease in the skin blood flow.     

Integration,coincidence detection and resonance in networks of spiking neurons expressing Gamma oscillations and asynchronous states

Gamma oscillations are widely seen in the awake and sleeping cerebral cortex, but the exact role of these oscillations is still debated. Here, we used biophysical models to examine how Gamma oscillations may participate to the processing of afferent stimuli. We constructed conductance-based network models of Gamma oscillations, based on different cell types found in cerebral cortex. The models were adjusted to extracellular unit recordings in humans, where Gamma oscillations always coexist with the asynchronous firing mode. We considered three different mechanisms to generate Gamma, first a mechanism based on the interaction between pyramidal neurons and interneurons (PING), second a mechanism in which Gamma is generated by interneuron networks (ING) and third, a mechanism which relies on Gamma oscillations generated by pacemaker chattering neurons (CHING). We find that all three mechanisms generate features consistent with human recordings, but that the ING mechanism is most consistent with the firing rate change inside Gamma bursts seen in the human data. We next evaluated the responsiveness and resonant properties of these networks, contrasting Gamma oscillations with the asynchronous mode. We find that for both slowly-varying stimuli and precisely-timed stimuli, the responsiveness is generally lower during Gamma compared to asynchronous states, while resonant properties are similar around the Gamma band. We could not find conditions where Gamma oscillations were more responsive. We therefore predict that asynchronous states provide the highest responsiveness to external stimuli, while Gamma oscillations tend to overall diminish responsiveness.     

Elongation and circumnutation oscillations of hypocotyl of pine seedlings (Pinus silvestris L.)

The elongation of pine seedlings (Pinus silvestris L.) is associated with nutation movements. Trajectories of these growth oscillations were recorded by film technique in horizontal and vertical projection during a three day period of growth. On the basis of these data the parameters of elongation and nutation oscillations,i.e. rate, amplitudes and frequency of oscillations, were calculated and their changes during plant development compared. Oscillation trajectories are circular or elliptic spirals the amplitudes of which are increasing with the age of hypocotyl from 1.5 mm to 7.5 mm. The frequencies of nutations are decreasing during the growth from 0.5 to 0.2 rev. h^?1. On the other hand, the growth rate of hypocotyl increased from values near to 10^?3 mm h^?1 at the beginning of the experiment to 4×10^?1 mm h^?1 recorded at the end of the third day. The zone of nutation curvature was slightly transferred from the middle of the hypocotyl toward the apex and its location has not been identical with that of elongation. This indicates that the system controlling nutation oscillations need not be identical with that controlling direction of elongation. At a certain stage of development behaviour of the decapitated pine hypocotyl is analogical to that of the root without the centre of georeception. A possibility of analogy of the system controlling direction of hypocotyl growth and of the system proposed for geotropical control of root growth is discussed.     

Self-organization and dynamics of an open futile cycle

The dynamic behaviour of an open futile cycle composed of two enzymes has been investigated in the vicinity of a steady-state. A necessary condition required for damped or sustained oscillations of the system is that enzyme E₂, which controls recycling of the substrate S₂, be inhibited by an excess of this substrate. In order for the system to be neutrally stable and therefore to exhibit sustained oscillations, it is not necessary for antagonist enzyme E₁ to be activated by its product S₂. If it is enzyme E₁ which is inhibited by an excess of its substrate S₁, the system has a saddle point. Other conditions for stability or instability of the system have been determined. If the enzyme E₁, which is not inhibited by the substrate, exhibits a slow conformational transition of the mnemonical type, this transition dramatically alters the stability behavior of the system. If the mnemonical enzyme E₁ were exhibiting a positive kinetic co-operativity, decreasing the rate of the conformational transition of the mnemonical enzyme will increase the stability of the whole system and will tend to damp the oscillations in the vicinity of the steady-state. If conversely the mnemonical enzyme E₁ were exhibiting a negative kinetic co-operativity, decreasing the rate of the enzyme conformational transition will decrease the stability of the system and will tend to create or amplify oscillations of the system taken as a whole. If these results may be extended to more complex metabolic cycles, involving more than two enzymes, it may be tentatively considered that positive co-operativity associated with slow transition has emerged in the course of evolution in order to limit temporal instabilities of metabolic cycles. Alternatively one may speculate that the “biological function” of negative co-operativity is to create or amplify these temporal instabilities.     

Towards nano-physiology of insects with atomic force microscopy

Little study of insects with modern nanotechnology tools has been done so far. Here we use one of such tool, atomic force microscopy (AFM) to study surface oscillations of the ladybird beetles (Hippodamia convergens) measured in different parts of the insect at picometer level. This allows us to record a much broader spectral range of possible surface vibrations (up to several kHz) than the previously studied oscillations due to breathing, heartbeat cycles, coelopulses, etc. (up to 5-10 Hz). Here we demonstrate three different ways with which one can identify the origins of the observed peaks - by physical positioning the probe near a specific organ, and by using biological or chemical stimuli. We report on identification of high frequency peaks associated with H. convergens heart, spiracular closer muscles, and oscillations associated with muscles activated while drinking. The method, being a relatively non-invasive technique providing a new type of information, may be useful in developing “nanophysiology” of insects.     

Simulation of Spontaneous Ca Oscillations in Astrocytes Mediated by Voltage-Gated Calcium Channels

The purpose of this computational study was to investigate the possible role of voltage-gated Ca²⁺ channels in spontaneous Ca²⁺ oscillations of astrocytes. By incorporating different types of voltage-gated Ca²⁺ channels and a previous model, this study reproduced typical Ca²⁺ oscillations in silico. Our model could mimic the oscillatory phenomenon under a wide range of experimental conditions, including resting membrane potential (−75 to −60 mV), extracellular Ca²⁺ concentration (0.1 to 1500 μM), temperature (20 to 37°C), and blocking specific Ca²⁺ channels. By varying the experimental conditions, the amplitude and duration of Ca²⁺ oscillations changed slightly (both <25%), while the frequency changed significantly (∼400%). This indicates that spontaneous Ca²⁺ oscillations in astrocytes might be an all-or-none process, which might be frequency-encoded in signaling. Moreover, the properties of Ca²⁺ oscillations were found to be related to the dynamics of Ca²⁺ influx, and not only to a constant influx. Therefore, calcium channels dynamics should be used in studying Ca²⁺ oscillations. This work provides a platform to explore the still unclear mechanism of spontaneous Ca²⁺ oscillations in astrocytes.     

Simple oscillations in photosynthesis of higher plants

Illumination of leaves of Vicia faba L. provoked oscillations in the rates of CO₂ uptake and O₂ evolution. The oscillations were marked under anaerobic conditions, but were absent at 20% O₂. The minimum CO₂ concentration required for the appearance of oscillations was 600 μl · l^?1. The higher the CO₂ concentration, the stronger the oscillations. The effect of CO₂ concentration was saturated at 1000 μl CO₂ · l^?1. The period of the oscillations was 5–6 min at a light intensity of 80 nE · cm^?2 · s^?1 and became longer on lowering of the intensity. No oscillations appeared at intensities below 12 nE · cm^?2 · s^?1. Oscillations could also be generated by increasing the CO₂ concentration in the atmosphere during strong illumination under anaerobic conditions. The chlorophyl a fluorescence yield showed oscillations, similar in shape and frequency to those of photosynthesis, after such an environmental change. Oscillations were also observed in photosynthesis of other C₃ plants, Lycopersicon esulentum Mill and Glycine max Merrill, under the same conditions as those required for V. faba, but were absent for the C₄ plants, Zea mays and Amaranthus retroflexus L.     

Up-regulation of circadian clock gene Period 2 in the prostate mesenchymal cells during flutamide-induced apoptosis

Androgen regulates the proper development and physiological function of the prostate. Here, we investigated the modulation of androgen and androgen receptor (AR) antagonist on circadian oscillations of a clock core gene Period 2 (Per2) in rat prostate mesenchymal cells (PMCs). Circadian oscillations were analyzed with the real-time monitoring system of gene expression using transgenic rats introduced with mouse Per2 promoter fused to a destabilized luciferase (Per2-dLuc) reporter gene. Analyses of circadian oscillations, immunofluorescence, and androgen response element (ARE)-luciferase reporter assay revealed that circadian clocks are operative and the AR protein is functional in PMCs in vitro. Androgen such as testosterone (T) and dihydrotestosterone (DHT) did not cause any changes in circadian Per2-dLuc oscillations of confluent cells. Conversely, flutamide (FL) up-regulated the amplitude of circadian Per2-dLuc oscillations in a dose-dependent manner, whereas T antagonized the action of FL. The PER2 protein was markedly accumulated by FL treatment and localized in both the nucleus and cytoplasm during the first peak period of circadian Per2-dLuc oscillations. Simultaneously, FL treatment increased apoptotic cell death. Collectively, the present study demonstrates that a clock gene Per2 is up-regulated in PMCs during FL-induced apoptotic cell death. Thus, circadian oscillations of Per2 gene expression may be closely linked to the cellular states of PMCs such as apoptotic cell death.