Diurnal brain temperature rhythms and mortality after brain injury: a prospective and retrospective cohort study

Objective: To determine the clinical relevance of brain temperature (TBr) variation in patients after traumatic brain injury (TBI). Design: Cohort study with prospective (healthy participant) and retrospective (TBI patient) arms. Setting: Single neuroimaging site in the UK (prospective arm); intensive care sites contributing to the Collaborative European NeuroTrauma Effectiveness Research in TBI (CENTER-TBI) High Resolution ICU (HR ICU) Sub-Study (retrospective arm). Participants: 40 healthy adults aged 20-40 years recruited for non-invasive brain thermometry and all patients up to May 2020 that had TBr measured directly and were not subjected to Targeted Temperature Management (TTM). Main outcome measures: A diurnal change in TBr (healthy participants); death in intensive care (patients). Results: In healthy participants, mean TBr (38.5 SD 0.4{degrees}C) was higher than oral temperature (36.0 SD 0.5{degrees}C), and 0.36{degrees}C higher in luteal females relative to follicular females and males (95% confidence interval 0.17 to 0.55, P=0.0006 and 0.23 to 0.49, P<0.0001, respectively). TBr increased with age, most notably in deep brain regions (0.6{degrees}C over 20 years; 0.11 to 1.07, P=0.0002). The mean maximal spatial TBr range was 2.41 (SD 0.46){degrees}C, with highest temperatures in the thalamus. TBr varied significantly by time of day, especially in deep brain regions (0.86{degrees}C ; 0.37 to 1.26, P=0.0001), and was lowest in the late evening. Diurnal TBr in cortical white matter across participants ranged from 37.0 to 40.3{degrees}C. In TBI patients (n=114), mean TBr (38.5 SD 0.8{degrees}C) was significantly higher than body temperature (TBo 37.5 SD 0.5{degrees}C ; P<0.0001) and ranged from 32.6 to 42.3{degrees}C. Only 25/110 patients displayed a diurnal temperature rhythm; TBr amplitude was reduced in older patients (P=0.018), and 25/113 patients died in intensive care. Lack of a daily TBr rhythm, or an age increase of 10 years, increased the odds of death 12-fold and 11-fold, respectively (OR for death with rhythm 0.09; 0.01 to 0.84, P=0.035 and for death with ageing by 1 year 1.10; 1.05 to 1.16, P=0.0002). Mean TBr was positively associated with survival (OR for death 0.45 for 1{degrees}C increase; 0.21 to 0.96, P=0.040). Conclusions: Healthy TBr exceeds TBo and varies by sex, age, menstrual cycle, brain region, and time of day. Our 4-dimensional reference resource for healthy TBr can guide interpretation of TBr data in multiple clinical settings. Daily temperature variation is frequently disrupted or absent in TBI patients, in which TBr variation is of greater prognostic use than absolute TBr. Older TBI patients lacking a daily TBr rhythm are at greatest risk of death in intensive care. Appropriately controlled trials are needed to confirm the predictive power of TBr rhythmicity in relation to patient outcome, as well as the clinical utility of TTM protocols in brain-injured patients. Registration: UK CRN NIHR CPMS 42644; ClinicalTrials.gov number, NCT02210221.

Establishing how TBr varies in health is critical; deviations from normal may have transformative diagnostic and/or prognostic value in neurological disease and injury, but only if these deviations can be distinguished from physiological variation over time. 22 With magnetic resonance spectroscopy (MRS), spatially resolved TBr data can now be obtained through non-invasive brain imaging. 13 Brain thermometry has proven to be a powerful research application of MRS but, with respect to healthy humans, it has only been used in studies that were poorly controlled for parameters that influence physiological temperature variation (Supplementary Table S1). We sought to establish the daily spatiotemporal variation of healthy TBr to enable evidence-driven appraisal of the clinical value of TBr monitoring in brain-injured patients. We hypothesized that healthy TBr would vary diurnally, and that disruption of diurnal temperature variation would be associated with outcome after TBI.

