Introduction

Maintaining a relatively constant core body temperature (Tb) is a vital homeostatic need. While excessive deviations (e.g. >2°C) from the typical mammalian Tb of approximately 37°C can be harmful (Nixdorf-Miller et al., 2006), Tb normally fluctuates within a narrow range (typically 0.5 -1°C) across sleep-wake cycles and behavioral states (Blessing and Ootsuka, 2016; Refinetti, 2020, 1999). These fluctuations often exhibit bimodal or multimodal distributions, indicating distinct thermal states. Additionally, Tb variations are increasingly understood not as passive outcomes of animal physiology, but as regulated, brain-initiated transitions between defended Tb “balance-points” (Romanovsky, 2018) (Blessing and Ootsuka, 2016; Morrison and Nakamura, 2019; Škop et al., 2024). For example, the transition from the rest to the active balance-point involves modulation of both autonomic and behavioral thermoeffector pathways (i.e., regulated heat loss and production) to meet the energetic demands of movement and arousal (Blessing and Ootsuka, 2016; Harding et al., 2020; Škop et al., 2024). However, the neural mechanisms that govern these physiological transitions remain incompletely understood (Harding et al., 2020, 2019).

Thermoregulation is orchestrated by the integration of autonomic and behavioral effectors. Autonomic outputs, such as brown adipose tissue (BAT) thermogenesis and vasomotor tone, are well-characterized at the circuit level. This pathway integrates sensory signals from peripheral thermosensors in the preoptic area of the hypothalamus (POA), then relays this information to the dorsomedial hypothalamus (DMH), and then to sympathetic efferent neurons of the rostral medullary raphe (rMR) (Morrison, 2011; Morrison and Nakamura, 2019). By contrast, the neuronal substrates of behavioral thermoregulation—such as huddling, nesting, or physical activity (Gordon, 2012; Landen et al., 2024; Škop et al., 2024, 2021)—and how they integrate with autonomic outputs are less well defined. Although recent studies have made progress in identifying brain regions underlying some of the associated neural pathways, (Almeida et al., 2015; Jung et al., 1998; Sotelo et al., 2022; Tan and Knight, 2018), much remains unknown about how animals behaviorally thermoregulate.

Behavioral thermoregulation strategies are naturally organized into alternating sequences across bouts of arousal and quiescence and are modulated by environmental conditions (including social context) and internal state (Endo et al., 2018; Gordon et al., 2014; Harding et al., 2020, 2019; Landen et al., 2024; Sotelo et al., 2024, 2022). Huddling is a cooperative behavior that provides thermal and energetic benefits through shared body heat and reduction of heat loss (Alberts, 1979; Fransson et al., 2005; Gilbert et al., 2010; Haig, 2008). When mice experience mild cold stress, huddling facilitates both entry into, and exit from, an energy-saving, low-Tb quiescent state (Landen et al., 2024). This pattern suggests active neural regulation of behavior and sympathetic arousal in tandem; however, the brain circuits that support this coordination remain unidentified.

Oxytocin (OT)-producing neurons in the paraventricular nucleus of the hypothalamus (i.e., PVN^OT) are a compelling candidate for such coordination (Rogers et al., 2024). PVN^OT neurons have been implicated in autonomic arousal (Chen et al., 2021), energy expenditure (Noble et al., 2014), and thermoregulation— including activation of BAT thermogenesis (Fukushima et al., 2022; Kasahara et al., 2013; Sutton et al., 2014). The oxytocin system is also classically associated with social behaviors (Froemke and Young, 2021) including physical touch (Tang et al., 2020)—a key element of huddling. Indeed, OT-mutant rat pups display deficits in both warm-seeking and huddling behavior (Harshaw et al., 2018), and OT neurons display bursts of pulsatile activity during bouts of lactation and direct physical contact with pups (Yaguchi et al., 2024). Yet, whether PVN^OT neural activity is functionally associated with the organization of thermoregulatory states is not known.

Here, we set out to identify the relationship between PVN^OT activity and daily patterns of behavioral thermoregulation in mice. We initially observed that PVN^OT activity is linked to huddling states. Then, we discovered that PVN^OT calcium dynamics during huddling are predictive of transitions to body warming and arousal. Intriguingly, in a parallel fashion, PVN^OT activity also tracks transitions from rest to arousal in solo animals. Last, using automated thermal feature tracking and optogenetics, we demonstrate that PVN^OT stimulation initiates state-dependent transitions that coordinate physiological warming with behavioral arousal. These findings reveal a previously unrecognized role for PVN^OT neurons in regulating dynamic, context-dependent thermo-behavioral strategies.

Results

PVN^OT neural activity is associated with huddling substates

Mice use huddling to facilitate thermoregulatory state transitions. We sought to identify huddling-associated brain regions in groups of adult females, which show stronger huddling-associated Tb changes than male groups (Landen et al., 2024). Using a home-cage surveillance system that eliminates direct human-animal interaction (Landen et al., 2024) in conjunction with FOS immunohistochemistry (IHC), we examined patterns of neural FOS activity in 18 brain regions. We focused on regions associated with thermoregulation, social behavior, and identified to be FOS-activated during huddling substates in preliminary studies (Fig. 1A-B). We examined animals engaged in active huddling (close physical contact while displaying some physical activity), quiescent huddling (close physical contact with no physical activity) and—as a control for physical touch in the absence of social contact—solo-housed animals engaged in self-grooming. For IHC, we scored cellular FOS activity after the animals had been engaged in state-specific behavior (i.e., active huddle, quiescent huddle, solo-groom) for at least 15 minutes.

Figure 1.

We initially focused on active huddling-associated FOS expression because this state requires that an animal actively maintain close physical contact with another individual. We identified three regions that were activated 90 minutes after the onset of active-huddling: the dorsal medial hypothalamus (DMH), the rostral lateral septum (LS), and the paraventricular nucleus of the hypothalamus (PVN) (Fig. 1C-F and S1A-C; all statistical results are reported in Table S1). The DMH and the PVN are areas associated with energy regulation and thermoregulation, with inputs to the rMR (Morrison et al., 2014; Morrison and Nakamura, 2019; Seoyoung Son et al., 2022). The LS is associated with a variety of functions, from motivation to fear responses (Besnard and Leroy, 2022; Rizzi-Wise and Wang, 2021) and is strongly associated with the control of social behavior (Menon et al., 2022). These results suggest active huddling recruits brain regions with known roles in thermoeffector and social interaction pathways.

The PVN is one of two primary locations of OT-producing neurons (Seoyoung Son et al., 2022). Given the established associations of PVN^OT neurons with thermoregulation and social behavior, we asked whether they were activated during active and quiescent huddling. To increase the temporal resolution of our assay from the timescale of protein translation to that of mRNA transcription, we used single-molecule fluorescence in situ hybridization (smFISH) to identify PVN^OT neurons that became Fos-positive 30 minutes after the onset of a behavioral state. We quantified neurons in the dorsal and ventral PVN, which have unique cyto-architectures and projection patterns (Biag et al., 2012) (Fig. 1G-I). PVN^OT neurons were Fos-positive during both active huddling and quiescent huddling. During active huddling, the dorsal population was more Fos-positive than the ventral population, and this dorsal population displayed greater Fos activity during active compared to quiescent huddling (Fig. 1I).

These results suggest the PVN is a huddling-associated brain region and that PVN^OT neurons are activated during huddling substates. Given the coarse timescale of Fos dynamics (≥30 minutes) compared to the faster timescale of rest-arousal cycles in mice (minutes; (Li et al., 2022; McShane et al., 2010)), we next employed in vivo calcium recordings to capture finescale population-level dynamics of PVN^OT neurons.

PVN^OT Ca²⁺ peaks track rest and arousal behavior states in social and non-social conditions

