Long-term potentiation prevents ketamine-induced aberrant neurophysiological dynamics in the hippocampus-prefrontal cortex pathway in vivo

N-methyl-D-aspartate receptor (NMDAr) antagonists such as ketamine (KET) produce psychotic-like behavior in both humans and animal models. NMDAr hypofunction affects normal oscillatory dynamics and synaptic plasticity in key brain regions related to schizophrenia, particularly in the hippocampus and the prefrontal cortex. It has been shown that prior long-term potentiation (LTP) occluded the increase of synaptic efficacy in the hippocampus-prefrontal cortex pathway induced by MK-801, a non-competitive NMDAr antagonist. However, it is not clear whether LTP could also modulate aberrant oscillations and short-term plasticity disruptions induced by NMDAr antagonists. Thus, we tested whether LTP could mitigate the electrophysiological changes promoted by KET. We recorded HPC-PFC local field potentials and evoked responses in urethane anesthetized rats, before and after KET administration, preceded or not by LTP induction. Our results show that KET promotes an aberrant delta-high-gamma cross-frequency coupling in the PFC and an enhancement in HPC-PFC evoked responses. LTP induction prior to KET attenuates changes in synaptic efficiency and prevents the increase in cortical gamma amplitude comodulation. These findings are consistent with evidence that increased efficiency of glutamatergic receptors attenuates cognitive impairment in animal models of psychosis. Therefore, high-frequency stimulation in HPC may be a useful tool to better understand how to prevent NMDAr hypofunction effects on synaptic plasticity and oscillatory coordination in cortico-limbic circuits.

