Brown adipose tissue thermogenesis among a small sample of reindeer herders from sub-Arctic Finland
Journal of Physiological Anthropology volume 41, Article number: 17 (2022)
Interest in human physiological responses to cold stress have seen a resurgence in recent years with a focus on brown adipose tissue (BAT), a mitochondria dense fat specialized for heat production. However, a majority of the work examining BAT has been conducted among temperate climate populations.
To expand our understanding of BAT thermogenesis in a cold climate population, we measured, using indirect calorimetry and thermal imaging, metabolic rate and body surface temperatures of BAT-positive and BAT-negative regions at room temperature, and mild cold exposure of resting participants from a small sample of reindeer herders (N = 22, 6 females) from sub-Arctic Finland.
We found that most herders experienced a significant mean 8.7% increase in metabolic rates, preferentially metabolized fatty acids, and maintained relatively warmer body surface temperatures at the supraclavicular region (known BAT location) compared to the sternum, which has no associated BAT. These results indicate that the herders in this sample exhibit active BAT thermogenesis in response to mild cold exposure.
This study adds to the rapidly growing body of work looking at the physiological and thermoregulatory significance of BAT and the important role it may play among cold stressed populations.
Much of biological anthropology in the past 50 years has focused on human health and disease, with less focus on how physical environments influence human biology. However, with the impacts of climate change growing ever more dire, scholars have begun to better integrate human biology, health, disease, and the physical environment into their research [1,2,3,4]. Though in many places climate change will lead to warming, other localities will experience an increase in extreme cold weather events, while drastically different day- and night-time temperatures will still exist elsewhere [5, 6]. As such, cold will always be a stress with which humans will need to cope . Though humans are a hot climate originating species , a variety of adaptations have evolved in order to survive and thrive in the relatively recent stressor—extreme cold . One such adaptation is brown adipose tissue (BAT), the main driver of non-shivering thermogenesis. Here, we present the results of BAT activity measurements conducted among a small sample of reindeer herders living in sub-Arctic Finland.
We are physiologically limited in our ability to maintain our core body temperature in the face of cold stress; thus, we also employ behavioral and cultural mechanisms that enable us to survive and thrive in extreme cold environments. Much of the work revolving around these mechanisms focuses on vascular and metabolic responses to cold stress. Vascular responses, which work to reduce heat loss, include peripheral vasoconstriction and countercurrent heat exchange [7, 9, 10]. Metabolic mechanisms, which increase body heat production, include shivering thermogenesis , increased resting metabolic rate (RMR, kcal/day) [11,12,13,14,15,16,17,18], and non-shivering thermogenesis via brown adipose tissue [19,20,21,22,23]. Shivering, produced by small muscle contractions in response to cold stress, does not substantially increase body heat production nor persist for a long period of time . RMR, the metabolic cost of maintenance, among cold climate populations is typically higher than expected based on predictive equations and higher relative to temperate climate populations [11, 12, 14, 16]. This is thought to be driven by high levels of thyroid hormone [11, 13, 14, 25], but there is also a high degree of interindividual RMR variation . Among the herders in the present study, after controlling for body and fat-free mass, females had significantly higher RMRs than predictive equations and significantly higher RMRs than males. Males, however, showed no distinct pattern relative to predictive equations , only partially matching the pattern seen among other cold climate populations.
Brown adipose tissue is a mitochondria dense tissue that interrupts the electron transport chain via mitochondrial uncoupling protein-1, resulting exclusively in heat rather than adenosine triphosphate production [19, 21, 26, 27]. BAT presence has been well known among hibernating and newborn mammals and human infants, but is, as we now know, also present in healthy adults in cold [19, 28] and temperate climate populations [20, 23, 26, 27, 29, 30]. BAT deposits (~ 100 g) are most prominent in supraclavicular and paracervical regions as well as along major deep blood vessels and around numerous internal organs . These deposits consist of brown and beige adipose tissue, with beige adipose tissue having a reduced thermogenic capacity and is also found among white adipose tissue deposits . Beige adipose tissue is not discussed nor analyzed here.
Those with more BAT thermogenesis typically experience a greater increase in metabolic rate, warmer body surface temperatures over BAT positive regions, and higher glucose utilization in response to mild cold stress (~ 10-15 °C ) [33,34,35]. This thermogenesis is induced by the sympathetic nervous system and thyroid hormones. Among individuals with active BAT, there is little decline or even an increase in body surface temperatures at the supraclavicular (known BAT location) region during mild cold exposure, whereas individuals without active BAT show significant surface temperature declines in this region [19, 20, 36]. While surface temperature patterns are typically consistent, there is currently no consensus on BAT substrate utilization and some have suggested a higher glucose utilization during BAT activity [19, 36, 37], while others found greater fatty acid utilization. Chondronikola et al.  found that the increase in metabolic rate associated with BAT was 70% fueled by free fatty acids and 30% by glucose, while others found an increase in low-density lipoprotein and high-density lipoprotein (HDL) cholesterol levels , indicating a high degree of variation in BAT substrate utilization [39, 40]. The amount of BAT is inversely associated with body-mass index, body adiposity, and age [29, 41].
Our understanding of inter- and intra-populational variation, morphological correlates, metabolic responses, genetic indicators, and potential therapeutic benefits of BAT are still limited. Among the Yakut of Siberia, males experienced a decrease in metabolic rate during mild cooling, though no decline in surface temperatures, while females did not . Potential explanations for this decline in metabolic rate may be habituation, increased vasoconstriction, or the Q10 effect, which is slowing of enzymatic activity in muscle due to cold exposure though the true cause is still unknown [19, 42, 43]. In the present study, we assessed BAT activity by measuring metabolic rate, substrate utilization, and surface temperatures during acute, mild cold exposure among a small sample of reindeer herders from sub-Arctic Finland. Though our sample size is small, we expected the reindeer herders in this study to have active BAT in response to mild cold exposure, as has been documented among another cold climate population. This work expands our current understanding of populational variation in BAT activity as a response to cold stress.
Methods and participant population
The results of this work are one part of a larger project, The metabolic cost of living among reindeer herders of sub-Arctic Finland, that took place during October of 2018 during the annual herd roundup and in January 2019. The BAT measurements discussed here took place in January of 2019. Reindeer herders (N = 22, 6 females, and 16 males) from herding districts within 180 km of Rovaniemi, Finland (66.5° N, 25.7° E) participated in this study. The herding districts included Palojärvi, Niemelä, Narkaus, Pyhä-Kallio, Oivanki, Isosydänmaa, and Poikajärvi; for a map, please refer to Ocobock et al. . The skewed sex ratio in this study reflects the actual sex ratio among reindeer herders in Finland, among whom 30% are female . Herders were recruited through existing contacts as well as advertising by the Reindeer Herders Association, and interested herders were contacted by the research team. This resulted in a fairly large age range (20–64 years) among the participating herders. For half of the participating herders, reindeer herding was their primary occupation. For the other half, herders supplemented their income with tourism, meat processing, land measurement, research, and construction. In Rovaniemi, the mean temperature in January of 2019, when BAT measurements were conducted, was − 16.4 °C (2.5 °F) .
The Arctic Centre of the University of Lapland in Rovaniemi served as the operation base for the larger study and was the location for much of the data collection. However, when herders were unable to come to the Arctic Centre, measurements were conducted at herders’ homes, cabins within the herding district, or hotel rooms near the herding district of interest. At the time of participation, all herders were active members of their herding district and self-reported being healthy at the time of measurement. All participating female herders were not pregnant or lactating; however, information on menstrual cycle phase was not documented during this study.
This study was conducted with Institutional Review Board approval from the University at Albany (17-E-165) and with the approval of the Ostrobothnian Health Care District from the University of Oulu, Finland (EETTMK: 4/2018). Participants were provided with an oral introduction and written information sheet about the study, and informed written consent (written in Finnish) was obtained from all participants .
Height and weight were measured following standard procedure using a portable stadiometer (Seca Corporation, Hanover, MD) and an electronic scale (AccuWeight, New York, NY), respectively. Height was recorded to the nearest 1 mm, and weight was recorded to the nearest 0.1 kg. Body composition was measured using bioelectrical impedance (RJL systems, Clinton, MI). Participants were in a supine position, with arms and legs abducted from the body, and small electrodes were placed on the right side of the body at the dorsum of the wrist, middle finger, ankle, and middle toe. Reactance and resistance were recorded and the NHANES-III equations were used to determine fat-free mass (FFM) and body fat percentage . All participants were asked to refrain from consuming alcohol in the 24 h before measurement, and they arrived to the data collection period 12 h fasted. Total cholesterol, glucose, and HDL cholesterol levels were measured immediately after the anthropometric measurements and before the metabolic measurements. Finger prick whole blood samples were analyzed using the CardioChek PA (Polymer Technology Systems, Indianapolis, IN) point of care devise with glucose-cholesterol test strips.
Brown adipose tissue and respiratory quotient
BAT thermogenesis was non-invasively inferred through the simultaneous measuring of metabolic rate and surface temperatures of herders, at an anticipated BAT and non-BAT region, during mild cold exposure following Levy [27, 48]. BAT thermogenesis is indicated by an increase in metabolic rate and higher surface temperatures at a BAT region, the supraclavicular region (TSC) compared to a non-BAT region at least among humans, the sternum (TST). In order to assess the impact BAT thermogenesis has on these two parameters, metabolic rate and surface temperatures were first measured in resting subjects at room temperature (20–27 °C) and then measured again at mild cold exposure (15–18 °C), which lasted for 30 min. To expose participants to mild cold, they wore a 3-piece liquid cooling garment (BCS4, Allen Vanguard; Ontario, Canada) fitted with internal tubing through which cold water flowed. For the room temperature condition, no water flowed through the suit.
Prior to the start of data collection, participants put on the 3-piece cooling suits. Participants were then asked to rest in the supine position at room temperature for 30 min prior to the start of measurements allowing them to come to a more complete rest, familiarize themselves with the environment, and minimize possible anxiety associated with testing. Room temperature metabolic rate (MRRT, kcal/day) and respiratory quotient (RQRT, VCO2/VO2) were measured using a K5 portable indirect calorimetry unit (Cosmed, Chicago, IL) for 30 min with the last 10 min of data used for analysis. This system utilizes a facemask with bi-lateral, unidirectional valves (allowing inspiration but not expiration), which is preferable to a metabolic hood as the hood interferes with thermal imaging of the key BAT locations. The K5 unit measures breath-by-breath O2 consumption and CO2 production. MRRT and RQRT values were calculated and recorded using the Cosmed Omnia software. In general, a RQ value closer to 1.0 indicates carbohydrates utilization, while a RQ closer to 0.70 signifies that fat is the primary metabolic substrate.
During the 30-min MRRT, thermal images of a BAT region and a non-BAT region were taken every 5 min using an E75 thermal imaging camera (FLIR, Wilsonville, OR) held 30 cm from the area of interest with an emissivity of 0.98. The supraclavicular region (a known BAT location) was defined medially by the sternocleidomastoid muscle, laterally by greater tubercle of the humerus, and inferiorly by the clavicle. The sternum served as a control as BAT is not found here. The images were uploaded to the FLIR Tools+ software, and maximum temperature of the supraclavicular and sternal regions were recorded.
Once the 30-min room temperature exposure was complete, cold water (7–10 °C) was pumped into the cooling suit which takes roughly 2–3 min to completely fill resulting in exposure temperatures of approximately 15–18 °C. The cold exposure similarly lasted 30-min, and cold metabolic rate (MRC), RQ (RQC), and surface temperatures were all measured exactly the same as described above for the room temperature exposure. Participants were monitored the entire time for potential shivering, and participants were instructed to alert the researcher if shivering started. Shivering was only observed in two participants, the time of shivering was noted so that data could be removed from analysis, and warm water was added to the cooling suit to keep the participants at mild cold exposure. These shivering episodes lasted no longer than 60 s.
We conducted all statistical analyses using SPSS version 26 (IBM, Armonk, NY), and we produced graphical figures using R Studio (version 2021.09.0). Data were normally distributed (p > 0.06 in all cases). One-way ANOVAs were used to compare each age, height, weight, sum of skin folds, percent body fat percentage, FFM, and biomarkers. We used paired Student’s t-tests to determine if TSC, TST, metabolic rates, and RQ were significantly different between room temperature and mild cold exposure within each individual. We calculated the difference (∆) and percent change for metabolic rate, RQ, TSC, TST, and the difference between cold exposure surface temperatures for the supraclavicular region (SCC) and sternal region (STC). We then performed multiple linear regression analysis to determine if body weight, fat free mass, body fat percentage, sex, age, and biomarkers (glucose, total cholesterol, and HDL cholesterol) were predictors for MRC, ∆MR, ∆RQ, ∆TSC, ∆TST, and SCC–STC. Tables display individual participant results, and pooled results are presented as the mean ± the standard deviation. Multivariate regression tables provide the adjusted R2, F statistics, p values, and β coefficient. We considered results significant at the p ≤ 0.05 level.
Descriptive statistics and a summary of the means ± the standard deviation for anthropometric, demographic, and biomarker data can be found in Table 1. For the overall sample, age ranged from 20 to 64 years, height ranged from 154.8 to 194.6 cm, body mass ranged from 46.7 to 129.6 kg, FFM ranged from 34.2 to 86.0 kg, and body fat percentage ranged from 20.4 to 40.5%. Females were significantly younger (20–37 years) than males 36–64 years (F = 26.557, p < 0.01). Among females, age ranged from 20 to 37 years, height ranged from 154.8 to 163.0 cm, body mass ranged from 46.7 to 78.9 kg, FFM ranged from 32.2 to 46.9 kg, and body fat percentage ranged from 26.8 to 40.5%. Among males, age ranged from 36 to 64 years, height ranged from 168.0 to 194.6 cm, body mass ranged from 62.9 to 129.6 kg, FFM ranged from 48.2 to 86.0 kg, and body fat percentage ranged from 20.4 to 36.8%. Females were significantly shorter (F = 33.168, p < 0.01) and had significantly lower body mass (F = 11.282, p < 0.01) and FFM (F = 31.043, p < 0.01) relative to males. Females had significantly more body fat than males (F = 16.370, p < 0.01). There was no significant difference between females and males for any of the biomarker measures (p > 0.2 in all cases).
Metabolic rate, ∆RQ, and surface temperatures at room temperature and cold exposure
A summary of the means ± the standard deviation for metabolic rate and surface temperatures at room temperature and during cold exposure can be found in Table 2. For a detailed description of room temperature resting metabolic rates among the reindeer herders, please see Ocobock et al. . Overall, MRRT ranged from 1004 to 2614 kcal/day. Among females, MRRT ranged from 1599 to 2324 kcal/day. Among males, MRRT ranged from 1004 to 2614 kcal/day. Females had a significantly higher MRRT than males, the reasons for which has been addressed previously . Overall MRC ranged from 1123 to 2787 kcal/day with females ranging from 1599 to 2324 kcal/day and males ranging from 1123 to 2787 kcal/day. MRC was significantly higher than MRRT (t = − 3.629, p = 0.002) (Fig. 1).
Room temperature supraclavicular surface temperatures (SCRT) (Table 2) ranged from 30.0 to 33.2 °C for the overall sample and 30.5–31.8 °C among females and 30.0–33.2 °C among males. Room temperature sternum surface temperatures (STRT) ranged from 27.2 to 32.6 °C for the overall sample and 29.7–31.2 °C among females and 27.2–32.6 °C among males. Cold exposure supraclavicular surface temperatures (SCC) ranged from 28.7 to 31.2 °C for the overall sample and 28.7–30.5 °C among females and 29.5–31.2 °C among males. Cold exposure sternum surface temperatures (STC) ranged from 24.2 to 29.8 °C for the overall sample and 24.2–28.0 °C among females and 24.8–29.8 °C among males. There was no significant difference in surface temperatures between females and males for SCRT (F = 1.723, p = 0.204), STRT (F = 0.259, p = 0.616), or STC (F = 0.522, p = 0.478). However, males had significantly warmer STC surface temperatures than females (F = 9.840, p = 0.005). Supraclavicular surface temperatures were significantly higher than sternum surface temperatures for both the room temperature condition (t = 5.274, p < 0.001) and cold condition (t = 10.675, p < 0.001). Room temperature surface temperatures were significantly higher than cold exposure temperatures (supraclavicular: t = 10.179, p < 0.001; sternum: t = 12.845, p < 0.001). SCC–STC ranged from 0 to 5.5 °C for the entire sample, 0.7–5.0 °C for females, and 0–5.5 °C for males (Table 2, Fig. 2). There was no significant difference in SCC–STC between males and females (F = 0.071, p = 0.792).
The mean difference ± the standard deviation in metabolic rate (∆MR), RQ (∆RQ), supraclavicular surface temperatures (∆TSC), and sternum surface temperatures (∆TST) are listed in Table 3. Overall, ∆MR ranged from − 320 to 517 kcal/day with females ranging from 27 to 279 kcal/day and males ranging from − 320 to 517 kcal/day. For ∆RQ, the overall sample ranged from − 0.07 to 0.13 with females ranging from − 0.04 to 0.13 and males ranging from − 0.07 to 0.07. For ∆TSC, the overall sample ranged from − 3.10 to − 0.40 °C, females ranged from − 3.0 to − 0.80 °C, and males from − 3.1 to − 0.4 °C. For ∆TST, the overall sample raged from − 6.20 to − 1.50 °C, females ranged from − 6.2 to − 2.60 °C, and males from − 5.7 to − 1.5 °C. For the whole sample, there was a significantly greater decrease in sternum surface temperatures relative to supraclavicular surface temperatures (t = 7.671, p < 0.001). There were no significant differences between the sexes for ∆MR (F = 0.627, p = 0.438), ∆RQ (F = 2.262, p = 0.148), ∆TSC (F = 0.914, p = 0.350), or ∆TST (F = 0.1.665, p = 0.212).
The relationship between changes in metabolic rate, RQ, surface temperatures, and anthropometric variables
Table 4 displays a summary of the multiple linear regression results for MRC, ∆MR, ∆RQ, ∆TSC, ∆TST, and SCC–STC assessing potential associations with age, sex, height, body mass, body fat percentage, FFM, and biomarkers. No significant regression equation was found for MRC or ∆MR, ∆RQ, or ∆TST. For ∆TSC, the only significant predictor was ∆RQ. Once the non-significant predictors were removed, the final model statistics were R2 = 0.304, F = 8.725, p = 0.008, and β = − 0.038 (Fig. 3). For SCC–STC, fat free mass was a significant predictor. Once the non-significant predictors were removed, the final model statistics were R2 = 0.527, F = 12.711, p < 0.01, and β=1.516. None of the biomarkers were significant predictors for any regression analysis.
We found that in response to mild cold exposure, herders experienced a significant increase in metabolic rate; however, there was a high degree of interindividual variation with percent change in metabolic rate ranging from − 13.7% to + 35.0%, with four of the 16 males experiencing a metabolic rate decrease. In response to cold exposure, surface temperatures decreased at both the supraclavicular region and sternum; however, the temperature decrease at the sternum was significantly greater than the supraclavicular region. The relatively warmer supraclavicular surface temperatures are indicative of BAT thermogenesis at this region as well as potentially greater heat loss through conduction with the cooling suit at the sternum. From the RQ measurements, there appeared to be no change in substrate utilization, with herders exhibiting a metabolic preference for metabolizing fatty acids in both room temperature and cold conditions. Herders did tend to have a lower RQ in both conditions, which is not uncommon among people who consume a high amount of dietary fat as the herders do . There was no correlation of MRC, ∆MR, or ∆RQ with any anatomical or physiological variables. We found a negative correlation between ∆RQ and ∆TSC, suggesting that a decrease in RQ during cold exposure was associated with smaller surface temperature changes at the supraclavicular region, indicating that fatty acids were overwhelmingly utilized during cold exposure. We also found a positive correlation between SCC–STC and fat free mass, indicating that those with greater FFM experienced a blunted drop in supraclavicular temperature during cold exposure.
Similar work has been conducted among temperate climate populations; however, there is not a great deal of consensus. A study done among a population in Albany, NY, looked at seasonal variation in BAT activity and found a greater increase in metabolic rate in response to cold exposure during the summer relative to the winter, 12% and 11%, respectively . In both seasons, though, glucose was the preferred metabolic substrate. Among Japanese males, the opposite was found, with a greater metabolic rate in response to cold exposure exhibited during the winter . Nirengi et al.  demonstrated a comparatively blunted BAT response to cold exposure. The Japanese participants in this study had a mean SCC–STC of 2.2 °C, similar to that of the Yakut , whereas the herders in the present study had a much greater SCC–STC–, this is indicative of greater BAT activity. Similarly, among a population from the Netherlands, there was a greater non-shivering thermogenic response after cold acclimation, an analogue to winter, than prior to acclimatization resulting in a higher metabolic rate associated with BAT activity . We have summarized the results of the studies that have employed the thermal imaging methodology for assessing BAT in Supplemental Table 1.
Compared to other cold climate populations, there are a few similarities and numerous differences between the herders of this study and the Yakut participants from the only other BAT study among a circumpolar population [19, 53]. Among the few similarities are that males (four in the present study) more frequently experienced a decrease in metabolic rate during cold exposure, that there was no correlation between any demographic and anatomical variables with MRC, and that both populations exhibited a correlation between ∆RQ and ∆TSC. The decrease in metabolic rate is not an unusual phenomenon; however, we do not yet fully understand what might be driving this pattern among some individuals. The decrease could be due to vasoconstriction, the Q10 effect, and/or habituation [19, 42, 43], though more work is needed to elucidate the cause and pathway.
Conclusive anatomical correlates with BAT thermogenesis remain elusive with little to no correlation with fat free mass. Some studies demonstrate a negative correlation between body fat and BAT thermogenesis [54,55,56], some show a positive correlation , and many demonstrate no correlation including the present study [19, 58,59,60,61]. The lack of anatomical correlation may be surprising given BAT is thought to derive from muscle progenitor cells and beige adipose tissue appears to have a mix of white fat and muscle progenitors [62, 63].
However, work by Levy et al.  conducted among the Yakut suggests that environmental exposures to cold stress during key periods of childhood development may determine BAT thermogenic capacity later in life. Investigators explored if there was any association between measurements of adult BAT thermogenesis  and, using retrospective surveys and weather station data, cold exposure during gestation, infancy, early childhood, middle childhood, and adolescence. They found three significant associations with adult BAT thermogenesis: (1) a negative association with mean ambient temperature during early childhood and adolescence, (2) a positive association with the number of below − 40 °C days during early childhood, and (3) a positive association with participation in outdoor winter activities for ages five to seven and 11 to 13. This work among the Yakut suggests that BAT is developmentally plastic and that key periods of development likely shape BAT thermogenesis in adulthood. This is an intriguing suggestion and works well within established theory on environmental signals altering developmental trajectories [64,65,66].
There are a number of differences between these two cold climate populations. First, the reindeer herders (Table 2), especially the females, have higher resting metabolic rates relative to the Yakut sample—female Yakut MRRT = 1129 ± 31 kcal/day, MRC = 1094 ± 43 kcal/day; male Yakut MRRT = 1575 ± 48 kcal/day, MRC = 1528 ± 49 kcal/day. The Yakut sample also had overall warmer surface temperatures than our herder sample at the supraclavicular region (~ 36 °C vs. 30 °C) and sternum (~ 34 °C vs. 26.7 °C), for both females and males. Finally, the Yakut sample had much higher RQ values during room temperature and cold exposure, and there was a correlation between BAT thermogenesis and blood glucose levels, demonstrating a metabolic preference for carbohydrate utilization. While both populations had a correlation between ∆RQ and ∆TSC, herders demonstrated a negative correlation while the Yakut demonstrated a positive one. Two potential reasons for this difference could be a negative energy balance or the high-fat diet consumed by the herders in this study, which can lead to lower RQ values .
Despite these differences, both populations exhibited evidence of BAT activity, which could be a critical adaptation to their respective harsh, cold climates. Furthermore, the differences, such as the substrate utilization mentioned above, should not be entirely surprising as one would not and should not expect two different populations to arrive at an environmental adaptation through the same physiological routes. Differences in adaptive routes are seen elsewhere, such as the wide array of differing adaptations observed among the high altitude populations of Tibet, Peru, and Ethiopia [67, 68]. It is likely that cold climate adaptations are similarly variable and evolutionarily subject to differential selective pressures as well as genetic drift and gene flow .
Besides potential differences in evolutionary trajectories, there are other reasons why the reindeer herders in the present study may differ from the Yakut population. First, there are likely different developmental exposures to cold stress between these two populations. For example, the Yakut have their infants nap indoors, whereas in Finland there is a tradition for infants (2 weeks to 2 years old) to nap outside in temperatures ranging from − 27 °C to + 5 °C [70, 71]. This is a tradition we were able to witness first hand when socializing with one of the herding families. The mother prepared the infant by dressing him warmly, placing him in a stroller that had a blanket at the base, draping a blanket over the top of the stroller, and then taking him outside for a 2-to-3-h nap. Infants napping outdoors could have a significant impact on adult BAT thermogenesis. Furthermore, given the need for outdoor physical activity among reindeer herders throughout winter and the strong cross country skiing culture in Finland, even among very young children , it is likely that children in this population have a high level of cold exposure during key developmental periods. The potential for BAT developmental plasticity among the herders in the present study is an exciting avenue for future research on adaptive responses to cold stress across the life course.
Seasonal acclimatization is another possible reason for the differences seen in this population, since in the present study, BAT thermogenesis was assessed during the winter, while the study among the Yakut was performed during early fall. There is a known increase in BAT thermogenesis during the winter among people living in temperate climates [52, 58, 73]. Furthermore, there appears to be a seasonal difference in BAT substrate utilization, at least among a population in Albany, NY, in which carbohydrate utilization during BAT thermogenesis is higher in the winter than it is in the summer . This indicates that BAT thermogenesis is not only plastic during key developmental stages but that it is also highly responsive to environmental cues in adulthood. Gathering more seasonal data on BAT thermogenesis among the herders and other cold climate populations will help elucidate seasonal and well as populational variation.
Finally, other confounding factors that may lead to variation in BAT thermogenesis, particularly substrate utilization, are the levels of physical activity and dietary intake during the days leading up to the measurement. The herders in this sample had an exceptionally high total energy expenditure during the annual herd roundup in the fall . Though physical activity levels in the winter are likely not as high as they are in the fall, herders typically provide supplemental feed to their herds—a demanding chore . At this time, herders also prefer to consume energy-rich foods. The RQ associated with BAT thermogenesis in this study could have been affected if in the days before the BAT measurements herders were highly active and consumed fat-rich foods. This situation could have left the herders, who were at least 12 h fasted at the time of measurement, with relatively depleted stores of glycogen  and more readily available fatty acids for BAT to utilize. Future work should account for diet and activity levels in the days prior to BAT thermogenesis assessment.
There are a number of limitations to the present study. First, there was a small number of participants reducing the overall statistical power. Second, males in this study outnumbered females. Though the sex-bias reflects the current reindeer herder demographics , the small number of females makes it difficult to determine if there may be any sex-based differences in BAT thermogenesis. Third, this data is cross-sectional in nature as the data collection took place only in winter. As mentioned above, BAT thermogenesis is sensitive to seasonal temperatures. To get a broader view of BAT variation in this population, seasonal measures are needed. Fourth, to better understand BAT substrate utilization, it would be beneficial to have whole blood measures of glucose and cholesterol before and after cold exposure in addition to the RQ measurement. Fifth, these measures took place in several different locations for which room conditions (temperature and humidity) could not be easily controlled. Future work should utilize one location for data collection. Finally, this study did not collect data on childhood cold stress nor physical activity and dietary intake in the days leading up to the BAT measurement. These factors likely alter BAT thermogenesis and should be included in any future studies.
This study assessed BAT thermogenesis among a small sample of Finnish reindeer herders from sub-Arctic Finland. We found that the herders do indeed have active BAT thermogenesis that resulted in a mean ~8.7% increase in metabolic rate in response to a mild cold stress. There were no statistically detectable anatomical correlates with BAT thermogenesis in this sample, though that is common among other BAT studies. Finally, this sample demonstrated a metabolic preference for fatty acid utilization during both room temperature and mild cold exposure, which was different from findings among the circumpolar Yakut. This study adds to the rapidly growing body of work looking at the physiological and thermoregulatory significance of brown adipose tissue and the key role it may play among cold stressed populations. The current state of brown adipose research is one in which we have far more questions than answers, making it rich topic for future research.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article, with the age and sex removed to protect participant confidentiality. The full dataset will be made available upon request.
Brown adipose tissue
Resting metabolic rate
Total energy expenditure
High density lipoprotein
Fat free mass
- MRRT :
Room temperature metabolic rate
- MRC :
Cold exposure metabolic rate
- RQRT :
Room temperature respiratory quotient
- RQC :
Cold exposure metabolic rate
- TSC :
Supraclavicular surface temperature
- TST :
Sternum surface temperature
- STRT :
Room temperature sternum surface temperature
- STC :
Cold exposure sternum surface temperature
- SCRT :
Room temperature supraclavicular surface temperature
- SCC :
Cold exposure supraclavicular surface temperature
Lock M. Recovering the body. Annu Rev Anthropol. 2017;46(1):1–14.
Lock M, Nguyen V-K. An Anthropology of Biomedicine. Oxford, UK: Wiley-Blackwell; 2018.
Rosinger AY. Household water insecurity after a historic flood: diarrhea and dehydration in the Bolivian Amazon. Soc Sci Med. 2018;197:192–202.
Markkula I, Turunen M, Rasmus S. A Review of climate change impacts on the ecosystem services in the Saami Homeland in Finland PALOMA View project Gradual changes and abrupt crises-changing operational environment of Finnish reindeer View project A review of climate change impacts on the. Sci Total Environ. 2019;692:1–16.
Cohen J, Agel L, Barlow M, Garfinkel CI, White I. Linking Arctic variability and change with extreme winter weather in the United States. Science. 2021;373(6559):1116–21.
Kivinen S, Rasmus S, Jylhä K, Laapas M. Long-Term climate trends and extreme events in Northern Fennoscandia (1914–2013). Clim. 2017;5(1):16.
Steegmann AT. Human cold adaptation: an unfinished agenda. Am J Hum Biol. 2007;19(2):218–27.
Conroy G, Pontzer H. Reconstructing human origins: a modern synthesis. New York: WW Norton & Company; 2012.
Stocks J, Taylor N, Tipton J, Greenleaf J. Human physiological responses to cold exposure. Aviat Sp Environ Med. 2004;75:444–57.
Moran E. Human adaptabilityan introduction to ecological anthropology. 3rd ed. Boulder, CO: Wesyview Press; 2008.
Galloway VA, Leonard WR, Ivakine E. Basal metabolic adaptation of the Evenki reindeer herders of Central Siberia. Am J Hum Biol. 2000;12(1):75–87.
Leonard WR, et al. Seasonal variation in basal metabolic rates among the yakut (Sakha) of Northeastern Siberia. Am J Hum Biol. 2014;26(4):437–45.
Leonard WR, et al. Climatic influences on basal metabolic rates among circumpolar populations. Am J Hum Biol. 2002;14(5):609–20.
Leonard WR, Snodgrass JJ, Sorensen MV. Metabolic adaptation in Indigenous Siberian populations. Annu Rev Anthropol. 2005;34(1):451–71.
Levy SB, et al. Seasonal and socioeconomic influences on thyroid function among the Yakut (Sakha) of Eastern Siberia. Am J Hum Biol. 2013;25(6):814–20.
Snodgrass JJ, Leonard WR, Tarskaia LA, Alekseev VP, Krivoshapkin VG. Basal metabolic rate in the Yakut (Sakha) of Siberia. Am J Hum Biol. 2005;17(2):155–72.
Rode A, Shephard RJ. Basal metabolic rate of inuit. Am J Hum Biol. 1995;7(6):723–9.
Ocobock C, Soppela P, Turunen MT, Stenbäck V, Herzig K-H. Elevated resting metabolic rates among female, but not male, reindeer herders from sub-arctic Finland. Am J Hum Biol. 2020;32(6):e23432.
Levy SB, et al. Brown adipose tissue, energy expenditure, and biomarkers of cardio-metabolic health among the Yakut (Sakha) of northeastern Siberia. Am J Hum Biol. 2018;30(6):e23175.
van der Lans AAJJ, Vosselman MJ, Hanssen MJW, Brans B, van Marken Lichtenbelt WD. Supraclavicular skin temperature and BAT activity in lean healthy adults. J Physiol Sci. 2016;66(1):77–83.
Cannon B, Nedergaard J. Yes, even human brown fat is on fire! J Clin Invest. 2012;122(2):486–9.
Ouellet V, et al. Brown adipose tissue oxidative metabolism contributles to energy expenditure during cold exposure in humans. J Clin Invest. 2012;122(2):545.
Niclou A, Ocobock C. Seasonal patterns of BAT activity imply metabolic buffering and role in human cold adaptation. Am J Hum Biol. 2020;32(S1):37.
Folk GJ. Introduction to environmental physiology. Philadelphia: Lea and Febiger; 1966.
Cepon TJ, et al. Circumpolar adaptation, social change, and the development of autoimmune thyroid disorders among the Yakut (Sakha) of Siberia. Am J Hum Biol. 2011;23(5):703–9.
Ouellet V, et al. Brown adipose tissue oxidative metabolism contributes to energy expenditure during acute cold exposure in humans. J Clin Invest. 2012;122(2):545–52.
Levy SB. Field and laboratory methods for quantifying brown adipose tissue thermogenesis. Am J Hum Biol. 2019;31(4):e23261.
Virtanen KA, et al. Functional brown adipose tissue in healthy adults. N Engl J Med. 2009;360(15):1518–25.
Cypess AM, et al. Identification and importance of brown adipose tissue in adult humans. N Engl J Med. 2009;360:1509–17.
Virtanen KA, et al. Functional brown adipose tissue in healthy adults. N Engl J Med. 2009;15(9):1518-25.
Sacks H, Symonds ME. Anatomical locations of human brown adipose tissue functional relevance and implications in obesity and type 2 diabetes. Diabetes. 2013;62:1783–90.
Roh HC, et al. Warming induces significant reprogramming of beige, but not brown, adipocyte cellular identity. Cell Metab. 2018;27(5):1121–1137.e5.
Alvarez-Dominguez JR, et al. De novo reconstruction of adipose tissue transcriptomes reveals long non-coding RNA regulators of brown adipocyte development. Cell Metab. 2015;21(5):764–76.
Ding C, et al. De novo reconstruction of human adipose transcriptome reveals conserved lncRNAs as regulators of brown adipogenesis. Nat Commun. 2018;9(1):1–14.
Hao Q, et al. Transcriptome profiling of brown adipose tissue during cold exposure reveals extensive regulation of glucose metabolism. Am J Physiol Endocrinol Metab. 2015;308(5):E380–92.
Chondronikola M, et al. Brown adipose tissue improves whole-body glucose homeostasis and insulin sensitivity in humans. Diabetes. 2014;63(12):4089–99.
Vallerand AL, Jacobs I. Rates of energy substrates utilization during human cold exposure. Eur J Appl Physiol Occup Physiol. 1989;58(8):873–8.
Hoeke G, et al. Short-term cooling increases serum triglycerides and small high-density lipoprotein levels in humans. J Clin Lipidol. 2017;11(4):920–928.e2.
Gagnon DD, et al. Cold exposure enhances fat utilization but not non-esterified fatty acids, glycerol or catecholamines availability during submaximal walking and running. Front Physiol. 2013;4 MAY:99.
McCue A, Munten S, Herzig KH, Gagnon DD. Metabolic flexibility is unimpaired during exercise in the cold following acute glucose ingestion in young healthy adults. J Therm Biol. 2021;98:102912.
Drubach LA, Palmer EL, Connolly LP, Baker A, Zurakowski D, Cypess AM. Pediatric brown adipose tissue: detection, epidemiology, and differences from adults. J Pediatr. 2011;159:939–44.
Wakabayashi H, Nishimura T, Wijayanto T, Watanuki S, Tochihara Y. Effect of repeated forearm muscle cooling on the adaptation of skeletal muscle metabolism in humans. Int J Biometeorol. 2017;61(7):1261–7.
Bennett AF. Temperature and muscle. J Exp Biol. 1985;115(1):333–44.
RHA, “Reindeer Herders Association Statistics,” 2021. [Online]. Available: https://paliskunnat.fi/.
FMI, “Finnish Meteorological Institute,” 2021. [Online]. Available: https://en.ilmatieteenlaitos.fi/. [Accessed: 11-Aug-2021].
Kohonen I, Kuula-Luumi A, Spoof S-K. The ethical principles of research with human participants and ethical review in the human sciences in Finland. Helsinki: Finnish National Board on Research Integrity; 2019.
Sun SS, et al. Development of bioelectrical impedance analysis prediction equations for body composition with the use of a multicomponent model for use in epidemiologic surveys. Am J Clin Nutr. 2003;77(2):331–40.
Nirengi S, et al. An optimal condition for the evaluation of human brown adipose tissue by infrared thermography. PLoS One. 2019;14(8):e0220574.
Ocobock C, et al. Reindeer herders from subarctic Finland exhibit high total energy expenditure and low energy intake during the autumn herd roundup. Am J Hum Biol. 2021;34(4):e23676.
Niclou A, Ocobock C. Weather permitting: increased seasonal efficiency of non-shivering thermogenesis through brown adipose tissue activation in the winter. Am J Hum Biol. 2021:e23716.
Yoneshiro T, et al. Brown adipose tissue is involved in the seasonal variation of cold-induced thermogenesis in humans. Am J Physiol Regul Integr CompPhysiol. May 2016;310(10):R999–R1009.
Van Der Lans AAJJ, et al. Cold acclimation recruits human brown fat and increases nonshivering thermogenesis. J Clin Invest. 2013;123(8):3395–403.
Levy SB, et al. Evidence for a sensitive period of plasticity in brown adipose tissue during early childhood among indigenous Siberians. Am J Phys Anthropol. 2021;175(4):834–46.
Saito M, et al. High Incidence of metabolically active brown adipose tissue in healthy adult humans: effects of cold exposure and adiposityorg/licenses/by-nc-nd/3.0/ for details. Diabetes. 2009;58:1526–31.
Matsushita M, Yoneshiro T, Aita S, Kameya T, Sugie H, Saito M. Impact of brown adipose tissue on body fatness and glucose metabolism in healthy humans. Int J Obes. 2014;38(6):812–7 Nov. 2013.
Hanssen MJW, et al. Serum FGF21 levels are associated with brown adipose tissue activity in humans. Sci Rep. 2015;5(1):1–8.
Lee P, et al. Cold-activated brown adipose tissue is an independent predictor of higher bone mineral density in women. Osteoporos Int. 2013;24(4):1513–8.
Bahler L, et al. Differences in sympathetic nervous stimulation of brown adipose tissue between the young and old, and the lean and obese. J Nucl Med. 2016;57(3):372–7.
Franssens BT, Hoogduin H, Leiner T, Van Der Graaf Y, Visseren FLJ. Relation between brown adipose tissue and measures of obesity and metabolic dysfunction in patients with cardiovascular disease. Int Soc Magn Reson Med. 2017;46:497–504.
Yoneshiro T, et al. Brown adipose tissue, whole-body energy expenditure, and thermogenesis in healthy adult men. Obesity. 2011;19(1):13–6.
Yoneshiro T, Aita S, Kawai Y, Iwanaga T, Saito M. Nonpungent capsaicin analogs (capsinoids) increase energy expenditure through the activation of brown adipose tissue in humans. Am J Clin Nutr. 2012;95(4):845–50.
Waldén TB, Hansen IR, Timmons JA, Cannon B, Nedergaard J. Recruited vs. nonrecruited molecular signatures of brown, ‘brite,’ and white adipose tissues. Am J Physiol Endocrinol Metab. 2012;302(1):19–31.
Long JZ, et al. A smooth muscle-like origin for beige adipocytes. Cell Metab. 2014;19(5):810–20.
Wells JC. Adaptive variability in the duration of critical windows of plasticityImplications for the programming of obesity. Evol Med Public Heal. 2014;2014(1):109–21.
Wells JCK. Flaws in the theory of predictive adaptive responses. Trends Endocrinol Metab. 2007;18(9):331–7.
Wells JCK. Developmental plasticity as adaptation: adjusting to the external environment under the imprint of maternal capital. Philos Trans R Soc B Biol Sci. 2019;374:1–8.
Beall CM, Decker MJ, Brittenham GM, Kushner I, Gebremedhin A, Strohl KP. An Ethiopian pattern of human adaptation to high-altitude hypoxia. Proc Natl Acad Sci U S A. 2002;99(26):17215–8.
Beall CM. Two routes to functional adaptation: tibetan and Andean high-altitude natives. Proc Natl Acad Sci U S A. 2007;104(S1):8655–60.
Roseman CC, Auerbach BM. Ecogeography, genetics, and the evolution of human body form. J Hum Evol. 2015;78:80–90.
Tourula M, Isola A, Hassi J. Children sleeping outdoors in winter: parents’ experiences of a culturally bound childcare practice. Int J Circumpolar Health. 2008;67(2–3):269–78.
Tourula M, Isola A, Hassi J, Bloigu R, Rintamäki H. Infants sleeping outdoors in a northern winter climate: skin temperature and duration of sleep. Acta Paediatr. 2010;99(9):1411–7.
Husu P, Paronen O, Suni J, Vasankari T. Suomalaisten fyysinen aktiivisuus ja kunto 2010. Terveyttä edistävän liikunnan nykytila ja muutokset [The physical activity and physical condition in 2010 among Finns. The current state and changes in health enhancing physical activity]. Opetus Ja Kultt Julk. 2011;15:1-87.
Bahler L, Deelen JW, Hoekstra JB, Holleman F, Verberne HJ. Seasonal influence on stimulated BAT activity in prospective trials: a retrospective analysis of BAT visualized on 18F-FDG PET-CTs and 123I-mIBG SPECT-CTs. J Appl Physiol. 2016;120(12):1418–23.
Turunen M, Soppela P, Ocobock C. How reindeer herders cope with harsh winter conditions in northern Finland: insights from an interview study. Arctic. 2021;74(2):188–205.
McArdle WD. Sports and exercise nutrition. Philadelphia: Lippincott Williams & Wilkins; 2018.
We thank Executive Director Anne Ollila and Advisors Anna-Leena Jänkälä and Maaren Angeli from the Reindeer Herder’s Association in Finland who helped us to recruit the reindeer herding districts and herders to this study. We are grateful to the participants and participating districts for not only taking part in this study but also providing valuable input along the way. Thank you to the University at Albany for supporting the early efforts of this research. Thank you also to incredibly helpful conversations with Alexandra Niclou. Thank you to Adam Gordon for his support and guidance. This study also would not have happened without the aid of Pertti Leinonen who provided fantastic logistical support at the Arctic Centre. This study also would not have been possible without the early advice and input from Tim Frandy, Bob Jarvenpa, Hannu Heikkinen, Hannu Rintamäki, and Terhi Vuojala-Magga.
Funding for this project was provided by the National Science Foundation High Risk Research in Biological Anthropology and Archaeology grant (1724819) and the American Scandinavian Foundation.
Ethics approval and consent to participate
This study was conducted with Institutional Review Board approval from the University at Albany (17-E-165) and with the approval of Ostrobothnian Health Care District from the University of Oulu, Finland (EETTMK: 4/2018).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ocobock, C., Soppela, P., Turunen, M. et al. Brown adipose tissue thermogenesis among a small sample of reindeer herders from sub-Arctic Finland. J Physiol Anthropol 41, 17 (2022). https://doi.org/10.1186/s40101-022-00290-4