Why circadian rhythms are circadian: competitive population dynamics of biological oscillators

Network Structure of the Master Clock Is Important for Its Primary Function

A master clock located in the suprachiasmatic nucleus (SCN) regulates the circadian rhythm of physiological and behavioral activities in mammals. The SCN has two main functions in the regulation: an endogenous clock produces the endogenous rhythmic signal in body rhythms, and a calibrator synchronizes the body rhythms to the external light-dark cycle. These two functions have been determined to depend on either the dynamic behaviors of individual neurons or the whole SCN neuronal network. In this review, we first introduce possible network structures for the SCN, as revealed by time series analysis from real experimental data. It was found that the SCN network is heterogeneous and sparse, that is, the average shortest path length is very short, some nodes are hubs with large node degrees but most nodes have small node degrees, and the average node degree of the network is small. Secondly, the effects of the SCN network structure on the SCN function are reviewed based on mathematical models of the SCN network. It was found that robust rhythms with large amplitudes, a high synchronization between SCN neurons and a large entrainment ability exists mainly in small-world and scale-free type networks, but not other types. We conclude that the SCN most probably is an efficient small-world type or scale-free type network, which drives SCN function.

Translational thresholds in a core circadian clock model

Organisms have evolved in a daily cyclic environment, developing circadian cell-autonomous clocks that temporally organize a wide range of biological processes. Translation is a highly regulated process mainly associated with the activity of microRNAs (miRNAs) at the translation initiation step that impacts on the molecular circadian clock dynamics. Recently, a molecular titration mechanism was proposed to explain the interactions between some miRNAs and their target mRNAs; new evidence also indicates that regulation by miRNA is a nonlinear process such that there is a threshold level of target mRNA below which protein production is drastically repressed. These observations led us to use a theoretical model of the circadian molecular clock to study the effect of miRNA-mediated translational thresholds on the molecular clock dynamics. We model the translational threshold by introducing a phenomenological Hill equation for the kinetics of PER translation and show how the parameters associated with translation kinetics affect the period, amplitude, and time delays between clock mRNA and clock protein expression. We show that our results are useful for analyzing experiments related to the translational regulation of negative elements of transcriptional-translational feedback loops. We also provide new elements for thinking about the translational threshold as a mechanism that favors the emergence of circadian rhythmicity, the tuning of the period-delay relationship and the cell capacity to control the protein oscillation amplitude with almost negligible changes in the mRNA amplitudes.

Heterogeneity in relaxation rate improves the synchronization of oscillatory neurons in a model of the SCN

The circadian rhythms in mammals, that are regulated by the suprachiasmatic nucleus (SCN) of the brain, have been observed even in the absence of a light-dark cycle. The SCN is composed of about 10 000 autonomous neuronal oscillators, which are heterogenous in many oscillatory properties, including the heterogeneity in relaxation rates. Although the relaxation rate affects the entrainability of the SCN as a whole, not much is known about the reasons why the heterogeneity in relaxation rate exists. In the present study, based on a Poincaré model, we examine whether the heterogeneity in the relaxation rate affects the synchronization of the SCN neuronal oscillators under constant darkness. Both our simulations and theoretical results show that the heterogeneity improves the synchronization. Our findings provide an alternative explanation for the existence of the heterogeneity in the SCN neurons and shed light on the effect of neuronal heterogeneity on the collective behavior of the SCN neurons.

Phenotypic plasticity in the invasive pest Drosophila suzukii: activity rhythms and gene expression in response to temperature

Phenotypic plasticity may contribute to the invasive success of an alien species in a new environment. A higher plastic species may survive and reproduce in more diverse environments, thereby supporting establishment and colonization. We focused on plasticity in the circadian rhythm of activity, which can favour species coexistence in invasion, for the invasive species Drosophila suzukii, which is expected to be a weaker direct competitor than other Drosophila species of the resident community. We compared between the invasive D. suzukii and the resident D. melanogaster the circadian rhythms of the locomotor activity in adults and the expression of clock genes in response to temperature. We showed that D. suzukii is active in a narrower range of temperatures than D. melanogaster and that the activities of both species overlap during the day, regardless of the temperature. Both species are diurnal and exhibit rhythmic activity at dawn and dusk, with a much lower activity at dawn for D. suzukii females. Our results showed that the timeless and clock genes are good candidates to explain the plastic response that is observed in relation to temperature. Overall, our results suggest that thermal phenotypic plasticity in D. suzukii activity is not sufficient to explain the invasive success of D. suzukii and calls for testing other hypotheses, such as the release of competitors and/or predators.