In rodents, ketamine and other NMDAr antagonists are widely used as a translatable pharmacological model capable of inducing psychotic-related behaviors, such as hyperlocomotion, working memory impairments, prepulse inhibition disruption, and abnormal social interaction 14,15 . Several neurophysiological features of the HPC-PFC pathway are associated with psychotic-like behaviors induced by NMDAr antagonists. In vivo experiments showed that ketamine increased gamma  power in the PFC and theta-gamma phase-amplitude coupling in the HPC 16 . Besides, it affected the synaptic transmission efficiency in the HPC-PFC pathway, and disrupted long-term potentiation (LTP) [17][18][19] . The role of NMDA receptors for LTP is well-described 20,21 . However, the underlying mechanisms of NMDA antagonists-induced aberrant gamma oscillations and facilitation of synaptic transmission are still not fully understood. A plausible explanation is that NMDAr antagonists may preferentially exert their effect on inhibitory interneurons, leading to a hyperglutamatergic state produced by the disinhibition of principal cells 22,23 . Consistently, it has been demonstrated that acute treatment of MK-801, a non-competitive NMDAr antagonist, induces long-lasting changes in the synaptic efficacy in the CA1-PFC pathway of rats in vivo 34 . This effect was accompanied by disruption of performance in a PFC-dependent task. Interestingly, the authors showed that MK-801-induced facilitation of basal synaptic efficacy was occluded by tetanus-induced LTP and vice versa, suggesting shared mechanisms between both processes. In other words, the massive activation of α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors (AMPAr) during a hyperglutamatergic state that follows the NMDAr antagonism seems to trigger an aberrant form of LTP. On the other hand, a prior potentiation of AMPAr transmission by LTP induction 24 somehow produced a state of reduced sensitivity to MK-801 25 . Consistently, prior studies showed that allosteric modulation of AMPAr can prevent NMDAr antagonist-induced impairments 27,29,30 . It has been hypothesized that a prior AMPAr activation can produce depolarization of membrane potentials, relieving the Mg 2+ block of colocalized NMDAR and, consequently, attenuating NMDAr hypofunction 26,27 . In another scenario, not mutually exclusive, a prior LTP could occlude the effects of NMDAr antagonists by competing mechanisms. These mechanistic hypotheses remain to be tested.
Despite its limitations, the urethane anesthesia has been regarded as a valid model to test mechanistic hypotheses involving brain oscillations and synaptic plasticity in normal [28][29][30][31] and pathological situations 32,33 . Particularly relevant for this study, rodents under urethane anesthesia have been useful to examine which are the effects induced by NMDAr antagonists that are independent of the typical hyperlocomotion observed after treatment 18 . Besides, this preparation can be used to investigate how NMDAr antagonists affect oscillations during different sleep-like brain states 34 . In this sense, Kiss et al. (2011) demonstrated that acute MK-801 produced a new brain state of persistent slow oscillations (0.5-1.5 Hz) associated with short-term synaptic plasticity impairment at subiculum-mPFC pathway in urethane-anesthetized rats. These MK-801 induced effects were reversed by treatment with the AMPA receptor potentiator (LY451395), reinforcing the idea that AMPAr transmission enhancement could be an interesting target to attenuate alterations produced by NMDAr hypofunction 35 . Considering this, we tested whether prior LTP was also able to occlude the ketamine effects on synaptic plasticity and oscillatory patterns in the HPC-PFC pathway. To test this hypothesis, we first characterized sleep-like brain state oscillations under urethane 35 and compared PFC field responses induced by HPC stimulation and HPC-PFC oscillatory coupling before and after ketamine treatment. Then we tested the modulatory effect of deep brain stimulation (LTP at HPC-PFC synapses) on ketamine-driven oscillatory patterns. Our findings show that ketamine induces a state of increased gamma frequency power and abnormal cross-frequency coupling in the PFC, followed by an enhancement of HPC-PFC evoked responses. Although LTP induction also increases gamma frequency activity, it prevented ketamine-induced aberrant oscillatory coupling and potentiation of HPC-PFC synaptic transmission. Our data suggest that modulation of synaptic plasticity in HPC-PFC circuits might be a useful tool to better understand how to prevent dysfunctions induced by NMDAr blockade.
HPC-PFC connectivity under urethane anesthesia. Before directly examining the ketamine effects on HPC-PFC neurophysiological dynamics, we sought to further our understanding of how urethane brain states per se can modulate functional connectivity between these regions in control conditions (Fig. 2). We show that while DEA epochs show synchrony in delta rhythms, spectral coherence during ACT state demonstrates a peak in theta oscillations (4 Hz, Fig. 2A). These patterns of synchronicity closely follow the alternation of brain states ( Fig. 2A). We further explored the directionality of this connectivity by performing cross-correlation and Granger causality analyses. In Fig. 2B we show that the peak of cross-correlation in delta frequency (predominant oscillation on DEA epochs) has a positive lag (2 ms), while theta oscillations indicate a negative lag (−22 ms). Granger analysis reveals that ACT periods show theta oscillations with the HPC leading the PFC (t (7) = 2.4909, p = 0.0415), while on DEA epochs we found a peak in delta oscillations (~1.5 Hz) with the PFC driving the HPC oscillations (t (7) = 7.7401, p = 0.0001). These results indicate that during ACT states hippocampal theta coordinates activity in the  (1), and raw LFP during deactivated (2) and activated (3) states in PFC (top) and HPC (bottom). Representative 210 min recording of a raw LFP (middle) and spectrograms (bottom) in the PFC and HPC showing spontaneous brain state alternation during urethane anesthesia. (D) Brain state classification. Scatter plots of RMS values and the number of zero crossing were used to classify brain states in deactivated (DEA) and activated (ACT) epochs (left). RMS and state classification in a representative recording (right). (E) Relative power spectral densities in the PFC and HPC (line plots) during different brain states (SAL, n = 8). In both regions, activated states present theta oscillation and an increase in gamma frequencies, while in deactivated state there is a predominance of delta oscillation (bar plots). *p < 0.05, ~p = 0.0506.
NMDAr blockade effects on HPC-PFC connectivity. Next, we investigated the effects of acute ketamine injection on oscillatory patterns and basal synaptic efficacy in the HPC-PFC pathway of rats under urethane. In Fig. 3B we illustrate the ketamine effects on the evoked potentials of the HPC and PFC. We show that NMDAr blockade leads to an increase in evoked field post-synaptic potential (fPSP) amplitude (F (20,260) = 1.611 significant interaction between groups p = 0.050; Bonferroni post-hoc test are bottom of the graphs) and a robust reduction of paired-pulse facilitation (PPF) (F (20,260) = 8.205 significant interaction between groups p = <0,001). Figure 3C shows a representative raw LFP from both regions indicating ketamine effects on brain oscillations. As described in other studies 18,40,41 , we verified that NMDAr blockade increased low-gamma and high-gamma power in the PFC independent of brain state (DEA: low-gamma: t (6) =2.327, p = 0.0589; high-gamma: t (6) = 3.423, p = 0.0141; ACT: low-gamma: t (5) = 4.108, p = 0.0093; high-gamma: t (5) = 3.258, p = 0.0225; Fig. 3D). In PFC, no significant effects were observed in other frequency bands, while in HPC there is a slight increase in gamma power ( Figure S1A). Interestingly, systemic ketamine appears to increase the probability of DEA state. In Fig. 3E we demonstrate the probability density function (PDF) of DEA epoch occurrence in all animals in the KET group compared with animals of the SAL group. While the probability of SAL group oscillates between low and high values during the recording, this oscillation pattern is reduced after ketamine injection, with the probability of a DEA epoch occurrence kept high for at least 30 min after systemic injection (see also Figure S2). This increase in DEA oscillatory activity is not accompanied by alterations in synchrony between the HPC-PFC LFP accounting showing delta (~1 Hz) and theta (~4 Hz) synchrony alternation during brain states. Averaged spectral coherence plot (right). Inset showing coherence peak ~4 Hz and ~1 Hz during ACT and DEA epochs, respectively. The dashed line represents 95% confidence interval of bootstrapped data. (B) Histogram of cross-correlogram lags (peak) of delta and theta oscillations in DEA an ACT (left). Granger causality in the ACT (middle) and DEA states (left). The dashed line represents 95% confidence interval of permuted data. Insets: box plots distribution of maximum Granger theta (left) and delta (right) power, showing significant causality in the HPC-PFC direction during activated state (left), and PFC-HPC direction during deactivated state. (C) Representative comodulation map showing delta-high-gamma coupling during deactivated state (left). Middle: Representative filtered LFP (top) and mean high-gamma amplitude as a function of delta phase in DEA state and surrogate data (bottom). Mean delta-high-gamma coupling strength (MI) during brain states (right) and surrogate data. MI is significantly increased during deactivated states. All data correspond to SAL group (n = 8). *(p < 0.05).
for epoch separation. As we show in Fig. 3F there are no significant changes between spectral coherence of DEA or ACT states comparing before or after ketamine injection. However, as shown in a representative coherogram, systemic ketamine leads to a persistent synchronization in delta, which follows the increase in DEA epochs. www.nature.com/scientificreports www.nature.com/scientificreports/ NMDAr blockade promotes a different oscillatory state in the PFC. Since ketamine produced a persistent and global state of delta oscillations (DEA state), we asked whether this effect was associated with alterations of information processing in local networks by calculating the phase-amplitude coupling 42 in the PFC. In Fig. 3G we represent high-gamma amplitude as a function of delta phase comparing DEA epochs before and after ketamine injection. We show here that NMDAr blockade increases cross-frequency coupling between delta and high-gamma activity (t (6) = 4.947, p = 0.0026). To investigate if the increase of the CFC was dependent on power changes promoted by ketamine we performed a linear correlation between DEA baseline-normalized MI www.nature.com/scientificreports www.nature.com/scientificreports/ and delta or high-gamma power values. Figure 3H shows a weak correlation between delta or high-gamma with MI values after ketamine injection (r = 0.09 for delta and r = 0.22 for high-gamma). To further investigate if the increase in CFC could be influenced by high-gamma power after ketamine we compared MI averaged values before and after the drug. For this, we only used post-ketamine epochs with high-gamma within 95% confidence interval of the pre-ketamine high-gamma distribution. We show that even epochs with no increase in high-gamma values compared to baseline present an average significant enhancement in CFC (Fig. 3H right; t (882) = 8.0693, p < 0.001), indicating that ketamine effects in delta-high-gamma coupling may be independent of power increase in these frequencies. We confirmed this result with bootstrap analysis using 1000 repetitions and controlling for the number of trials. These data indicate that ketamine promotes an oscillatory state that is different from the traditional DEA state in the PFC. Figure 3I illustrates all the 20 s DEA epochs from the KET group in a 3D plot showing how ketamine modifies brain oscillatory activity in DEA epochs increasing gamma frequency band and delta-high-gamma coupling.
LTP increases gamma activity in the PFC. Before testing for a possible modulation of KET effects by prior LTP induction, we sought to examine the electrophysiological effects produced by high-frequency stimulation (HFS) alone. In Fig. 4 we show the effects of HFS on HPC-PFC synaptic efficacy and oscillatory patterns. As shown previously by our group 30,32,42 , HFS protocol induces a stable LTP in the HPC-PFC pathway for at least 120 min (fPSP1; F (20,240) = 13.761, p < 0.0001). This LTP can be divided into the transient decaying phase, termed short-term potentiation (STP, Fig. 4B) and the more stable LTP per se. LTP effects were also seen in the amplitude of the fPSP2 (paired-pulse stimulation; F (20,240) = 11.470, p < 0.0001). However, the increase in the fPSP2 is relatively lower than the one in fPSP1, which is seen as a reduction in the PPF ratio at least for 60 min after HFS (F (20,240) = 5.621, p < 0.0001; Fig. 4B). Figure 4C demonstrates the effects of HFS on LFP of the PFC. LTP induction produced an increase in low and high-gamma on the PFC restricted to DEA epochs (low-gamma: Mann-Whitney test, U = 30, p = 0.0006; high-gamma: Mann-Whitney test, U = 27, p = 0.0003). No significant alterations were observed at gamma frequencies in the HPC of LTP groups ( Figure S1B). Interestingly, the increase in high-gamma power was not related to an enhancement of delta-high-gamma coupling. As shown in Fig. 4D the modulation index did not differ from before or after LTP induction (p > 0.05).

