Indian Journal of PsychiatryIndian Journal of Psychiatry
Home | About us | Current Issue | Archives | Submission | Instructions | Subscribe | Advertise | Contact | Login 
    Users online: 19 Small font sizeDefault font sizeIncrease font size Print this article Email this article Bookmark this page
 


 

 
     
    Advanced search
 

 
 
     
  
    Similar in PUBMED
   Search Pubmed for
   Search in Google Scholar for
 Related articles
    Email Alert *
    Add to My List *
* Registration required (free)  


    Abstract
    Introduction
    Patients and Methods
    Results
    Discussion
    Conclusion
    References
    Article Figures
    Article Tables

 Article Access Statistics
    Viewed2953    
    Printed95    
    Emailed3    
    PDF Downloaded192    
    Comments [Add]    
    Cited by others 2    

Recommend this journal

 


 
ORIGINAL ARTICLE Table of Contents   
Year : 2010  |  Volume : 52  |  Issue : 1  |  Page : 42-47
The influence of negative mood on heart rate complexity measures and baroreflex sensitivity in healthy subjects


1 Department of Psychiatry, Friedrich-Schiller-University, Jena, Germany,
2 Medical Engineering and Biotechnology, University of Applied Sciences, Jena, Germany,
3 Medical School, University of Alberta, Edmonton, Canada,
4 Rajiv Gandhi University of Health Sciences, India,
5 Psychiatry and Behavioral Neurosciences, Wayne State University School of Medicine, Detroit, USA and University of Alberta, Edmonton, Canada,

Click here for correspondence address and email

Date of Web Publication13-Jan-2010
 

   Abstract 

Background: Decreased cardiac vagal function is linked with increased cardiac mortality and depression is associated with decreased heart rate variability. We have previously shown that the Mood Induction Procedure (MIP) in healthy subjects alters pain perception and thalamic activity during pain perception.
Aim: To study the effect of negative emotion on heart rate variability and complexity measures as well as on baroreceptor sensitivity, as these parameters reflect cardiac autonomic function.
Patients and Methods: We studied 20 healthy female controls before and after neutral MIP and 20 healthy female subjects before and after negative MIP. We investigated measures of valence of mood, heart rate variability and complexity and the baroreceptor sensitivity index.
Results: While there was a significant difference in the valence of mood between the neutral and the negative effect condition, there were no significant differences in any of the heart rate or baroreceptor sensitivity measures between the two groups.
Conclusions: Our findings did not show any significant influence of acute negative MIP on heart rate variability and complexity measures and baroreceptor sensitivity, even though depressive disorder and stress are associated with decreased heart rate variability. These findings are discussed in the context of clinical depression and anxiety and the increased risk for cardiac mortality. In contrast to the presented results here, we have previously shown that MIP in healthy subjects alters pain perception and thalamic activity.

Keywords: Autonomic nervous system, heart rate variability, major depressive disorder, parasympathetic nervous system, psychophysiology, psychosis, risk of cardiac mortality, vagal

How to cite this article:
Kobele R, Koschke M, Schulz S, Wagner G, Yeragani S, Ramachandraiah CT, Voss A, Yeragani VK, Bar KJ. The influence of negative mood on heart rate complexity measures and baroreflex sensitivity in healthy subjects. Indian J Psychiatry 2010;52:42-7

How to cite this URL:
Kobele R, Koschke M, Schulz S, Wagner G, Yeragani S, Ramachandraiah CT, Voss A, Yeragani VK, Bar KJ. The influence of negative mood on heart rate complexity measures and baroreflex sensitivity in healthy subjects. Indian J Psychiatry [serial online] 2010 [cited 2014 Sep 22];52:42-7. Available from: http://www.indianjpsychiatry.org/text.asp?2010/52/1/42/58894



   Introduction Top