Dispersion of the intrinsic neuronal periods affects the relationship of the entrainment range to the coupling strength in the suprachiasmatic nucleus

Living beings on the Earth are subjected to and entrained (synchronized) to the natural 24-h light-dark cycle. Interestingly, they can also be entrained to an external artificial cycle of non-24-h periods. The range of these periods is called the entrainment range and it differs among species. In mammals, the entrainment range is regulated by a main clock located in the suprachiasmatic nucleus (SCN) which is composed of 10 000 neurons in the brain. Previous works have found that the entrainment range depends on the cellular coupling strength in the SCN. In particular, the entrainment range decreases with the increase of the cellular coupling strength, provided that all the neuronal oscillators are identical. However, the SCN neurons differ in the intrinsic periods that follow a normal distribution in a range from 22 to 28 h. In the present study, taking the dispersion of the intrinsic neuronal periods into account, we examined the relationship between the entrainment range and the coupling strength. Results from numerical simulations and theoretical analyses both show that the relationship is altered to be paraboliclike if the intrinsic neuronal periods are nonidentical, and the maximal entrainment range is obtained with a suitable coupling strength. Our results shed light on the role of the cellular coupling in the entrainment ability of the SCN network.

Differences in intrinsic amplitudes of neuronal oscillators improve synchronization in the suprachiasmatic nucleus

In mammals, a main clock located in the suprachiasmatic nucleus (SCN) regulates the ∼24 h rhythms of behavioral and physiological activities exposed to a natural 24 light-dark cycle or even under constant darkness. The rhythms originate from self-sustained oscillations of the SCN neurons, which differ in both intrinsic periods and intrinsic amplitudes. The intrinsic periods and the intrinsic amplitudes were found to be bound to specific regions in the previous experiments. In particular, neurons of smaller amplitudes and larger periods are located in a ventrolateral part, and neurons of larger amplitudes and smaller periods are in a dorsomedial part. In the present study, we examined the effects of the differences in the intrinsic frequencies and the differences in the intrinsic amplitudes of neuronal oscillators on the synchronization, respectively. We found that the differences in the intrinsic frequencies weaken the synchronization, whereas the differences in the intrinsic amplitudes strengthen the synchronization. Our finding may shed light on the effects of the heterogenous properties of individual neurons on the collective behaviors of the SCN network and provide a way to enhance the synchronization.

The asymmetry of the entrainment range induced by the difference in intrinsic frequencies between two subgroups within the suprachiasmatic nucleus

The rhythms of physiological and behavioral activities in mammals, which are regulated by the main clock suprachiasmatic nucleus (SCN) in the brain, can not be only synchronized to the natural 24 h light-dark cycle, but also to cycles with artificial periods. The range of the artificial periods that the animal can be synchronized to is called entrainment range. In the absence of the light-dark cycle, the animal can also maintain the circadian rhythm with an endogenous period close to 24 h. Experiments found that the entrainment range is not symmetrical with respect to the endogenous period. In the present study, an explanation is given for the asymmetry based on a Kuramoto model which describes the neuronal network of the SCN. Our numerical simulations and theoretical analysis show that the asymmetry results from the difference in the intrinsic frequencies between two subgroups of the SCN, as well as the entrainment range is affected by the difference.

Dissociation between two subgroups of the suprachiasmatic nucleus affected by the number of damped oscillated neurons

In mammals, the main clock located in the suprachiasmatic nucleus (SCN) of the brain synchronizes the body rhythms to the environmental light-dark cycle. The SCN is composed of about 2×10^{4} neurons which can be classified into three oscillatory phenotypes: self-sustained oscillators, damped oscillators, and arrhythmic neurons. Exposed to an artificial external light-dark cycle with a period of 22h instead of 24h, two subgroups of the SCN can become desynchronized (dissociated). The ventrolateral (VL) subgroup receives photic input and is entrained to the external cycle and a dorsomedial (DM) subgroup oscillates with its endogenous (i.e., free running) period and is synchronized to the external light-dark cycle through coupling from the VL. In the present study, we examined the effects of damped oscillatory neurons on the dissociation between VL and DM under an external 22h cycle. We found that, with increasing numbers of damped oscillatory neurons located in the VL, the dissociation between the VL and DM emerges, but if these neurons are increasingly present in the DM the dissociation disappears. Hence, the damped oscillatory neurons in different subregions of the SCN play distinct roles in the dissociation between the two subregions of the SCN. This shows that synchrony between SCN subregions is affected by the number of damped oscillatory neurons and the location of these cells. We suggest that more knowledge on the number and the location of these cells may explain why some species do show a dissociation between the subregions and others do not, as the distribution of oscillatory types of neurons offers a plausible and novel candidate mechanism to explain heterogeneity.