Ketamine effects are attenuated by prior LTP. Because allosteric modulation of glutamatergic receptors
prevents psychotic-like effects 26,[43][44][45][46][47] , and induction of LTP enhances both AMPAr and NMDAr activity 48-50 , we asked whether previous LTP induction can prevent ketamine effects on HPC-PFC pathway. To test this hypothesis, we applied HFS 30 min prior to KET (see Methods) or vehicle administration (LTP-KET and LTP-SAL groups) and compared the results. Figure 4F indicates that ketamine treatment following LTP did not affect the www.nature.com/scientificreports www.nature.com/scientificreports/ evoked potentials on the PFC either on fPSP1, fPSP2 or PPF comparing with control group (p > 0.05). Similar to what we observed in the LTP-SAL group, LTP-KET also presented an increase in low and high-gamma activity in the DEA epochs of the PFC (low-gamma: t (7) = 2.229, p = 0.0611; high-gamma: t (7) = 2.931, p = 0.0220). While these effects are also seen in the LTP group LTP-KET shows an increase in gamma frequency also during ACT states (low-gamma: t (8) = 7.166, p < 0.0001; high-gamma: t (8) = 7.709, p < 0.001; Fig. 4G). These alterations appear to be related to the ketamine administration since it was observed in the KET group as well (Fig. 3D). Despite the high-gamma power increase, no significant alterations were observed in the delta-high-gamma coupling in DEA states (Fig. 4H), indicating that LTP induction prevents ketamine effects on PFC CFC. LTP prevents ketamine-induced aberrant oscillatory brain state. Finally, we explored the possibility that prior LTP induction attenuates the aberrant oscillatory patterns induced by KET. Figure 5A shows the absolute MI value across the entire recording in the KET and LTP-KET groups. We demonstrate that LTP induction attenuates ketamine effects on delta-high-gamma coupling especially in the first 30 minutes after drug injection. We next compared the baseline-normalized MI values for only DEA epochs in the initial 30 min after ketamine. Our data indicate that LTP attenuates the effects of ketamine on enhancing CFC coupling in the PFC (t (12) = 2.589, p = 0.0237). Figure 5B shows that KET and LTP-KET epochs after ketamine injection can be distinguished based on their electrophysiological features. We used principal component analysis for dimensionality reduction of electrophysiological features, and then applied a quadratic discriminant analysis to classify epochs in the two groups. For the quadratic function, we used the first three principal components, which have an explained variance of 79.21% (Fig. 5C). Our classification model with training data (filled circles) and then cross-validated with test data (open triangles) gave a classification accuracy of 85.56% accuracy. This data shows that ketamine effects can be distinguished when preceded by LTP based on their electrophysiological characteristics. We next compared the epochs mean score from each group for the first three principal components (Fig. 5D). Our results revealed that KET has a higher score on PC1 (p < 0.0001), while LTP-KET has a higher score in PC2 (p = 0.0069) and PC3 (p = 0.0001). Interpreting the correlation coefficients (loadings) of the original variables with the first three principal components we observe that PC1 has a high correlation with CFC values (r = 0.7696, p < 0.0001), high-gamma (r = 0.5975, p < 0.0001) and delta activity (r = 0.8963, p = 0.0001). PC2 and PC3 have high correlation with theta (r = 0.8368, p < 0.0001) and high-gamma (r = 0.6750, p < 0.0001), respectively, and negative correlation with fPSP (r = −0.5214, p < 0.0001 and r = −0.5375, p < 0.0001, respectively).

