Abstract
Mental stress-induced biological responses considerably differ across animals, which may be explained by intrinsic brain activity patterns. To address this hypothesis, we recorded local field potential signals from six cortical areas, electrocardiograms, and electromyograms from freely moving rats. Based on their stress-induced changes in cardiac signals, individual defeated rats were classified into stress susceptible and resilient groups. Rats with lower correlations in theta power across wide ranges of cortical regions before the stress challenge had higher probability to be stress-susceptible rats as defined based on the irregularity of heartbeat signals. A combination of principal component analysis and the support vector machine algorithm revealed that functional connections across cortical regions could be predictive factors accounting for individual differences in future stress susceptibility. These results suggest that individual differences in cortical activity may be a mechanism that causes abnormal activity of peripheral organs in response to mental stress episodes. This evidence will advance the understanding of the neurophysiological correlates of mind-body associations during mental stress exposure.
Similar content being viewed by others
Introduction
Stressful experiences trigger a series of complex biological responses in the whole body. Although many early studies revealed stress-induced biological factors at the molecular and cellular levels1,2, relatively few systemic physiological studies provided insights into stress responses, especially for functional interactions across the central and peripheral organs. At the macroscopic level, a number of behavioral studies have assessed how stressful experiences alter subsequent behavioral patterns, such as changes in anxiety-like and depressive-like behavior3,4,5. To connect existing evidence from the molecular to the behavioral levels, further studies are required to reveal the effects of stress loads on ongoing physiological organ dynamics.
Notably, stress-induced responses considerably differ across individual animals. A subset of animals exhibit large stress-induced reactions (so-called stress susceptibility), whereas other animals show no apparent changes (so-called stress resiliency)6,7,8. Since an initial step in most biological reactions induced by stress, especially mental stress, is information processing in the brain, the different degree of stress susceptibility may be accounted for by brain activity patterns. To support this hypothesis, recent studies have demonstrated that correlational activity across multiple brain regions before experiencing stressful episodes significantly differs and is associated with subsequent stress-induced responses, such as depression-like behavior9,10.
To address the hypothesis that brain activity is a crucial determinant of stress susceptibility, we performed multisite recordings of local field potential (LFP) signals from rat cortical regions. The rats were exposed to social defeat (SD) stress, which was shown to evoke profound changes in rat’s behavioral patterns, such as increased anxiety-like and depressive-like behavior6,7,8,11. Importantly, the stress loads applied to the rats in these studies yielded pronounced individual differences in their behavioral results. In this study, to quantify detailed stress responses in addition to behavioral phenotypes, we integrated recordings of electrocardiogram (ECG) signals, which represented the cardiac rhythm, and electromyogram (EMG) signals, which represented awake-related muscle contraction, into the recording system for our cortical-wide LFP recordings12,13,14. Based on the classification of stress-induced peripheral physiological activity, we analyzed whether multidimensional cortical LFP signals before experiencing SD stress were related to subsequent stress susceptibility.
Results
Simultaneous recording of ECG and EMG signals with cortical LFP signals
We recently established a recording system that enabled simultaneous monitoring of LFPs, ECG, and EMG in a freely moving rat using a custom-made plastic device (Fig. 1A). In this study, we created a plastic cover to mount on the rat’s head and protect the recording device against physical attacks from the other rat (Fig. 1A, middle). For the LFP recordings, six cortical regions were targeted, including the prelimbic cortex (PL), primary somatosensory cortex (S1), posterior parietal cortex (PPC), hippocampus (HPC), retrosplenial cortex (RSC), and primary visual cortex (V1) (Fig. 1B), and verified by histological inspections. These brain areas were selected so that they cover the wide range of cortex from anterior to posterior and from medial and lateral parts as much as possible. Figure 1C shows representative simultaneous recordings of all bioelectrical signals from the central and peripheral organs in a freely moving rat.
Classification of stress susceptibility based on heartbeat changes
In total, 19 rats were recorded and loaded with SD stress by exposing them to a resident rat for up to 10 min (termed a SD session). Before and after the SD session, the rat was placed in the rest box for 30 min and 120 min (termed the pre and post sessions, respectively), during which time the electrophysiological recordings of the LFP, ECG, and EMG signals were continuously monitored (Fig. 2A, top). The signals obtained from the SD session were not analyzed because these signals included a massive amount of electrical noise due to intense physical movement. During the pre and post session recordings, the heartbeat time was detected to quantify stress-induced changes in the irregularity of cardiac cycles. As demonstrated by a previous study15, some rats exposed to acute SD stress exhibited unstable cardiac cycles, including arrhythmia and increased HRV, in the post session. A susceptible rat was defined as a rat that showed either one of following criteria computed between the 10 min pre and 30 min post session: Δarrhythmia rate more than 0.2/min or ΔHRV more than 0.02. The classification results from all 19 defeated rats tested are summarized in Fig. 2B, showing a significant positive correlation between the two factors (r = 0.73, p = 3.5 × 10−4). Figure 2A (bottom) shows representative susceptible (magenta) and resilient (green) rats. As shown in Fig. 2C, irregular cardiac responses specifically occurred during the first 30 min of the post session in the susceptible group. The total duration of the SD and the duration of all attack patterns (i.e., severe attack, mild attack, sniffing, and no interactions) did not significantly differ between the susceptible and resilient groups (Fig. 2D; total duration: t17 = 0.91, p > 0.99; severe attack: t17 = 1.0, p > 0.99; mild attack: t17 = 0.31, p > 0.99; sniffing: t17 = 1.6, p = 0.65; no interactions, t17 = 0.32, p > 0.99, Student’s t-test with Bonferroni correction). These results suggest that the variation in stress susceptibility is not simply due to individual differences in the intensity and duration of the physical attacks from the resident rats. The running speed and the sleep duration in the pre session did not significantly differ between the two groups (Fig. 2E; speed: t17 = 0.054, p = 0.95; sleep duration: t17 = 0.74, p = 0.47, Student’s t-test). These results suggest that baseline locomotor activity before applying stress loads is not likely to be a determinant of stress susceptibility. In addition, we found no significant correlation between the duration of severe attack intensity and ΔHRV (r = 0.30, p = 0.21) or Δarrhythmia rate (r = 0.28, p = 0.27) (Fig. 2F), demonstrating that the variability of attack intensity during the SD session alone did not account for resultant stress susceptibility.
LFP signals at single cortical regions do not account for stress susceptibility
Next, we tested the possibility that cortical activity states during rest in the pre session might be associated with subsequent stress susceptibility. For each cortical region, LFP power in the pre session was computed by Fourier transformation. Figure 3A shows typical temporal changes in the root-mean-square (rms) of the EMG, time of arrhythmia, instantaneous R-R intervals, and hippocampal LFP power at the delta (1–4 Hz), 4–6 Hz (a middle band between delta and theta bands), theta (6–10 Hz), beta (10–25 Hz), and gamma (25–45 Hz) bands. For all cortical regions tested, we found no significant differences in the LFP power at all frequency bands from 1 to 45 Hz in the pre session between the susceptible and resilient groups (Fig. 3B; q > 0.05 at all frequency bands, FDR corrected). These results demonstrate that LFP power in single cortical regions alone cannot account for subsequent susceptibility against SD stress.
Correlational LFP power changes across multiple cortical regions are related to subsequent stress susceptibility
Then, we examined whether regional correlations of LFP power changes in the cortical network (so-called functional connectivity) might be associated with SD stress susceptibility. To obtain an overview of functional cortical organizations, first we computed correlation coefficients of LFP power changes in the 30-min pre session (n = 5 susceptible and 6 resilient rats). The averaged correlation coefficients obtained from all cortical region pairs (6C2 = 15) within each group are summarized in Fig. 4A,E,I,M, and 4Q for the delta, 4–6 Hz, theta, beta, and gamma bands, respectively. The cumulative distributions of these coefficients revealed significant differences in regional correlations between the susceptible and resilient groups at the lower ( < 10 Hz) bands (Fig. 4B,F,J; n = 60 and 66 pairs of cortical areas from 5 susceptible and 6 resilient rats; delta: Dmax = 0.27, p = 0.019; 4–6 Hz: Dmax = 0.32, p = 0.0019; theta: Dmax = 0.33, p = 0.0017, Kolmogorov-Smirnov test) but not the higher ( > 10 Hz) bands (Fig. 4N,R; beta: Dmax = 0.20, p = 0.16; gamma: Dmax = 0.23, p = 0.052, Kolmogorov-Smirnov test). These results demonstrate that the resilient rats show have higher temporal correlations of LFP power at the theta band across cortical-wide regions prior to loading SD stress, compared with the susceptible rats. To further examine whether such correlational changes in LFP power is associated with the strength of stress susceptibility, we computed correlations between averaged correlation coefficients of LFP power changes at each frequency band for each animal and its or ΔHRV. However, this analysis failed to detect significant correlations except for the gamma band (Δarrhythmia rate, delta: r = 0.51, p = 0.11; 4–6 Hz: r = 0.62, p = 0.068; theta: r = −0.17, p = 0.62; beta: r = −0.58, p = 0.059; gamma: r = 0.64, p = 0.034; ΔHRV, delta: r = 0.46, p = 0.15; 4–6 Hz: r = 0.051, p = 0.88; theta: r = −0.27, p = 0.42; beta: r = −0.50, p = 0.12; gamma: r = 0.63, p = 0.039), showing that the strength of stress-induced cardiac dysfunction is not simply accounted for by correlational LFP power changes across brain regions alone.
Based on these statistical results, we next addressed population analyses by applying the principal component analysis (PCA) and the support vector machine (SVM) algorithm to the same correlation coefficient datasets. First, the correlation coefficients of all cortical region pairs for each frequency band were decomposed by the PCA, and their first three principal components were plotted in a three-dimensional space (Fig. 4C,G,K,O,S). Next, a linear SVM, which is a supervised machine learning algorithm, was applied to the same datasets to determine the cutoff line that best separated the susceptible and resilient groups. Figure 4D,H,L,P,T shows all weights for individual cortical region pairs to define by the SVM analysis. The failure percentages calculated from the leave-one-out cross validation analysis were 63.7%, 9.1%, 9.1%, 18.2%, and 9.1%, for the delta, 4–6 Hz, theta, beta, and gamma bands, respectively. In other words, the SD stress susceptibility of more than 90% of the rats could be correctly predicted from cortical regional LFP power correlations, especially for > 4 Hz band. These results further support that functional connectivity across cortical-wide regions at resting states is a possible determinant of a rat’s susceptibility in response to subsequent stress episodes.
Post-stress LFP power changes in stress susceptible and resilient rats
The results described above focused on how cortical LFP power before applying stress was related to stress susceptibility. Next, we examined how the cortical LFP power changed after SD stress using the same analyses shown in Fig. 3B (Fig. 5). Fourier transformation detected no differences in the LFP power in individual cortical areas during the 30-min periods over a frequency band of 1–45 Hz between the resilient and susceptible groups (Fig. 5A; q > 0.05 at all frequency bands, FDR corrected). Cumulative distribution analyses revealed that the resilient group showed significantly higher regional correlation coefficients at the delta and theta band than the susceptible group (Fig. 5B; delta: Dmax = 0.29, p = 0.0075; 4–6 Hz: Dmax = 0.35, p = 0.00057; theta: Dmax = 0.35, p = 0.00071; beta: Dmax = 0.22, p = 0.072; gamma: Dmax = 0.080, p = 0.98).
Discussion
In this study, we obtained LFP signals simultaneously from six cortical regions together with ECG and EMG signals from a single rat exposed to acute SD stress. ECG recordings demonstrated that rats subjected to SD stress exhibited heterogeneous responses in cardiac signals, resulting in stress susceptible and resilient groups. A multidimensional analysis with the support vector machine showed that correlational power changes across cortical regions were predictive for such individual differences in stress susceptibility.
In the majority of stress studies, stress susceptibility has been quantified by changes in behavioral patterns, such as decreased social interactions and increased anxiety-like behavior in the social interaction and open field tests, respectively4,5,6,7,8. However, it is impossible to perfectly estimate animals’ mental states from specific behavioral signs alone. On the other hand, it has been well known that stress loads cause abnormal changes in peripheral activity patterns such as heartbeat signals and respiratory signals in both animal models and human patients, possibly via abnormal changes in autonomic transmission. Therefore, animal’s mental states can be more accurately evaluated by measuring autonomic physiological responses of internal organs. Based on this assumption, we recorded ECG and EMG signals during pre- and post-stress periods, allowing an ECG-based classification of the stress susceptible and resilient groups. These analyses confirmed that peripheral physiological signals can be an useful parameter to estimate animals’ psychiatric states in addition to simple behavioral patterns.
The rodent brain forms functional network organizations at a resting state (termed the default mode-like network)16,17 similar to the human brain. We showed that LFP power arising from individual cortical regions in the pre session, which was a period that possibly included a resting state, was not significantly different between the stress susceptible and resilient groups. Conversely, regional LFP power correlations in the pre session showed significant differences between the two groups. These results are consistent with previous studies showing that functional connectivity across certain brain regions, rather than activity levels of single brain areas, is crucial for accounting for mental disorders in rodents9,10,18,19 and humans20,21. Moreover, we demonstrated that individual differences in such functional organizations of cortical-wide networks observed before loading stress could predict the heterogeneity of future stress susceptibility. These insights are consistent with recent studies that have shown that monitoring cortical neuronal activity is a useful method to predict susceptibility to mood disorders, such as major depression disorder,10 and the therapeutic effects of drugs against mood disorders20.
Cortical LFP activity reflects populational neuronal activity in the brain and is linked with a variety of information processing, including cognition, learning and memory22. In the neocortex, generally delta LFP oscillations are larger during low arousal states, such as sleep and quiet states23,24. Our results showed that rats with higher theta power correlations, which are generally more dominant during movement or memory acquisition25,26, showed stable heartbeat signals against SD stress. These results suggest that rats showing stronger movement-related cortical activity, possibly representing higher arousal states, are more resilient to stress loads.
In traditional views, physiological organ activity has been described by the subdivisions of research area describing a single organ or single function, such as one brain region, a cardiac system, and a respiratory system. However, all organs continuously transmit information each other. To address this issue, a systemic viewpoint is further required to examine interactions across multiple organs. In this view, our study provides a typical example implying that peripheral dysfunctions may arise from cortical functional organizations with highlighting the importance of large-scale physiological data analyses in which both the central and peripheral organs are combined. In future studies, these research ideas will help deepen our understanding of the brain-body associations in health and disease.
Materials and Methods
Ethical approval
This study was carried out in accordance with NIH guidelines for the care and use of rats. The protocol was approved by the Experimental Animal Ethics Committee of the University of Tokyo (approval number: P29–7).
Animals
A total of 19 male Long-Evans rats (10–15 weeks old) with preoperative weights of 300–350 g were recorded and exposed to SD stress in this study. To apply SD stress to these rats, an additional 2 male Long-Evans rats (more than 4 months old) with weights of 450–550 g and 2 female Long-Evans rats (more than 4 months old) with weights of 300–400 g were used as the resident rats. The rats were housed individually and maintained on a 12-h light/12-h dark schedule with lights off at 7:00 AM unless otherwise specified. All rats were purchased from SLC (Shizuoka, Japan). After at least 1 week of adaptation to the laboratory, the rats underwent surgery.
Surgery
Surgery was performed to implant electrodes for recording ECG, EMG and LFP signals from a single rat. The detailed surgical procedures have been described elsewhere12,13,14. Briefly, before surgery, an electrode comprising a core body and a custom-made electrical interface board (EIB) accommodating at most 6 LFP channels, 2 ECG channels, 2 EMG channels, and 2 ground/reference channels was assembled. For the surgery, the rats were anesthetized with 1–2% isoflurane gas in air. For 19 rats, two incisions (∼1 cm) were made on both sides of the upper chest, and 2 stainless steel ECG electrodes with a tip dimeter of 0.147 mm diameter (AS633, Cooner Wire Company) in which the PTFE coating at the tip was peeled off at a length of ~5.0 mm were sutured to the tissue underneath the skin of the upper chest. Then, the rat was fixed in a stereotaxic instrument with two ear bars and a nose clamp. A midline incision was made from the area between the eyes to the incised neck area, and 2 stainless steel EMG electrodes with a tip dimeter of 0.147 mm diameter (AS633, Cooner Wire Company) in which the PTFE coating at the tip was peeled off at a length of ~5.0 mm were sutured to the dorsal neck muscles. For 11 of the 19 rats, LFP electrodes were implanted into the brain. Circular craniotomies 0.9 mm in diameter were made using a high-speed drill at the following coordinates: 3.0 mm anterior and 0.5 mm lateral to the bregma for the prelimbic cortex, 2.0 mm posterior and 2.8 mm lateral to the bregma for the primary somatosensory cortex, 3.8 mm posterior and 2.8 mm lateral to the bregma for the posterior parietal cortex and the hippocampus, 5.5 mm posterior and 0.5 mm lateral to the bregma for the retrosplenial cortex, and 5.5 mm posterior and 4.0 mm lateral to the bregma for the primary visual cortex. For the prelimbic cortex and hippocampus, electrodes were implanted at a depth of 2.5 mm. For the retrosplenial cortex, an electrode was implanted at a depth of 1.5 mm. For the primary somatosensory cortex, visual cortex, and posterior parietal cortex, electrodes were implanted at a depth of 1.2 mm. For the cerebellum, stainless steel screws were implanted on the skull attached to the brain surface to serve as ground/reference electrodes. The open edges of all electrodes were soldered to the corresponding channels on the EIB. All wires and the electrode assembly were secured to the skull using dental cement. After completing all surgical procedures, the anesthesia was removed, and the rats were allowed to awaken from the anesthesia spontaneously. Following surgery, each rat was housed in a transparent Plexiglas cage with free access to water and food.
Electrophysiological recording
At least 7 days after surgery, each rat with implanted electrodes was connected to recording equipment via Cereplex M (Blackrock), which is a digitally programmable amplifier, close to the rat’s head. The headstage output was conducted via a lightweight multiwire tether and a commutator to the Cereplex Direct recording system (Blackrock), which is a data acquisition system. LFP recordings were sampled at 2 kHz and filtered between 1 and 500 Hz.
Acute SD stress
A recording day for each rat consisted of three sessions: a 30-min pre session in a rest box (24 cm × 40 cm × 35 cm), an up to 10 min SD session in an open field (50 cm × 50 cm × 60 cm), and a 120-min post session in the same rest box. Acute SD stress was applied to the recorded rat similar to rat resident-intruder tests15,27. To prepare an experimental environment for the SD session, a male and a female rat were housed as resident rats in the open field at least one hour per day for 7 days before the SD session. When starting the SD session, the female resident rat was removed from the field, and the recorded rat was placed as an intruder rat into the open field. During the SD session, the intruder recorded rat was physically attacked by the resident male rat for 10 min. The SD session was immediately terminated if the intruder rat had a wound and bleeding due to the attack. The resident male rats were selected for their attack latencies shorter than 60 s when they were exposed to several rats that were not used for the subsequent experiments. The physical attack patterns of the resident rats were categorized into four types: (i) severe attack, including biting and chasing, (ii) mild attack, including mounting and touching, (iii) sniffing, and (iv) no interactions. The rat’s moment-to-moment position was tracked using a video camera attached to the ceiling with a frame rate of 3 Hz.
Histological analysis to confirm electrode locations
The rats were overdosed with urethane/α-chloralose, perfused intracardially with 4% paraformaldehyde (PFA) in phosphate-buffered saline (PBS; pH 7.4) and then decapitated. After dissection, the brains were fixed overnight in 4% PFA and then equilibrated with 30% sucrose in PBS. Frozen coronal sections (100 µm) were cut using a microtome, and serial sections were mounted and processed for cresyl violet staining. For cresyl violet staining, the slices were rinsed in water, counterstained with cresyl violet, and coverslipped with Permount. The positions of all electrodes were confirmed by identifying electrode tracks in the brain slices.
Analysis of electrophysiological traces
All analyses were performed using Matlab (Mathworks). ECG traces were bandpass filtered at 20–200 Hz, and beat-to-beat intervals (R-R interval) were calculated from the timestamp of the R-wave peaks. The root-mean-square (rms) of the EMG signals was computed in each bin with a bin size of 500 ms, which represented the absolute changes in EMG amplitude relative to the average. The LFP traces were downsampled to 200 Hz, and the LFP power was calculated by fast Fourier transform in each 5-s time window. Based on the time changes in LFP power, Pearson’s correlation coefficient was computed from a pair of cortical regions. Although our recordings targeted six cortical regions, which yielded complete datasets, one region was missed in 3 susceptible rats and 3 resilient rat, and two regions were missed in 1 resilient rat. In these rats, the data points at the missed areas were filled with average values of the corresponding areas computed from the other animals classified in the same group. The instantaneous speed of each frame was calculated based on the distance traveled within a frame (~333 ms). A sleep period was detected when the rms amplitude of the EMG was less than the mean–SD and lasted for 10 s, where the mean and SD were the average and standard deviation of the top 20–80% of the EMG rms in individual rats.
Definition of stress susceptibility by ECG signals
The defeated rats were classified into two groups based on two variables related to cardiac activity: the arrhythmia rate and heart rate variability (HRV). The time of arrhythmia (i.e., skipping an R-peak) was automatically detected and then manually scrutinized by eye in the original ECG traces. Arrhythmia rate was counted in each 1 min. HRV is computed as the coefficient of variation of all R-R intervals observed in each 1-min bin. In Fig. 2, the Δarrhythmia rate and ΔHRV were computed as differences in the averaged arrhythmia rates and averaged HRV between the first 30 min of the post session and the last 10 min of the pre session, respectively. Rats with a Δarrhythmia rate more than 0.2 or ΔHRV more than 0.02 were classified as the susceptible group, whereas the other rats were classified as the resilient group (Fig. 2B). These criteria were set based on the observation that all 13 control rats without experiencing SD stress never exhibited such higher values, Δarrhythmia rate more than 0.2 and ΔHRV more than 0.02.
PCA
Using PCA, three principal component dimensions were extracted from the set of correlation coefficients of LFP power changes at a specific frequency band computed from all pairs of cortical areas (e.g., Fig. 4A).
SVM
Linear SVM was applied to define the most appropriate separation line between the susceptible and resilient groups on the dataset of correlation coefficients of LFP power changes computed from all pairs of cortical areas. To assess the degree of separation by SVM, a leave-one-out cross-validation method was applied for all points (i.e., all rats). A failure percentage was defined as a probability in which a removed point was incorrectly assigned to the group to which the point belonged. All routines were written in Matlab (MathWorks).
Statistics
In Fig. 2D, the duration of individual behavioral patterns in the susceptible and resilient groups was compared by Student’s t-test, followed by post hoc Bonferroni correction. In Fig. 2E, the two variables in the susceptible and resilient groups were compared by Student’s t-test. In Figs 3B and 5A, LFP power at each frequency band in the susceptible and resilient groups was compared by Student’s t-test, followed by FDR correction. In Figs 4B,F,J,N,R and 5B, the difference between a pair of distributions was assessed by the Kolmogorov-Smirnov test. Correlational time changes in LFP power between groups were assessed by computing Pearson’s correlation coefficients. The null hypothesis was rejected at the p < 0.05 level unless otherwise specified. All data are presented as the mean ± standard error of the mean (SEM).
References
McEwen, B. S. et al. Mechanisms of stress in the brain. Nature neuroscience 18, 1353–1363, https://doi.org/10.1038/nn.4086 (2015).
Krishnan, V. & Nestler, E. J. The molecular neurobiology of depression. Nature 455, 894–902, https://doi.org/10.1038/nature07455 (2008).
Otabi, H., Goto, T., Okayama, T., Kohari, D. & Toyoda, A. The acute social defeat stress and nest-building test paradigm: A potential new method to screen drugs for depressive-like symptoms. Behavioural processes 135, 71–75, https://doi.org/10.1016/j.beproc.2016.12.003 (2017).
Bosch-Bouju, C., Larrieu, T., Linders, L., Manzoni, O. J. & Laye, S. Endocannabinoid-mediated plasticity in nucleus accumbens controls vulnerability to anxiety after social defeat stress. Cell reports 16, 1237–1242, https://doi.org/10.1016/j.celrep.2016.06.082 (2016).
Nie, X. et al. The innate immune receptors tlr2/4 mediate repeated social defeat stress-induced social avoidance through prefrontal microglial activation. Neuron 99, 464–479 e467, https://doi.org/10.1016/j.neuron.2018.06.035 (2018).
Golden, S. A., Covington, H. E. 3rd., Berton, O. & Russo, S. J. A standardized protocol for repeated social defeat stress in mice. Nature protocols 6, 1183–1191, https://doi.org/10.1038/nprot.2011.361 (2011).
Krishnan, V. et al. Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions. Cell 131, 391–404, https://doi.org/10.1016/j.cell.2007.09.018 (2007).
Berton, O. et al. Essential role of bdnf in the mesolimbic dopamine pathway in social defeat stress. Science 311, 864–868, https://doi.org/10.1126/science.1120972 (2006).
Kumar, S. et al. Prefrontal cortex reactivity underlies trait vulnerability to chronic social defeat stress. Nature communications 5, 4537, https://doi.org/10.1038/ncomms5537 (2014).
Hultman, R. et al. Brain-wide electrical spatiotemporal dynamics encode depression vulnerability. Cell 173, 166–180 e114, https://doi.org/10.1016/j.cell.2018.02.012 (2018).
Isingrini, E. et al. Resilience to chronic stress is mediated by noradrenergic regulation of dopamine neurons. Nature neuroscience 19, 560–563, https://doi.org/10.1038/nn.4245 (2016).
Okada, S., Igata, H., Sakaguchi, T., Sasaki, T. & Ikegaya, Y. A new device for the simultaneous recording of cerebral, cardiac, and muscular electrical activity in freely moving rodents. J Pharmacol Sci 132, 105–108, https://doi.org/10.1016/j.jphs.2016.06.001 (2016).
Sasaki, T., Nishimura, Y. & Ikegaya, Y. Simultaneous recordings of central and peripheral bioelectrical signals in a freely moving rodent. Biol Pharm Bull 40, 711–715, https://doi.org/10.1248/bpb.b17-00070 (2017).
Shikano, Y., Sasaki, T. & Ikegaya, Y. Simultaneous recordings of cortical local field potentials, electrocardiogram, electromyogram, and breathing rhythm from a freely moving rat. Journal of visualized experiments: JoVE, https://doi.org/10.3791/56980 (2018).
Sgoifo, A. et al. Incidence of arrhythmias and heart rate variability in wild-type rats exposed to social stress. The American journal of physiology 273, H1754–1760 (1997).
Upadhyay, J. et al. Default-mode-like network activation in awake rodents. PloS one 6, e27839, https://doi.org/10.1371/journal.pone.0027839 (2011).
Lu, H. et al. Rat brains also have a default mode network. Proceedings of the National Academy of Sciences of the United States of America 109, 3979–3984, https://doi.org/10.1073/pnas.1200506109 (2012).
Tamura, M., Mukai, J., Gordon, J. A. & Gogos, J. A. Developmental inhibition of gsk3 rescues behavioral and neurophysiological deficits in a mouse model of schizophrenia predisposition. Neuron 89, 1100–1109, https://doi.org/10.1016/j.neuron.2016.01.025 (2016).
Henckens, M. J. et al. Stress-induced alterations in large-scale functional networks of the rodent brain. NeuroImage 105, 312–322, https://doi.org/10.1016/j.neuroimage.2014.10.037 (2015).
Dunlop, B. W. et al. Functional connectivity of the subcallosal cingulate cortex and differential outcomes to treatment with cognitive-behavioral therapy or antidepressant medication for major depressive disorder. The American journal of psychiatry 174, 533–545, https://doi.org/10.1176/appi.ajp.2016.16050518 (2017).
Drysdale, A. T. et al. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nature medicine 23, 28–38, https://doi.org/10.1038/nm.4246 (2017).
Buzsaki, G. Rhythms of the brain. New York: Oxford University Press (2006).
Vyazovskiy, V. V. et al. Cortical firing and sleep homeostasis. Neuron 63, 865–878, https://doi.org/10.1016/j.neuron.2009.08.024 (2009).
Watson, B. O., Levenstein, D., Greene, J. P., Gelinas, J. N. & Buzsaki, G. Network homeostasis and state dynamics of neocortical sleep. Neuron 90, 839–852, https://doi.org/10.1016/j.neuron.2016.03.036 (2016).
Buzsaki, G. Theta oscillations in the hippocampus. Neuron 33, 325–340 (2002).
Colgin, L. L. Rhythms of the hippocampal network. Nat Rev Neurosci 17, 239–249, https://doi.org/10.1038/nrn.2016.21 (2016).
Marini, F. et al. Single exposure to social defeat increases corticotropin-releasing factor and glucocorticoid receptor mrna expression in rat hippocampus. Brain research 1067, 25–35, https://doi.org/10.1016/j.brainres.2005.10.002 (2006).
Acknowledgements
This work was supported by Kaken-hi (17H05939 and 17H05551), a grant-in-aid for Precursory Research for Embryonic Science and Technology, the Exploratory Research for Advanced Technology, and the Senri Life Science Foundation.
Author information
Authors and Affiliations
Contributions
R.N. and T.S. designed the study. R.N. acquired the electrophysiological data and performed the analyses. Y.I. supervised the project. R.N. and T.S. prepared all figures and wrote the main manuscript text. All authors reviewed the main manuscript text.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Nakayama, R., Ikegaya, Y. & Sasaki, T. Cortical-wide functional correlations are associated with stress-induced cardiac dysfunctions in individual rats. Sci Rep 9, 10581 (2019). https://doi.org/10.1038/s41598-019-47171-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-47171-y
This article is cited by
-
Mesolimbic dopamine release precedes actively sought aversive stimuli in mice
Nature Communications (2023)
-
Selective effects of acute and chronic stress on slow and alpha-theta cortical functional connectivity and reversal with subanesthetic ketamine
Neuropsychopharmacology (2023)
-
Weak representation of awake/sleep states by local field potentials in aged mice
Scientific Reports (2022)
-
Is depression a disorder of electrical brain networks?
Neuropsychopharmacology (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.