Entrainment range of the suprachiasmatic nucleus affected by the difference in the neuronal amplitudes between the light-sensitive and light-insensitive regions

Mammals not only can be synchronized to the natural 24-h light-dark cycle, but also to a cycle with a non-24-h period. The range of the period of the external cycle, for which the animals can be entrained to, is called the entrainment range, which differs among species. The entrainment range as a characteristic of the animal is determined by the main circadian clock, i.e., the suprachiasmatic nucleus (SCN) in the brain. The SCN is composed of ∼10000 heterogeneous neurons, which can be divided into two subgroups, i.e., the ventrolateral subgroup (VL) directly receiving the light information from the retina and relaying the information to the dorsomedial subgroup (DM). Among the SCN neurons, the amplitudes are different; however, it is unclear that the amplitude is related to the location of the neurons in experiments. In the present study, we examined the effect of the difference in the neuronal amplitude between the VL and the DM on the entrainment range of the SCN, based on a mathematical model, i.e., the Poincaré model, which is used to describe the circadian clock. We find that the maximal entrainment range is obtained when the difference is equal to a critical point. If the difference of the amplitudes of the VL neurons to the amplitudes of the DM neurons is smaller than a critical point, with the increase of the difference, the entrainment range of the SCN increases, while if the difference is larger than the critical point, the entrainment range decreases with the increase of the difference. Our finding may give a potential explanation for the diversity of the entrainment range among species.

The effects of non-self-sustained oscillators on the en-trainment ability of the suprachiasmatic nucleus

In mammals, the circadian rhythms of behavioral and physiological activities are regulated by an endogenous clock located in the suprachiasmatic nucleus (SCN). The SCN is composed of ~20,000 neurons, of which some are capable of self-sustained oscillations, while the others do not oscillate in a self-sustainable manner, but show arrhythmic patterns or damped oscillations. Thus far, the effects of these non-self-sustained oscillatory neurons are not fully explored. Here, we examined how the proportion of the non-self-sustained oscillators affects the free running period under constant darkness and the ability to entrain to the light-dark cycle. We find that the proportion does not affect the free running period, but plays a significant role in the range of entrainment. We also find that its effect on the entrainment range depends on the region where the non-self-sustained oscillators are located. If the non-self-sustained oscillatory neurons are situated in the light-sensitive subregion, the entrainment range narrows when the proportion increases. If they are situated in the light-insensitive subregion, however, the entrainment range broadens with the increase of the proportion. We suggest that the heterogeneity within the light-sensitive and light-insensitive subregions of the SCN has important consequences for how the clock works.

The circadian rhythm induced by the heterogeneous network structure of the suprachiasmatic nucleus

In mammals, the master clock is located in the suprachiasmatic nucleus (SCN), which is composed of about 20 000 nonidentical neuronal oscillators expressing different intrinsic periods. These neurons are coupled through neurotransmitters to form a network consisting of two subgroups, i.e., a ventrolateral (VL) subgroup and a dorsomedial (DM) subgroup. The VL contains about 25% SCN neurons that receive photic input from the retina, and the DM comprises the remaining 75% SCN neurons which are coupled to the VL. The synapses from the VL to the DM are evidently denser than that from the DM to the VL, in which the VL dominates the DM. Therefore, the SCN is a heterogeneous network where the neurons of the VL are linked with a large number of SCN neurons. In the present study, we mimicked the SCN network based on Goodwin model considering four types of networks including an all-to-all network, a Newman-Watts (NW) small world network, an Erdös-Rényi (ER) random network, and a Barabási-Albert (BA) scale free network. We found that the circadian rhythm was induced in the BA, ER, and NW networks, while the circadian rhythm was absent in the all-to-all network with weak cellular coupling, where the amplitude of the circadian rhythm is largest in the BA network which is most heterogeneous in the network structure. Our finding provides an alternative explanation for the induction or enhancement of circadian rhythm by the heterogeneity of the network structure.