Prospective study design and recruitment
We conducted a prospective, single-site, cohort study in healthy adults, controlled for age, sex, body mass index (BMI), menstrual cycle phase, seasonal variation, and individual chronotype. Our primary objective was to determine whether healthy TBr varies by time of day. Our secondary objectives were to compare variability in brain and oral temperatures, to test for differences between males and luteal-phase females, and explore brain-regional changes in TBr with time. We hypothesized that TBr would (i) exceed and vary more than oral temperature across the day, (ii) be higher in luteal females relative to males, and (iii) increase with increasing brain tissue depth.
Sample size was estimated for achieving the primary outcome (a change in mean global TBr between time points) using a linear mixed model, considering published data on the reliability of MRS brain thermometry in healthy men. 13 With 36 subjects, and a conservative true mean TBr difference of . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint 0.5°C, we estimated 80% power to detect a statistically significant difference between time points at the 5% significance level. A health-related finding (HRF) on MRI was the key exclusion criterion.
Completion of a feedback pathway for two volunteers was expected (based on 5% prevalence of HRFs using high-resolution MRI). 23 We aimed to scan 40 eligible participants (20 females) to account for potential withdrawal, exclusion, and/or technical scan failure.
Recruitment for our Circadian Brain Temperature (CiBraT) Study was based on meeting criteria for our primary outcome (Supplementary Table S2), and was conducted locally using mailshots to University of Edinburgh and NHS staff, social media posts, and posters displayed at University of Edinburgh campuses and NHS Lothian hospitals. By completing an online eligibility questionnaire, all prospective participants provided written informed consent for their personal data to be used to schedule consenting interviews, and to notify general practitioners of their intention to participate.
The questionnaire provided access to inclusion and exclusion criteria, the Participant Information Sheet and Consent to Participate Form (Supplementary Appendix 1), and Data Protection Information sheet. All participants provided written, informed consent to participate during face-toface interview conducted by the Chief Investigator (NMR) at the University of Edinburgh.
Additional written informed consent was obtained for publication of individual data which, by nature of its distinctive features, could potentially be recognized by participants as their own data. The Study Protocol is presented in Supplementary Appendix 2.

Prospective data collection
During a consenting interview at the study site, one week in advance of scanning, each participant was given a wrist-worn actimeter (ActTrust2, Condor Instruments, Sao Paulo, Brazil). Each participant then underwent three identical brain scans in the morning (9-10am), afternoon (4-5pm), and late evening (11pm-midnight) of their scheduled scanning day. Multiple time points spanning >12 hours were selected because the human circadian rhythm (body clock) impacts almost every . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint aspect of physiology (Supplementary Text 'Internal rhythms and health'). [24][25][26][27] The exact alignment, or phase relationship, between the body clock and the day-night cycle is dictated by individual chronotype, which is determined by genetic and lifestyle factors, and can be derived from longitudinal monitoring of locomotor activity. 28 To assign scan times to the appropriate part of each participant's circadian cycle, we determined individual chronotypes using wrist actigraphy to extract the sleep-corrected midpoint of sleep on free (non-work) days (MSFsc) (Supplementary Methods). 28 Height and weight were measured immediately before the morning scan to calculate BMI. Oral temperature was recorded before each scanning session using a digital Clinical Thermometer Limited, Middlesex, UK) and placed sublingually. For females, hormonal influences were controlled through urine-based ovulation testing (ClearBlue®), or documenting hormonal contraception type.
We aimed to scan females during the luteal phase of their natural menstrual cycle, or on a day when an active combined pill would be taken, or combined patch worn. Females using other forms of contraception (implant or intrauterine device) were excluded. On the day of scanning, food consumption was restricted to 6am-8am, 12noon to 2pm, and 6pm-8pm, and caffeine consumption was restricted to 6am-8am and 12noon to 2pm. Alcohol was strictly prohibited at all times.
Participants were asked not to participate in excessive physical activity on the day of scanning. Data collection was limited to a 14-week period between July and October 2019 to avoid daylight savings clock changes and large seasonal variation in environmental light and temperature conditions. Data management procedures are described in Supplementary Methods.

Brain imaging
All brain imaging was conducted at the Edinburgh Imaging (Royal Infirmary of Edinburgh) Facility using a 3-T MAGNETOM Prisma scanner (Siemens Healthcare, Erlangen, Germany) with a 32channel head coil. All participants were screened for MRI contraindications and changed into . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint hospital scrubs for each scanning session-conducted in a temperature-controlled room (target 21.5°C). Room lights were off and the scanner lighting and fan were maintained on their lowest setting. Ear protection was provided and a mirror was attached to the head coil so participants had the choice of closing their eyes or viewing the MRI control room; no visual or acoustic entertainment was provided. Participants were permitted to sleep during scans and were asked to report on this event at the end of each session. At each time point, after whole-brain structural MRI, MRS data were collected from 82 brain locations (voxels). The scanning protocol was well tolerated, with no serious adverse events reported during 7-day follow-up. Further details on the scanning protocol and MRS data processing are provided in Supplementary Methods and Appendix 3; the dedicated Study Participant Data Form (Case Report Form) is provided in Appendix 4.
Calculation of MRS-derived brain temperature MRS brain thermometry exploits the fact that the chemical shift of water is exquisitely temperaturedependent (-0.01ppm/°C), whilst that of the reference metabolite NAA is not. 29 The chemical shift difference between water and NAA can estimate absolute TBr in healthy men with a short-term precision of 0.14°C at 3-T. 13 TBr for each brain tissue voxel in this study was calculated using the following relationship: TBr = 100*[NAA frequency -H2O frequency +2.665] + 37 where frequency is in ppm and temperature is in °C The reliability and accuracy of TBr determination using this MRS protocol was thoroughly tested using in vivo human and in vitro phantom measurements; the latter validated with an MR-compatible industrial thermometer that meets international standards. 13 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