Emotions have diverse effects on autonomic nervous function as illustrated by symptoms such as palpitations and hyperventilation. [1] It is well documented that emotional processes result in changes in heart rate (HR), heart rate variability (HRV) and contractility. Heart rate variability and complexity are useful noninvasive techniques to assess the vagal component of cardiac autonomic function in response to emotional changes. Dishman and colleagues [2] have shown reduced cardiac vagal function among men and women who perceived more stress in the week before the assessment. Moreover, reducing the perceived stress by attending to an emotionally self-management training program increases vagal tone. [3] Similarly, anxiety is often accompanied by somatic manifestations that suggest marked changes in autonomic nervous system (ANS) activity. Reduced HRV has been demonstrated in subjects with various anxiety disorders including patients suffering from panic disorder. [4] This is also true of patients suffering from depression. [5] This is especially important due to the high incidence of sudden unexplained cardiac death in patients with anxiety as well as depression. [6],[7] The underlying mechanisms for cardiac vulnerability in depression remain unclear. An altered autonomic neuro-cardiac regulation might be one important pathophysiological factor. [8],[9] However, conflicting results were reported in studies that used HRV [10],[11],[12],[13] and complexity measures. [5] Baroreflex sensitivity (BRS) and blood pressure variability have also been shown to be decreased in patients with depression, [14],[15] thus suggesting decreased parasympathetic function in the disease.

Negative affect is an important symptom of anxiety as well as depression. To elucidate the influence of mood on cardiac autonomic function, we assessed parameters of HRV and complexity in healthy subjects before and after negative mood induction and in a control condition.


   Patients and Methods Top


Subjects

We investigated 40 female students in this study. Twenty students (age: 24.5 ± 5.1 years) participated in the negative mood induction protocol and the others (age: 25.7± 4.3 years) in the control condition. All of them were completely healthy (clinical examination and investigation, ECG, routine blood chemistry; BMI < 26 kg/ m 2 ), non-smokers and none was receiving medication that would confound the results of this study. The scores on the Hamilton Depression Rating Scale (HAM-D; [16] ) and Beck's Depression Inventory (BDI, [17] ) were unremarkable.

This study complies with the declaration of Helsinki. All participants gave written informed consent to a protocol approved by the local Ethics Committee of the Medical Faculty of the Friedrich-Schiller University, Jena.

Mood induction

To induce sad and neutral emotional states, the modified Velten Mood Induction Procedure (VMIP) was used. The VMIP is amongst the most widely used techniques for studying effective influences upon cognition and behavior, [18] and it has demonstrated effectiveness in altering subjective emotional states. [19] During the VMIP, participants were exposed to a series of 21 self-referent sad-mood statements, which were presented twice and the subjects had to read the statements aloud. While reading the statements, participants were asked to attempt to experience the mood suggested by the statements (e.g. 'Life is a heavy burden'). Additionally, to facilitate the sad MIP, participants were presented with individually tailored music during VMIP. [20] The procedure is effective and has been used by our group in previous studies. [21]

In contrast, during the neutral procedure, participants were exposed to a series of 21 neutral statements (e.g. 'An orange is a citrus fruit'), which was similarly presented as the sad statements. Mozart's Piano concerto No. 21 in C Major was chosen for all subjects as neutral music. The whole mood- induction procedure lasted approximately 12 min.

To assess effective changes during the experiment, participants were asked to rate the amount of VMIP on the dimensions of valence and arousal using the Self-Assessment Manikin (SAM), an effective rating system devised by Lang. [22] In this system, ratings of valence are indicated by five graphical representations of facial expressions ranging from a severe frown (most sad = 2 4) to a broad smile (most positive = 14). For arousal, the Manikin varies from a state of low to high agitation (9 represents a high rating and 1 represents a low rating). The assessment of mood was rated before the mood-induction (test-1), five minutes after induction (test-2) and at the end of the experiment (test-3).

Data acquisition and preprocessing

Examinations were performed between 3 and 6 PM and subjects were asked to relax. Respiratory rate was obtained for all patients. The electrocardiogram (high resolution, 1000 Hz) and continuous blood pressure was recorded (CNSystems® , Medizintechnik GmbH, Austria).

Heart rate variability

We computed measures of heart rate variability (HRV) in the time and frequency domain. [23] In particular, we obtained the standard deviation of the averages of N-N (normal-to-normal beat) intervals (SDNN, in milliseconds), the square root of the mean squared differences of successive NN intervals (RMSSD). Furthermore, the quotient of low- frequency (LF: 0.04-0.15 Hz) and high-frequency (HF: 0.15-0.4 Hz) powers (LF/HF) was calculated as a measure of relative sympathovagal balance.

Compression entropy