Behaviorally, mice modulate body temperature through physical activity, nesting, and huddling (Gordon, 2012; Landen et al., 2024; Škop et al., 2024, 2021). In a longitudinal study accounting for individual variance in physiology and behavior, we next asked how calcium activity in PVN^OT neurons was associated with real-time behavior during social thermoregulation in different thermal conditions (Fig. 2A-C). As a control for social context, we also examined solo-housed animals. Experiments were conducted at three floor temperatures: cold (15°C), cool (23°C; room temperature), and a temperature near the “lower critical temperature” (29°C; the transition point between energy-requiring and energy-neutral thermoregulatory mechanisms (Gordon, 2012; Maloney et al., 2014; Škop et al., 2020)). We tracked six behaviors in the solo condition and 10 behaviors in the paired condition (Fig. 2B).

Figure 2.

We used an AAV encoding the calcium indicator GCaMP (pGP-AAV9-syn-FLEX-jGCaMP8s-WPRE) (Zhang et al., 2023) and fiber photometry using an optical fiber implanted above the PVN (Fig. 2A-D). GCaMP AAVs have been validated to transfect OT neurons in OT-Cre mice (Resendez et al., 2020; Yukinaga et al., 2022; Zhan et al., 2024) and, in this study, showed co-labelling with OT-immunoreactive neurons (Fig. S2A-B). Neural-photometric recordings were conducted during the light/rest phase, when the distinction between different thermoregulatory states is more discernable (Landen et al., 2024). We observed large calcium transients (hereafter “peaks”) that appeared to be similar to those previously reported in lactating female mice (Yaguchi et al., 2024; Yukinaga et al., 2022) (Fig. 2D). Further analyses focused on peaks that were at least six standard deviations above baseline.

We first examined the effect of social context (i.e., paired-vs. solo-housed) on calcium peak amplitude and frequency across the three floor temperatures. While there was no effect of social context on peak amplitude (Fig. S2C), there was an increase in peak frequency in paired animals (Fig. S2C-D). To our surprise, in both social and nonsocial contexts, PVN^OT peaks were strongly linked to specific behaviors. In the solo context, peak amplitude and frequency were associated with nesting and quiescence (i.e., motionless rest) (Fig. 2E-G). In the paired context, peaks were associated with active and quiescent huddling (Video S1) and, to a lesser extent, nesting (Fig. 2H-J). We next examined the behavioral associations of PVN^OT peaks within each floor temperature. The behavior-state associations with calcium peaks were largely preserved across all three floor temperatures (Fig. S2G-J). Linear mixed model analysis of calcium peak amplitude and frequency revealed main effects of quiescence and nesting in the solo context (Fig. S3G,I), and main effects of active huddling, quiescent huddling, and nesting in the paired context (Fig. S3H,J). Thus, PVN^OT peaks are associated with distinct resting and waking behavioral states in social and non-social contexts.

Next, we analyzed the relationship between floor temperature and peak frequency across the two social contexts (Fig. 2K-P). Peaks were more frequent 29°C for both solo (Fig. 2K) and paired (Fig. 2N) animals. In mice, sleep is more common at warmer ambient temperatures (Harding et al., 2020). To disambiguate whether PVN^OT peaks display warm sensitivity vs. arousal-state dependency, we performed a multiple regression analysis relating the number of peaks to floor temperature and the duration of rest bouts (i.e., quiescence quiescent huddling) (Fig. 2F-G and I-J). In both the solo and paired conditions, calcium peak counts were strongly related to rest bout duration (Fig. 2L,O). Moreover, quiescence and quiescent huddling bouts appeared longer at 29°C (Fig. 2L,O, red dots). We therefore examined the contribution of bout length and floor temperature to peak counts. In the solo and paired contexts, increased peak counts at 29°C could be explained by longer rest bout lengths (Fig. 2M and 2P).

Together, these results suggest PVN^OT peaks track distinct resting and active behaviors in social and non-social settings, and are slightly enhanced in the social context. In addition, PVN^OT peaks are largely unaffected by the ambient temperature conditions tested; we therefore pooled data from the three floor temperatures in subsequent analyses.

PVN^OT neurons predict transitions towards thermogenesis and behavioral arousal in social and non-social contexts

PVN^OT peaks track two rest states (quiescence and quiescent huddling), and two active states (i.e., nesting and active huddling) (Fig. 2). Notably, each of these states are in part thermoregulatory and are naturally organized into alternating sequences of bouts across rest/active cycles (Endo et al., 2018; Gordon et al., 2014; Harding et al., 2020, 2019; Landen et al., 2024; Sotelo et al., 2024, 2022). To address how PVN^OT dynamics align with these sequences, we quantified peaks relative to the onset and offset of behavior bouts, as well as to changes in physical activity and core body temperature (Fig. 3A).

Figure 3.

We first examined physical activity (as measured by frame-to-frame pixel changes) before and after PVN^OT calcium peaks in the home-cage. In both solo and paired conditions, physical activity increased following PVN^OT peaks. This increase began at the approximate time of the peak and persisted through the subsequent 400 seconds (6.6 min) (Fig. 3B-E). We next examined the density (i.e., frequency of peaks per time bin) of calcium peaks with respect to the onset and offset of the two resting states: quiescence and quiescent huddling. For both behaviors, there was a higher density of peaks in the minutes prior to bout offset compared to onset (Fig. S3A-C). Moreover, there was a positive relationship between bout length and the number of peaks per bout, indicating that longer rest bouts are more likely to have peaks (Fig. S3B,D). To better understand the relationship between PVN^OT dynamics and bouts of quiescence and quiescence huddling, we examined (1) the per-second probability of observing a peak relative to bout onset/offset using logistic regression, and (2) changes in calcium baseline (Fig. 3F-M).

For bouts of quiescence, peak probability was approximately 2.7% near the time of onset and 14.0% near the time of offset (Fig. 3F). For bouts of quiescent huddling, peak probability was approximately 2.1% near the onset and around 25.6% near the time of offset (Fig. 3H). In accordance with the peak density and probability data, baseline calcium in PVN^OT neurons was higher in the minutes leading to quiescence offset compared to onset in all animals (Fig. 3G,I and S3I-J). Thus, PVN^OT neural activity strongly predicts quiescence/quiescent huddling offset compared to onset by at least five-fold, and signals an increase in physical activity—a correlate of behavioral arousal (De Lecea et al., 2012) and a means of increasing metabolic rate and Tb (Gordon, 1993).

We next examined PVN^OT peak density and probability during the onset and offset of the two active states: nesting and active huddling (Fig. 3J-M, and S3E-H). In contrast to the two resting states, peak density in both states was higher near bout onset compared to offset (Fig. S3E,G). Moreover, there was only a weak relationship between peak count and bout length for both states, indicating that longer bouts are not more likely to have peaks (Fig. S3F,H). Logistic regression analysis of nesting bouts showed that PVN^OT peak probability was approximately 7.7% near the time of onset and 3.3% near the time of offset (Fig. 3J). For active huddling bouts, peak probability was approximately 15.3% near the time of onset and 5.2% near the time of offset (Fig. 3L). Thus, for nesting and active huddling, PVN^OT peaks are two-to three-fold more likely to predict bout onset than offset. Notably, within the paired condition, there was a higher peak probability for active huddling compared to nesting (Poisson regression, nesting vs. active huddling peaks counts, P < 0.001), suggesting that social interaction enhances the probability of PVN^OT activity.

Nesting and active huddling often occur immediately before or after (or in between) bouts of quiescence (i.e., the peri-quiescent phase) (Gordon et al., 2014; Landen et al., 2024; Sotelo et al., 2024, 2022). We next examined the quiescence phase for PVN^OT peaks associated with nesting and active huddling. PVN^OT peaks were more likely to occur during the post-quiescence phase for both nesting (Fig. 3K) and active huddling (Fig. 3M). Notably, although nesting and active huddling are commonly expressed during the pre-quiescent phase, these two behaviors are rare during the post-quiescent phase (Fig. S3K-L). Thus, even though nesting and active huddling are relatively rare during the post-quiescent phase, those are precisely the active states during which PVN^OT activity increases. These results suggest PVN^OT neurons distinguish nesting/active huddling in the post-quiescent phase from these same behaviors in the pre-quiescent phase.

Together these results suggest that PVN^OT activity dynamics predict within approximately 100 seconds the offset of two rest states (quiescence and quiescent huddling) and within around 20 seconds the onset of two post-quiescence active states (nesting and active huddling) in solo and paired mice, respectively.

