Contrary to common belief, T_b is not fixed but varies with arousal state, motor activity, and environmental temperature (Kiyatkin EA, 2019; Hankenson et al., 2018; Wang et al., 2014). Using a sterile IPTT-300 implantable temperature transponder (Bio Medic Data Systems, LLC) with a 0.025-second time constant and 0.1°C accuracy, T_b measurements were recorded every five-minutes over six hours from 2–3-week-old wildtype mice in an open-field arena at room temperature (Figure 1A). This analysis revealed a stable median temperature of 36.5°C with a maximum and minimum fluctuation of 0.5°C (Figure 1A). However, within-subject analysis showed brief T_b elevations into the fever range (shaded red bars) during increased activity like cage climbing or digging (Figure 1B-C). Exposure to infrared light elevates ambient temperatures and induces prolonged T_b elevations into the fever range, with an onset within minutes and a duration of several hours (Figure 1D).

To investigate how postnatal cortical neurons respond to fever-range temperature, we used whole-cell current clamp electrophysiology to record natural spiking, evoked synaptically, as brain slice temperature gradually rose from 30°C to 36°C to 39°C in PNs (Figure 2A-B). Acute S1 cortex slices (350 μm thick) were prepared using standard methods (Antoine et al., 2019). PNs were visually identified via infrared DIC optics, and physiological verification for regular spiking was done in current clamp. Neuronal activity in PNs was evoked by stimulating in the barrel center of cortical layer (L) 4 using a bipolar electrode (0.2 ms pulses) at a specific stimulation intensity i.e. 1.4X Eθ (Figure 2A). Eθ is defined as the minimal intensity evoking an excitatory postsynaptic current (EPSC) during more than 3 of 5 consecutive sweeps with 10 s inter-sweep interval. Eθ was determined in voltage clamp for each recorded cell, prior to recording post-synaptic potentials (PSPs) and spiking in current clamp. L4-evoked PSPs and spiking, were recorded from L2/3 PNs, with resting membrane potential (V_m) set just below spike threshold (-50mV) to simulate in vivo conditions (Yamashita et al., 2013).

Recordings of synaptically-evoked activity in L2/3 wildtype cortical excitatory PNs from P12-14 mice showed significant effects of temperature elevations from 30°C to 39°C. These changes resulted in an increased depolarization required to reach spike threshold (ST) (defined as the minimal V_m that just elicits spiking) (Figure 2B-C). Consistent with previous findings, measurements of input resistance (R_in) which reflect the extent to which membrane channels are open, showed significant reductions across all recorded neurons with increasing temperature (Figure 2D). Similar temperature-induced effects on R_in have been reported in L2/3 PNs in the rat visual cortex (Hardingham and Larkman, 1997) and in hippocampal CA1 and CA3 PNs from P13-P16 mice (Kim and Connors, 2012). The depolarized ST and reduced R_in in PNs would likely decrease cell excitability and spiking. Indeed, the percentage of spiking cells decreased during temperature elevations from 30°C to 39°C (Figure 2E). However, contrary to expectations, average spiking remained remarkably stable (Figure 2F-G). We also examined the temperature-induced effects on L2/3 PNs in the absence of excitatory and inhibitory synaptic inputs. Intrinsic spiking was measured in the presence of synaptic blockers (100 μM APV, 10 μM NBQX, and 3 μM gabazine). PNs were held at -80 mV, followed by positive current injection (500 ms) steps to determine rheobase (the minimum current injection needed to produce one spike). At 39°C, rheobase was increased compared to 30°C (Figure 2H). Current was injected up to 250 pA above rheobase at 30°C, 36°C, and 39°C, using a new brain slice for each recording. Although rheobase was decreased at 39°C compared to 30°C, we found that the mean number of spikes produced did not differ across current steps (Figure 2I).

This stability in average spiking at 39°C, despite the loss of some previously active neurons, suggests the presence of temperature-adaptive changes in PNs that help to maintain normal circuit activity levels. These changes may either initiate firing in previously non-spiking neurons or either maintain or elevate spiking levels in currently active cells. Interestingly, in P7–8 mice, temperature elevations from 30 °C to 39 °C produced effects on ST and R_in similar to those observed in P12–14 mice, with ST increasing and R_in decreasing (Figure 3A–D). However, in P7-8 PNs, during temperature elevations from 36 to 39°C, ST did not increase, and R_in decreased (Figure 3C-D). Additionally, the percentage of spiking PNs decreased during temperature elevations from 30°C to 39°C (Figure 3E). In contrast to P12-14 PNs, the average spiking levels in P7-8 PNs significantly decreased during these temperature elevations (Figure 3F-G). This suggests that temperature-adaptive changes, possibly involving R_in and ST, which help maintain average spiking levels at fever temperatures despite the loss of some previously active neurons, may be underdeveloped or absent at this earlier age.