Impact of dispersed coupling strength on the free running periods of circadian rhythms

The dominant endogenous clock, named the suprachiasmatic nucleus (SCN), regulates circadian rhythms of behavioral and physiological activity in mammals. One of the main characteristics of the SCN is that the animal maintains a circadian rhythm with a period close to 24 h in the absence of a daily light-dark cycle (called the free running period). The free running period varies among species due to heterogeneity of the SCN network. Previous studies have shown that the heterogeneity in cellular coupling as well as in intrinsic neuronal periods shortens the free running period. Furthermore, as derived from experiments, one neuron's coupling strength is negatively associated with its period. It is unknown what the effects of this association between coupling strength and period are on the free running period and how the heterogeneity in coupling strength influences this free running period. In the present study we found that in the presence of a negative relationship between one neuron's coupling strength and its period, surprisingly, the dispersion of coupling strengths increases the free running period. Our present finding may shed new light on the understanding of the heterogeneous SCN network and provides an alternative explanation for the diversity of free running periods between species.

Heterogeneity induces rhythms of weakly coupled circadian neurons

The main clock located in the suprachiasmatic nucleus (SCN) regulates circadian rhythms in mammals. The SCN is composed of approximately twenty thousand heterogeneous self-oscillating neurons, that have intrinsic periods varying from 22 h to 28 h. They are coupled through neurotransmitters and neuropeptides to form a network and output a uniform periodic rhythm. Previous studies found that the heterogeneity of the neurons leads to attenuation of the circadian rhythm with strong cellular coupling. In the present study, we investigate the heterogeneity of the neurons and of the network in the condition of constant darkness. Interestingly, we found that the heterogeneity of weakly coupled neurons enables them to oscillate and strengthen the circadian rhythm. In addition, we found that the period of the SCN network increases with the increase of the degree of heterogeneity. As the network heterogeneity does not change the dynamics of the rhythm, our study shows that the heterogeneity of the neurons is vitally important for rhythm generation in weakly coupled systems, such as the SCN, and it provides a new method to strengthen the circadian rhythm, as well as an alternative explanation for differences in free running periods between species in the absence of the daily cycle.

Generation of slow phase-locked oscillation and variability of the interspike intervals in globally coupled neuronal oscillators

To elucidate how a biological rhythm is regulated, the extended (three-dimensional) Bonhoeffer-van der Pol or FitzHugh-Nagumo equations are employed to investigate the dynamics of a population of neuronal oscillators globally coupled through a common buffer (mean field). Interesting phenomena, such as extraordinarily slow phase-locked oscillations (compared to the natural period of each neuronal oscillator) and the death of all oscillations, are observed. We demonstrate that the slow synchronization is due mainly to the existence of "fast" oscillators. Additionally, we examine the effect of noise on the synchronization and variability of the interspike intervals. Peculiar phenomena, such as noise-induced acceleration and deceleration, are observed. The results herein suggest that very small noise may significantly influence a biological rhythm.

Entrainment range of nonidentical circadian oscillators by a light-dark cycle

The suprachiasmatic nucleus (SCN) is a principal circadian clock in mammals, which controls physiological and behavioral daily rhythms. The SCN has two main features: Maintaining a rhythmic cycle of approximately 24 h in the absence of a light-dark cycle (free-running period) and the ability to entrain to external light-dark cycles. Both free-running period and range of entrainment vary from one species to another. To understand this phenomenon, we investigated the diversity of a free-running period by the distribution of coupling strengths in our previous work [Phys. Rev. E 80, 030904(R) (2009)]. In this paper we numerically found that the dispersion of intrinsic periods among SCN neurons influence the entrainment range of the SCN, but has little influence on the free-running periods under constant darkness. This indicates that the dispersion of coupling strengths determines the diversity in free-running periods, while the dispersion of intrinsic periods determines the diversity in the entrainment range. A theoretical analysis based on two coupled neurons is presented to explain the results of numerical simulations.

Entrainment of the suprachiasmatic nucleus network by a light-dark cycle