PVN^OT peaks occur during low Tb (∼36.2°C) and prior to body warming

In rodents, transitions from rest to arousal are thought to result, in part, from an ultradian (i.e., occurring on cycles < 24 hr), brain-driven switch of the defended body temperature from a rest balance-point (i.e., lower bound) to an awake balance-point (i.e., upper bound) (Blessing and Ootsuka, 2016; Škop et al., 2024). Because PVN^OT calcium peaks aligned with transitions from resting to active states, we next examined individual core Tb dynamics before, during, and after PVN^OT peaks (Fig. 3N-Q). In both solo and paired conditions, PVN^OT peaks occurred when Tb was around 0.5°C lower than baseline (baseline-Tb: 36.7°C; peri-peak-Tb: 36.2°C) (Fig. 3N,P). Moreover, PVN^OT peaks aligned with the low point of a U-shaped body temperature profile: on average, Tb decreased before, and increased after, the time of the calcium peak in both solo and paired conditions (Fig. 3O,Q). Together, these results suggest that PVN^OT peaks occur during a low Tb and predict a forthcoming increase in Tb.

Validation of PVN^OT neural calcium peaks

PVN^OT neurons display pulsatile activity bursts during lactation and direct contact with pups, with heightened intensity and broadening waveforms during late-stage lactation (PPD 8 to 14) compared to early-stage lactation (PPD 2-7) (Yaguchi et al., 2024; Yukinaga et al., 2022; Yukinaga and Miyamichi, 2025). However, the extent to which this type of activity in PVN^OT neurons is exclusive to lactation and nursing events remains a matter of debate (Perkinson et al., 2021; Resendez et al., 2020; Yaguchi et al., 2024). To address whether peaks occurring during transitions from quiescence to active states resemble those during lactation, we recorded from virgin females as described above and then allowed them to reproduce and nurse their young (Fig. 4A-C). As documented in other studies (Yaguchi et al., 2024; Yukinaga et al., 2022), PVN^OT neurons displayed peaks that increased in frequency and amplitude from early-stage to late-stage lactation—an observation that provides physiological verification of the oxytocinergic identity of the recorded neurons in our system (Fig. 4D-E and Video S2). PVN^OT peaks in virgin females appeared similar in amplitude in early-stage compared to late-stage lactation and appeared narrower overall (Fig. F-H). Systematic comparisons of peak kinetics revealed that virgin female calcium peaks displayed narrower full-width half maximum (FWHM) values than either early-stage or late-stage lactation (Fig. 4H-I). Next, while virgin and early-stage lactation peak amplitudes were equivalent, both were less than that of late-stage lactation (Fig. 4J). Finally, while virgin interpeak intervals (IPIs) were equivalent to early-stage lactation, they were longer than those of late-stage lactation (Fig. 4K). Together, these results suggest that PVN^OT peaks observed in the virgin females of this study during naturalistic home-cage behaviors share similar calcium kinetics to those during early-stage, but not late-stage lactation. As such, this data uncovers a behavioral and physiological context in which PVN^OT neurons display calcium burst-like peaks.

Figure 4.

PVN^OT peaks signal increases in BAT thermogenesis and vasomotor tone

Our data suggest PVN^OT peaks signal transitions towards thermogenesis and arousal as measured by core Tb and physical activity. However, core Tb represents the integration of several different thermoeffector pathways like BAT thermogenesis and vasomotor responses. We therefore sought to automate a computational strategy that uses video thermography recordings of awake, behaving animals to approximate BAT thermogenesis and vasomotor processes at sub-second resolution—a timescale relevant to sympathetic activity.

We developed a computer vision system called Skeleton-Guided Bodypart Segmentation (SGBS) that uses deep learning to record the surface temperature of three thermal features of an animal: the surfaces over (1) interscapular brown fat (BAT), (2) the rump, and (3) the entire dorsal surface (Fig. 5A). Our pipeline integrates DeepLabCut-based animal tracking (Mathis et al., 2018) and Mask R-CNN technologies (He et al., 2017) in a Python framework (TorchVision maintainers and contributors, 2016) to process high-resolution thermography (FLIR) videos. An attention layer integrates 10 anatomical skeleton keypoints (e.g., head, nose, shoulders and tail) with a model trained on manual annotations of the three thermal features. Leveraging these keypoints, SGBS can precisely localize and delineate thermal features for mask prediction (Fig. 5B). SGBS detects above-average temperatures in the BAT region and below-average temperatures rump region without the need for shaving the back—a manipulation alters heat loss and can alter behavioral and autonomic thermoregulation (Meyer et al., 2017)(Fig. 5C). Moreover, SGBS demonstrates 100-fold lower loss compared to an unmodified Mask R-CNN using identical training data (Fig. 5D). We used SGBS to segment BAT, rump, and dorsal surface features on a frame-by-frame basis, extracting the mean temperature for each feature (Fig. 5A and Video S3). SGBS data were synchronized with timestamps and integrated with other data streams, including calcium traces and behavior tracking. For these analyses, solo-housed females were used to eliminate the confounding effects of heat exchange with social partners.

Figure 5.

To test the performance of SGBS, we examined the temporal dynamics of pairwise combinations of the three thermal features with respect to core Tb and identified the lag time that gives the best correlation (i.e., cross-correlation analysis) (Fig. 5E). Core and dorsal surface temperatures were very similar with no lag. Next, changes in core Tb lagged changes in BAT by approximately 90 seconds, a result consistent with a previous report (Ootsuka et al., 2009). Finally, rump temperature lagged core Tb by 20 seconds. Thus, SGBS captures the dynamics of thermoregulatory features at the resolution of seconds. These results support the notions that (1) core and dorsal surface temperatures are highly correlated because they both integrate multiple thermoeffector pathways, (2) BAT thermogenesis can drive temperature increases in the mouse body on the order of minutes, and (3) vasomotor processes in the rump regulate Tb on a relatively faster timescale (i.e., seconds) (Blessing et al., 2016).

We next used SGBS to align thermal feature data with PVN^OT calcium peaks from fiber photometry recordings (Fig. 5F-M). Consistent with our previous core Tb measurements (Fig. 3N-Q), BAT and dorsal surface temperatures trended downward before the calcium peak but increased afterward, a shift confirmed by a significant difference in slope before vs. after the peak (Fig. 5H,J). Moreover, BAT and dorsal surface temperatures were significantly greater after the peak compared to before (Fig. 5I,K). By contrast, although rump temperature showed no difference in slope before vs. after peak, the mean temperature was lower after the peak compared to before (Fig. 5L-M). Together, these results suggest that PVN^OT activity is associated with a thermogenic state characterized by the onset of BAT thermogenesis and, to a lesser degree, increased rump vasoconstriction, following PVN^OT activation.

Optogenetic activation of PVN^OT neurons initiates thermogenic responses

To elucidate the functional link between PVN^OT activity during quiescence and thermoeffector outputs, channelrhodopsin (AAVDJ-EF1a-DIO-hChR2(H134R)-EYFP-WPRE-pA) was expressed in the PVN of Oxytocin-Cre mice (ChR2+) to allow for optogenetic stimulation of PVN^OT neurons (Fig. 6A). To control for potential warm-sensitive physiological responses due to heat from light stimulation, a “light-only” cohort of Oxytocin-Cre mice received optic fiber implants over the PVN. These experiments were conducted at room temperature (23°C) in solo-housed females to focus on the precise physiology of PVN^OT neurons and to eliminate confounding effects of social-partners. Blue light stimulation (10 Hz; 20 ms pulse width (20% duty cycle); 10 mW power; 10 or 20 sec pulse train) occurred after solo animals had been quiescent for at least two minutes (a condition similar to when PVN^OT calcium peaks occur).

Figure 6.