The temperature-induced increases in ST at P12-14, would make it harder for a cell to spike, as larger depolarization levels would be required to maintain spiking (Figure 4A). Therefore, we hypothesized that in cells with increased ST that remain spiking, a temperature-adaptive change might involve a proportional increase in depolarization. To determine if larger depolarizations were associated with higher (less negative) STs, we analyzed the relationship between peak postsynaptic potential (PSP), which indicates maximum depolarization, and ST (Figure 4B).

Correlation analysis revealed a strong positive correlation between peak PSP and spike threshold at P12-14, which was absent at P7-8 and P20-23 (Figure 4C). This correlation was specific to these parameters, as it did not occur for comparisons involving input resistance, an intrinsic property which is affected by temperature changes (Figure 4D-E).

One mechanism that could promote enhanced depolarization is a reduction in the inhibitory PSP (IPSP) (Figure 4F). In addition to the early component of the PSP (the peak), where excitation is maximal, we also quantified the late component, where inhibition levels are greatest and thus largely represents the IPSP (Bhatia et al., 2019; Antoine et al., 2019) (Figure 4B). This analysis revealed that the late PSP (LPSP) was significantly reduced with temperature elevations from 30°C to 39°C (Figure 4G). Further analysis indicated that there was no correlation between the ST and LPSP at any age (P7-8: r=0.17, P=0.41, XY pairs=25; P12-14: r=0.24, P=0.31, XY pairs=19; P20-23: r=0.33, P=0.28, XY pairs=13 Pearson’s test). However, a significant, moderate relationship occurred between the magnitude of the LPSP and the PSP threshold (Figure 4H). Thus, PNs with increased STs and larger PSPs were part of a circuit where synaptic inhibition was reduced across all or most PNs, regardless of whether STs were elevated. These results suggest that at P12-14, ST is sensitive to temperature increases and spiking is maintained by larger depolarization levels in these neurons.

At a later age (P24-P26), we applied external heat to elevate normal T_b (36°C) into the fever range and used in vivo high-density, extracellular electrophysiological recordings to capture single-unit activity in the S1 cortex of ketamine/xylazine anesthetized wildtype mice (Figure 5A). Unlike in vitro whole-cell patch-clamp recordings, extracellular recordings allow analysis only from cells that spike throughout the recording period during temperature elevations from 36°C (baseline) to fever-range. This is largely due to the inability of spike sorting algorithms to reliably track neurons that stop firing during the recording.

Color maps (with yellow and red indicating spiking activity) illustrate multiunit activity obtained in the left and right cortices during baseline temperature (left), warming to fever Tb (middle), and cooling back to 36°C (recovery upon heat removal) are shown in Figure 5B. Figure 5C presents example traces of spiking activity recorded from three cortical positions (marked by black, gray, and white asterisks on the right side of panels in Figure 5B). Following manual curation, single units (n = 633) were classified as putative inhibitory interneurons or excitatory pyramidal neurons (PNs) based on their spike waveform duration (Barthó et al., 2004; see Methods). The spike duration was calculated as the time difference between the trough and the subsequent waveform peak of the mean filtered (300 – 6000 Hz bandpassed) spike waveform. Durations of extracellularly recorded spikes showed a bimodal distribution (Hartigan’s dip test; p < 0.001) characteristic of the neocortex with shorter durations corresponding to putative interneurons (narrow spikes) and longer durations to putative principal cells (wide spikes). Next, k-means clustering was used to separate the single units into these two groups, which resulted in 140 interneurons (spike duration < 0.6 ms) and 493 principal cells (spike duration > 0.6 ms), corresponding to a typical 22% - 78% (interneuron – principal) cell ratio. Figure 5D displays representative single units isolated via spike sorting, with unit 3 representing a putative interneuron.

Finally, the firing rates of neurons were computed separately during the three distinguished periods (baseline, thermal fever (39°C) and recovery back to 36°C). To decrease the effect of transient changes (e.g., tissue recovery during the baseline period or the short-term effect of heating during the thermal fever period), for each 45-minute-long period, we used only the last 25-minutes to calculate the mean firing rates of neurons. The normalized average firing rate for each animal (n= 5) was calculated for these putative inhibitory interneurons and excitatory PNs across the recording session. We found that in vivo, putative excitatory PNs that remained active increased their spike rates at fever temperatures (Figure 5E). The normalized firing rates of five representative units are shown in Figure 5F during the baseline, warming (red shaded area), and recovery periods. Additionally, we observed that putative inhibitory interneurons that remained active also increased their spike rates at fever temperatures (Figure 5E).