The synchronization of biological activity with the alternation of day and night (circadian rhythm) is performed in the brain by a group of neurons, constituting the suprachiasmatic nucleus (SCN). The SCN is divided into two subgroups of oscillating cells: the ventrolateral (VL) neurons, which are exposed to light (photic signal), and the dorsomedial (DM) neurons, which are coupled to the VL cells. When the coupling between these neurons is strong enough, the system synchronizes with the photic period. Upon increasing the cell coupling, the entrainment of the DM cells has been recently shown to occur via a very sharp (jumping) transition when the period of the photic input is larger than the intrinsic period of the cells. Here, we characterize this transition with a simple realistic model. We show that two bifurcations possibly lead to the disappearance of the endogenous mode. Using a mean-field model, we show that the jumping transition results from a supercritical Hopf-like bifurcation. This finding implies that both the period and strength of the stimulating photic signal, and the relative fraction of cells in the VL and DM compartments, are crucial in determining the synchronization of the system.

Mechanism of phase splitting in two coupled groups of suprachiasmatic-nucleus neurons

The phase-splitting behavior of coupled suprachiasmatic-nucleus neurons has been observed in many mammals, and its mechanism is still not completely understood. Based on our previous work [C. Gu, J. Wang, and Z. Liu, Phys. Rev. E 80, 030904(R) (2009)] on the free-running periods of neurons in the suprachiasmatic nucleus, we present here a modified Goodwin oscillator model to explain the mechanism of phase splitting. In contrast to the previous phase model, the modified Goodwin oscillator model contains the information on both the phase and amplitude and, thus, can show more features than the purely phase model, including all three behaviors of synchronization, phase splitting, and amplitude death and the distributed periodicity in the regions of synchronization and phase splitting, etc. An analytic phase model is extracted from the modified Goodwin oscillator model to explain the dependence of periodicity on the parameters. Moreover, both the modified Goodwin oscillator model and the analytic phase model show that the ensemble frequency can be enhanced or reduced by the time delay.

Free-running period of neurons in the suprachiasmatic nucleus: Its dependence on the distribution of neuronal coupling strengths

The suprachiasmatic nucleus (SCN) pacemaker shows a free-running period ranging from 20 to 28 h for different species, which was usually explained from the angle of coupling strength. Based on the assumption of nonidentical coupling strengths in SCN, we find an alternative mechanism that the diversity of free-running period can be also caused by the distribution of coupling strengths. The free-running period is proportional to the average coupling strength and inverse proportional to the dispersion of couplings. Moreover, we present an analytic phase model to confirm our finding, which shows a solid foundation of our finding and opens a window to study the collective behaviors of SCN oscillators.

Basic Science Review on Circadian Rhythm Biology and Circadian Sleep Disorders

The sleep-wake cycle displays a characteristic 24-hour periodicity, providing an opportunity to dissect the endogenous circadian clock through the study of aberrant behaviour. This article surveys the properties of circadian clocks, with emphasis on mammals. Information was obtained from searches of peer-reviewed literature in the PUBMED database. Features that are highlighted include the known molecular components of clocks, their entrainment by external time cues and the output pathways used by clocks to regulate metabolism and behaviour. A review of human circadian rhythm sleep disorders follows, including recent discoveries of their genetic basis. The article concludes with a discussion of future approaches to the study of human circadian biology and sleep-wake behaviour. Key words: Circadian clocks, Entrainment, Human circadian sleep-wake disorders

Synchronization of excitatory neurons with strongly heterogeneous phase responses

In many real-world oscillator systems, the phase response curves are highly heterogeneous. However, the dynamics of heterogeneous oscillator networks has not been seriously addressed. We propose a theoretical framework to analyze such a system by dealing explicitly with the heterogeneous phase response curves. We develop a method to solve the self-consistent equations for order parameters by using formal complex-valued phase variables, and apply our theory to networks of in vitro cortical neurons. We find a novel state transition that is not observed in previous oscillator network models.

Inhibitory neurons can facilitate rhythmic activity in a neural network

Synchrony emerges in a population of oscillators interacting through in-phase couplings. We addressed a question of whether inhibitory neurons simply hinder the emergence of the synchronous activity among excitatory neurons, or facilitate it. An analysis of a simple phase model revealed that both cases may take place. Numerical simulations of the more realistic models revealed that inhibitory neurons rather facilitate rhythmic activity.