We first examined the effect of blue light stimulation of PVN^OT neurons on core Tb. Approximately 400 sec (6.6 min) after a 20 sec, but not a 10 sec, pulse train, Tb increased in ChR2+ but not light-only controls (Fig. S4A). Further examination of Tb data revealed that some blue light stimulations occurred during a relatively high Tb— in contrast to the low Tb context in which PVN^OT calcium peaks naturally occur (Fig. 3N-Q). To address a possible mismatch in internal physiology during optogenetic stimulations vs. endogenous calcium peaks, we stratified the data according to whether Tb was above or below average (per-mouse) at the time of stimulation (i.e., Tb-T0). The low Tb-T0 distribution was similar to the Tb distribution observed during PVN^OT calcium peaks in solo animals (Fig. 6B blue line, cf. Fig. 3N). We then examined the effect of 20 sec optogenetic stimulation in the stratified data and found ChR2+ stimulation caused increased Tb in the low Tb-T0 context but not the high Tb-T0 context; this increased Tb was observable 420 sec after stimulation and persisted for several minutes (Fig. 6C). Thus, precisely timed, context-dependent stimulation of PVN^OT neurons is sufficient to increase core Tb.

Next, we examined the effect of PVN^OT stimulation on physical activity and behavior (Fig. 6D, S4B-C). While no significant differences were detected between ChR2+ and controls, ChR2+ mice showed a trend toward increased physical activity over time in the low Tb-T0 context (Fig. 6D). A separate slope-based analysis of post-stimulation activity revealed that ChR2+ animals exhibited a steeper increase in physical activity compared to light-only controls in the low (but not high) Tb-T0 context. (Fig. S4C). Analysis of manually-annotated behaviors suggested that PVN^OT stimulation did not activate a specific motor pattern output but instead resulted in combined increases in the time spent in nesting (linear mixed model estimate coefficient of ChR2+ stimulation: +38.0 sec), locomotion (+54.0 sec), and grooming (+14.5 sec), but not in eating/drinking (- 0.4 sec) (Fig. S4B).

We next examined the effect of light stimulation on thermal features using SGBS. Compared to light-only controls, 20 sec pulse trains triggered increased BAT and rump surface temperatures in the low (but not high) Tb-T0 context approximately 360 sec post stimulus (Fig. 6E-F). Moreover, in the low Tb-T0 context, blue light stimulation in ChR2+ animals increased the slopes of the surface temperatures of BAT and rump surface temperatures compared to light-only control animals (Fig. S4D-F). In contrast, in the high Tb-T0 context, blue light stimulation did not affect the BAT or rump surface temperatures and caused a slight decrease in the dorsal surface temperature (Fig. S4D-F).

These data suggest that 20 seconds of PVN^OT stimulation during the low Tb-T0 resting states is sufficient to trigger a sustained increase in Tb, and to a lesser extent, an increase in physical activity. Thus, PVN^OT neurons signal state dependent transitions in thermogenesis and behavioral arousal.

PVN^OT cellular projections to the rMR

We next asked what pathways could PVN^OT neurons promote transitions towards thermogenesis and physical activity. rMR efferent neurons regulate non-shivering thermogenesis, cutaneous vasoconstriction, and heart rate (Morrison and Nakamura, 2019), and the surrounding gigantocellular reticular nucleus (GiA) regulates locomotion and movement (Brownstone and Chopek, 2018; Engmann et al., 2020). In rats, PVN^OT to rMR projections (PVN^OT→rMR) stimulate thermogenesis and cardiovascular responses in an oxytocin receptor (OXTR)-dependent fashion (Fukushima et al., 2022). To investigate possible PVN^OT→rMR cell types in mice, we used tracing tools to disambiguate magno-vs. parvocellular PVN^OT projections (Althammer and Grinevich, 2018) (Fig. S5A-C). The majority of PVN^OT→rMR projections were non-magnocellular cells located in the caudal and lateral PVN, suggesting they are parvocellular (i.e., PVN^OT(parvo)) (Fig. S5C), consistent with previous studies (Althammer and Grinevich, 2018; Fukushima et al., 2022; Li et al., 2024). In addition, we observed OXTR-positive neurons in the rMR and GiA (i.e., rMR^OXTR) using smFISH (Fig. S5D). Thus, PVN^OT(Parvo) →rMR^OXTR neurons are a candidate pathway underlying transitions towards thermogenesis and behavioral arousal in mice.

Conclusions

This study uncovers a distinct functional pattern of state-dependent PVN^OT activity during natural patterns of thermoregulatory and arousal behaviors. We initially set out to examine the role of PVN^OT neurons in social thermoregulation. Leveraging long-duration in vivo recordings, we identified large PVN^OT calcium peaks that emerged approximately one hour into the home-cage recording. To our surprise, these peaks were observed in both social and non-social contexts. PVN^OT peaks were temporally aligned with transitions from resting states (quiescence and quiescent huddling) to active states (nesting, active huddling, and physical activity) and increased thermogenesis. Notably, the PVN^OT activity patterns and kinetics we observed in virgin female mice resemble the well-characterized bursts observed during the early stages of lactation (Yaguchi et al., 2024, 2023; Yukinaga et al., 2022), suggesting a broader physiological role for oxytocin neurons in coordinating transitions between thermoregulatory and behavioral states.

PVN^OT dynamics were temporally aligned to the offset of a low-Tb resting state. In the minutes prior to the end of a rest bout, baseline calcium increased and large peaks occurred in a probabilistic fashion: there was an approximately 20% chance of peaks occurring within around 100 seconds prior to quiescence offset, as opposed to ∼2.5% during the first 100 seconds of rest onset. In addition, PVN^OT calcium peaks occurred more consistently when core Tb was relatively low—approximately 36.2°C, which is about 0.5°C below the average rest-phase Tb for C57B/6 animals in our colony (Landen et al., 2024). In mice, Tb displays an ultradian rhythm characterized by limit-cycle-like oscillations between biologically constrained upper and lower bounds (Blessing and Ootsuka, 2016; Škop et al., 2024). A Tb threshold of 36.2°C may represent an approximate lower-bound of the resting state and a point at which a transition to the active-state upper-bound is initiated (Škop et al., 2024). PVN^OT neurons might therefore signal the transition toward thermogenesis and arousal during ultradian cycles of rest-to-arousal during the light phase. Supporting this interpretation, PVN^OT calcium activity and optogenetic stimulation during a low Tb—but not a high Tb—rest state signaled increases in Tb and locomotor activity within a few minutes. This state-dependency indicates that these neurons integrate information to modulate sympathetic activation (Fig. 7).

Figure 7.

Our deep-learning pipeline (SGBS) further allowed us to track thermal features (dorsal surface, brown fat, and rump) at milli-second resolution. We detected increases in BAT temperature following calcium peaks and optogenetic stimulation, supporting a role of PVN^OT neurons in driving state-dependent transitions towards thermogenesis. Several neural populations that promote the transition from sleep to wake have been described (Sulaman et al., 2023), some of which are directly involved in thermoregulation (Harding et al., 2019). In conjunction with another report showing that chemogenetic activation of PVN^OT neurons increases wakefulness (Chen et al., 2021), our data now place PVN^OT neurons as a key node in the initiation of thermogenic and behavioral arousal from rest. Here, we suggest a model, congruent with two other reports, whereby the transition from rest to arousal begins with a change in neural activity leading to increased thermogenesis that is either coincident with, or followed by, increased physical activity (Ootsuka et al., 2009; Škop et al., 2024) (Fig. 7). Notably, we found no evidence that decreases in PVN^OT activity signal transitions towards quiescence and heat loss, indicating that distinct systems drive entry to, and exit from, resting states.

Some caveats of this study should be considered. First, we did not link PVN^OT in vivo calcium activity or optogenetic stimulation with endocrinological release of OT. Future studies should address how the neuropeptide oxytocin (e.g., in addition to glutamate release (Hrabovszky and Liposits, 2008)) drives transitions in behavioral arousal and thermogenesis. Second, although some of our experiments were conducted with only a handful of animals, we nevertheless detected robust associations between PVN^OT activity and thermo-behavioral changes within individuals. Future work, incorporating more animals, could better capture variation across individuals and internal physiological states. Additionally, here we only studied females based on their discernable thermo-behavioral states; future work could focus on determining the role of PVN^OT neurons in males.