Consistent with these in vivo findings, ex vivo recordings at a similar age (P20-23) showed that 11.1% of recorded excitatory PNs (n=27 cells) spiked at both 36°C and 39°C (Figure 5G). However, similar to the in vivo condition, neurons that remained spiking at 36°C and 39°C increased their spike/firing rate (Figure 5H). Along with the 11.1% of recorded excitatory PNs that stayed spiking during temperature increases from 36°C to 39°C, roughly equal proportions of neurons stopped spiking and started spiking (Figure 5G). Thus, neurons that maintain spiking at fever temperatures play a crucial role in determining the overall activity levels of the circuit. Additionally, we found that the proportion of PNs that remain spiking from 36°C to 39°C is higher at earlier ages (P12-14) compared to later ages (P20-23) (Figure 5I).

As neurons that remain spiking with temperature elevations may be crucial in setting overall cortical circuit activity levels, we sought to uncover additional information about the intrinsic characteristics of cortical PNs at P12-14 that facilitate their activity at 36°C to 39°C. We refer to neurons that retained spiking at all three temperatures (30°C, 36°C, and 39°C) as STAY neurons, while those cells that stopped spiking upon temperature transitions from 30°C to 36°C or 36°C to 39°C are denoted as STOP neurons. Firstly, in STOP and STAY neurons, ST increased during temperature elevations from 30°C to 39°C (Figure 6A). However, in STOP PNs, the increases in ST were greater than in those neurons that continued spiking (Figure 6A). The ST of neurons that continued spiking at 39°C was similar in magnitude to that of neurons that stopped spiking at 36°C.

There was an interaction of temperature and whether the PNs were STOP or STAY cells in the PSP (P=0.048, two-way ANOVA, mixed-effects model), with a trend toward larger PSPs in the STAY PNs (Figure 6B). STAY cells showed a strong positive correlation between PSP and ST (Figure 6C-D). At 30°C, STAY neurons had larger LPSPs (Figure 6E). However, inhibition levels in these PNs were unique in that inhibition levels significantly decreased with temperature elevations (Figure 6E). Moreover, we found a significant moderate inverse correlation between ST and LPSP in the STAY cells that remain spiking at 39°C (r=-0.57, P=0.046, XY pairs=13, Spearman). Thus, in PNs that stay spiking, the excitation and inhibition levels are matched to facilitate greater depolarization and maintain spiking activity. Temperature-induced alterations in spike height and spike afterhyperpolarization (AHP) were not different between STAY and STOP cells, while R_in was largely insensitive to temperature changes in STAY PNs but not in STOP PNs (Figure 6F-H).

Next, we sought to identify molecular determinants that might underlie the properties of STAY cells. We focused on the transient receptor potential vanilloid 1-4 (TRPV1-4) channels, a family of thermosensitive, Ca²⁺-permeable, nonselective cation channels. Of these, only TRPV3 and TRPV4 are activated within the innocuous warm temperature range (31-39°C), whereas TRPV1 and TRPV2 respond to noxious heat above 42°C and 52°C, respectively (Kashio and Tominaga, 2022). TRPV3, in particular, has a reported temperature threshold for activation between 31– 39°C, with Smith et al., 2002 noting activation at 39°C.

We examined TRPV3 protein expression by immunohistochemical staining at P7, P14, and P21 (Figure 7A-B). Staining specificity was confirmed at all ages by co-incubating the TRPV3 antibody with a TRPV3 blocking peptide, with an example shown for a P7 brain section (Figure 7A, right panel). TRPV3 expression was detected in S1 cortex at all ages, with strong expression in L2/3 at P14 but not at P7 or P21 (Figure 7A). Outside S1 cortex, TRPV3 was also present in the hippocampus, thalamus, and striatum (Figure 7A-B). Co-immunostaining for TRPV4 revealed several brain regions where TRPV4 (red color) and TRPV3 (green color) were non-overlapping (Figure 7B, top-left panel). In the S1 cortex at P14, TRPV3 was detected in L4 barrels and showed a patchy distribution in L2/3, with some patches displaying more intense staining (Figure 7B, bottom-right panel). This pattern was not observed for TRPV4 (Figure 7B, bottom-left panel).