Retrospective study design and patient data sources
To determine the clinical relevance of TBr variation, we conducted a multicentre, retrospective cohort study of TBI patients that had high temporal-resolution TBr data collected directly from the brain.
Data for all eligible patients were extracted using version 2.0 of the CENTER-TBI dataset, compiled between 2015 and 2017. Additional eligible patients monitored at one of the contributing sites (the Intensive Care Unit, Western General Hospital, Edinburgh, UK) were included up to May 2020 and comprised 109 of the 134 eligible patients screened. The Western General Hospital is the tertiary referral centre in South East Scotland for neurosurgical emergencies. Patients with moderate to severe TBI admitted to intensive care requiring intubation, sedation, and intracranial pressure (ICP) management also received brain oxygen tension and temperature monitoring using the Integra Licox system (Integra, France). Patients were managed in accordance with Brain Trauma Foundation guidelines. 30 Patients were either admitted directly to intensive care or following surgical intervention for mass lesions. TBr was measured via a thermistor, inserted into the brain parenchyma via a dedicated bolt placed via a burr hole (Integra Neurosciences, Andover, UK). The bolt was placed so that the thermistor inserted into frontal white matter; for diffuse injuries this was into the non-dominant hemisphere. When the main injury was focal, the bolt was placed on the side of maximal injury, unless this would place the monitors into non-viable tissue. High temporalresolution physiological data were recorded at a minimum of 1-minute intervals to either a bedside computer running ICU Pilot software (CMA, Sweden) or to a Moberg neuromonitoring system (Moberg Research Inc., USA). Data were collected continuously (except for interruptions due to computed tomography scanning or surgical intervention) and until ICP monitoring was no longer required, or the patient died. Data for the CENTER-TBI study were collected through the Quesgen e- is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint system (Moberg Research Inc., USA) were used. For TBo, the primary method of measurement was documented in 26 of 134 screened patients and included tympanic (21), bladder (3), external axillary (1), and nasopharyngeal (1). Secondary sites included rectal, external axillary, oesophageal, and skin.

Patient temperature data processing
Four inclusion criteria levels were applied to ensure that sufficient temperature data were available to assess for a diurnal rhythm ( Table 1), and that any data potentially affected by TTM protocols were excluded. Analysis of patient temperature data was blinded to outcome. Data from the first 2 hours of monitoring were excluded from the analysis to ensure the results were not influenced by the time required for the electrode to stabilise. Raw data processing was performed in Excel to exclude artefactual data, identify any gaps in the time series, and define the analysis window. Temperature data were visualized in GraphPad Prism version 8.2 and assessed for the presence of daily rhythmicity. Visual analyses were validated with a combination of rhythm-detection algorithms using GraphPad Prism, BioDare2 (biodare2.ed.ac.uk) 31 and the Harmonic Regression package in R. 32 To be categorized as diurnally rhythmic, the patient's temperature pattern need not be conventionally aligned with the day-night cycle, but it had to meet both of the following criteria: (1) a period length of ~22-26h in at least part (but not necessarily all) of the time series as determined by blinded visual analysis of the raw data in GraphPad Prism and (2) a period length of 22-26h as determined by period analysis in (i) cosinor analysis in GraphPad Prism and/or (ii) statistically significant output from Harmonic Regression in R and/or (iii) BioDare2.
In GraphPad Prism, period results were only considered valid if a cosinor curve fit was significantly preferred over a straight line. When using the Harmonic Regression package, the period length term (Tau) of the model to test for was set to 24 hours. In BioDare2, period analysis was performed using . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint 6 different algorithms. A full description of these algorithms can be found at https://biodare2.ed.ac.uk/documents/period-methods.
Level A: Criteria for extracting maximum and minimum daily brain and/or body temperatures • Known sex • Known age • Minimum 24 hours of temperature data collection under 'constant' conditions. The first data point recorded in the intensive care setting that exceeds the minimum recorded temperature from that patient in the absence of TTM will be taken as the start point (to exclude low temperature points surrounding insertion of probe or those relating to patient hypothermia on arrival in intensive care) • For data where only the maximum and minimum daily TBo (and in some cases TBr) are recorded with their respective times, a minimum of two days' worth of data is needed • If TTM was applied, only data relating to time preceding TTM or after the first inflection of data after cessation of TTM can be used and must meet the above requirements for minimum time length in the absence of TTM • When extracting the time of the minimum and maximum temperature point, the first occurrence of that specific temperature point under intensive care 'constant' conditions will be selected Level B: Additional criteria for performing diurnal rhythmic temperature analyses • Minimum hourly TBo or TBr data with TBr data extracted via intracranial probe (standard depth and positioning in cortical white matter) recorded continuously over a minimum of 36 hours. The same rules as above apply in relation to TTM. Ideally injury type (focal/diffuse; from CT scoring) and/or severity (IMPACT imputed GCS) on admission to Study Hospital and/or TIL score (including individual components of this) • Ideally site of probe insertion for focal injury (ipsilateral or contralateral to injury-to be determined using CT/MRI images if available) is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

