Transition of phase response properties and singularity in the circadian limit cycle of cultured cells

The circadian system has been regarded as a limit cycle oscillator constructed by the integrated interaction of clock genes and proteins. Here, we investigated a mammalian circadian oscillation geometrically before and after a perturbation. We detected the singular point and transition from a type 1 to type 0 phase response curve (PRC) and determined the embedding dimension to show how many variables are needed to describe the limit cycle oscillation and relaxation process after a perturbation. As a perturbation, forskolin (FK) was administered to Rat-1 cells expressing the Per2::luc gene. By broadly and finely changing the phase and strength of the perturbation, we detected the transition of the PRC from type 1 to type 0 and a possible singular transition point, the property of which agreed quite well with our numerical simulation of the noisy Goodwin model, a simple yet canonical model for the transcription-translation feedback loop of the core clock genes. Furthermore, we estimated the embedding dimension of the limit cycle before and after the perturbation. The trajectory of the limit cycle was embedded in two dimensions but with the perturbation of the state point moved out of the trajectory, the relaxation process was generally embedded in higher dimensions. The average number of embedding dimensions at each dose of FK increased as the FK dose increased but most of the relaxation process was generally embedded within four dimensions. These findings support the existence of a circadian limit cycle oscillator in mammalian cells and suggest that a small number of variables determine the relaxation process after a perturbation.