Using whole-cell voltage clamp, we recorded TRPV3 currents at different voltages in cortical excitatory PNs with bath application of camphor (5 mM), a TRPV3 agonist (Moqrich et al., 2005) (Figure 8A). We observed a significant increase in TRPV3 currents at 39°C (Figure 8B-C, two-way repeated measures ANOVA, with Fisher’s LCD test). We also measured the net TRPV3 currents by comparing TRPV3 current densities recorded in the absence and presence of forsythoside B, a natural TRPV3 inhibitor found in the plant Lamiophlomis rotate with an IC 50 of 6.7 μmol/L and no obvious inhibitory effects on the other TRPV channels like TRPV1 or 4 (Zhang et al., 2019). Consistent with Smith et al., 2002, who noted TRPV3 activation at 39°C, we also observed a net positive TRPV3 current at 39°C (Figure 8D-G).

To investigate the potential role of TRPV3 and TRPV4 in regulating STAY PN phenotype in the S1 cortex during temperature increases from 30-39°C, we applied the following blockers: forsythoside B (50μM) or RN1734 (10μM, TRPV4 blocker) intracellularly in our recordings of L4-evoked spiking at 1.4X Eθ in L2/3 PNs (Figure 9A). RN1734 is a highly selective antagonist for TRPV4 with an IC 50 of 5.9 μmol/L and minimal to negligible affinity for the other TRPV channels (Vincent et al., 2009). These blockers phenocopy electrophysiological effects observed in Trpv4 ^-/- KO and Trpv3 ^-/- KO mice (Chen et al., 2022). For instance, Trpv3 ^-/- KO mice exhibit a hyperpolarizing effect on the resting membrane potential (RMP) in striatal MSNs and L2 stellate neurons in the entorhinal cortex (Chen et al., 2022). This hyperpolarization reflects reduced TRPV3 activity, which normally promotes cation influx (Shibasaki et al., 2007; Chen et al., 2022). Although increases in temperature alone can depolarize the RMP, consistent with prior reports in Trpv3^⁻/⁻ mice, our recordings of L4-evoked spiking in L2/3 PNs of S1 cortex showed a similar trend toward RMP hyperpolarization with intracellular TRPV3 block via forsythoside at 30 °C and 36 °C, and a significant enhancement at 39°C. Specifically, in the absence of block (n = 13 cells) versus with TRPV3 blocker (n = 15 cells), the RMP values were: 30 °C: −75.1 ± 1.1 vs. −77.1 ± 1.0 mV (P = 0.273); 36 °C: −73.5 ± 1.2 vs. −76.3 ± 1.1 mV (P = 0.273); 39 °C: −71.3 ± 1.0 vs. −74.4 ± 1.1 mV (P = 0.048); two-way repeated-measures ANOVA with Fisher’s LSD test.

Furthermore, we found that intracellular blockade of TRPV3 caused a significant reduction in the percentage of STAY neurons (No block: 27% (37 recorded cells), TRPV3 block: 8% (n=26 recorded cells), TRPV4 block: 17% (n=24 recorded cells), two-tailed binomial test) (Figure 9B). Consistent with the STAY cells being largely responsible for maintaining stable average spiking during temperature elevations, spiking was reduced in PNs with intracellular TRPV3 blockade during the temperature transitions from 30°C to 39°C (Figure 9C). TRPV3 block did not occlude the increases in the depolarization required to reach ST (TRPV3 blocker: 30°C: 7.15 ± 0.57, 36°C: 8.39 ± 0.26, 39°C: 10.0 ± 0.94, 30°C vs 39°C: P=0.007, One-way repeated measures ANOVA, mixed-effects model with Tukey’s test). However, it did abolish the temperature induced increases in the PSP (TRPV3 blocker: 30°C: 2.37 ± 0.36, 36°C: 2.45 ± 0.38, 39°C: 2.22 ± 0.25; No blocker: 30°C: 2.47 ± 0.32, 36°C: 2.90 ± 0.32, 39°C: 3.32 ± 0.52, 30 vs 39°C:

P=0.030, One-way repeated measures ANOVA, mixed-effects model with Dunnett’s test). Consistently, TRPV3 block abolished the positive correlation between ST and PSP (No block: XY pairs=27, TRPV3 blocker: XY pairs=18, Pearson’s test) (Figure 9D).