Discussion
In this study, we demonstrated that a prior enhancement of synaptic efficacy at hippocampal-prefrontal projections is sufficient to attenuate the disturbing effects of ketamine on oscillatory coupling and basal synaptic transmission in the prefrontal cortex. Ketamine was shown to increase gamma frequency power and delta-gamma phase-amplitude coupling in the PFC and boosted HPC-PFC synaptic plasticity with PPF disruption. We also observed that LTP induction was associated with an increase of gamma power in DEA states in the PFC and no alteration in delta-gamma coupling in the PFC.
Under urethane anesthesia, rodents show a spontaneous alternation of brain states characterized by activated and deactivated periods 34,36 . We expanded previous work demonstrating a 4 Hz synchronization from the HPC to the PFC during the activated state and a 1 Hz synchronization from the PFC to the hippocampus in the deactivated. Furthermore, we describe a ~2 Hz oscillation concomitant with theta activity during ACT epochs. Our finding possibly agrees with Roy et al. 39 , which indicates that this oscillatory activity can be important to HPC-PFC coupling and that the nucleus reuniens of thalamus is crucial for the coupling in this rhythm 39 . We also show that DEA states displays a prominent coupling in delta-gamma oscillatory activity. NMDAr antagonism by S + Ketamine (12,5 mg/Kg) changed the brain state dynamics by (1) increasing the numbers of DEA states, affecting oscillatory brain state alternation; and (2) inducing a distinct DEA state with high gamma power and abnormal cross-frequency coupling between delta phase and high-gamma amplitude. The urethane preparation has been regarded as a valid model to understand sleep physiology 36 and so our data suggest that NMDAr antagonism could interfere with sleep architecture. Interestingly, rodent and human data indicate that subanesthetic doses of ketamine diminishes the REM period of sleep 51 and increases sleep slow-wave activity and NREM sleep 52,53 . However, as indicated by our data, ketamine induces a distinct NREM-like state, marked by an increase in gamma power and aberrant delta-gamma coupling. An interesting question is if ketamine, as seen in our state-space, could also modulate gamma power and frequency coupling during sleep and what are the possible cognitive and behavioural outcomes.
The induction of aberrant gamma oscillations in cortical and subcortical regions is a typical effect of acute ketamine treatment 54 . A variety of in vivo studies have demonstrated that non-competitive NMDAr antagonism increases the firing rate of PFC pyramidal neurons reducing their synchrony, and enhancing broadband gamma activity 40,[55][56][57][58] . This broadband gamma enhancement is believed to be an aberrant and diffuse noise at the network level that can cause dysfunction in cognitive and sensory-motor integration 18,59 . Furthermore, gamma activity increased by ketamine has been observed regardless of the level of consciousness, which indicates that the gamma power increase observed in our study are not particularly of the urethane preparation. In our study, however, gamma increase after KET injection was robust in the PFC, with subtle effects in the HPC (see also Figure S1A). Robust gamma increases have been previously reported in the HPC following NMDAr antagonism in freely moving 16,60 and anesthetized animals 18 as well. To our knowledge, none of these studies tested the effects of S + ketamine, which has a higher affinity to NMDAr compared to racemic ketamine 61 . Indeed, different NMDAr antagonists are known to produce different effects on gamma oscillatory dynamics 55,62,63 .
Interestingly, we show that in DEA states, the increase in high-gamma activity induced by ketamine is strongly coordinated by slow oscillations, rather than by a generalized increase of activity. We suggest that this rhythmic coupling is independent of enhancement in gamma activity since there was low correlation between gamma Scientific RepoRtS | (2020) 10:7167 | https://doi.org/10.1038/s41598-020-63979-5 www.nature.com/scientificreports www.nature.com/scientificreports/ power and CFC. Post-ketamine epochs of similar gamma power to pre-ketamine showed significantly higher delta-gamma coupling. In freely moving animals, it was shown that acute ketamine can increase cross-frequency coupling between gamma activity and theta oscillations in the hippocampus of rats 16 . Chronic ketamine treatment also causes long-term alterations on theta-gamma CFC depending on the behavioral state 64 . Theta oscillations are prominent during exploratory activity 65 and theta-gamma CFC has been related to working memory and memory consolidation mechanisms [66][67][68] .
In contrast, delta-gamma CFC has been described in cortico-striatal network and could reflect cortico-mesolimbic connectivity 63 , and be modulated by dopaminergic activation 69,70 . If the aberrant coupling observed in our experiment is a summation of a possible anesthetic effect of KET under a urethane preparation it remains to be elucidated. Increase of delta-low gamma and theta-low gamma in the frontal cortex has been described as a common neurophysiological feature of several anesthetic agents, including ketamine 71 . However, Pal et al. (2017) observed a lack of effect in the delta-high-gamma CFC and a reduction of high-gamma power using anesthetic doses of ketamine, which suggests that our results are independent of the sedative properties of ketamine, since we used a lower dose than the mentioned study 71 . Furthermore, an increase in delta-high frequency oscillations CFC has been observed with sub-anesthetic doses of ketamine in the striatum of awake animals 60 . One possibility is that the observed aberrant coupling induced by NMDAr antagonism could be associated with an overprocessing of information in the state-space generated by urethane anesthesia (i.e. slow oscillatory activity). However, little is known about the function of delta-gamma coupling in natural physiological states. To further explore a possible relationship with psychotic symptoms, future studies should consider the investigation of ketamine effects on CFC in neocortical circuits during sleep and relevant behavior.
In agreement with other studies, we observed that S + KET administration leads to a sustained increase in HPC-PFC synaptic efficiency and a reduction in PPF. Blot et al. showed that under urethane anesthesia, systemic injection of MK-801 produced an increase of PFC evoked responses induced by ventral HPC stimulation 25 . Several studies suggest that in low-doses NMDAr antagonist acts preferentially in the fast-spiking GABAergic interneurons of the PFC 22,72,73 . Thus, the potentiation of the HPC-PFC evoked potential could be a consequence of synaptic disinhibition on PFC principal neurons. However, Blot et al. (2015) demonstrated by occlusion experiments, that MK-801-induced potentiation shares the expression mechanism with LTP 25 . In this case, the disinhibition of PFC principal neurons would be a trigger for long-lasting synaptic plasticity modifications 25 . Here, we reinforced this hypothesis by showing that prior induction of LTP precluded the S + KET-induced enhancement of the HPC-PFC synaptic efficiency. Consistent with our findings on PPF, Kiss et al. observed a stronger but shorter reduction of HPC-PFC PPF after systemic MK-801 35 . Given the ability of short-term synaptic plasticity to influence information processing, it is plausible that dysfunctional sensory/cognitive processing in schizophrenia may arise from modified short-term synaptic plasticity [74][75][76] . Indeed, short-term synaptic plasticity alterations are a consistent finding in genetic mouse models of schizophrenia (for a review: Ruggiero et al. 77 ; Sigurdsson et al. 78 ).
The main finding of our study is that LTP induction prevents the neurophysiological effects of ketamine. LTP induction attenuated the ketamine effects on cortical fPSP similarly as observed with MK-801 25 . One possibility is that the LTP induction promoted a ceiling effect on the amplitude of the fPSP and the effects of KET could not be observed under this saturation. This does not seem to be the case since in our protocol KET was injected after an initial decreased in the evoked fPSP, and it was showed that NMDAr antagonist effects were blocked after LTP even when the stimulus intensity was lowered to avoid saturation 25 . Since ketamine blocks active NMDAr, and LTP in the HPC-PFC pathway depends on NMDAr activation 19,20 , it is conceivable that LTP could prevent KET effects by direct competitive mechanisms. We also demonstrated that LTP induction in the HPC-PFC projections abolished ketamine reduction of PPF. This is in agreement with Kiss et al. (2011), showing that treatment with the AMPAr allosteric modulators (LY451395) reverts the increase in slow oscillations and PPF deficits produced by MK-801. Classically, PPF arises from calcium accumulation in the presynaptic terminal due to a first stimulation pulse that results in increased neurotransmitter release in response to a second stimulus followed by a short interval 79 . However, short-term synaptic plasticity in the HPC-PFC pathway may encompass a more complex interaction, involving GABAergic interneurons terminating in the PFC pyramidal cells, since CA1 collaterals innervate both pyramidal and GABAergic neurons 80 . Electrical stimulation of CA1 induces a relevant burst activity in PFC interneurons in contrast to the few spikes elicited in pyramidal neurons 81 . This feed-forward inhibition is proposed to constrain the excitatory influence of the HPC on PFC pyramidal cells, supporting rhythmic synchronization between the hippocampus and cortical activity 81 . Thus, we speculate that PPF deficits induced by NMDAr antagonism could contribute to reduced synchrony between HPC and PFC and LTP induction could attenuate this effect. Also, the decrease in HPC-PFC pathway PPF response is related to an increase of delta activity in the 0.5-2 Hz band 35 . Interestingly, we show that ketamine produced an increase in the number of deactivated epochs that are characterized by ~1 Hz power. In contrast, LTP induction prevented this effect on brain state alternation, which could explain the reduction of KET effects on PPF.
Extending the effects of deep brain stimulation, we observed that cortical LTP induction increased gamma activity early after HFS (0-30 min) specifically in DEA epochs. This effect is consistent with previous reports that show gamma increase in the posterior HPC-PFC pathway following LTP, but no LTD 82 . However, we did not observe enhanced CFC following gamma increase. These findings suggest that in contrast to KET, LTP promotes less coordinated gamma oscillations during slow delta activity. Furthermore, HFS attenuated aberrant delta-high-gamma CFC during DEA epochs. NMDAr dependent LTP is divided into different components that do not rely on gene transcriptions: short-term potentiation (STP) and LTP 83 . Importantly, in our protocol KET was injected 30 min after HFS, which corresponds to the time course of STP. STP is induced by presynaptic mechanisms and has been shown in HPC-PFC pathway 84 . STP involves an increase in the probability of neurotransmitter release and the activation of GluN2B and GluN2D containing NMDAr 83 . In contrast, LTP requires mainly activation of GluN2A and GluN2B containing NMDAr. It has been shown that ketamine inhibits STP more potently than LTP 85 and it is possible that is the induction of STP that prevents ketamine effects. This hypothesis Scientific RepoRtS | (2020) 10:7167 | https://doi.org/10.1038/s41598-020-63979-5 www.nature.com/scientificreports www.nature.com/scientificreports/ needs to be addressed in future studies since STP has been speculated as a physiological mechanism supporting working memory 86 . Another interesting hypothesis to be tested is if tetanization would mitigate aberrant oscillatory activity after ketamine injection. Several studies demonstrated that NMDAr antagonists blocked LTP induction (Blot et al., 2015;Rame et al., 2017). Indirect evidence has been provided by Kiss et al. 2011 by showing that LY451395 reverts the increase in slow oscillations and PPF deficits produced by MK-801. These results suggest that AMPAr transmission enhancement could attenuate the effects of ketamine on brain oscillatory activity. However, we were unable to test whether LTP could revert KET effects since the dual HFS protocol spanning ~45 min would surpass the duration of KET effect on oscillatory brain activity as shown to be 30-40 min (Fig. 5).
As a demonstration that prior deep brain stimulation significantly changes the neural dynamics induced by KET, we were able to separate these two brain states (KET and LTP-KET) using an unsupervised algorithm solely based on the electrophysiological features of each state. It remains an open question, however, whether LTP induction would also attenuate functional and behavioral deficits induced by ketamine. Direct evidence supporting hypothesis come from experimental deep brain stimulation (DBS) studies. High-frequency stimulation of the ventral HPC was shown to normalize auditory evoked responses in the MAM model of schizophrenia 87 . Using the same animal model, Perez et al. 2013 showed that ventral hippocampal DBS: (1) normalized aberrant dopamine neuron activity, (2) decreased locomotor response to amphetamine, and (3) restored deficits of cognitive flexibility 88 . Furthermore, the application of DBS to other cerebral regions, such as PFC, nucleus accumbens 89,90 , and medial septum 89 showed promising results for alleviating behavioral deficits in animal models of schizophrenia. Taken together, these data suggest that HFS of limbic circuits should be further investigated as a possible treatment for drug-resistant schizophrenia.