Statistical analysis
To determine healthy temperature variation, we applied a linear mixed modelling approach. The fixed effects (predictors) were specified a priori based on published literature describing factors that were most likely to affect body and/or brain temperature in humans and other mammals. 22 In each case, the upper limit for the number of fixed effects was set to a maximum of five to avoid overfitting each model within the confines of our sample size. 33  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

BrainTempij = [intercept (β0) + Time (β1) + Sex (β2) + BrainRegion (β3) + Age (β4) + Sleep (β5)] + εij (residuals for subject i at time j) + U1i (intercept for subject i) + U2i (slope for subject i in relation to time)
where fixed effects include: Survival in intensive care or 'alive' was specified as a miss, and death or 'dead' as a hit. The model incorporated fixed effects and random effect for intercept and was built as follows: . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

BrainRange (TBr range across analysis window)
Diurnal (presence or absence of a daily temperature rhythm within analysis window-categorized as 'yes' or 'no'; see above for details on how tests for diurnal rhythmicity were performed) with random effects for intercept by subject The final choice of fixed effects (predictors) to include in the model was based on our core study objectives, avoiding redundant terms, and optimising the model fit (Appendix 5). Missing data values for any of the model components were input as 'NA', and thus patients with values missing for one or more of the components were excluded from the model output. The most conservative approach was taken i.e. multiple imputation was not performed since the random nature of missing data could not be assumed.
Statistical modelling and other circular analyses were performed using R version 3.6.3 (R Core is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

Patient and public involvement
Edinburgh Imaging staff were actively involved in the design of the prospective study. Pre-existing patient temperature data helped to inform the selection of most appropriate time points for scanning in healthy participants. TBI patients were not actively involved in study design since their data was anonymised prior to extraction and analysis. Healthy participants were invited to outline reasons for non-willingness to participate; where appropriate, this information was used to refine the recruitment approach. The findings of this study will be disseminated via Open Access publication. Conference presentations and public engagement activities will be used to explain the purpose of the research and its potential future impact. End of Study Information Sheets will be emailed to individual healthy participants who gave consent to receive these, and a lay summary of the study results will be posted online after publication.

Spatiotemporal measurements of healthy brain temperature
Of the 77 volunteers screened for eligibility, we recruited 20 males and 20 females (aged 20-40 years) between July and September 2019 (Fig.1A). Participants represented 15 nationalities across five continents, and the last participant was scanned on October 8 th 2019. One male attended only for morning scanning and another male volunteer missed afternoon scanning; available data from both of these participants was included in the analysis. Of the females scanned, 11 had natural menstrual cycles, eight were taking a combined contraceptive pill and seven of these took an 'active' pill on the day of scanning. The female subject on a 'pill break' reported day one of menstruation at their afternoon scan; their TBr data was included in the luteal group. Of the females with natural cycles, six were confirmed luteal (urine test), two were in menstruation, and three were in non-menstrual follicular phase at scanning. Five females thus formed a non-luteal group. One female wore a . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021.  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint temperature, total light exposure, and activity by sex during actigraphy week (left). Females n=20, males n=20. Associated rose plots with circular means (acrophases±SD) displayed (right). For each data type, radial uniformity was rejected for both groups (Rayleigh uniformity test p<0.0001) and there were no significant differences in circular mean between them (Watson's two-sample test for homogeneity, p>0.1). (D) Linear mixed modelling results for oral temperature. Solid blue lines represent model fits, shaded areas and double-ended error bars represent 95% confidence intervals, dark grey circles display residuals (single temperature data points), and smoothed dashed yellow lines represent partial residuals. The x-axis for Time summarizes the continuous variable of time distance since the participant's MSFsc (proportion of a linearized unit circle, where 0=MSFsc and 1=24 hours). Note time-dependent trend but lack of significant diurnal variation in oral temperature, likely reflecting inherent practical challenges of obtaining accurate oral temperature readings in human subjects. 34 All subjects exhibited diurnal variation in wrist skin temperature, which was anti-phasic with their rhythm in activity and light exposure, in the week preceding their scans (Fig.1B-C and fig.S1-S2).
BMI was marginally higher in males (P=0.014; Table 2). Oral temperature was 0.29°C higher in luteal females relative to males (95% confidence interval 0.03 to 0.58, P=0.029), and 0.04°C higher for a unit increase in BMI (0.005 to 0.083, P=0.024; Fig.1D). There were no differences in oral temperature by age or time of day however, despite daily changes in environmental temperature ( Fig.1D, fig.S3). Brain locations for MRS data sampling are shown in Fig.2A. MRS data from one female were excluded due to a HRF; 24 TBr data points from a total of 9434 (0.25%) were excluded because they did not meet quality control criteria for MRS spectral fitting (Supplementary Methods, fig.S4, Table S3). The data points that failed quality control derived from 15 of the 40 subjects scanned. Together, these data confirmed that our cohort was representative of healthy adult men and women with respect to basic physiological parameters, chronotype distribution, and sleep patterns. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. difference between sleep onset on free and work days when average sleep duration was longer on free than work days; if average sleep duration was longer on work days than free days, SJLsc was calculated as the absolute difference between sleep offset on free and work days). Note that this parameter was calculated only for participants that reported at least one of each 'day type' (free or scheduled) during data collection. CFI = circadian function index; this parameter ranged from 0.43-0.73 in an age-matched group of healthy volunteers. 35

Regional variation in brain temperature by sex and age
Global TBr was higher than oral temperature (38.5 SD 0.4°C versus 36.0 SD 0.5°C), and was 0.36°C higher in luteal females relative to follicular females and males (95% confidence interval 0.17 to 0.55, P=0.0006 and 0.23 to 0.49, P<0.0001, respectively). This sex difference appeared to be driven by menstrual cycle phase ( fig.S5). Despite age-selective recruitment, we captured an age-dependent increase in TBr, most notably in deep brain regions (thalamus and hypothalamus; 0.6°C over 20 years; 0.11 to 1.07; P=0.0002). Sex, age, and spatial effects on TBr are summarized in Fig.2B and . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  TBr substantially exceeds oral temperature and varies by sex, age, menstrual cycle, and brain region. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

Diurnal variation in brain temperature
Absolute TBr is ultimately determined by a balance between the rate of heat generated by the brain, and its rate of heat loss, mediated principally by CBF. 36,37 Since blood arrives to the brain from the body at a lower temperature, this temperature gradient should enable effective brain heat removal, as long as cerebral perfusion is maintained. 38 It follows that TBr must be partially determined by TBo.
Since TBo and CBF both show clear diurnal regulation in humans, with lower temperature and higher CBF at night, 20,21 we reasoned that human TBr should drop in the evening. Our linear mixed model Diurnal temperature variation was significantly greater in deep brain regions than in the cerebrum or the body (oral temperature; Fig.3C and fig.S6B-C), and for all brain regions, TBr was lowest in the late evening. Robust, approximately sinusoidal, daily TBo rhythms are a very well-characterised aspect of human physiology, and similar temperature rhythms have been extensively documented in other diurnal mammals in the brain and body. 39,40 Since TBr is expected to depend (at least in part) on TBo, we used the simplest and most appropriate mathematical model (cosinor analysis) to predict diurnal human TBr in a continuous fashion. We interpolated a sinusoidal time series for TBr in six brain regions of interest (Fig.3D). The predicted average minimum (anticipated around MSFsc, ~3am) was 38.4°C in luteal females and 38.0°C in males. Importantly, the diurnal TBr range across individuals was ~37.0-40.3°C in healthy cortical white matter-the location measured in patients with moderate-to-severe brain injury. In summary, these data show that normal human TBr varies substantially over the day, in a sex-and brain region-dependent fashion. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint healthy participant (n=39). Temperature varied more by time of day in the hypothalamus than orally (repeated measures one-way ANOVA with Sidak's multiple comparisons test ****P<0.0001; see fig.S6B for other brain regions). (D) 24-hour temperature rhythms of the healthy brain, double plotted. Interpolated average TBr rhythms in healthy luteal females (n=14) and males (n=20) based on sinusoidal fit using temperatures measured at three time points. Note higher temperatures in all regions in luteal females relative to males and marked variation in deep brain temperatures in males.

HEATWAVE-a 4D map of human brain temperature
Combining our spatial and temporal observations, we built HEATWAVE-a 4-dimensional map to model human TBr at hourly resolution (Movies S1 and S2). HEATWAVE can be dynamically explored at (https://www2.mrc-lmb.cam.ac.uk/groups/oneill/research/heatwave/). These comparisons highlight the relatively hot deep brain regions and their greater diurnal variation in males than females. The HEATWAVE videos complement the voxel maps in Fig.3A, which represent a reference resource for interpreting human TBr at each of the time points tested. Since each data point in each map is an average of data from multiple individuals, it incorporates the range of ages, BMIs, and chronotypes expected for each sex in the demographic tested. Our data collection points also cater for the times (morning and afternoon) when most patients would present for MR-based neuroimaging in the non-acute setting. In addition to modelling diurnal human TBr in a continuous fashion, HEATWAVE thus provides the first comprehensive spatially-resolved description of normal human TBr at three clinically-relevant time points; a rich reference dataset for future studies in different age groups and patient cohorts. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint

Daily brain temperature rhythms predict patient survival
The characterisation of physiological TBr variation in humans allows its dysregulation by brain injury to be understood in context for the first time. Of 134 eligible patient records screened, 114 had at least 24h of temperature data recorded (criteria level A). Of these, 110 patients had sufficient temperature data (≥36h) for diurnal rhythm analysis (criteria levels A and B; for eight of these patients, sufficient data was available for TBo only). Outcome in intensive care was available for 113/114 patients (criteria levels A and C), and a complete set of injury severity scores (PLR, GCS, and GCSM) was available for 109/114 patients (criteria levels A, C, D). A total of 107 patients met all criteria levels, and 100 patients had sufficient data to test for an association between diurnal TBr rhythmicity and outcome (mortality). Summary data are shown in Table 3. As in our healthy cohort, mean TBr (38.5 SD 0.8°C) was significantly higher than mean TBo (37.5 SD 0.5°C; P<0.0001, Fig.4A), but the range was much wider (32.6 to 42.3°C). TBr was not affected by the site of intracranial probe placement relative to focal injury ( fig.S7). We found an approximately daily temperature rhythm in 25/110 patients, of which 23 had a daily TBr rhythm (Fig.4B, fig.S8).
However, across the cohort, the timings of temperature maxima and minima were poorly aligned with the external day-night cycle. This desynchronization of internal timing from the external solar cycle is a hallmark of circadian rhythms when external timing cues are diminished, 41 and lies in stark contrast to rectal temperature data from healthy individuals maintained in the presence of daily light/dark and feed/fast cues (Fig.4C). is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint data from publicly-available database (right). 42 (D) Linear regression of patient temperature ranges with age; reduction in temperature amplitude significant for brain (slope of -0.016 significantly different from zero; 95% confidence interval -0.029 to 0.003). Shaded areas represent 95% confidence intervals for lines of best fit. Max, maximum; Min, minimum; NS, not significant.
As for healthy adults, there was a relationship between TBr and age; TBr amplitude was reduced in older patients (P=0.018), dominated by an upward trend in minimum temperature (Fig.4D, fig.S9).
Twenty-five patients died in intensive care. Applying a GLMM (Fig.5), we found that lack of a daily TBr rhythm, or an age increase of 10 years, increased the odds of death in intensive care 12-fold and 11-fold, respectively (OR for death with rhythm 0.09; 95% confidence interval 0.01 to 0.84, P=0.035 and OR for death with ageing by 1 year 1.10; 1.05 to 1.16, P=0.0002). These relationships could not be explained by a general elevation in TBr, since mean TBr was positively associated with survival (OR for death 0.45 for 1°C increase, 0.21 to 0.96, P=0.040). The presence of a diurnal TBr rhythm did not correlate with either age or mean TBr (Appendix 5). Together, these data show that daily temperature variation is frequently disrupted or absent in TBI patients and that TBr variation is of greater prognostic use than absolute TBr. Older TBI patients lacking a daily TBr rhythm are at greatest risk of death in intensive care, and presence of a daily TBr rhythm appears to be the strongest single predictor of survival after TBI. 43 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint in intensive care transformed from the data in (A-C); the results for these three predictors are significant since the 95% confidence intervals (double-ended error bars) do not include 1. Note also that confidence intervals become numerically asymmetric once transformed from log odds to regular odds. Only factors that demonstrated a statistically significant relationship with mortality are shown.
Note logarithmic scale on x-axis and large effect size for presence of a daily rhythm in TBr in (D).

See Methods and Supplementary
Methods for further details on the GLMM, and Appendix 5 for all numerical outputs and related code.

Principal findings
We have established a 4-dimensional map of human TBr, and shown how this parameter varies with time of day, brain region, age, and sex in healthy adults. These data provide clinicians with an urgently-needed and readily-accessible reference resource for evidence-based interpretation of TBr data in patients. Furthermore, we have found a relationship between the presence of a daily TBr rhythm and survival of TBI patients. Our findings demonstrate the high prognostic value of timeresolved TBr measurements in neurocritical care, thus empowering a temperature-based prediction of mortality. 22 Overall, this work reveals marked heterogeneity and dynamism of human TBr that must influence neural cell activity, and represents an important correlate of brain health.

Strengths and limitations
Although time-based human neuroimaging studies are sparse, some morning/afternoon comparisons are consistent with diurnal regulation of brain morphometry, 44,45 as well as diurnal variation in neural activity and metabolism. [46][47][48] However, prior studies were underpowered without consideration of chronotype and a late evening time point, which provides greater insight into healthy brain . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint physiology by incorporating the period approaching habitual sleep time (Fig.3). It was neither practical nor clinically-relevant to deprive our participants of all external timing cues (to derive circadian TBr variation), but this diurnal TBr variation is almost identical to direct measurements obtained in healthy non-human primates under stringent conditions. 18 Though unlikely, it is conceivable that gross regional differences in the activity of cellular water impact upon the apparent spatial TBr variation we observe within an individual at a given time point. However, we are unaware of any supporting evidence for this, nor can it be attributed to simple grey versus white matter distribution. 49 Moreover, such differences cannot influence the TBr variation we have found in relation to time of day, sex, age, or menstrual cycle stage. This is illustrated well when we limit our model to a subset of deep brain regions of more homogenous tissue structure, where TBr variation persists with respect to all of the aforementioned fixed effects (Fig.2B, 3B, fig.S6). Crucially, our robust statistical approach caters for multiple physiologically-relevant confounders within and between individuals that would have prevented the detection of significant TBr variation in previous studies. 49,50 Alongside the patient data (Fig.5), and multiple parallel methods of temperature measurement in healthy subjects by us and others, 22,51-53 our results offer compelling evidence of a daily temperature rhythm throughout the normal human brain (Supplementary Text 'Temperature rhythms and sleep').
The within-brain temperature gradient is remarkable (Fig.2B). As an 'open' thermodynamic system performing no mechanical work, aerobic metabolism of the brain releases heat at ~0.66 J/min/g of tissue which is primarily removed by CBF. 38,[54][55][56] It is therefore highly likely that regional variation in neurovascular anatomy plays the chief role in creating spatial TBr gradients (Supplementary Text 'Temperature gradients'). Although we cannot completely exclude a contribution from regional differences in water content, 49,57 these are unlikely to explain the temperature difference between the thalamus and hypothalamus (both grey matter structures devoid of cerebrospinal fluid). We suggest . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint that the lower temperature of the hypothalamus might reflect its closer proximity to major vascular networks such as the Circle of Willis. In principle, technical limitations (Supplementary Text 'Technical limitations of brain thermometry') could potentially exaggerate MRS-derived temperature differences at the extreme edges of regions of interest (cerebral layer Sup4). The spatial distribution we have found is however very similar to non-human primates, excepting a larger gradient magnitude that is entirely consistent with the difference in brain volume between humans and rhesus monkeys. 19 Unlike previous studies, 49,50 we make no baseline assumption that temperature should be homogeneous across brain regions, nor between different tissue types within the brain. Importantly, we did not apply a post-acquisition correction to our MRS data to equalize temperatures between grey and white matter, 49,50 since this would perpetuate the above assumption, and overlooks the clear tissue temperature differences observed in non-human primates and normothermic human patient brains. 14, 16,19 Indeed, higher temperatures in white matter-rich areas concur with predictions based on modelling perfusion, blood volume fraction, and heat generation in different brain tissues. 13,58-60

Possible mechanisms and implications
An increase in mean TBr (Fig.2B) and a trend upwards in minimum TBr ( fig.S9) with age suggests that overnight brain cooling becomes less efficient in older people, leading to a damped TBr rhythm.
This age-dependent reduction in TBr amplitude is consistent with studies of TBo and may contribute to the disrupted sleep patterns and 'sundowning' symptoms of dementia patients. 24,61-63 Cerebral blood supply is considered so efficient that heat removal is achieved without the need for other mechanisms under most circumstances, 36,37,64 which seems intuitive for the young, healthy brain.
However, the vast literature linking neurodegeneration to cerebrovascular compromise indicates that our key brain cooling mechanism progressively deteriorates with age (Supplementary Text 'Internal rhythms and health'). 65,66 Neuronal activity is highly sensitive to temperature change, with a Q10 of ~2.3, although this is generally considered to be most problematic in the acute setting. 38 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint study of 1130 epilepsy patients, 80-92% showed a 24-hour cycle of seizure rates, with events most common at ~8am, when TBr should increase most steeply (Fig.3D). 68 Given that cooling can terminate epileptic discharges, 69 diurnal changes in TBr may well contribute to diurnal variation in the incidence of seizures and cluster headache. 70,71 TBo increases, and its overnight drop is blunted, in luteal versus follicular phase women. 52,72 This menstrual variation predicts that TBo is ~0.4°C higher in the early luteal phase. 72 3D). 76 It is widely accepted that BMI positively correlates with TBo as found here (Fig.1). 22,77 Since BMI was slightly higher in males relative to females in our healthy cohort, a difference in BMI cannot explain the higher TBr observed in luteal females and, notably, there was no relationship between BMI and TBr overall ( Fig.2B; Appendix 5). This supports our conclusion that TBr cannot be solely dependent on, nor predicted from, TBo since brain heat removal also occurs through routes that are unaffected by adipose deposition. 36,37 Essential to clinical diagnostics is the comparison of patient data with reference ranges from healthy individuals; MRS-thermometry now makes this possible for TBr. We have validated our core MRS findings using multiple complementary methods of temperature measurement. This is most pressing for TBr, where such methods are effectively mutually exclusive in healthy individuals and . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint neurocritical care patients. TTM is the mainstay of neuroprotection subsequent to out-of-hospital cardiac arrest. 78 Here, the objective is to reduce TBr, which is rarely measured directly in trials that test the therapeutic value of TTM in the context of brain injury. Cooling adults at the 'wrong' biological time or fixing patient temperatures at a constant target value may further compromise thermoregulation by abolishing physiologically-important, health-critical temperature variation. The highest temperature we observed in any healthy individual was 40.9°C in the thalamus of a luteal female in the afternoon; whilst the perception exists that a TBr of this value would cause brain damage, there is no direct evidence for this, and similar deep brain temperatures are observed physiologically in other mammalian species. 79 Furthermore, the TBr range in our volunteers raises doubt over whether TBr was abnormally high in some patient reports. 15 Current temperature management guidelines do not consider physiological differences by sex or time of day, 80 and whether adults should be cooled at all in neurocritical care remains controversial. A clear understanding of how and why TBr varies in health and disease is thus imperative. Here we report a healthy cortical white matter maximum TBr of 40.3°C, but we caution strongly against overinterpreting single TBr values or transitory trends. Rather, we recognize the need for technological solutions that allow individualized target temperature ranges to be determined, facilitating decision making that incorporates chronotype, age, sex, menstrual cycle, and time of day.

Unanswered questions and future research
Prospective controlled trials are needed to confirm the predictive power of TBr rhythmicity in relation to patient outcome, as well as the clinical utility of TTM protocols in brain-injured patients. There may be high clinical value in exploiting TBr variation to detect or monitor focal pathologic processes such as neoplasia, trauma, vascular insults, and epileptogenesis, but also more distributed inflammatory, metabolic, and neuropsychiatric diseases. 54,[81][82][83][84][85][86][87] In particular, future work should be directed to address whether abnormal TBr rhythmicity may serve as an early biomarker of . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint  88 We have shown that more sophisticated analyses can better exploit temperature as a clinical tool. Wearable devices now permit easy and convenient recording of daily rhythms in many physiological parameters. Algorithm-based temperature profiling will help accomplish the goals of precision medicine, not just for individuals, 89 but at scale. For example, in an infectious disease outbreak, real-time screening for fever development could rapidly identify high-risk individuals by deviation from their own temperature rhythm, rather than a population 'mean' or by random testing.
Personalized, digital, round-the-clock temperature monitoring would thus advance remote health tracking and evidence-based enforcement of global health policy in the context of emerging disease.
Whilst providing excellent spatial resolution, MRS brain thermometry is clearly impractical for routine use in most clinical settings. Since core TBo is not a faithful proxy for TBr, 90  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint Contributors: NMR conceived the idea for the work, designed and orchestrated the study, conducted recruitment and data collection for the prospective study, analyzed and interpreted all data, prepared all tables and figures, and wrote and revised the manuscript. JSON made conceptual contributions to the prospective study and contributions to the overall study design, interpreted data, contributed to figure preparation, and revised the manuscript. IM and MJT collected, analyzed, and interpreted spectroscopy data, contributed to figure preparation, and revised the manuscript. FMC contributed to study design, devised the statistical analysis plan for the prospective study, and revised the manuscript. AE and MC contributed to study design, CENTER-TBI High Resolution ICU data curation and extraction, statistical analysis, and revised the manuscript. GM analyzed, interpreted, and reported on prospective structural MRI data, and revised the manuscript. JR is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint Research and Development Office. This research was independent from funders. Funders had no role in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the article for publication.
Competing interests: all authors have completed the Unified Competing Interest form (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years, no other relationships or activities that could appear to have influenced the submitted work. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint CRF. Ethical approval was obtained for each recruiting site. The list of sites, Ethical Committees, approval numbers and approval dates can be found on the website: https://www.centertbi.eu/project/ethical-approval.
Data sharing: Individual patient data contained within the CENTER-TBI database are not publicly available but permissions for access can be requested at https://www.center-tbi.eu/data. We are committed to sharing all other anonymised individual participant and patient data that would support the clinical community. All shareable items are available immediately upon publication and indefinitely, or ending 5 years following article publication, by reasonable request from the Lead Author at ninar@mrc-lmb.cam.ac.uk. Shareable items will be available to anyone who wishes to access them and for any purpose. Code for statistical modelling is provided in Supplementary Appendix 5.
Transparency: The lead author (NMR) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.

Copyright:
The Corresponding Author has the right to grant on behalf of all authors and does grant on behalf of all authors, an exclusive licence (or non-exclusive for government employees) on a worldwide basis to the BMJ Publishing Group Ltd to permit this article (if accepted) to be published in BMJ editions and any other BMJPGL products and sublicences such use and exploit all subsidiary rights, as set out in our licence.
. CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted January 27, 2021. ; https://doi.org/10.1101/2021.01.23.21250327 doi: medRxiv preprint