Amplitude of circadian oscillations entrained by 24-h light-dark cycles

An intriguing property of circadian clocks is that their free-running period is not exactly 24h. Using models for circadian rhythms in Neurospora and Drosophila, we determine how the entrainment of these rhythms is affected by the free-running period and by the amplitude of the external light-dark cycle. We first consider the model for Neurospora, in which light acts by inducing the expression of a clock gene. We show that the amplitude of the oscillations of the clock protein entrained by light-dark cycles is maximized when the free-running period is smaller than 24h. Moreover, if the amplitude of the light-dark cycle is very strong, complex oscillations occur when the free-running period is close to 24h. In the model for circadian rhythms in Drosophila, light acts by enhancing the degradation of a clock protein. We show that while the amplitude of circadian oscillations entrained by light-dark cycles is also maximized if the free-running period is smaller than 24h, the range of entrainment is centered around 24h in this model. We discuss the physiological relevance of these results in regard to the setting of the free-running period of the circadian clock.

Engaging space: extraterrestrial architecture and the human psyche

The human fascination with exploring and inhabiting the space that lies beyond Earth's atmosphere continues to grow. Nevertheless, 40 years of experience to date have clearly established that humans in outer space routinely suffer significant psychological impairment arising from their stressful extraterrestrial living conditions. This paper explores those extraterrestrial conditions through the interactions between the extraordinarily harsh environment of outer space, the sensations that humans encounter in space, and the qualities of a habitat that physically interposes itself between the two. The objective of this paper is to develop a habitat that expresses the extraterrestrial condition while supporting the mental health of its inhabitants, so as to augment the success of prolonged extraterrestrial residence and interplanetary travel.

Population dynamics of clock-controlled biological species: models and why circadian rhythms are circadian

Conventional models of population dynamics of biological species do not take into account the fact that activity of most living organisms is under control of biological clocks, which are oscillators with periods near 1 day, 1 year, and so on. Here proposed are new categories of population dynamics models which are combinations of Lotka-Volterra type of equations (as a simplest example; in fact, an arbitrary type may be used) and equations of the biological oscillators. Such models may be particularly useful in addressing to those problems which are inaccessible with conventional models. A simple example of such models with distributed clock periods is proposed and studied numerically and analytically to consider why periods of biological rhythms are not precisely one day or one year, etc., as is often questioned. For this purpose, a notion of survival index is introduced to measure the degree of success in survival. It turns out that clock periods equal or close to that of the milieu are not necessarily advantageous for survival and can even lead to extinction; survival is then most successful at the edge of entrainment. This may provide a clue to answer the question. It is also found that survival is difficult for those species which cannot be entrained by the environmental cycle.

Quantification of heart rate variability by discrete nonstationary non-Markov stochastic processes

We develop the statistical theory of discrete nonstationary non-Markov random processes in complex systems. The objective of this paper is to find the chain of finite-difference non-Markov kinetic equations for time correlation functions (TCF) in terms of nonstationary effects. The developed theory starts from careful analysis of time correlation through nonstationary dynamics of vectors of initial and final states and nonstationary normalized TCF. Using the projection operators technique we find the chain of finite-difference non-Markov kinetic equations for discrete nonstationary TCF and for the set of nonstationary discrete memory functions (MF's). The last one contains supplementary information about nonstationary properties of the complex system on the whole. Another relevant result of our theory is the construction of the set of dynamic parameters of nonstationarity, which contains some information of the nonstationarity effects. The full set of dynamic, spectral and kinetic parameters, and kinetic functions (TCF, short MF's statistical spectra of non-Markovity parameter, and statistical spectra of nonstationarity parameter) has made it possible to acquire the in-depth information about discreteness, non-Markov effects, long-range memory, and nonstationarity of the underlying processes. The developed theory is applied to analyze the long-time (Holter) series of RR intervals of human ECG's. We had two groups of patients: the healthy ones and the patients after myocardial infarction. In both groups we observed effects of fractality, standard and restricted self-organized criticality, and also a certain specific arrangement of spectral lines. The received results demonstrate that the power spectra of all orders (n=1,2, ...) MF m(n)(t) exhibit the neatly expressed fractal features. We have found out that the full sets of non-Markov, discrete and nonstationary parameters can serve as reliable and powerful means of diagnosis of the cardiovascular system states and can be used to distinguish healthy data from pathologic data.