3.1 Electrical performance characterization of MEAs
Figure 2a and b showed the packaged MEA and an enlarged view of the modified MEA tip under the microscope. Black modification materials were uniformly distributed on each microelectrode site. The scanning electron microscopy (SEM) image indicated that PtNPs/PEDOT:PSS nanocomposite adhered closely to the surface of the microelectrode (Figure 2c and d). The PtNPs improved the roughness of the microelectrode, and PEDOT:PSS structure covered the entire surface of the electrode, which was conducive to improving the stability of the electrode. The modification of PEDOT also improved the biocompatibility and long-term detection capability of microelectrodes in vivo (Figure S2).
We compared the electrical performance of bare microelectrodes and modified microelectrodes in the frequency range from 10 Hz to 1 MHz. The phase delay and impedance of the bare microelectrode were both large in the frequency domain from 10 Hz to 1 MHz (Figure 3a and b), which was not conducive to the detection of electrophysiological signals in vivo. The phase delay and impedance of the PtNPs/PEDOT:PSS-modified microelectrodes were significantly better than the bare microelectrodes and the PEDOT:PSS-modified microelectrode at 10 Hz to 1 MHz. The central frequency of the neural activities is around 1 kHz, at which the impedance and phase delay of the three were compared in this work. The microelectrodes impedance decreased from 847.07 ± 334.3 kΩ (bare) to 156.18 ± 64.14 kΩ (PEDOT:PSS) and then to 25.55 ± 10.33 kΩ (PtNPs/PEDOT:PSS) (Figure 3d). There was no significant difference in phase between PEDOT:PSS-modified microelectrodes (-31.41 ± 6.28°) and PtNPs/PEDOT:PSS-modified microelectrodes (-28.89 ± 1.15°), which indicated that PEDOT played a major role in phase delay (Figure 3e). In addition, the Nyquist plots of the impedance spectra showed that the PtNPs/PEDOT:PSS-modified microelectrode had intense diffusion characteristics (Figure 3c). Our results meant that the PtNPs/PEDOT:PSS-modified microelectrodes had better electrical performance, which was helpful for the detection of neuronal information transmission.
3.2 Characterization of behavioral and electrophysiological signals during 2MT-induced defensive behavior
To interrogate the role of dPAG and vPAG in the defense behaviors related to innate fear, we combined the electrophysiological activities of dPAG and vPAG neurons recorded by MEA with the behavior changes recorded by the digital camera before and during the defense behavior evoked by exposure to 2MT. We associated the changes in motor output related to the defense behaviors with their neural activity patterns, which allowed us to evaluate the endogenous neural mechanism of dPAG and vPAG under different defense behaviors. After exposure to 2MT, we observed two different innate defense behaviors in rats (Figure 4a and d). Velocity and freezing rate are important indexes of defense behavior (flight and freezing) in rats. The rats showed the obvious flight avoidance reaction, which was characterized by an extremely low freezing rate and high moving velocity. We define this state as flight. Then, after the rats found that there was no hope of escape, they stopped running and stayed in place almost motionless and began to enter the freezing state. This freezing mechanism made the prey more likely to avoid being captured when threatened. As shown in Figure 4e, we also counted the temporal heatmap of the spatial position of the rat. It could be found that the rats in the control moved evenly throughout the space, and the rats in the flight circled the space. When rats entered the freezing stage, they stayed far away from the filter paper (the source of predator odor) for a long time.
Our fabricated MEA was a brain-computer interface (BCI) with a high spatiotemporal resolution that provided a powerful strategy for monitoring real-time neuronal electrophysiological activities. The recording sites of MEA were in close contact with the neurons in the target area, so the neural information recorded in each channel could accurately reflect the neuronal activities. Figure 4b and c showed the neural spikes and LFP recorded by the same rat during the experiment. Compared with the control, the firing density of spikes was significantly more intensive and LFP had higher amplitude and frequent fluctuation during the flight in the dPAG, while the spike firing and LFP amplitude during the freezing were between the control and flight. The difference was that in the vPAG, the release density of spikes was the most intensive in the freezing, and the fluctuation and amplitude of LFPs increased to the most frequent. The above results indicate that dPAG and vPAG have different activation levels in flight and freezing.
In conclusion, the above results suggested that the highly efficient predator odor analog 2MT induced innate fear (defense) behaviors in rats. Both dPAG and vPAG were significantly activated during the defensive behavior. dPAG was highly activated during the flight, while vPAG was highly activated during freezing.
3.3 Information Transmission and Functional Analysis of vPAG and dPAG neurons in Defense State
To explore the activities of dPAG and vPAG neurons in the whole process of defensive behavior, we counted the change in the average spike firing rate during the experiment (Figure 5a). After exposure to 2MT, the spike firing rate of dPAG neurons reached a peak during the flight, and gradually decreased and stabilized after the peak. It indicated that dPAG was fastly activated to a high degree after exposure to 2MT, and then the activation degree was gradually stable during the freezing. The spike firing rate of vPAG neurons began to increase after exposure to 2MT, which indicated that the activation degree of vPAG neurons gradually increased after exposure to 2MT. LFP reflects the real-time transmission of information across the neural network and is an important supplement to measure neural information activity. The energy change of the dPAG and vPAG LFP spectrogram is consistent with the changing trend of the spike firing rate in the defense state (Figure 5b). To analyze the correlation and spread of defense behavior regulation in PAG subregions, the time when peaks of LFP and spike-firing occurred in the dPAG and vPAG were compared after 2MT-induced innate fear (Figure 5c). LFP peak and spike-firing first appeared in dPAG and gradually transmitted to the vPAG. In addition, the LFP waveforms of dPAG and vPAG had roughly the same components. The transmission was found in multiple channels and accompanied by the release of spikes. The transmission time between typical channels was about 80 ms.
As shown in Figure 5d, the difference in the average spike firing rate of dPAG and vPAG was not significant during the control. During the flight, the average spike firing rate of dPAG and vPAG neurons increased significantly, but dPAG neurons were significantly higher than that vPAG neurons. During the freezing, the average firing rate of vPAG neurons continued to increase, while that of dPAG neurons decreased compared with the flight and was lower than that of vPAG neurons. The average power of LFP in dPAG and vPAG was calculated during different behaviors. Compared with the control state, the average LFP power of dPAG in the flight and freezing state increased significantly, the LFP power in the flight state increased to 0.49 ± 0.07mW, and the LFP power in the freezing decreased to 0.36 ± 0.05 (Figure 5e). In the vPAG (Figure 5f), the average LFP power in the flight and freezing increased significantly compared with the control, and the LFP power in the freezing state (0.34 ± 0.03 mW) was slightly higher, but there was no significant difference compared with the flight (0.32 ± 0.04 mW). Then, we normalized the power of LFPs in three states and analyzed their relative characteristics in the 0-20 Hz range. We divided it into the delta frequency band (1-4 Hz), theta frequency band (4-8 Hz), alpha frequency band (8-13 Hz), and beta frequency band (13-30 Hz) (Figure S4b). After exposure to 2MT, the energy of dPAG was more concentrated in the delta band, and the proportion of the delta band in the flight was significantly higher than freezing and control, while the change of vPAG was reflected in the increase of the proportion of theta band.
The results showed that the rapid activation of dPAG was consistent with the flight behavior after exposure to 2MT, and accompanied by the transmission of neuronal defense information from dPAG to vPAG. When entering the freezing state, the activation degree of vPAG neurons gradually rose, and the activation degree of dPAG neurons began to decline. And the spike firing rate of vPAG neurons was highly correlated with the freezing rate (Figure S6). In other words, the activities of dPAG may be more related to sports and active defense response, while vPAG was consistent with the passive defensive.
3.4 Comparison and Analysis of dPAG and vPAG Neurons Spike Characteristics
As we detected that dPAG and vPAG neurons had different activation degrees and exhibited different regulation modes during flight and freezing, this work compared and analyzed the diverse firing pattern of the two types of neurons in defense information transmission. The signal of 84 channels recorded was separated into spike units using principal component analysis (PCA) and valley searching methods (VSM). We calculated the average waveform of individual spike units and analyzed their autocorrelation. As shown in Figure 6a, the average waveform of dPAG neurons had shorter peak intervals and a narrower waveform than that of vPAG neurons. The spike-time autocorrelogram of dPAG neurons had a shorter latency to firing than that of vPAG neurons (Figure 6b). Pooling all dPAG and vPAG neurons, we calculated the trough to right peak latency and peak amplitude asymmetry of the average waveform (see Statistical methods for the calculation method). Figure 6c showed that there was a clear distinction between the waveform types of dPAG and vPAG neurons. The trough-to-right peak latency of the dPAG neurons' average waveform was generally shorter than that of vPAG neurons, and the spike asymmetry index was higher.
In summary, The firing characteristics of dPAG and vPAG neurons were obviously different. The waveform of dPAG neurons was narrow and the spike-time autocorrelograms latency period was shorter, which also confirms that dPAG neurons were more likely to exhibit burst or burst-like firings, thus providing physiological conditions for controlling the flight behavior that requires the instant response.
3.5 Changes in spike characteristics of neurons in different defense states
The spike-firing characteristics of neurons provided the physiological basis for neurons in the regulation of defense behavior. To study the characteristic changes of neuronal spikes under different defense states, we also calculated the average spike waveform under different behavior states (Figure 6d). After exposure to 2MT, the average amplitude of spikes dPAG-neurons spikes during control was significantly lower than during flight and freezing. The amplitude was the highest during flight, and the amplitude during freezing was between flight and control. For the vPAG average spike waveform, the amplitude of flight and freezing state was similar, which was slightly greater than the control state. The above changes confirmed that the waveform of PAG neurons changed significantly under the defensive state. It was noteworthy that dPAG neurons showed more significant waveform changes, which seemed that dPAG neurons were more susceptible to being affected by innate fear.
The autocorrelograms of PAG neurons displayed behavior-associated firing changes under different defensive. dPAG neurons during flight had a higher probability of firing with short inter-spike intervals than during control and freezing (Figure 7a). The distribution of vPAG neuron inter-spike intervals was relatively uniform in freezing than in flight and control (Figure 7b). Pooling all dPAG and vPAG neurons, we also calculated the average autocorrelogram of neurons before and after exposure to 2MT. As shown in Figure 7c, compared with vPAG neurons, dPAG neurons had a higher probability of firing with short spike-inter intervals. We linearly fit the spike firing rate of dPAG and vPAG neurons to the amplitude and latency of their autocorrelogram peak respectively. We found that the spike firing rate was positively correlated with the amplitude of the peak value of the autocorrelogram (Figure S5a), and negatively correlated with the latency of the peak value of the autocorrelogram (Figure S5b).
These results indicated that a higher amplitude of the autocorrelogram peak and a shorter peak latency predicted a higher neuron firing rate. dPAG and vPAG neurons had changes in firing characteristics under different defense states. While dPAG was more prone to burst or burst-like firing in the flight. The spike firing time of vPAG neurons was more evenly distributed in the freezing state so that vPAG can provide stable and continuous neuroregulatory information to match the long-term freezing behavior.