This measure calculates the ratio of an original time-series length to a compressed version. An approach to describe the entropy of a text was introduced in the framework of algorithmic information theory. [24] Here, the entropy (complexity) of a given text is defined as the smallest algorithm that is capable of generating the text. In this study, we applied the LZ77 algorithm for loss-less data-compression introduced by Lempel and Ziv. [25] This algorithm is widely used and implemented in many file compressors such as Winzip® . Its application in R-R time-series has been described in detail elsewhere. [9],[26] The ratio of the compressed to the original time-series length represents an index of entropy and is referred to as compression entropy (Hc).

Symbolic dynamics for high and low variability

For analysis of symbolic dynamics, every single heart beat was compared to the preceding beat. Whenever the beat differed within a specific time limit, this was encoded with the letter '0'. In contrast, when beats were different exceeding the given time limits (10 ms), the letter '1' was attributed. [27] From these data, letter sequences were analyzed. For low variability parameters (plvar), the occurrence of sequences containing six consecutive '0s were assessed, whereas for high variability parameters (phvar) six consecutive '1's were relevant. The time limit applied is always indicated by adding the respective number to the parameter assessed, e.g. plvar10 for low variability with a time limit of 10 ms.

Baroreflex sensitivity (b-slope; t-slope)

The baroreflex sensitivity (BRS) was assessed using the sequence method. [28] A detailed description has been published previously. [29] In brief, spontaneous sequences of at least three consecutive beats were analysed, when an increased systolic blood pressure (SBP) of at least 1 mmHg caused an increased BBI (beat-to-beat interval) of at least 5 ms (bradycardic sequence; b-slope) or a decreased SBP caused a decreased BBI (tachycardic sequence; t-slope). Furthermore, the numbers of bradycardic and tachycardic baroreflex sequences (b-count; t-count) were computed as indicators of baroreflex activation.

Nonlinear joint symbolic dynamics (JSDsym; JSDdiam)

To assess heart rate and blood pressure dynamics in a more complex way, an analysis based on joint symbolic dynamics (JSD) was applied, [30] which has been described in detail previously. [31] Here, the beat-to-beat changes of R-R interval and SBP were each coarse-grained to two different symbols: Increasing values are coded as '1,' whereas decreasing and unchanged values were coded as '0,' respectively. Subsequently, short patterns of symbol sequences (words) were formed, and their distribution properties are analysed [probability of symmetric baroreflex-related words (JSDsym), probability of diametric non-baroreflex-related words (JSDdiam), etc.]. Based on the considerations mentioned above, words of a three-letter length are feasible for short- term recordings. [Figure 1]

Statistical analysis

For statistical analyses, SPSS for Windows (version 14.0) was used. Independent t-tests were used to examine differences of dimension of valence and arousal between groups at each point in time (t1, t2, t3).

To compare different effects of both mood-induction procedures on autonomic function, a multivariate analysis of variance for repeated (t1 vs. t2) measures (MANOVA) was performed. In particular, HR and the natural logarithm of RMSSD, HC, sdNN, LF/HF, systolic BP, phvar10, diastolic BP, t-slope, b-slope, t-count, b-count, JSDsys and JSDdiam were investigated. Follow-up univariate ANOVAs were performed for each parameter. We used Pearson's correlations to explore the bivariate relationship between the questionnaire scores and each physiological parameters (t1, t2, t3). Since groups were matched according to age and sex, no covariates were used.


   Results Top


Inferential statistics

There were highly significant differences in mood-valence between the neutral and negative group [Figure 1]a after mood-induction (t2) as shown by an independent t-test. (t = 9.03, df = 39, P > 0.001). As described in [Figure 1]a, the average dimension of valence was significantly decreased after induction in the negative effect group only, and subjects with sad mood induction reported a significantly decreased sad mood after VMIP. For the dimension of valence before induction and at the end of the test (t1, t3), we found no significant differences between the two groups. We have not found any significant differences for the dimension of arousal.

Multivariate statistics