Furthermore, in keeping with the correlation between the magnitude of the LPSP and the PSP threshold, we found that TRPV3 blockade also prevented the temperature induced decreases in the LPSP (Figure 9E). Moreover, these reductions in spiking with TRPV3 blockade were not dependent on temperature-induced decreases in R_in, as changes to PNs were similar in the TRPV3 block and no-block conditions (Figure 9F). Altogether, these results suggest that TRPV3 expression or its functional effects—such as mediating Ca²⁺ or other cation (Na⁺, K⁺) influx to promote increased depolarization with spike threshold increases—may define the subset of STAY neurons. Given this role, the loss of TRPV3 is expected to promote temperature-induced hypoactivity in cortical circuits rather than exacerbate seizure activity.

Recordings of L4-evoked PSPs and spiking in L2/3 cortical PNs from wildtype (WT) and homozygous Trpv3 knockout (Trpv3^-/-) (Moussaieff A et al., 2008) mice revealed similar degrees of temperature sensitivity in ST and R_in. Febrile temperatures increased ST and decreased R_in (Figure 10A-C). Consistent with findings that TRPV3 blockade using 50μM forsythoside B reduces spiking in cortical L2/3 PNs, we observed significantly reduced spiking in Trpv3^-/- mice as well (Figure 10D). Analysis of postsynaptic potentials in these neurons, showed that, in WT mice, PSP amplitude increased with temperature elevation into the febrile range, whereas this temperature-dependent depolarization was absent in Trpv3^-/- mice (Figure 10E). This reduction in spiking is consistent with the idea that, during increases in STs, concurrent increases in PSP amplitude help sustain neuronal spiking. These findings support the conclusion that TRPV3 contributes to enhanced depolarization and excitability at febrile temperatures.

We next examined the effects of reduced Trpv3 expression on hyperthermia-induced seizures in WT(Trpv3^+/+), heterozygous (Trpv3^+/-), and homozygous knockout (Trpv3^-/-) P12 pups. As described in Figure 1D, mice were exposed to infrared light to raise their core T_b into the range of hyperthermia (>40°C; Figure 10F). To enable wireless, non-invasive T_b recordings, all mice were implanted with temperature transponders. Baseline T_b measurements prior to infrared exposure revealed no significant differences across genotypes (Trpv3^+/+ : 33.7 ± 0.44°C, Trpv3^+/-: 33.7 ± 0.33°C, Trpv3^-/- : 32.8 ± 0.10°C; One-way ANOVA with Tukey’s multiple comparison test, F (2,11) = 2.41, P = 0.136)). As T_b increased during infrared exposure, P12 mice typically transitioned from rest to active walking and exploratory behavior. The T_b at which mice became ambulatory did not significantly differ across genotypes (Trpv3^+/+ : 37.5 ± 0.55°C, Trpv3^{+/ -} : 36.3 ± 0.43°C, Trpv3^-/- : 36.6 ± 0.34°C; One-way ANOVA with Tukey’s multiple comparison test, F(2,11) = 1.34, P = 0.301).

As T_b reached 41–42°C, mice exhibited loss of postural control (LPC), defined as collapse and failure to maintain upright posture. LPC is a key behavioral marker of cortical and motor circuit dysfunction. The T_b at which LPC occurred was significantly lower in Trpv3^-/- mice compared to wildtype controls (Figure 10G). LPC was followed by the rapid onset of generalized seizures. While the T_b at seizure onset was not significantly different among genotypes (Figure 10H), the time to seizure onset was significantly prolonged in both Trpv3^+/-and Trpv3^-/- mice. Specifically, Trpv3^+/- mice exhibited a 49.4% increase, and Trpv3^-/- mice a 91% increase in latency to seizure compared to WT mice (Figure 10I). Importantly, the rate of T_b increase leading to seizure did not differ across genotypes (Trpv3^+/+: 0.011 ± 0.003°C/s, Trpv3^{+/ -}: 0.013 ± 0.002°C/s, Trpv3^-/- :0.012 ± 0.001°C/s; One-way ANOVA with Tukey’s multiple comparison test, F(2,11) = 0.1055, P = 0.9008), suggesting the delayed seizure onset was not due to slower heating. Together, these findings indicate that reduced TRPV3 function increases resistance to seizure initiation and/or propagation under febrile conditions, likely by decreasing neuronal excitability.

