1. Introduction
Contextual fear conditioning (CFC) is a form of Pavlovian learning in which an association is made between non-threatening environmental stimuli and painful, dangerous, or threatening stimulation [
1,
2]. Studies of Pavlovian fear conditioning and extinction in rodents and humans have suggested that a neural circuit including the hippocampus, amygdala, and medial prefrontal cortex is involved in the learning and memory processes that enable context-dependent behavior [
3]. Recently, the roles of other parts of the brain in fear learning, responses, and extinction have been more extensively studied [
4,
5,
6,
7]. Although an early study showed that thalamic lesions did not prevent the acquisition of conditioned auditory fears [
8], later studies found that large lesions in the posterior thalamus (including posterior intralaminar nuclei) prevented the acquisition of fear-potentiated startle [
9,
10]. Recently, the thalamus (midline and intralaminar nuclei) and its role in conditioned fear behaviors have received more attention. A part of the midline thalamus, mediodorsal thalamic nucleus (MD), has been implicated in the control of fear memory processes [
11,
12]. Specifically, it has been shown that the tonic firing in MD neurons contributes to extinction learning, whereas increased burst firing has a suppressive effect on fear extinction [
11,
12]. The paraventricular thalamus (PV, part of the midline thalamus) has been established as having an important role in the consolidation and retrieval of fear memories [
13,
14,
15,
16]. The PV–amygdala pathway constitutes a circuit that is essential for both the establishment of fear memories and the expression of fear responses [
13,
14]. The posterior intralaminar nucleus of the thalamus has been studied in the context of pain induced by unconditioned stimuli, mainly because it projects to the basolateral amygdala [
17].
The central medial nucleus of the thalamus (CMT), as a part of the intralaminar nucleus of the thalamus, is well-positioned to modulate fear-conditioning responses. The CMT has connections with two key structures that are important for fear behaviors—the prefrontal cortex and amygdala [
18]. Additionally, the CMT is rich in Ca
V3.1 T-type calcium channels (T-channels) and can modulate different behaviors through dual tonic and burst firing modes [
19,
20]. In comparison to high voltage-activated (HVA) calcium channels, T-channels need weaker depolarization to open and can form “window” currents around the resting neuronal membrane potential [
21,
22]. T-channels are de-inactivated during neuronal hyperpolarization, can be opened after depolarization, and are known to elicit low threshold spikes (LTSs) with characteristic rebound burst firing patterns [
23]. Due to these unique properties, T-channels have an important role in regulating neuronal excitability, synaptic plasticity, and oscillatory behaviors [
24,
25]. T-channel dysfunction has been implicated in sleep disorders, anesthetic effects, the absence of epilepsy, pain, and neurological, neuropsychiatric, and cognitive disorders [
22,
26,
27,
28,
29,
30].
Recent studies have shown that pharmacological targeting of T-channels can also regulate the acquisition and recall of conditioned fear in animals [
31,
32]. In a genetic model with impaired fear extinction, enhanced burst firing in the thalamus was accompanied by increased T-currents, and pharmacological blockers of T-channels rescued the fear extinction deficiency [
12]. These results suggest that T-channels may contribute to the neural systems that mediate the learning and memory of conditioned fear. However, the specific roles of the different thalamic T-channel isoforms in conditioned fear behavior are not clear. As we previously reported that Ca
V3.1 T-channels are essential for CMT excitability [
19,
20], in this study, we set out to investigate the role of these channels in fear behavior using global knock-out (KO) and CMT-specific knock-down (KD) mice, thereby establishing thalamic T-channels as a potential novel target for the treatment of affective disorders.
3. Discussion
In this study, we demonstrated for the first time that by silencing thalamic CMT CaV3.1 T-channels, we can modulate fear expression, thus suggesting that CaV3.1 T-channels could be a potential cellular target for disorders where there is an altered fear response such as post-traumatic stress disorder, schizophrenia, anxiety, or substance abuse.
In the context of fear conditioning, the thalamus is considered to play a crucial role as a sensory relay station, receiving sensory information from the periphery and transmitting it to the amygdala, which is the primary brain region responsible for processing fear. In addition, recent studies have implicated non-sensory parts of the thalamus (midline and intralaminar nuclei of the thalamus) in fear-conditioning behaviors [
13,
15,
16,
17]. The midline and intralaminar thalamic nuclei are well-placed to contribute to fear behaviors due to the extensive projections that they receive from the periaqueductal gray while projecting to the prefrontal cortex and amygdala [
17,
34]. In addition, most of the midline and intralaminar thalamic nuclei are part of the thalamo-hippocampal circuitry and can therefore influence memory-related behaviors [
35]. Specifically, animal studies have shown that the PV regulates fear processing through selective inactivation of the lateral division of the central amygdala [
13] and by tuning the inhibitory functions in the prefrontal cortex [
15,
16]. Additionally, conditioned stimulus–unconditioned stimulus presentations caused increased c-Fos expression in the prefrontal cortex, midline, and intralaminar thalamus (including the PV and CMT), the lateral amygdala, and retrograde-labeled midline thalamic afferents to the prefrontal cortex [
13,
36].
Fear conditioning is one of the most frequently used experimental procedures not only for modeling learning and memory but also for anxiety disorders. The conditional freezing response is used as the output for most fear-conditioning experiments in rodents. Here, we showed that global Ca
V3.1 KO animals froze more at baseline and immediately after each tone/shock pairing in the acquisition–training phase. The higher baseline freezing in Ca
V3.1 KO mice implies that these mutant animals have different emotional responses to novel and non-aversive environments, at least in the first few minutes. This baseline difference in freezing can partially explain the stronger freezing behavior related to tones/shocks observed in the mutant animals. Additionally, there is a possibility that the brief increase in freezing of global KO animals after each tone/shock pairing was due to the known role of thalamic Ca
V3.1 T-channels in auditory processing [
33]. Ca
V3.1 T-type calcium channels play a significant role in auditory perception. They contribute to the electrical activity of neurons in the auditory pathway, particularly in the thalamus, which is crucial for regulating the firing pattern of neurons, impacting how sound information is transmitted to the auditory cortex. Hence, they can potentially influence aspects like sound localization and the perception of sound intensity [
33]. It has been shown that the switch between firing modes depends on thalamic Ca
v3.1 T-type calcium channels and that pharmacologic and/or genetic inhibition of these channels in the auditory thalamus substantially influences auditory processing [
33]. On the other hand, based on our findings, we propose that an increase in freezing behavior after each tone/shock pairing was not due to different pain thresholds as our data showed no pain threshold differences between WT and Ca
V3.1 KO animals, but it could be present due to an altered emotional processing of pain.
As Ca
V3.1 T-channels are not only abundantly expressed in the thalamus but also in the cortex and other brain regions important for pain processing and modulation [
37], it is not surprising to see different fear responses in global Ca
V3.1 KO animals after aversive stimuli, i.e., control (WT) mice exhibited normal freezing behaviors when re-exposed to the context associated with the aversive stimulus, whereas Ca
V3.1 KO mice had stronger freezing responses, suggesting the existence of altered fear expression in these mutant animals. Similarly, CMT-specific Ca
V3.1 KD animals spent more time freezing than the control animals during the testing phase, emphasizing the important role of thalamic Ca
V3.1 T-channels in fear expression.
Previously, we showed that the Ca
V3.1 T-channel subtype is the dominant isoform in the CMT and that CMT neurons from Ca
V3.1 KO animals have a reduced ability to burst fire [
20]. Ca
V3.1 T-channels are important in the burst-firing mode of thalamic neurons [
19,
20], and here, we were able to modulate the fear-related responses in Ca
V3.1 KD animals by reducing T-currents and diminishing the ability of CMT neurons to fire in the burst mode. Interestingly, CMT-specific Ca
V3.1 T-channel silencing did not affect the fear behaviors during training, as was seen with the global Ca
V3.1 KO animals, suggesting a more selective effect of thalamic Ca
V3.1 T-channels on fear expression. Our data are supported by previous findings reporting that the injection of a T-channel blocker before the acquisition of fear conditioning significantly increased freezing in rats [
31]. These authors commented that the enhanced fear behavior observed prior to acquisition and extinction, as well as during acquisition, could reflect an increase in anxiety after the systemic T-channel blocker injection [
31]. This could also explain the enhanced fear behavior we observed in Ca
V3.1 KO mice during acquisition and fear expression. However, although the animals with global deletion of Ca
V3.1 T-channels may have hippocampal-dependent learning and memory deficits, the general locomotor activity and anxiety-related behaviors of Ca
V3.1 KO male mice were unchanged [
38]. Additionally, in our experiments, we did not observe changes in the locomotor activity and anxiety-related parameters of the thalamic Ca
V3.1 KD animals in comparison to the control mice, suggesting that the enhanced freezing behaviors during testing cannot be attributed to an increase in anxiety in these animals.
As thalamic Ca
V3.1 T-channels are important for auditory processing [
33], we used a protocol with a reduced number of tone/shock pairings in CFC experiments. The CFC behavioral paradigms with five or more tone/shock pairings typically have stronger and persistent fear memory, often accompanied by enhanced generalization to similar cues and contexts [
39]. In contrast, two pairings, as in our experiment, usually produce weaker fear responses, which facilitate the assessment of subtler cognitive processes such as discrimination, extinction, or precise neural circuit involvement [
39]. Thus, the number of pairings directly influences the intensity and stability of conditioned fear responses and may impact both behavioral outcomes and neurobiological mechanisms studied. Although in our experiments all the animals had a higher freezing percentage after each tone in a different context (Day 2), we did not observe differences in freezing between the groups related to the tone. This suggests that Ca
V3.1 T-channels may be more important for contextual conditioning to an unpleasant stimulus in comparison to the tone-related fear or that our conditioning to the cue was too subtle. As this could be a limitation of our study, further experiments are needed to fully investigate the role of CMT Ca
V3.1 T-channels on cue-related fear conditioning and the processing of fear extinction.
It is known that the shRNA approach, while effective for gene knock-down, may have potential off-target effects [
40]. However, in our experiments, knocking down Ca
V3.1 T-channels only affected rebound burst firing but not tonic firing. This strongly suggests that, as we used the shRNA approach, we did not significantly affect other major targets that control neuronal excitability. Additionally, our experimental method is justified as the use of shRNA to manipulate genes in different tissues is a commonly used approach in animal studies and shRNA has enormous potential as a precision-based therapy to treat numerous neurodegenerative disorders [
41,
42].
Using CMT-specific Ca
V3.1 KD mice, we found that fear expression can be modulated by reducing the ability of thalamic CMT neurons to fire in burst mode. Similarly, it has been reported that the firing mode of the MD thalamus is critical for the modulation of fear extinction [
12]. These authors showed that the mutant mice exhibited impaired fear extinction, accompanied by enhanced burst firing and increased T-currents in MD neurons. Importantly, the abnormal fear extinction was rescued by pharmacologically blocking T-channels in vivo [
12]. Together with our data, these findings highlight the key role of midline and intralaminar thalamic T-channels in modulating fear expression and extinction. We conclude that, through regulating T-channel activity, we can modulate the firing mode of thalamic neurons and consequently affect fear-mediated responses.
4. Materials and Methods
4.1. Animals
The experimental procedures with animals were performed according to the guidelines approved by the University of Colorado Anschutz Medical Campus. All efforts were made to minimize animal suffering and to only use the number of animals necessary to produce reliable scientific data. All experiments were conducted during the light cycle in male adult (2–5-month-old) C57BL/6J wild-type (WT) and Ca
V3.1 knock-out (KO; Ricken BioResources Centre, Japan) mice, as it is well-documented that the estrous cycle can have a significant effect on associative learning [
43]. All researchers were blinded during the behavioral and ex vivo electrophysiological experiments.
4.2. Generation of CMT-Specific CaV3.1 Knock-Down (KD) in Mice
The use of global Ca
V3.1 knock-out (KO) mice is useful for proof-of-principle confirmation, but the widespread distribution of these channels and possible compensatory modulations is a major disadvantage of this approach. To specifically investigate the contribution of Ca
V3.1 T-channels in the thalamus, we used mice with CMT-specific silencing of the Cacna1g gene, which encodes the Ca
V3.1 channel. Our method consisted of knocking down Cacna1g by injecting short hairpin RNA (shRNA) using a procedure described elsewhere [
38,
44]. Adult male C57BL/6J mice were anesthetized with 1–2% isoflurane, given 2.5 L/min of oxygen, and transferred to a standard stereotaxic frame (Kopf Instruments, CA, USA). The mice were randomized into two treatment groups: scrambled control (AAV2-GFP-U6-scrmb-shRNA; titer: 1.1 × 10
13 GC/mL; Vector Biolabs, Malvern, PA) or Cacna1g shRNA targeting Ca
V3.1 (AAV2-GFP-U6-mCACNA1G-shRNA; titer: 6.8 × 10
12 GC/mL; Vector Biolabs) groups. The shRNAs were injected into the CMT [anteroposterior (AP): −1.35 mm; mediolateral (MD): 0; and dorsoventral (DV): −3.6 mm] using high-titer AAV2 (0.5 μL, using a 5 μL Hamilton syringe at a rate of 0.1 µL/min). The animals were allowed to recover for at least three weeks before the experiments (behavioral or ex vivo brain slice recordings) were performed to ensure adequate virus expression. Although we aimed to target just CMT neurons in our experiments, we are aware that other parts of the thalamus could have been partially transfected (mostly the IMD, PV, and MD;
Figure 1B). We excluded all animals that did not have the highest Cacna1g shRNA or scrambled shRNA expression in GFP-positive neurons in the CMT in comparison to other areas. The WT animals were randomly assigned to the different groups such that animals from the same litter were used for the generation of both control and Ca
V3.1 KD animals. Animals from at least 3 different litters were used for the experiments.
4.3. Ex Vivo Brain Slice Preparation
The mice were briefly anesthetized with 5% isoflurane and decapitated. Their brains were quickly removed and placed in a cold (4 °C) oxygenated (95 vol% O2 and 5 vol% CO2) solution. Live 250 μm thick coronal brain slices were sectioned at 4 °C using a vibrating microslicer in a cold solution containing 260 mM sucrose, 10 mM D-glucose, 26 mM NaHCO3, 1.25 mM NaH2PO4, 3 mM KCl, 2 mM CaCl2, and 2 mM MgCl2 (Laica VT 1200S). The brain slices were immediately incubated for 30 min in a solution containing the following at 37 °C: 124 M NaCl, 10 mM D-glucose, 26 mM NaHCO3, 1.25 mM NaH2PO4, 4 mM KCl, 2 mM CaCl2, and 2 mM MgCl2. The electrophysiology experiments were carried out at room temperature. During the incubation, the slices were constantly perfused with a gas mixture of 95 vol% O2 and 5 vol% CO2.
Whole-cell recordings of CMT neurons were performed using Zeiss optics (Zeiss AXIO Examiner D1, ×40 objective). Glass microelectrodes (Sutter Instruments, borosilicate glass with filament OD of 1.2 mm) were pulled using a Sutter Instruments P-1000 model and fabricated to maintain an initial resistance of 4–6 mΩ. Neuronal membrane responses were recorded using a Multiclamp 700 B amplifier (Molecular Devices, Foster City, CA, USA). The voltage–current commands and the digitization of the resulting voltages and currents were performed using Clampex 8.3 software (Molecular Devices) running on a PC-compatible computer. The resulting current traces were analyzed using Clampfit 10.5 (Molecular Devices). The statistical and graphical analyses were performed using GraphPad Prism 9.0 (GraphPad Software) or Origin 7.0 (OriginLab) software. The results are presented as the mean ± SEM unless stated otherwise.
4.4. Electrophysiology Experiments Using Brain Slices
The external solution for current-clamp and voltage-clamp electrophysiology experiments consisted of the following: 125 mM NaCl, 25 mM D-glucose, 25 mM NaHCO
3, 1.25 mM NaH
2PO
4, 2.5 mM KCl, 1 mM MgCl
2, and 2 mM CaCl
2. For the current-clamp experiments, the external solution contained the synaptic blockers picrotoxin (20 μM), D-2-amino-5-phosphonovalerate (D-AP5; 50 µM), and 2,3-dihydroxy-6-nitro-7-sulfamoyl-benzo [f]quinoxaline-2,3-dione (NBQX; 5 µM). For the voltage-clamp experiments measuring T-type calcium channels, tetrodotoxin (TTX; 1μM) was added to the extracellular medium as a blocker of voltage-dependent sodium currents. The internal solution for the current-clamp recordings consisted of the following: 130 mM potassium-D-gluconate, 5 mM ethylene-glycol-bis (β-aminoethylether)
N,N,N′,N′-tetra acetic acid (EGTA), 4 mM NaCl, 0.5 mM CaCl
2, 10 mM HEPES, 2 mM Mg-ATP, and 0.5 mM Tris-GTP, pH 7.2. For the T-channel recordings, we used an intracellular cesium-based internal solution containing the following: 110 mM Cs-methanesulfonate, 14 mM phosphocreatine, 10 mM HEPES, 9 mM EGTA, 5 mM Mg-ATP, and 0.3 mM Tris-GTP, with the pH adjusted to 7.15–7.20 using CsOH (standard osmolarity: 300 mOsm) [
19]. For the voltage-clamp experiments, maximal inactivation curves were generated using a standard double-pulse protocol with 3.6 s long prepulses (voltage: −120 mV) and test potentials (−50 mV). For additional voltage-clamp experiments, activation (IV) curves were generated using a standard protocol where the holding potential (Vh) of −90 mV was increased to depolarized test potentials (Vt) from −80 to −40 mV in 5 mV increments. The voltage dependencies of activation were described using a single Boltzmann distribution: G (V) = G
max/ (1 + exp [-(V − V
50)
/k]), where G
max is the maximal conductance (calculated by dividing the current amplitude by the estimated reversal potential), V
50 is the voltage at which half of the current is activated, and
k represents the voltage dependence (slope) of the distribution. The average current amplitudes were analyzed using the appropriate statistics.
In the current clamp experiments, both the stimulated tonic and burst-firing properties of the CMT neurons were characterized using multistep protocols. To investigate the stimulated tonic firing patterns in CMT cells, we injected a depolarizing current pulse through the recording pipette with a 400 ms duration in 25 pA increments, starting from 50 pA. To investigate the rebound burst-firing patterns, the neurons were injected with hyperpolarizing currents in 25 pA intervals from 0 to −225 pA. The subsequently stimulated tonic action potential (AP) frequencies, rebound burst-firing thresholds, and APs in the rebound burst were determined. The resting membrane potential was measured at the beginning of each recording and the liquid junction potential was not corrected. All recordings were obtained using a −60 mV membrane potential.
4.5. Quantitative Real-Time PCR
After anesthetizing the animals with isoflurane and removing the brain, each brain was sliced on a vibratome (Leica VT 1200S) in ice-cold PBS, creating 500 μm thick horizontal slices. CMT tissue punches with GFP-positive labeled neurons were collected under a dissecting fluorescence microscope. Samples from 3 animals from the Cacna1g shRNA group and 4 animals from the control (scrambled shRNA) group were pooled and the samples were stored at -80 °C until further use. RNA was isolated using the RNeasy Microarray Tissue Mini Kit with QiAzol (QIAGEN), and quantitative real-time reverse transcription (qRT)–PCR was performed on a BioRad Icycler using the RT2 First Strand Kit and RT2 SYBR Green qPCR Mastermix (QIAGEN) according to the manufacturer’s protocols. The primers for CaV3.1 and cyclofilin (internal control) were purchased from QIAGEN (CACNA1G NM_031601 catalog number: PPR52633A-200; cyclophilin NM_001004279 catalog number: PPR59729A). The qRT-PCR data were analyzed as previously described in [
45].
4.6. Behavioral Experiments
The investigators were blinded to the genotype in all the behavioral experiments. All tests were performed by trained researchers and monitored and analyzed using specialized software (Ethovision XT system, Noldus, Wageningen, The Netherlands). The behavioral experiments involving virus-injected mice began at least three weeks after injection to allow for adequate expression levels. The animals were habituated to the testing room for at least 30 min before all behavioral experiments.
For the contextual fear conditioning test, the mice were placed into conditioning boxes that were 30.5 cm × 24.1 cm × 29.2 cm in size. Inside the box was a metal grid floor, a house light, and an odor cue (70% ethanol, which was also used to clean the apparatus). On the first day of training, the mice were placed in the chamber for 7 min, during which 2 tones (30 s in duration) paired with shocks (1 s, 0.7 mA) were presented. On the second day, the mice were placed in a different environment for 7 min, and just the tone (30 s, 2 times) was present without a shock. On the third day, the mice were placed into the same box as Day 1 (same box, but no shock or tone applied), and the mouse’s behavior was monitored for 7 min. Fear behavior was analyzed and presented as the percent of time spent in a freezing state. The conditioning boxes were cleaned with ethanol after each animal, making sure that the ethanol had evaporated completely before placing a new mouse in the chamber. The cleaning was performed the same way during the experiments (Days 1–3) so that the odor cue was always constant.
For the open-field test, the mice were placed in the center of the arena (gray box, 44 cm × 44 cm × 30 cm) and tested for 10 min without prior habituation. This test provided important information about the general locomotor activity (total distance traveled), as well as the anxiety-like behaviors (number of entries and time spent in the center of the arena) of the mice.
The elevated zero maze is a modification of the elevated Plus Maze test and it is used to investigate anxiety-like behaviors in mice. It has an elevated ring-shaped runway (50 cm diameter, 5 cm width, and 20 cm wall height) with equal areas for the adjacent open and closed quadrants. The animals were placed in the center of the open quadrant and tested for 10 min. Similar to the open-field test, the test gives information about general locomotor activity (total distance traveled), as well as anxiety-like behaviors (number of entries and time spent in the center quadrat) of mice.
The open-field and zero-maze apparatuses were cleaned with ethanol after each animal, making sure that the ethanol had evaporated completely before placing a new mouse in the chamber. Cleaning was performed the same way during all the behavioral experiments.
4.7. Data Analysis
In every ex vivo electrophysiological experiment, we attempted to obtain as many recordings as possible from the neurons to minimize the number of animals used. The statistical analysis was performed using two-way repeated measure (RM) ANOVA, as well as Student’s unpaired two-tailed t-test where appropriate. For significant interactions between factors according to two-way RM ANOVA, Sidak’s post hoc comparisons were used. p values < 0.05 were considered significant. Additionally, Cohen’s d and effect size r were reported. The statistical and graphical analyses (presented as the mean ± standard error of the mean) were performed using the GraphPad Prism 10 software (GraphPad Software, La Jolla, CA, USA) and the traces were prepared using Origin 2018 (OriginLab, Northampton, MA, USA). The schematic presentations of the CFC, open-field, and zero-maze tests were created using bioRender (BioRender.com).