Repeated measures MANOVA (t1 vs. t2) for autonomic parameters revealed no significant effect of TIME (Wilks' Lambda = 0.46; F(13,11) = 0.98; P = 0.52). The analysis also revealed that there was no multivariate difference between both groups (Wilks' Lambda = 0.45, F(13,11) = 0.77, P = 0.68) or the interaction of TIME 3 GROUP (Wilks' Lambda = 0.52, F(11,24) = 0.18, P = 0.1).

Univariate statistics

Follow-up univariate ANOVAs showed no significant main effects for any variable of the autonomic nervous system (for mean and standard deviation of values, see [Table 1] and [Figure 1]b-d.

Correlations

There were no significant correlations between the effect scores and the HRV parameters.


   Discussion Top


Emotions involve a complex mix of cognitive, affective, behavioral and physiological responses. [32] The influence of negative emotions on autonomic cardiac function is important for psychiatric research due to the association of anxiety and depression with increased cardiac mortality. [8],[33] Sad mood is a negative affect and in this study, we investigated changes after induction of sad mood in healthy female students and assessed the dimensions of valence and arousal. Significant differences in the dimension of valence are not mirrored in autonomic function, although we performed an analysis applying linear as well as nonlinear parameters of HRV. Besides linear HRV parameters describing the variance of heart beats, nonlinear complexity parameters have been developed to describe the regularity of HR time-series. The application of these novel analyses has led to a higher sensitivity for detecting autonomic dysfunction [5],[34] as well as patients at risk for sudden death in different diseases. [35] These results indicate that negative mood per se does not change cardiac autonomic function significantly in healthy subjects as assessed here in the acute condition. [36] This result has important consequences and adds to an ongoing debate on how autonomic function is changed in depression [5] as it might differ between acute versus chronic dysphoric conditions. In addition, decreased BRS has been reported in depression. [14] We did not find any significant changes in BRS during negative VMIP neither by means of linear nor nonlinear parameters for BRS. This probably goes along with our findings of no significant changes in HRV in this study.

One important difference between negative affect in normal controls and patients with clinical depression is one of immediate mood change versus sustained. Our study might indicate that autonomic changes due to negative mood in healthy control subjects are not directly comparable with the depressed mood state in major depression. [5] It is also possible that repeated bouts of negative affect may decrease cardiac vagal function, and thus our acute induction of negative effect is a limitation in our study. The physiological state in depressed patients might be changed due to additional factors, such as hormonal changes in the pituitary gland or even a lack of physical fitness, which probably occurs gradually over the course of time.

On the other hand, duration of the disease and the influence of the disease on other biological systems such as the HPA (hypothalmo-pituitary-adrenal) axis might contribute to autonomic dysfunction. [37] The duration of the disease might be influential since long-term adaptation might occur in depressed patients. The degree of physical fitness is an important issue, as changes of autonomic function in depression might potentially reflect a low degree of physical fitness in such patients. This is of great relevance for cardiac function as well. It is unlikely that the amount of mood induction was not sufficient in our experiment, since we were able to demonstrate significant differences for the dimension of valence and we have found in a previous study differences for pain perception and thalamic activity during pain perception applying the same method. [21] However, the degree of arousal was not different between the neutral and negative condition. It is likely that the degree of arousal is more important than the actual dimension of valence. This degree of arousal should be investigated in future studies and should be applied in psychopathological rating scales used in clinical studies.


   Conclusion Top


In contrast to our previous findings in depressed patients, [5],[11] the induced sad mood in healthy subjects in this study revealed opposite findings. Experimentally, induced sad mood in healthy controls highlighted the difference between clinically relevant depressed mood and 'normal' variation of sad mood.

 
   References Top

1.Sakuragi S, Sugiyama Y, Takeuchi K. Effects of laughing and weeping on mood and heart rate variability. J Physiol Anthropol Appl Human Sci 2002;21:159-65.  Back to cited text no. 1      
2.Dishman RK, Nakamura Y, Garcia ME, Thompson RW, Dunn AL, Blair SN. Heart rate variability, trait anxiety, and perceived stress among physically fit men and women. Int J Psychophysiol 2000;37:121-33.  Back to cited text no. 2      
3.McCraty R, Atkinson M, Tiller WA, Rein G, Watkins AD. The effects of emotions on short-term power spectrum analysis of heart rate variability. Am J Cardiol 1995;76:1089-93.  Back to cited text no. 3      
4.Yeragani VK, Rao R. Effect of nortriptyline and paroxetine on measures of chaos of heart rate time series in patients with panic disorder. J Psychosom Res 2003;55:507-13.  Back to cited text no. 4      
5.Boettger S, Hoyer D, Falkenhahn K, Kaatz M, Yeragani VK, Bär KJ. Nonlinear broad band dynamics are less complex in major depression. Bipol Disorders 2008;10:276-84.  Back to cited text no. 5      
6.Zellweger MJ, Osterwalder RH, Langewitz W, Pfisterer ME. Coronary artery disease and depression. Eur Heart J 2004;25:3-9.  Back to cited text no. 6      
7.Ariyo AA, Haan M, Tangen CM, Rutledge JC, Cushman M, Dobs A. Depressive symptoms and risks of coronary heart disease and mortality in elderly Americans: Cardiovascular Health Study Collaborative Research Group. Circulation 2000;102:1773-9.  Back to cited text no. 7      
8.Carney RM, Freedland KE, Miller GE, Jaffe AS. Depression as a risk factor for cardiac mortality and morbidity: A review of potential mechanisms. J Psychosom Res 2002;53:897-902.  Back to cited text no. 8      
9.Voss A, Baier V, Schulz S, Bär K. Linear and nonlinear methods for analyses of cardiovascular variability in bipolar disorders. Bipolar Disord 2006;8:441-52.  Back to cited text no. 9      
10.Guinjoan SM, Bernabo JL, Cardinali DP. Cardiovascular tests of autonomic function and sympathetic skin responses in patients with major depression. J Neurol Neurosurg Psychiatry 1995;59:299-302.  Back to cited text no. 10      
11.Rechlin T, Weis M, Spitzer A, Kaschka WP. Are affective disorders associated with alterations of heart rate variability? J Affect Disord 1994;32:271-5.  Back to cited text no. 11      
12.Bär KJ, Greiner W, Jochum T, Friedrich M, Wagner G, Sauer H. The influence of major depression and its treatment on heart rate variability and pupillary light reflex parameters. J Affect Disord 2004;82:245-52.  Back to cited text no. 12      
13.Moser M, Lehofer M, Hoehn-Saric R, McLeod DR, Hildebrandt G, Steinbrenner B. Increased heart rate in depressed subjects in spite of unchanged autonomic balance? J Affect Disord 1998;48:115-24.  Back to cited text no. 13      
14.Broadley AJ, Frenneaux MP, Moskvina V, Jones CJ, Korszun A. Baroreflex sensitivity is reduced in depression. Psychosom Med 2005;67:648-51.  Back to cited text no. 14      
15.Tulen JH, Bruijn JA, de Man KJ, Pepplinkhuizen L, van den Meiracker AH, Man int Veld AJ. Cardiovascular variability in major depressive disorder and effects of imipramine or mirtazapine. J Clin Psychopharmacol 1996;16:135-45.  Back to cited text no. 15      
16.Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry 1960;23:56-62.  Back to cited text no. 16      
17.Beck AT, Ward CH, Mendelson M, Mock J, Erbaugh J. An inventory for measuring depression. Arch Gen Psychiatry 1961;4:561-71.  Back to cited text no. 17      
18.Velten E Jr. A laboratory task for induction of mood states. Behav Res Ther 1968;6:473-82.  Back to cited text no. 18      
19.Frost RO, Teghtsoonian R. Use of psychophysical scaling procedures in the assessment of phobias: A case example. J Behav Ther Exp Psychiatry 1982;13:191-3.  Back to cited text no. 19      
20.Sutherland G, Newman B, Rachman S. Experimental investigations of the relations between mood and intrusive unwanted cognitions. Br J Med Psychol 1982;55:127-38.  Back to cited text no. 20      
21.Wagner G, Koschke M, Leuf T, Schlösser R, Bär KJ. Reduced heat pain thresholds after sad-mood induction are associated with changes in thalamic activity. Neuropsychologia 2009;47:980-7.  Back to cited text no. 21      
22.Lang PJ. Behavioral treatment and bio-behavioral assessment. In: Sidowski JB, Johnson JH, Williams TA editors. Technology in mental health care delivery systems. Norwood, NJ: Ablex; 1980. p. 119-37.  Back to cited text no. 22      
23.Ziegler D, Laux G, Dannehl K, Spuler M, Muhlen H, Mayer P. Assessment of cardiovascular autonomic function: Age-related normal ranges and reproducibility of spectral analysis, vector analysis, and standard tests of heart rate variation and blood pressure responses. Diabetes Med 1992;9:166-75.  Back to cited text no. 23      
24.Li M, Vitnyi P. An introduction to Kolmogorov complexity and its applications, 2.Aufl. ed. Springer Verlag; 1997.  Back to cited text no. 24      
25.Lempel A, Ziv A. Universal algorithm for sequential data-compression. IEEE Trans Inf Th 1977;20:337-43.  Back to cited text no. 25      
26.Baumert M, Baier V, Haueisen J, Wessel N, Meyerfeldt U, Schirdewan A, et al. Forecasting of life threatening arrhythmias using the compression entropy of heart rate. Methods Inf Med 2004;43:202-6.  Back to cited text no. 26      
27.Bär KJ, Boettger MK, Koschke M, Schulz S, Chokka P, Yeragani VK. Non-linear complexity measures of heart rate variability in acute schizophrenia. Clin Neurophysiol 2007;118:2009-15.  Back to cited text no. 27      
28.Bertinieri G, Di Rienzo M, Cavallazzi A, Ferrari AU, Pedotti A, Mancia G. A new approach to analysis of the arterial baroreflex. J Hypertens Suppl 1985;3: S79-81.  Back to cited text no. 28      
29.Bär KJ, Boettger MK, Boettger S, Groteluschen M, Neubauer R, Jochum T, et al. Reduced baroreflex sensitivity in acute alcohol withdrawal syndrome and in abstained alcoholics. Drug Alcohol Depend 2006;85:66-74.  Back to cited text no. 29      
30.Baumert M, Walther T, Hopfe J, Stepan H, Faber R, Voss A. Joint symbolic dynamic analysis of beat-to-beat interactions of heart rate and systolic blood pressure in normal pregnancy. Med Biol Eng Comput 2002;40:241-5.  Back to cited text no. 30      
31.Bär KJ, Boettger MK, Berger S, Baier V, Sauer H, Yeragani VK. Decreased baroreflex sensitivity in acute schizophrenia. J Appl Physiol 2007;102:1051-6.  Back to cited text no. 31      
32.Oatley K, Jenkins JM. Human emotions: Function and dysfunction. Annu Rev Psychol 1992;43:55-85.  Back to cited text no. 32      
33.Lesperance F, Frasure-Smith N. Depression and heart disease. Cleve Clin J Med 2007;74: S63-6.  Back to cited text no. 33      
34.Baumert M, Baier V, Haueisen J, Wessel N, Meyerfeldt U, Schirdewan A, et al. Forecasting of life threatening arrhythmias using the compression entropy of heart rate. Methods Inf Med 2004;43:202-6.  Back to cited text no. 34      
35.Voss A, Kurths J, Kleiner HJ, Witt A, Wessel N, Saparin P. The application of methods of non-linear dynamics for the improved and predictive recognition of patients threatened by sudden cardiac death. Cardiovasc Res 1996;31:419-33.  Back to cited text no. 35      
36.Koelsch S, Remppis A, Sammler D, Jentschke S, Mietchen D, Fritz T, et al. A cardiac signature of emotionality. Eur J Neurosci 2007;26:3328-38.  Back to cited text no. 36      
37.Swaab DF, Bao AM, Lucassen PJ. The stress system in the human brain in depression and neurodegeneration. Ageing Res Rev 2005;4: 141-94.  Back to cited text no. 37      

Top
Correspondence Address:
Karl-Jurgen Bar
Department of Psychiatry and Psychotherapy, University Hospital Jena, Philosophenweg 3, 07743, Jena, Germany

Login to access the Email id


DOI: 10.4103/0019-5545.58894

PMID: 20174517

Get Permissions



    Figures

  [Figure 1]
 
 
    Tables

  [Table 1]

This article has been cited by
1 Self-referential processing influences functional activation during cognitive control: an fMRI study
G. Wagner,K. Koch,C. Schachtzabel,G. Peikert,C. C. Schultz,J. R. Reichenbach,H. Sauer,R. G. Schlosser
Social Cognitive and Affective Neuroscience. 2013; 8(7): 828
[Pubmed]
2 Autonomic nervous system in euthymic patients with bipolar affective disorder
Latalova, K., Prasko, J., Diveky, T., Grambal, A., Kamaradova, D., Velartova, H., Salinger, J., Opavsky, J.
Neuroendocrinology Letters. 2010; 31(6): 829-836
[Pubmed]



 

Top