All experimental procedures were performed in accordance with the National Institutes of Health guidelines for care and use of laboratory animals and the EC Council Directive of September 22, 2010 (2010/63/EU). All experiments were approved by the Animal Care and Use Committee of the National Institute on Alcohol Abuse and Alcoholism (NIAAA) (Protocol # LIN-MA-1) and the Animal Care Committee of the HUN-REN Research Centre for Natural Sciences (RCNS) and by the National Food Chain Safety Office of Hungary (license number: PE/EA/1004- 5/2021). Mice were housed under standard conditions with ad libitum access to food and water, under a 12-h light/dark cycle. The following mouse strains were obtained from The Jackson Laboratory (JAX):C57BL/6J (JAX 000664), FVB/N mice (JAX 004828), and Trpv3 (JAX 010773) (Moussaieff A et al., 2008; Moqrich et al., 2005). The day of birth was denoted as postnatal day (P)1. Both male and female mice were used for body temperature (T_b) recordings and electrophysiological experiments.

All mice used for body temperature (T_b) recordings were implanted with sterile IPTT-300 Implantable Programmable Temperature transponders (Bio Medic Data Systems, LLC) at 2 weeks of age. The injection site was pre-cleaned with a betadine solution, and the sterile transponders were injected subcutaneously using a pre-loaded 12-gauge syringe into animals anesthetized with 4-5% isoflurane. The injection site was on the left side, approximately 11 mm from the base of the hip, with the animal in the prone position. Post-implantation, the mice were able to move, eat, and drink autonomously, with no adverse phenotypes noted. T_b recording commenced at least 5 days after injection. Prior to recording, mice were brought into the procedure room and allowed to habituate for 1 hour. T_b was recorded every 5 minutes for 6 hours using an IPTT-300 Implantable Programmable Temperature reader (DAS-8027) in a 37L x 16W x 13H cm cage. T_b was recorded in the presence (Figure 1D) or absence (Figure 1A-C) of infrared light. Infrared light exposure was provided via a 250-W temperature-controlled infrared lamp (catalog #50320, Stoelting), set to 32.5-33.5 °C. The lamp was positioned ∼15 cm above the cage, and the ambient temperature was monitored every 5 minutes (Figure 1D).

Using standard methods, acute primary somatosensory cortex (S1) slices (350 µm thick) from P7-8, P12-14, or P20-23 mice were cut in the "across-row" plane, oriented 35° toward coronal from midsagittal (Antoine et al., 2019). The cutting solution contained (in mM): 85 NaCl, 75 sucrose, 25 D-(+)-glucose, 4 MgSO₄, 2.5 KCI, 1.25 NaH₂PO₄, 0.5 ascorbic acid, 25 NaHCO₃, 0.5 CaCl₂. Once cut, slices were transferred to a submerged-style holding chamber containing standard Ringer’s solution (in mM: 119 NaCl, 2.5 KCI, 1.3 MgSO₄, 1 NaH₂PO₄, 26.2 NaHCO₃, 11 D-(+)-glucose and 2.5 CaCl₂) and then incubated at 32°C for 30 minutes. Both solutions were at neural pH (i.e., 7.3), 300 mOsm, and saturated with 95% O₂ and 5% CO₂. Slices were kept at room temperature for at least 30 minutes before being transferred to a submerged recording chamber.

Whole-cell current-clamp recordings were made from pyramidal neurons (PNs) visually identified via infrared DIC optics. Physiological verification for regular spiking was done in current clamp. Recordings were made using 3-6 MΩ micropipettes containing (in mM):116 K gluconate, 20 HEPES, 6 KCI, 2 NaCl, 0.5 EGTA, 4MgATP, 0.3 NaGTP, 10 Na phosphocreatine, with a Multiclamp 700B amplifier (Molecular Devices, Sunnyvale, CA). Where applicable, forsythoside B (50 µM, TRPV3 blocker, Millipore Sigma Cat #: PHL83313) or RN1734 (10 µM, TRPV4 blocker, Tocris Bioscience Cat #: 3746/25) was added to the internal solution. Camphor (MedChem, Cat. No: HY-N0808) was added to the bath solution for experiments in Figure 8.

Signals were filtered (2-6 kHz) and digitized (10-20 kHz). Perfusate temperature for all in vitro experiments was regulated by a PC-connected Peltier heater (SM-4600, Scientifica, UK) and temperature controller (Linlab2, Scientifica, UK), with a high-accuracy, low-noise temperature control system, containing both a built-in temperature sensor and bath sensor for accurate feedback. ACSF flow rates of ∼3 ml/min were used to facilitate efficient heating/cooling.

To determine the effect of temperature on an activated cortical network, layer 4-evoked spiking was quantified during perfusate temperature increases from 30°C to 36°C, and then to fever range (∼39°C) while recording membrane potential. The slicing plane facilitated clear identification of whisker barrel columns, and neuronal activity was evoked by stimulating the barrel center using a bipolar electrode (0.2 ms pulses) at a stimulation of 1.4 Eθ. Eθ is defined as the minimal intensity that evoked a consistent excitatory postsynaptic current (EPSC) during more than 3 of 5 consecutive sweeps with 10 s ISI (Antoine et al., 2019). Eθ was determined in voltage clamp for each recorded cell prior to recording PSPs and spiking in current clamp. A new brain slice was used for each recording. L4-evoked feedforward post synaptic potentials (PSPs) and spiking were recorded in single L2/3 PNs from a pre-stimulus baseline membrane potential (V_m) of -50 mV. This V_m, just below spike threshold, was selected to mimic in vivo conditions, as PNs in vivo can reside within this V_m range during whisker exploration of objects or surfaces (Yamashita et al., 2013). Evoked spikes were analyzed over a 150-ms interval, starting 2-3 ms post-stimulus, for 11 sweeps at 10 s inter-sweep interval (ISI). Spike threshold was defined as the membrane potential (V_m) at which the second derivative of V_m was >5 SDs above the pre-stimulus period. Intrinsic spiking excitability was measured in the presence of glutamate and GABA-A receptor blockers (in μM: 100 APV, 10 NBQX, 3 gabazine). F-I curves were obtained from PNs held at –80 mV by applying incremental current injections relative to this V_m.

In vivo experiments were done under general anesthesia where five mice (age, P24-26; body weight, 7-10 g; both genders) received an intraperitoneal injection of ketamine (100 mg/kg) and xylazine (10 mg/kg). Regular doses of the ketamine/xylazine cocktail were given intramuscularly to maintain the depth of anesthesia during the experimental sessions. Up until the thermal fever protocol started, the T_b of the animals was kept at 36°C using a homeothermic heating pad connected to a temperature controller (Supertech, Pécs, Hungary). To measure the internal body temperature, a Type T thermocouple microprobe (MT-29/5; Physitemp Instruments, Clifton, NJ, USA) was placed in the rectum of the animal. The microprobe has a shaft diameter of 330 µm, a time constant of 0.025 second and 0.1°C accuracy.

To perform high-density extracellular electrophysiological recordings, mice were placed in a stereotaxic frame (David Kopf Instruments, Tujunga, CA, USA), then two circular craniotomies (∼1.5 mm in diameter) were made with a dental drill above the left and right barrel cortices. The craniotomies were centered at the following stereotaxic coordinates: anterior-posterior (AP): −1.0 mm; medial-lateral (ML): 3.5 mm (with respect to the bregma; Paxinos and Franklin 2001). Two commercially available Neuropixels 1.0 silicon probes (imec, Leuven, Belgium; Jun et al., 2017) mounted on two motorized stereotaxic micromanipulators (Neurostar, Tubingen, Germany) were implanted into the brain, one into the left and the other into the right barrel cortex, to a depth of 1.5 mm. An insertion speed of 2 µm/s was used to decrease the mechanical trauma caused by the probe insertion (Fiáth et al. 2019).

In order for the probe tracks to be perpendicular to cortical layers, the probes were inserted at an angle of 20 degrees from vertical. The dura mater in the craniotomy was left intact, except when it was too thick and thus the silicon probe could not pierce through this layer (which was indicated by significant probe buckling; n = 3 insertions). In these cases a 36 gauge, slightly bent needle was used to carefully cut the dura above the targeted cortical area. After the probe reached its final insertion depth, to allow the brain tissue to settle, we waited at least 10 minutes before electrophysiological recording was started. Spiking activity of cortical neurons (action potential band, 300–10.000 Hz) was recorded on 384 channels (768 channels in total for the two probes), with a sampling rate of 30 kHz/channel and with a gain of 500 (yielding a resolution of 2.34 μV per bit). Data were acquired using the SpikeGLX open-source software (http://billkarsh.github.io/SpikeGLX/). A common stainless steel wire inserted into the neck muscle of the animal served as the external reference and ground electrode. To avoid the dehydration of the cortex, the skull and the craniotomy was kept moist during the whole experiment using body temperature, sterile physiological saline solution and Gelaspon. Before starting the experimental protocol, manual whisker stimulation (by repetitively touching the whiskers of the animal with a cotton swab) was used to verify the recording position. In all cases (n = 10 probe insertions), we detected strong whisker-evoked neuronal activity.

First, we recorded cortical activity for 45 minutes at ∼36°C body temperature ("baseline " period), then the body temperature of the animal was elevated with the aid of the heating pad and a power bank with warming capability which was placed next to the mouse. The body temperature was increased from 36°C to 39°C in about 5 minutes (∼0.01°C/s). The elevated body temperature was maintained for 45 minutes ("thermal fever " period). Next, we turned off the heating until the T_b reached physiological temperatures (36°C; ∼5 min; ∼0.01°C/s). Cortical activity was recorded for another 45 minutes at 36°C body temperature ("recovery " period).

Continuous recordings with a total duration of 145 minutes were obtained for each mouse.

To extract cortical single-unit activity, spike sorting was performed with Kilosort2 (https://github.com/MouseLand/Kilosort; Pachitariu et al., 2016a, 2016b) using the default parameter set (available in the StandardConfig.m file). Channels containing activity acquired by recording sites located outside the cortex were removed before spike sorting (usually ∼130 channels/probe recorded from the barrel cortex). The list of single unit clusters generated by Kilosort2 was visually inspected to remove units considered as noise (e.g., units with abnormal spike waveform shapes) or multi-unit activity (e.g., clusters with a contaminated refractory period). Manual curation of the Kilosort2 results was done with the Phy Python library, which provides a graphical user interface for interactive visualization of high-density data and supply operations for merging, splitting and marking of clusters (https://github.com/cortex-lab/phy; Rossant et al., 2016). In this dataset, we aimed to keep only those single units which had at least 900 spikes (>∼0.1 Hz firing rate), a clear refractory period, a consistent waveform shape and whose spikes were present throughout the 145-minute-long recording.

Following manual curation, based on their spike waveform duration, the selected single units (n = 633) were separated into putative inhibitory interneurons and excitatory principal cells (Barthó et al., 2004). The spike duration was calculated as the time difference between the trough and the subsequent waveform peak of the mean filtered (300 – 6000 Hz bandpassed) spike waveform. Durations of extracellularly recorded spikes showed a bimodal distribution (Hartigan’s dip test; p < 0.001) characteristic of the neocortex with shorter durations corresponding to putative interneurons (narrow spikes) and longer durations to putative principal cells (wide spikes). Next, k-means clustering was used to separate the single units into these two groups, which resulted in 140 interneurons (spike duration < 0.6 ms) and 493 principal cells (spike duration > 0.6 ms), corresponding to a typical 22% - 78% (interneuron – principal) cell ratio. Finally, the firing rates of neurons were computed separately during the three distinguished periods (baseline, thermal fever and recovery). To decrease the effect of transient changes (e.g., tissue recovery during the baseline period or the short-term effect of heating during the thermal fever period), for each 45-minute-long period, we used only the last 25-minutes to calculate the mean firing rates of neurons.

Mice were anesthetized and transcardially perfused using standard methods. Brains were post-fixed overnight in 4% PFA at 4 °C with shaking, rinsed in 1X PBS, and cryoprotected in 30% sucrose. They were embedded in optimal cutting temperature (OCT) compound, frozen on dry ice, and stored at –80 °C until sectioning. Tissue was cut at 30 µm on a cryostat, mounted onto slides. Sections were permeabilized in 0.3% Triton X-100 (in 1X PBS) for 10 min at room temperature. Endogenous peroxidase activity was quenched by incubation in 3% H₂O₂ (in 1X PBS) at room temperature. Slides were blocked in 10% normal goat serum for 1 h at room temperature, then incubated overnight at 4 °C with anti-TRPV3-Biotin antibody (#ACC-033-B), TRPV3 blocking peptide (#BLP-CC033), and/or anti-TRPV4 antibody(#ACC-034) (Alomone Labs, Israel). For antibody–peptide controls, the anti-TRPV3-Biotin antibody was pre-incubated with its blocking peptide at a 1:2.5 ratio, using 1 mg/ml of each. The following day, sections were incubated for 1 h at room temperature with fluorescently labeled secondary antibodies (streptavidin-conjugated or unconjugated). Nuclei were counterstained with DAPI Fluoromount-G (Electron Microscopy Sciences, USA). Finally, sections were mounted onto slides, cover-slipped, and stored at 4 °C until imaging. Images were acquired using a Keyence BZ-X810 microscope (BZ-X800 Viewer, Version 1.1.1.8) with BZ-X filter sets: DAPI (Ex 360/40 nm; OP-87762), GFP (Ex 470/40 nm; OP-87763), and TRITC (Ex 545/25 nm; OP-87764).

Statistical analyses were conducted using Prism 9.0 (GraphPad). For experiments, at least 3 mice from a minimum of 2 separate litters were used, with littermates typically included across the age groups P7-8, P12-14, and P20-23. Non-Gaussian data were either log-transformed for parametric testing or subjected to nonparametric tests, as specified in the results and figure legends. Significance was consistently reported at α = 0.05. Details on the statistical tests used are provided in the figure captions.