conclusion
Our findings expand previous studies showing that systemic treatment with S + KET produces complex changes in connectivity, synaptic plasticity, and oscillatory patterns in the HPC-PFC pathway in vivo. The prevention of most of these electrophysiological effects through LTP induction supports the idea that NMDAr antagonist effects share common mechanisms with LTP induction. Additional studies are needed to clarify the underlying molecular mechanisms of this interference induced by this form of deep brain stimulation. Additionally, our results suggest that HFS applied to the hippocampus could be a useful strategy to test the attenuation of cognitive impairments in animal models of schizophrenia. We hope these results will contribute to the development of non-pharmacological treatments aimed at preventing or mitigating cognitive deficits associated with psychiatric disorders.  (Fig. 1A). Manufactured electrodes were made of single Teflon-coated tungsten wires (60 μm, AM-Systems) for recording and two twisted wires for bipolar stimulation (~500 μm inter-pole distance). An epidural screw placed in the right parietal bone was used for reference and ground. Temperature (37 ± 0.5 °C) was kept constant during all the procedure by a heating pad.

Subjects.
HPC recording electrode positioning was adjusted by monitoring typical LFP and audio-monitor signals from the hippocampus (i.e., prominent spikes and theta oscillation). The stimulus electrode in HPC was adjusted by applying low-intensity test-pulses (square monophasic pulses, ~150 μA 200 μs, 0.05 Hz) aiming for a consistent fPSP in the PFC (i.e., the latency of first negative peak of 14-17 ms and amplitude>0.25 mV 91,92 . After electrode adjustment, an input-output curve (I/O curve; 60-500 μA) was used to establish the current intensity necessary to evoke 70% of the maximal fPSP amplitude. This current was used to apply paired monophasic pulses (same parameters as in test pulses and 80 ms of inter-pulse interval; S88; Grass Instruments) during the entire experiment. Recorded signals consisted of evoked fPSP and concomitant LFPs of PFC and HPC (Fig. 1C). Signal was amplified 100×, band-pass filtered (0.3-1000 Hz; Grass), and digitized at 10 kHz (ADInstruments). LTP induction was induced by applying a high-frequency stimulation (HFS): two series of 10 trains (50 pulses at 250 Hz every 10 s) separated by 10 min. Figure 1B illustrates the experimental design. In Experiment 1 we investigated the ketamine effects on HPC-PFC connectivity. fPSPs and LFPs were monitored for 120 min after ketamine (S(+)-ketamine; 12.5 mg/Kg ip) or saline (0.9%) injection and compared with the 90 min baseline (groups KET, n = 7 and SAL, n = 8, respectively). We choose S( + )-ketamine given its high affinity for NMDAr 61 and because it reproduces the metabolic effects observed in psychotic patients 93 . In Experiment 2 we explored the hypothesis that LTP induction could prevent ketamine effects on the HPC-PFC connectivity. After a 30 min baseline, two HFS protocols Scientific RepoRtS | (2020) 10:7167 | https://doi.org/10.1038/s41598-020-63979-5

Experimental design.
www.nature.com/scientificreports www.nature.com/scientificreports/ were applied at 30 and 60 min. Following HFS, ketamine or saline was injected and field potentials were monitored for an additional 120 min (groups LTP-KET, n = 9 and LTP-SAL, n = 6, respectively). All experiments were conducted during urethane anesthesia. Data analysis. All data processing was performed using customized scripts in Matlab (Mathworks). The amplitude of evoked fPSP1 and fPSP2 (Fig. 1C) were normalized as a percentage of baseline mean (90 and 30 min for experiments I and II, respectively). PPF was calculated as the ratio of fPSP2 and fPSP1 as an indication of short-term synaptic plasticity. All fPSP measures were averaged in blocks of 10 min. LFP signal was re-sampled to 1000 Hz and then high-pass filtered at 0.5 Hz. The whole recording data was epoched in a 20 s period following HPC electrical stimulation. Time windows of 0.5 s containing the evoked fPSP and electrical stimulation artifact (in both regions) were eliminated from all epochs.
Brain state classification. We classified deactivated and activated periods by plotting PFC epoch values for RMS (root mean square) and the number of zero-crossing in the LFP (Fig. 1D). These measures directly reflect amplitude and presence of faster rhythms in the signal (Fig. 1C,D). We used a k-means algorithm (squared Euclidean distance for three groups) for an initial clustering and manually refined the classification eliminating epochs at the cluster edge. Spectrum content of each state was analyzed to confirm classification. Not classified epochs were not analyzed in this work. Spectral analysis. Power spectral density estimates (PSD) were calculated using Welch's method in which Discrete Fourier Transform (FFT Matlab algorithm) is applied in overlapping windows and the periodogram is calculated for each segment individually and the magnitude squared result of the FFT is averaged. We used 3 s Hamming tapered windows, with 50% overlap and a 2 12 points FFT. PSDs estimates were then averaged over trials and animals. For representative spectrograms, we used a short-time Fourier Transform in the whole recording using a 60 s window with 2 16 points FFT and 50% overlap. For statistics, power was integrated into specific frequency bands (delta: 0.5-2 Hz, theta: 3-5 Hz, low-gamma: 30-55 Hz, and high-gamma: 65-100 Hz). Relative power was obtained by dividing PSD estimation by the integrated power over all frequencies.
Spectral coherence was estimated using Welch's periodogram method to compute the cross-PSD of PFC and HIPO (Pxy) and the PSD of both region (Pxx and Pyy). The magnitude squared coherence was calculated as: Cxy(f) = |Pxy(f)| 2 /Pyy(f) Pxx(f). The parameters used were the same as described for the PSDs estimate. Coherograms were calculated using the same approach while using a moving window of 90 s and 50% overlap. Power and coherence values from 58-62 Hz (line noise contamination) were removed from the analysis. To evaluate HFS effects on spectral parameters we combined all LTP animals (LTP-SAL and LTP-KET groups) into one group (LTP group) since HFS was the only manipulation from the second electrical stimulation until drug injection.
Phase-amplitude coupling. Cross-frequency coupling was estimated by the modulation index as described by Tort and adapted in previous work from our group 29,94 . Briefly, comodulation maps were constructed applying Hilbert transform to the signal filtered in bins of 0.5 Hz from 1 to 20 Hz on steps of 1 Hz for the phase modulating signal, while the amplitude modulated signal was filtered in bins of 1 Hz from 10 to 120 Hz on steps of 5 Hz. Shannon entropy of the distribution of mean amplitudes per phase (divided into 18 bins) in each frequency bin was calculated to obtain the cross-frequency modulation index (MI) for each period. The MI between delta oscillations (1-2 Hz) and high gamma band (65-100 Hz) was calculated for comparisons. Directionality analysis. We inferred directionality by cross-correlation and Granger causality.
Cross-correlation measures the similarity of two time series by performing the sliding dot product between the signals. We calculated the peak lag of the cross-correlation of delta and theta filtered signals for deactivated and activated epochs, respectively. Wiener-Granger causality spectra were performed using the MVGC toolbox developed by Barnett and Seth 95 , which is freely available online (http://users.sussex.ac.uk/~lionelb/MVGC/). This algorithm uses vector auto-regressive models to estimate prediction of a time series A based on another time series B compared with the prediction obtained by using the past values of time series A alone. We used pairs of HPC and PFC LFP separating for DEA or ACT epochs. Initially, raw LFP was decimated to 200 Hz, and the model order was estimated by Akaike Information Criterion for each animal separately (model order range: [36][37][38]. We fixed the model order of 40, which gave an adequate frequency resolution for the slow oscillations that predominate on our signals with reasonable computation cost. For statistical significance, we calculated the 95% confidence interval (CI) from the empirical null distribution of the frequency-domain Granger estimates, based on randomly permuting one LFP in bins the size of the model order. For directionality comparison (HPC → PFC vs. PFC → HPC; in DEA or ACT) we used a Bonferroni-corrected paired t test for the peak frequency in the delta (0.5-2 Hz) or theta (3-5 Hz) bands.
Principal component and discriminant analysis. Principal component analysis, using singular-value decomposition (pca Matlab function), was used for dimensionality reduction in order to find patterns of variance among multivariate data. Variables analyzed were MI, high-gamma, delta and theta power and fPSP amplitude. We used these features since they were associated with ketamine effects (Fig. 3). LFP data were normalized by the mean value of baseline DEA epochs, while for fPSP we used the normalized mean value for the 10 min before drug injection (all epochs regardless of state classification to account for LTP induction). The normalized values were extracted in 5 min epochs from 10-40 min after drug injection (KET and LTP-KET groups) in order to capture stable effects of ketamine. Data were z-scored for each variable. Following PCA analysis, data was projected against principal components (PCs) and the mean score was compared between conditions (KET vs. LTP-KET). Correlation coefficients between original variables and scaled components were obtained Scientific RepoRtS | (2020) 10:7167 | https://doi.org/10.1038/s41598-020-63979-5 www.nature.com/scientificreports www.nature.com/scientificreports/ by multiplying eigenvectors by the square root of the eigenvalues. For interpretation, we used correlation coefficients value > 0.5 96 . A discriminant analysis classifier was used based on the first three principal components extracted (explained variance of 79.21%). A quadratic function fit was heuristically determined and a 50 fold cross-validation was performed with 81 epochs for training data and 9 epochs for data test. Cross-validation of the quadratic discriminant model using PC1-PC3 resulted in 85.56% accuracy, which was better than the 83.33% cross-validation accuracy of the discriminant fit using all the 5 dimensions of the original data.
Statistics. Normal distribution was evaluated in all data sets using Kolmogorov-Smirnov test. fPSP data were analyzed by two-way ANOVA with repeated measures and Bonferroni post-hoc test to compare treatment over time. For LFP (power and CFC) we used paired t tests for within-group comparisons, and unpaired t tests for between-group comparison or Wilcoxon matched-pairs signed rank and Mann-Whitney test, respectively, as non-parametric equivalents. For comparison of more than two conditions we used one-way ANOVA with tukey-kramer post-hoc or Kruskal-Wallis test with Dunn's post-hoc for non-Gaussian distribution. Significance of spectral coherence and MI were estimated calculating a CI for a surrogate data using a bootstrapped shuffled data (8000 iterations). Probability density function after drug injection was estimated calculating the presence of DEA classified epochs for all the animals in Saline or ketamine groups. Pearson's correlation coefficient was calculated to investigate linear dependency between CFC and LFP power. Data are expressed as the mean ± standard error of the mean (SEM) for bar and line plots, and for Box-plot data are expressed as 1 st quartiles, medians, and 3 rd quartiles, with whiskers representing minimum and maximum values. The significance level was set to 0.05.

Histology.
To confirm electrode positioning we performed an electrolytic lesion (1 mA, 1 s) at the end of the trial. After an additional dose of anesthesia animals were decapitated and had their brains removed and placed in solutions for fixation (10% formaldehyde in phosphate-buffered saline, PBS) and cryoprotection (20% sucrose in PBS). Coronal sections (30 μm) stained with cresyl violet were evaluated through a bright-field microscope (Fig. 1A).

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.