Overall, these findings extend the role of in vivo PVN^OT neural activity dynamics beyond social behavior (Hegoburu et al., 2024; Resendez et al., 2020; Sun et al., 2025; Tang et al., 2020), parturition/lactation (Perkinson et al., 2021), maternal behavior (Carcea et al., 2021), sexual behavior (Qian et al., 2023), and stress responses (Zhan et al., 2024), and highlight their contribution to daily patterns of rest/wake transitions, thermoregulation, and energy balance. In our model PVN^OT neurons respond to internal physiological conditions, such as core body temperature, to coordinate autonomic and behavioral outputs necessary for adaptive thermoregulation. These insights open new directions for understanding the neural circuits that integrate social context, internal state, and metabolic demands to regulate mammalian behavior.

Methods

Animals

OT and FOS histology

C57B6/J (Strain #: 000664) female mice were purchased from Jackson Laboratory and were 8 to 11 weeks old during the experimental period. Virgin, group-housed females were used for all trials unless otherwise specified.

PVN^OT manipulation and recording

Oxytocin-Ires-Cre mice were purchased from Jackson Laboratory (Strain #: 024234) (Wu et al., 2012). Animals were bred and maintained at the University of Wyoming. Experimental female mice ranged from 9-20 weeks old during the experimental period. All experimental mice were group-housed with ad libitum access to food and water in a room maintained at 23±2°C with a standard 12 h light / dark cycle (light: 6:00-18:00; dark: 18:00-6:00) until they were used for surgery or experiments. Animals were housed in Optimice® IVC cage systems; cages were enriched with red huts, nestlets, and a mix of cob and pine bedding. Mice used were 18-26g during the experimental period. Animal genotyping was outsourced to Transnetyx using ear tissue biopsies. All procedures were approved by the University of Wyoming Institute Animal Care and Use Committee (IACUC) in compliance with NIH guidelines for the care and use of laboratory animals.

FOS/Fos behavior assays

To measure FOS protein expression across brain regions, mice were placed into the home cage recording suite as described (Landen et al., 2024) for at least 24 hours prior to behavior being monitored, with the cage’s camera being connected to a Zoom meeting. After the 24 hours, beginning at the start of the light cycle when huddling is most frequent, a human observer monitored behavior from another room, without interrupting natural behavior. Once the target behavior was observed for a minimum of 15 consecutive minutes, 90 minutes were counted from the onset of the behavior, and animals were pulled from the recording suite for histology.

To determine Fos transcription in PVN^OT neurons during different huddling states, C57B6/J female mice that were 10 weeks old were housed in home cages in trios. Mice were habituated to the conditions, handling, and the cage lid camera set up for ∼5 days prior to recording. Animals were split into two groups: active huddle (N=6) and quiescent huddle (N=3). To capture active huddle Fos transcription, mice were recorded using the cage lid camera connected to a Zoom meeting. The time of active huddle onset was documented, and mice were sacrificed 30 minutes later. To determine quiescent huddle Fos transcription, mice were recorded as described above. Mice were sacrificed 30 minutes after the onset of quiescent or active huddle. All brains were dissected, embedded in OCT (Fisher Scientific, #23-730-571) and stored in -80°C.

Animal Surgery

Pre- and post-surgical procedures

Prior to surgery, animals were given an analgesic dose of carprofen (20 mg/kg, s.c.). Mice were deeply anesthetized using a mixture of ketamine and xylazine (100 mg/kg and 10 mg/kg, respectively, administered i.p.). Post-surgery, the mice were returned to a warm and clean home cage. For post-operative care, antibiotic ointment (Neomycin and Polymyxin B Sulfates and Dexamethasone Ophthalmic Ointment, Bausch + Lomb) was topically applied to the incision, animals received 20 mL/Kg of lactated ringer’s solution s.c., and closely observed until active. Upon emergence from anesthesia, 1 mg/kg of buprenorphine extended-release solution was administered subcutaneously between the shoulder blades. Animals were housed with cagemates and monitored daily. Carprofen was administered for three days post-surgery.

Temperature logger implants

To continuously measure internal body temperature, laparotomies were performed to implant mice with Star-Oddi DST nano-T temperature loggers. The loggers were sterilized with glutaraldehyde and rinsed with sterile PBS prior to being implanted. Mice received pre-operative care and were anesthetized as described above. A one-centimeter incision was made in the skin and underlying peritoneum, and the temperature logger was placed into the abdominal cavity. The peritoneum was sutured with PGA absorbable sutures, while the skin was closed with nylon (nonabsorbable) sutures. Mice received post-surgical care as described above.

Stereotaxic GCaMP injections and fiber implants

For fiber photometry recordings, oxytocin-Ires-Cre heterozygous, pair-housed, female mice at 8-12 weeks of age underwent stereotaxic surgery in which pGP-AAV9-syn-FLEX-jGCaMP8s-WPRE (#162377, Addgene, titer: 2.7 x 10¹³ gc/mL) was delivered to the PVN. Mice were head-fixed on a digital mouse stereotax (Stoelting, 51730D). For added analgesia, a local anesthetic (5% lidocaine) was delivered subcutaneously along the incision site before cutting. A midline incision was made along the top of the skull, the periosteum was cleared with 3% hydrogen peroxide, and the skull was dried off with canned air. Vertical measurements of bregma and lambda were taken and aligned to the same plane (DV of ± 0.05 mm). The AAV was delivered using a programmable nanoliter injector, (Nanoject III, Drummond Scientific) (300nL per side: 6 cycles of 50nL delivered at 2nL/s with 10 second delays). To allow for proper diffusion of solution into brain tissue, the needle was left in for 10 minutes before retracting. PVN coordinates relative to bregma and DV measured from the top of the brain: AP = -0.76mm, ML= ±0.19mm, DV=-4.62mm.

A fiber-optic cannula (Ø1.25mm ferrule, Ø200 µm core, NA = 0.37, L= 6.5 mm, Neurophotometrics Ltd.) was then lowered into the injection site until the tip was located ∼ 0.02 mm above the viral injection. The cannula was anchored in place with dental cement (C&B Metabond) such that a dome-shaped cap was halfway up the black ceramic ferrule and covered the entire incision. Post-operative care and recovery took place as described above. Two weeks post recovery, animals were tested for GCaMP signal using the Neurophotometrics photometry system (as described in “Fiber photometry recordings”). If a signal was detected, animals underwent a laparotomy to implant a temperature logger as described above. Animals with no detectable signal were excluded from recordings.

ChR2 injection & optic fiber implant

For optogenetic stimulation of PVN^OT neurons, Oxytocin-Ires-Cre heterozygous, pair-housed, female mice at 8-12 weeks of age were split into two different groups: ChR2+ and control (ChR2-) animals. ChR2+ animals underwent stereotaxic injections with AAVDJ-EF1a-DIO-hChR2(H134R)-EYFP-WPRE-pA (Salk, titer: 4.03×10¹³ gc/mL) using the stereotaxic surgical procedure described above. Viral vectors were delivered using a programmable nanoliter injector, (Nanoject III, Drummond Scientific) (200nL per side: 4 cycles of 50nL delivered at 1nL/s with 10 second delays). Both ChR2+ and control animals were implanted with a fiber optic cannula (Ø1.25mm ferrule, Ø200 µm core, NA = 0.37, L= 6.5 mm, Neurophotometrics Ltd.). All optogenetic animals were given a minimum of seven days to recover before undergoing a laparotomy to implant a temperature logger.

Stereotaxic DREADD, CTB and tail vein injections

To identify PVN^OT projections to rMR, a subset of mice received two stereotaxic injections: one of AAV5-hSyn-DIO-hM3D(Gq)-mCherry (#44361, Addgene, titer: 2.6 x 10¹³ gc/mL) in the PVN and one of cholera toxin subunit B conjugated with Alexa Fluor 647 (Invitrogen, C34778). Using the same surgical procedure and PVN coordinates as mentioned above, 300 nL of AAV was delivered per side (6 cycles of 50nL delivered at 2nL/s with 10 second delays). The incision was closed using absorbable PGA suture threads and antibiotic ointment was applied topically along the incision. Animals recovered for two weeks to maximize viral expression. A second stereotaxic surgery was performed to deliver a volume of 10 nL of 0.5% CTB into the rMR. RMR coordinates relative to bregma and DV measured from the top of the brain: AP = -5.80mm, ML= 0.00mm, DV= -4.95mm. Mice received post-operative care as described above and recovered for 3-4 days. For intravenous labeling of magnocellular and parvocellular OT neurons, mice received 15 uL of 4% Fluoro-Gold (Fluorochrome) diluted in 100 uL of sterile saline. Prior to injection, mice were given an analgesic dose of carprofen (20 mg/kg, s.c.). Mice were briefly restrained using a modified 50 mL conical tube, in which holes were drilled to allow for proper air flow and respiration. Mouse tails were interposed between two heating pads to enhance visibility of the tail vein. Tails were wiped down with 70% ethanol and FG was administered via either right or left lateral tail vein using a 0.5 mL 28G syringe. Mice were sacrificed 24-48 hours post-FG administration.

Histology

FOS Immunohistochemistry

For quantitative analysis of FOS expression, mice were given a lethal dose of isoflurane after removal from the home cage recording suite. Animals were then transcardially perfused with 1X PBS and 4% paraformaldehyde (PFA) in PBS. Brains were removed and kept in 4% PFA overnight at 4°C. Brains were sliced into 60 um sections using a vibratome (Leica VT1000S) and every other slice was saved in 48-well plates, wrapped in parafilm in 1X PBS, overnight at 4°C. The next day, 0.3 mL of block solution (0.20% triton, 10% natural goat serum (NGS) (Jackson ImmunoResearch 005-000-121), in 1X PBS) was placed into each well and incubated for one hour at room temperature on a plate shaker. Block was removed and replaced with freshly made primary antibody solution (c-Fos (9F6) Rabbit mAb, Cell Signaling, 14609, 1:1000 dilution in block solution) and incubated for 48 hours on a shaker at 4°C. Slices were washed three times with 1X PBS for 15 minutes. Slices were then incubated in secondary antibody solution (Goat anti-Rabbit IgG, Alexa Fluor Plus 555, Thermo Scientific, A32732, 1:500 dilution in block) for 2 hours on a shaker at room temperature and protected from light. Wash steps were repeated with 1X PBS. Slices were then mounted on Superfrost Plus slides using a paintbrush and stained with DAPI mounting medium (Vectashield H-2000). Slides were imaged with a Zeiss Axio Scan.Z1 fluorescent microscope. Regions were identified by finding all locations that displayed FOS expression and cropped using ImageJ alongside the Allen Brain Atlas. The quantity of cells and overlap was calculated using a custom CellProfiler 4.2.1 script.

Fluorescent in situ hybridization (RNAscope)

For quantitative analysis of Oxt and Fos overlap in PVN^OT mice were anesthetized with a lethal amount of isoflurane. Brains were extracted, embedded in OCT and stored at -80°C. Brains were sectioned at 16 µm using a cryostat and subsequently stored at -80°C. For in situ, slides containing brain sections were placed in 4% PFA for 45 minutes and then were dehydrated in 50%, 70% and 100% ethanol for 5 minutes each at RT. Slides dried for 5 minutes after dehydration and hydrogen peroxide was added for 10 minutes at RT. Slides were then treated with Protease III for 30 minutes at RT to allow for antigen accessibility. The probe cocktail contained Oxt (ACD Bio, #493171) and Fos (ACD Bio, #316921) and was diluted at 1:2 ratio with probe diluent. An aliquot of 40µL of probe cocktail was added to each section and incubated in a 40°C oven for two hours. For amplification, 40µL of AMP 1 was added each slice and placed in the 40°C oven for 15 minutes. AMP 2 was then added to each slice and slides were again incubated at 40°C for 15 minutes. Finally, 2-3 drops of AMP 3 were added to each slice and slides were incubated at 40°C for 30 minutes.

Fluorophores were added to each channel and diluted at 1:1500. Each fluorophore was added by first opening the channel using 2-3 drops of HRP on each slice and incubating at 40°C for 15 minutes. The fluorophore was added at a volume of 40 µL to each slice and slides were incubated at 40°C for 30 minutes. To close the channel, 2-3 drops of blocker were added to each slice and incubated for 15 minutes. Lastly, a drop of DAPI (ACD Bio, #320858) was added to each slice for 30-45 seconds. DAPI was allowed to dry for approximately 3 minutes before immediately coverslipping with 1-2 drops of ProLong Gold (Invitrogen, P36930). Slides were protected from any surrounding light, left to dry overnight and imaged with a Zeiss Axio Scan.Z1.

AAV expression

To verify AAV injection sites, animals who had received stereotaxic DREADD, GCaMP8s, or CHR2 injections were sacrificed at ≤18 weeks. Mice were deeply anesthetized with a lethal amount of isoflurane and perfused with 4% PFA. Extracted brains were then stored in 4% PFA overnight at 4°C before being sliced at 60 µm using a vibratome (Leica VT1000S). PVN regions were identified via the Allen Brain Atlas. PVN slices were stained with a DAPI mounting medium (200 µl/slide, Vectashield H-2000) and scanned using a Zeiss Axio Scan.Z1. Images were enhanced and/or pseudo-colored using ImageJ.

Fiber photometry recordings

Recording setup

The Neurophotometrics fiber photometry system (FP3001), operating with Bonsai (Version 2.72) software, was used to acquire photometric data with concurrent video recordings (Arducam, OV9281). Patchcords (MBF Bioscience, NPM-BPC-4) were attached to the fiber ferrules, and held in place by either ceramic sleeves or by a quick-release interconnecter (ThorLabs, ADAL3). Calcium dependent GCaMP fluorescence signals were recorded with light from a 470-nm LED that was band-pass filtered, narrowed with a collimator, reflected by a dichroic mirror, and focused on the sensor of a CCD camera by a 20x objective. To account for auto-fluorescence and possible motion artifacts, a 415-nm LED for calcium-independent isosbestic fluorescence signal was used. The recorded signals were acquired with a sampling rate of 30 Hz for each channel at 100% duty cycle. LED light power was calibrated prior to every recording using a power meter (ThorLabs, PM100USB and S120C). LED light was delivered at a power that resulted in 60µW of 415-nm light, and 85 µW of 470-nm light at the tip of the patchcord.

Recording conditions

Photometric recordings occurred in the home cage, with the cage-lid replaced by a custom-fit metal sleeve that extended the cage wall height to 16”. Optic fiber-implanted mice and their pair-housed sister were habituated to the connector and patchcord prior to commencement of experimental schedule. Habituation periods ranged from several days to a week depending on individual animal reaction to the patchcord and to sufficiently aversion train the cagemate. Aversion training was performed during habituation periods by lightly coating the patchcord with diluted capsaicin (EMD Millipore, 211275). The cagemate and implant mouse were observed, and training was repeated until the cagemate associated the unpleasant sensation of capsaicin with chewing of the patchcord. Recordings were approximately 2-hours and were conducted at three different floor temperatures (15C°, 23C° and 29°C), with and without the sister cage-mate. Floor temperature was controlled by placing the home-cage onto a custom aluminum chassis that attached to the Hot/Cold Plate (Ugo Basile, 35300), temperature was then validated by a temperature logging iButton (DS1922L-F5, Embedded Data Systems) placed inside the home-cage. A portion of the nesting material was removed before experiments to prevent occlusion of the animals during the recording (see Supplementary Video 1 for typical amount of bedding left in cage).

Validation trials

To verify the bursting activity of PVN oxytocin neurons observed in photometry recordings, virgin females were crossed with stud males after completing at least two photometry trials (as described above). The dams and pups were housed in their home cages throughout the trial period. Dams were recorded during nursing and lactation. Maternal oxytocin neural activity was recorded at the following timepoints: postpartum day (PPD) 2, 5, 7, 8, 9, and 14. For this experiment the minimum peak amplitude was set to be > 7 standard deviations (sds) above the mean. This threshold was obtained by examining peaks in mothers that clearly displayed pulsatile activity bursts during lactation and closely resembled the peak dynamics that have been previously described during lactation (Yaguchi et al., 2024, 2023; Yukinaga et al., 2022).

Behavioral characterization

Videos were uploaded to Noldus Ethovision XT 16, where the physical activity of the mice (calculated total percentage of pixel change between frames) and behaviors were identified. Activity in the paired context was defined as the summed activity for both animals.

The behaviors for solo photometry animals were manually scored as follows: Locomotor Activity (LMA), Eating or Drinking (EaDr), Grooming (Groom), Nesting or Nest Building Nest (Nest) where the animal engaged in nest-building behavior or activity within its nest other than grooming, Quiescence (Quies) where the animal is inactive inside its nest, and Stationary (Sta) which was defined as animal inactivity outside of the nest.

The behavior of the photometry animal for paired assays was manually scored as follows: Locomotor Activity (LMA), Eating or Drinking (EaDr), contact initiated (ConI), contact received (ConR), Grooming (Groom), Nesting or Building Nest (Nest) characterized by the photometry animal engaging in nest-building behavior or its activity inside the nest while its cagemate engaged in nest-building behavior, Active huddle (AHud) defined as physical contact and activity between both animals in the nest, Quiescent Huddle (QHud) referred to physical contact and inactivity between both animals in the nest, Quiescence (Quies) where the photometry animal was engaging in inactivity alone in the nest, and Stationary (Sta).

Behavior and activity data was exported at 10 frames per second as a .csv file and imported to R for analysis.

Calcium data preprocessing

Fluorescent signal data preprocessing was done with a custom R script as follows: 415 and 470 channels de-interleaved; 470 channel filtered to remove high-frequency noise using a 4th-order Butterworth filter (Van Boxtel et al., 2021) modified to remove end-effect transients; trimmed to remove approximately the first six minutes of data (to mitigate the effect of photobleaching); if necessary, further trimmed to remove effects of disruptions to the patchcord (chewing or twisting); smoothed with a sliding window function (rollapply; width = 9) (Zeileis et al., 2004).

To correct the calcium-dependent signal from effects of photobleaching and heat-mediated LED decay, we first extracted the “fitted values” from the linear relationship between the 470 (calcium-dependent) and 415 (calcium independent) fluorescent signals. Next, to calculate dF/F we divided the 470 signal by these fitted values. Because we observed occasional artifactual shifts in baseline dF/F values, we baseline corrected dF/F using the iterative least squares method (Liland and Mevik, 2011). The Z-scored dF/F was calculated as Z = (χ-μ)/ σ, where μ is the mean dF/F and σ is the standard deviation of dF/F. To define peaks, we used the findpeaks function (Borchers, 2011), with the minimum peak distance set at two seconds and the minimum peak height set at six standard deviations above the mean dF/F. Photometric data were aligned to behavioral data (exported from Noldus Ethovision) and Tb data (exported from Mercury) using common centisecond-resolution timestamps and reduced to a sampling frequency of 10 frames per second.

Peri-event data oriented on specific events (i.e., calcium peaks or start/stop frames of behavioral bouts) were extracted using a custom function that (1) extracts the indices of the events and (2) generates a list of vectors for a specified number of rows before and after the event.

To compare longitudinal physiological patterns of activity in PVN^OT neurons in virgin females (primary subjects of this study) compared to nursing and lactation, we made the following modifications to the protocol as described above. First, to quantify PVN^OT peak characteristics (amplitude, full width half maximum, interpeak interval) before linear normalization and baseline correction, (i.e., to better describe the dynamics of their native characteristics), we corrected the calcium-dependent signal from effects of photobleaching and heat-mediated LED decay using an exponential decay model. Specifically, we used a self-starting nls (nonlinear least squares) function in R to fit the isosbestic data with the following model: y(t) ∼ SSasymp(timeS, yf, y0, log_alpha), where the measured isosbestic value y starts at y0 and decays towards yf at a rate α. SSasymp is a shortcut that guesses its own parameters; instead of fitting the rate constant α directly, it searches for the logarithm of α: y(t)∼yf+(y0−yf)e^−exp(logα)t. We then linearly scale this fitted decay to the calcium-dependent data using robust regression (MASS::rlm with the psi function set to bisquare); finally, we divided the calcium-dependent data by this scaled fit to get a corrected signal.

In accordance with previous work, we examined three epochs in the same animals: virgin females; early-stage lactation (PPD 2-7); and late-stage lactation (PPD 8-14). To calculate the full-width half maximum (FWHM) of PVN^OT peaks, we used a custom function that performs the following steps: finds each peak; determines the half-max value; interpolates where the signal crosses the half-max value; computes the FWHM as the difference between these crossing points. To calculate the interpeak interval (IPI), we used the lead function in R to compute the time (in seconds) from one peak to the next on a per-recording basis.

Computer vision analysis

Processing Flir .seq files

RGB videos, timestamps, and thermal data were extracted from Flir thermographic .seq files using a Python-based adaptation of ThermImageJ ³. Raw images were obtained via Flirpy ⁴, converted to RGB, and merged into AVI video files. Metadata, including timestamps and frame rates, was extracted using Exiftool ⁵. This workflow provided thermal images for model training, AVI videos for analysis, and timestamps and temperature data for integration with behavioral information. The processing script is available at https://github.com/j-landen/seq_process.

Custom Mask R-CNN for body part segmentation

To extract surface temperature from thermographic videos, we developed a computer vision pipeline (SGBS) integrating keypoints from DeepLabCut ⁶ with a custom Mask R-CNN ⁷ model for segmenting specific body regions. A DeepLabCut deep learning model was first trained on manually labeled data from 1800 images to detect ten keypoints along the dorsal surface of the mice (nose, implant, head tip, neck base, center back, tail base, mid tail, patchcord, and right & left shoulder). These keypoints were then incorporated into an attention layer alongside the thermal images to enhance texture representation and improve segmentation accuracy. The custom Mask R-CNN was trained using 457 images on Darwin (V7 Labs) to segment the BAT region, the rump region, and the entire dorsal surface of the mouse. Finally, segmented regions were mapped to their corresponding pixel-wise temperature data, providing per-frame temperature measurements (mean, min, max, SD) for each target region. The SGBS code is available at: https://github.com/j-landen/SGBS.

To determine the efficacy of the custom segmentation software, it was compared with an unmodified Mask R-CNN model, using the same pre-trained weights (ResNet50_Weights IMAGENET1K_V2) and the same 450 manually-trained images (80% training, 20% validation). Both models ran over 200 epochs, and loss was measured after each epoch. Loss is calculated as the sum of three components: (1) a classification loss (cross-entropy) that assesses how accurately each body part is identified, (2) a bounding box regression loss that measures how close the predicted box is to the ground truth, and (3) a mask loss (binary cross-entropy) that evaluates pixel-by-pixel agreement between the predicted mask and the ground truth mask ^7,8.

Optogenetic stimulation

Recording setup

Blue light (450 nm) stimulation was delivered using the Neurophotometrics fiber photometry system (FP3001), with simultaneous video and thermal recording. Mono-filament patchcords (Doric Lenses, MFP_200/220/900-0.37_3m_FC-MF1.25) were attached to the fiber ferrules and held in place by a quick-release interconnecter (ThorLabs, ADAL3). Optogenetic laser power was measured prior to every recording using a power meter (ThorLabs, PM100USB and S120C). Optogenetic activation used a frequency of 10 Hz, with a 20 ms pulse width (20% duty cycle) and 10 mW of power measured at the tip for either 10- or 20-second trains.

Optogenetic recordings occurred in the home cage, with the same setup as fiber photometry recordings, described above. Mice were pair-housed until recording. One animal was tested at a time, while the other animal was placed in a holding cage for the duration of the experiment. During experiments, animals were observed remotely from a different room. The human observer turned on blue light stimulation after the animal had been inactive for a minimum of two minutes. Stimulations were a minimum of ten minutes apart from one another. Experiments were two hours long, averaging four to seven stimulations. Every animal was tested at least twice using each stimulation time (10s and 20s), until a minimum of 10 total stimulations occurred per animal, per stimulation time.

Optogenetic data preprocessing

For each stimulation, data spanning 5 min (–300 s) before to 10 min (+600 s) after light onset were extracted, with the onset designated as time 0. Core body temperature (Tb) data was recorded once every five seconds for the duration of the experiment. We therefore subset all other data to the same rate for core Tb analyses.

Activity levels were measured at a rate of 10 Hz, so all data were subset to the same rate for activity analyses. As activity is measured using percent pixel movement over the entire image, activity was first set to a logarithmic scale, to measure small differences in change, then normalized using a Z score across mouse IDs to account for differences in baseline activity levels.

Feature-specific surface temperatures using SGBS were measured at a rate of around 20-30 Hz, depending on the Flir A50’s output. As previously described, all surface temperature data was subset to 10 Hz and aligned to surface temperature data based on the nearest timestamps when the laser was turned on (< 0.1s deviation). To eliminate the effect of outliers originating from an incorrect segmentation output, all surface temperature data were smoothed with a rolling average over one half of a second. To eliminate differences in baselines per animal, every value was baseline-normalized using ΔT/T, where each temperature value was expressed as a change from the pre-stimulation mean (3-minutes prior to light on) relative to that baseline.

To assess whether optogenetic effects varied with initial thermal state, core body temperature at the onset of each blue-light pulse (Tb-T0) was recorded. Because Tb-T0 values were generally higher than those observed during spontaneous PVNOT calcium peaks, stimulations were stratified for each mouse into “low Tb-T0” and “high Tb-T0” trials using that animal’s mean Tb-T0 as the cutoff. The two strata were then analyzed separately.

Statistics

Statistical results are reported in Table S1

Effect of behavior state on FOS expression

To determine the effect of behavior on FOS expression, brain ROIs (divided into left and right sides, where appropriate) were scanned by CellProfiler 4.2.1 to determine the total number of cells, as marked by DAPI. DAPI ROIs were shrunk to a point and then expanded by a factor of 6, and classified as being either FOS-positive or FOS-negative, as marked by presence of the secondary antibody (Alexa Fluor 555), using the same parameters for every slice. These per-ROI portions were then standardized using a Z-score between different histology days, to account for any batch effects of antibody binding or fluorescence. A linear regression model was used on a per-ROI basis, and a post-hoc tukey test was used to determine the effect of behavior on FOS expression, for each brain region. P-values were adjusted for multiple comparisions using the Holm method.

Relationship between PVN^OT calcium activity and behavior state, activity, and body temperature

To determine the effect of different behavioral states on calcium peak amplitude and frequency, we fit LMMs with amplitude or frequency as the dependent variable, behavioral state as the independent variable, and mouseID as a random effect. These models were run for both solo and paired conditions.

To determine the effect of social context and floor temperature on calcium peak amplitude and frequency, we fit LMMs with amplitude or frequency as the dependent variable, social context or floor temperature as the independent variable, and mouseID as a random effect.

To determine how calcium peak frequency is affected by behavioral state bout length (i.e., for quiescence and quiescent-huddle) and floor temperature, we fit a LMM with per-bout peak counts as the dependent variable; bout-length, floor temperature, and their interactions as the independent variables; and mouseID as the random effect. We subsequently extracted the standardized beta coefficients using the sjPlot package (Lüdecke, 2013).

To determine how locomotor activity changed during a 300 sec. interval of time surrounding calcium peaks (i.e., peri-peak time), we first z-scored the activity data (estimated at 10 Hz), then smoothed it according to the mean using the rollapply function (Zeileis et al., 2004), and finally calculated the per-second mean for each individual. Next, we calculated normalized activity data across the 600 sec interval, including the mean normalized activity before and after each peak (i.e., “before-after”). Last, we fit a LMM with normalized activity as the dependent variable, “before-after” as the independent variable, and mouse-ID as a random effect.

To determine peak probability during defined behaviors we used logistic regression. First, we used custom R code to (1) extract the onset and offset of behavioral bouts, and (2) label rows of data (at 10 frames per second) according to whether they occurred within a minute containing (or not containing) a calcium peak (i.e., PeakMinute). Next, we fit a binomial generalized linear model with PeakMinute as the dependent variable; the independent variable was the time span before, during, and after the behavioral epoch of interest (e.g., quiescence offset). Next, we extracted the per-frame predicted values (including standard error) from this model.

To compare the statistical distributions of Tb during calcium peaks vs. baseline, we plotted Tb according to the PeakMinute designation described above. We then used a LMM with Tb as the dependent variable; the independent variables were PeakMinute and mouseID as a random effect.

PVN^OT peak characteristics: before mating vs. early/late stages of lactation

To determine the effect of life-history phase (virgin, early-stage and late-stage lactation) on peak amplitude, FWHM, and interpeak interval we used LMMs with epoch (PPD 2-7; PPD 8-14; virgin) as a fixed effect and mouseID as a random effect to control for repeated measures. We performed posthoc tukey comparisons for each pairwise combination of fixed effect level using the emmeans package.

Effect of optogenetic stimulation on thermal profile

To determine the effect of light stimulation on Tb in ChR2+ and light-only controls separately, Tb data were binned per minute, with each value representing the mean of the 60 seconds following. We used a linear mixed-effects model (LMM) with Tb as the dependent variable, the interaction between light-status and minute as independent variables, and mouseID as a random effect. A post-hoc tukey test was used to determine per-minute differences.

To determine the effect of light stimulation on physical activity in ChR2+ vs. light-only controls separately, activity was binned per minute, and we used the same LMM as above, but with activity as the dependent variable.

To determine the effect of light stimulation on surface temperatures (BAT, rump, or dorsal surface) in ChR2+ vs. light-only controls separately, each surface temperature was analyzed separately, first by binning temps per minute, and then using the same LMM as above, with surface temps as the dependent variable.

To further determine the effect of light stimulation on direction of temperature change, slopes before and after light onset were analyzed with a second LMM that included fixed effects of time, proximity to stimulation (pre vs. post), and their interaction, plus a random intercept for mouse ID; frame-level weights were set to the inverse of each frame’s temperature SD. A significant interaction term indicated that temperature trajectories differed between ChR2+ and control groups across the stimulation boundary.

To assess the effect of light stimulation on behavior in ChR2+ versus light-only controls, we quantified cumulative time spent in each behavioral state following stimulation. First, linear regression models were applied to individual arousal behaviors to evaluate their contribution to group differences. As no single behavior showed a significant effect, all arousal behaviors were grouped, and a second linear regression model was used to assess whether total arousal time differed between groups.

Data and code availability

Raw data and code will be made available on a public repository.

Acknowledgements

We thank University of Wyoming Sensory Biology Center for input; Robert Carrol for assistance with animal husbandry; members of the Nelson and Bedford laboratory for input; Sean Harrington for coding assistance; the UW Engineering Shop for designing and building hardware.

Additional information

Resource availability

Lead contact

For further information or requests, contact should be directed towards the lead contact, Adam Nelson (anelso74@uwyo.edu).

Materials availability

This study did not generate any new reagents.

Author contributions

M.V., J.G.L, J.F.R., N.L.B., and A.C.N. designed the study. M.V. performed the fiber photometry studies. M.V, J.G.L, and A.C.N. performed the optogenetic studies. M.V., J.G.L., B.A., C.P., and A.C.N. performed the histology studies. J.G.L, S.K. and G.J.T. performed code development. J.F.R. contributed to a previous draft of this manuscript. N.L.B. and A.C.N secured funding and provided expertise and feedback.

Additional files

VideoS1. Photometry recording of virgin mouse in the homecage. A large-amplitude PVN^OT peak occurs at 0:00:03s. Video is sped up to 4x.

VideoS2. Photometry recording of a lactating dam in the homecage with pups on PPD 14. The first large-amplitude PVN^OT peaks occur at approximately 0:00:04s. An even larger peak typical of late-stage lactation occurs at 0:00:15s. Video is sped up to 4x.

VideoS3. SGBS enables real-time segmentation of thermally defined anatomical features in freely moving mice. Example video shows the output of the SGBS model, which accurately identifies and tracks three anatomical regions – BAT, rump, and dorsal surface – in a freely moving mouse in the homecage using thermal imaging. Video is sped up to 4x.

Note

This reviewed preprint has been updated to correct level headings and formatting of figure supplement titles.