Abstract

Per‐ and polyfluoroalkyl substances (PFAS) are associated with several disrupted physiological and endocrine parameters. Regarding endocrine mechanisms, laboratory studies suggest that PFAS could disrupt the thyroid hormone system and alter circulating thyroid hormone concentrations. Thyroid hormones play a ubiquitous role—controlling thermoregulation, metabolism, and reproduction. However, evidence for disruption of thyroid hormones by PFAS remains scarce in wildlife. The present study investigated the associations between concentrations of PFAS, thyroid hormones, and body condition in an arctic seabird, the black‐legged kittiwake (Rissa tridactyla). We collected blood from kittiwakes sampled in Svalbard, Norway (2013 and 2014). Plasma samples were analyzed for total thyroxine (TT4) and total triiodothyronine (TT3) concentrations; detected PFAS included branched and linear (lin) C8 perfluoroalkyl sulfonates (i.e., perfluoroctane sulfonate [PFOS]) and C9‐C14 perfluoroalkyl carboxylates (PFCAs). The dominant PFAS in the kittiwakes were linPFOS and C11‐ and C13‐PFCAs. Generally, male kittiwakes had higher concentrations of PFAS than females. We observed positive correlations between linPFOS, C10‐PFCA, and TT4 in males, whereas in females C12‐14‐PFCAs were positively correlated to TT3. Interestingly, we observed contrasted correlations between PFAS and body condition; the direction of the relationship was sex‐dependent. Although these results show relationships between PFAS and circulating thyroid hormone concentrations in kittiwakes, the study design does not allow for concluding on causal relationships related to effects of PFAS on the thyroid hormone system. Future experimental research is required to quantify this impact of PFAS on the biology of kittiwakes. The apparently different associations among PFAS and body condition for males and females are puzzling, and more research is required. Environ Toxicol Chem 2021;40:820–831. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Abstract

The thyroid hormones (THs) triiodothyronine and thyroxine and per‐ and polyfluoroalkyl substances (PFAS) were examined in plasma from black‐legged kittiwakes (Rissa tridactyla). We found a positive association between several PFAS and the THs. T3 = triiodothyronine; T4 = thyroxine.

INTRODUCTION

Per‐ and polyfluoroalkyl substances (PFAS) are synthetically produced chemicals used as surfactants and repellents. Once released into the environment they are able to undergo long‐range transport from areas of high use to remote areas with few local sources such as the Arctic (Young et al. 2007). Following transport to the Arctic, PFAS are available for uptake by local wildlife and have been detected in a range of species such as mammals, fish, and birds (Letcher et al. 2010).

A number of studies—laboratory and field studies—have investigated the potential effects of PFAS on organisms. These substances seem to affect lipid storage; for instance, a laboratory‐based study has demonstrated that PFAS can affect fatty acid and lipid synthesis, storage, and metabolism in a fish model (Cheng et al. 2016). Tartu et al. (2017) investigated the effects of organohalogenated contaminants on biomarkers of energy metabolism in polar bears (Ursus maritimus) and found a positive association between PFAS and the transcript levels of the sterol regulatory element–binding protein 1 gene (involved in synthesis of fatty acids) and plasma levels of cholesterol, high‐density lipoprotein, and triglycerides.

Furthermore, PFAS are endocrine‐disrupting chemicals, targeting, among other endocrine systems, the hypothalamus–pituitary–thyroid axis and the secretion and regulation of thyroid hormones. The thyroid hormones thyroxine (T4) and triiodothyronine (T3) are released from the thyroid gland into the bloodstream, where the majority of the hormones bind to transthyretin and albumin (Ishihara et al. 2003). The active form of the hormone is T3, and T4 is converted to T3 by iodothyronine deiodinases (Ishihara et al. 2003). Thyroid hormones are required for proper metabolism and thermoregulation, body weight control, growth, reproduction, molting, brain development, and cell differentiation (Merryman and Buckles 1998).

Laboratory exposure experiments on fish indicate effects of PFAS on thyroid hormone production (Shi et al. 2009; Liu et al. 2011; Zhang et al. 2018). But less is known for wildlife, and there are contrasting reports of the potential effects of PFAS on the thyroid hormone system in wildlife species in the scientific literature. Positive correlations between various PFAS and circulating concentrations of T4 and/or T3 have been reported for seabirds (Braune et al. 2011; Nøst et al. 2012; Melnes et al. 2017; Hovden 2018). In contrast, negative correlations have been found between PFAS and thyroid hormones in polar bears (U. maritimus; Bourgeon et al. 2017), suggesting that the impact of PFAS on this endocrine axis may depend on the degree of PFAS contamination, the species, or the environmental context.

In the Arctic, wildlife faces extreme environmental challenges because of the low temperatures and demanding climate. Even for migratory species, such as seabirds, which only remain in the Arctic for the summer season, the conditions during summer are characterized by low ambient temperatures and frequently inclement weather. Reproduction is a demanding event (Drent and Daan 1980), which when combined with challenging environmental conditions necessitates proper thyroid function, in particular because of the function of thyroid hormones in thermoregulation, reproduction, and energetics (Merryman and Buckles 1998). If the thyroid system is disrupted, it may impact the ability of the species to cope with environmental challenges, store energy, grow, and reproduce, thus potentially affecting avian population demography (Ruuskanen et al. 2021).

The black‐legged kittiwake (Rissa tridactyla, henceforth “kittiwake”), an arctic and boreal seabird, is currently experiencing population declines along the coast of mainland Norway (Anker‐Nilssen et al. 2019). The causes are still unknown, and although factors related to climate change and prey availability are likely partially responsible, further research on the significance of endocrine‐disrupting anthropogenic pollutants, such as PFAS, is required to understand these population declines.

Studies have reported PFAS in blood and liver of kittiwakes (e.g., Letcher et al. 2010; Nøst et al. 2012; Tartu et al. 2014; Blévin et al. 2017a). Because PFAS are a relatively new class of toxicants to be investigated, there have been no extensive studies on how they affect thyroid hormones in kittiwakes. To our knowledge, only 2 studies have investigated this so far: Nøst et al. (2012) reported positive associations between PFAS and thyroid hormone concentrations in kittiwake chicks in Svalbard; Blévin et al. (2017b), on the other hand, did not find any significant associations between PFAS and thyroid hormone concentrations in adult kittiwakes. It is possible that the relatively low sample size used in that particular study (n = 35) may have prevented the detection of significant relationships. There are differences in thyroid hormone status and function between chicks and adults (McNabb and Hughes 1983). The role of thyroid hormones in chicks is important for the development of endothermy and metabolism and for growth and development, where the hormones act to stimulate growth and cell differentiation (McNabb 2007). There is a large body of evidence demonstrating effects of thyroid hormones on tissue oxygen consumption and a strong link with basal or resting metabolism (Hulbert 2000; Kim 2008; Welcker et al. 2013), and basal metabolic rate (BMR) has a role in short‐term thermal acclimation and endothermy (McKechnie 2008; Sirsat and Dzialowski 2020). In addition, thyroid hormones are important factors in regulating reproduction (Nishiwaki‐Ohkawa and Yoshimura 2016) and thus reproductive success. Therefore, a disruption of thyroid hormones caused by exposure to PFAS in arctic breeding seabirds may impair their capability to deal with cold temperatures, interfere with lipid stores and lipid metabolism, and affect reproductive success, especially in a context of environmental changes (Ruuskanen et al. 2021).

The aim of the present study was therefore to further examine the relationships between plasma concentrations of PFAS and thyroid hormone concentrations (total T4 [TT4] and total T3 [TT3]) in chick‐rearing kittiwakes breeding in Kongsfjorden, Svalbard, Norway. Furthermore, the study investigated whether PFAS affected the body condition of the kittiwakes, which was used as a proxy for lipid storage.

MATERIALS AND METHODS

Sample collection and field procedures

Kittiwakes were all sampled from one colony (“Krykkjefjellet” located 6 km southeast of Ny‐Ålesund, 78°55′N, 11°56′E) in Kongsfjorden, during the chick‐rearing period of their breeding cycle. The sampling periods were between 3 and 6 July in 2013 and between 14 and 20 July in 2014. Each year 15 male and 15 female kittiwakes were sampled in the colony, resulting in 60 samples. The kittiwakes were caught on their nests using a telescopic rod with a nylon noose attached at the end. A maximum of 2.5 mL of blood was collected from the brachial vein. The blood was stored on ice during fieldwork and centrifuged on returning to the laboratory in Ny‐Ålesund. After centrifugation the plasma was transferred into another Eppendorf tube and frozen at –20 °C. The remaining red blood cells were kept in the original Eppendorf tubes and frozen at –20 °C to be used for sex determination.

Birds were weighed to the nearest 2 g using a Pesola spring balance, and their skull length (head and bill), wing length, and tarsus length were measured to the nearest 0.5 mm with a sliding caliper (data summarized in Supplemental Data, Table S1). The birds were monitored after release to ensure that they returned to the nest, which they all did.

The study was approved by the governor of Svalbard and the Norwegian Animal Research Authority. All handling and sampling of the birds were in accordance with the regulations of the Norwegian Animal Welfare Act.

Analyses of PFAS

The PFAS analyses were performed at the Norwegian Institute for Air Research in Tromsø, Norway, on plasma samples. The PFAS analytes are listed in Table 1 together with their chemical structures. The samples were analyzed for PFAS as described by Hanssen et al. (2013). More details on how the analysis was performed can be found in the Supplemental Data.

Table 1.

The per‐ and polyfluoroalkyl substances analyzed in chick‐rearing black‐legged kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014, divided into perfluoroalkyl carboxylic acids and perfluoroalkyl sulfonatesa

GroupAbbreviationAnalyteChemical structure
PFCAsPFHxAPerfluorohexanoateC5F11COO
PFHpAPerfluoroheptanoateC6F13COO
PFOAPerfluorooctanoateC7F15COO
PFNAPerfluorononanoateC8F17COO
PFDAPerfluorodecanoateC9F19COO
PFUnDAPerfluoroundecanoateC10F21COO
PFDoDAPerfluorododecanoateC11F23COO
PFTrDAPerfluorotridecanoateC12F25COO
PFTeDAPerfluorotetradecanoateC13F27COO
PFSAsPFHxSPerfluorohexane sulfonate anionC6F13SO3
brPFOSBranched perfluoroctane sulfonate anionC8F17SO3
linPFOSLinear perfluoroctane sulfonate anionC8F17SO3
PFDSPerfluorodecane sulfonate anionC10F21SO3
GroupAbbreviationAnalyteChemical structure
PFCAsPFHxAPerfluorohexanoateC5F11COO
PFHpAPerfluoroheptanoateC6F13COO
PFOAPerfluorooctanoateC7F15COO
PFNAPerfluorononanoateC8F17COO
PFDAPerfluorodecanoateC9F19COO
PFUnDAPerfluoroundecanoateC10F21COO
PFDoDAPerfluorododecanoateC11F23COO
PFTrDAPerfluorotridecanoateC12F25COO
PFTeDAPerfluorotetradecanoateC13F27COO
PFSAsPFHxSPerfluorohexane sulfonate anionC6F13SO3
brPFOSBranched perfluoroctane sulfonate anionC8F17SO3
linPFOSLinear perfluoroctane sulfonate anionC8F17SO3
PFDSPerfluorodecane sulfonate anionC10F21SO3
a

The limit of detection was <0.05 ng/mL for all compounds.

PFCA = perfluoroalkyl carboxylic acid; PFSA = perfluoroalkyl sulfonate.

Table 1.

The per‐ and polyfluoroalkyl substances analyzed in chick‐rearing black‐legged kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014, divided into perfluoroalkyl carboxylic acids and perfluoroalkyl sulfonatesa

GroupAbbreviationAnalyteChemical structure
PFCAsPFHxAPerfluorohexanoateC5F11COO
PFHpAPerfluoroheptanoateC6F13COO
PFOAPerfluorooctanoateC7F15COO
PFNAPerfluorononanoateC8F17COO
PFDAPerfluorodecanoateC9F19COO
PFUnDAPerfluoroundecanoateC10F21COO
PFDoDAPerfluorododecanoateC11F23COO
PFTrDAPerfluorotridecanoateC12F25COO
PFTeDAPerfluorotetradecanoateC13F27COO
PFSAsPFHxSPerfluorohexane sulfonate anionC6F13SO3
brPFOSBranched perfluoroctane sulfonate anionC8F17SO3
linPFOSLinear perfluoroctane sulfonate anionC8F17SO3
PFDSPerfluorodecane sulfonate anionC10F21SO3
GroupAbbreviationAnalyteChemical structure
PFCAsPFHxAPerfluorohexanoateC5F11COO
PFHpAPerfluoroheptanoateC6F13COO
PFOAPerfluorooctanoateC7F15COO
PFNAPerfluorononanoateC8F17COO
PFDAPerfluorodecanoateC9F19COO
PFUnDAPerfluoroundecanoateC10F21COO
PFDoDAPerfluorododecanoateC11F23COO
PFTrDAPerfluorotridecanoateC12F25COO
PFTeDAPerfluorotetradecanoateC13F27COO
PFSAsPFHxSPerfluorohexane sulfonate anionC6F13SO3
brPFOSBranched perfluoroctane sulfonate anionC8F17SO3
linPFOSLinear perfluoroctane sulfonate anionC8F17SO3
PFDSPerfluorodecane sulfonate anionC10F21SO3
a

The limit of detection was <0.05 ng/mL for all compounds.

PFCA = perfluoroalkyl carboxylic acid; PFSA = perfluoroalkyl sulfonate.

For every batch of 10 samples, one tube with a blank and one with a standard reference material (SRM 1957; National Institute of Standards and Technology) were included. Together with an internal and a recovery standard, this formed the basis of the quality assurance. The internal standard contained known concentrations of 13C‐labeled PFAS. This standard was added at the beginning of the extraction with the assumption of an equal loss of the internal standard and the 12C PFAS in the sample through the extraction and cleanup processes. The recovery standard was added just before running the samples on the instrument for the purpose of determining the performance of the analytical method used. The reference samples were within the acceptable range, and the blank samples were below the limit of acceptable contamination determined by the laboratory.

Analyses of thyroid hormones

The hormone assays for TT3 and TT4 were performed at the Centre d'Etudes Biologiques de Chizé, France, using radioimmunoassay, as described by Blévin et al. (2017b).

For quality control, reference materials were used, which were within the acceptable range determined by the laboratory. Furthermore, because every plasma sample was analyzed in duplicate, only samples with a coefficient of variance (CV) <15% were used for statistical analyses. The mean CV calculated for the TT4 assays was 3.82% in 2013 and 3.34% in 2014. The mean CV for the TT3 assays was 4.76% in 2013 and 4.50% in 2014. The mean interassay variation for TT3 was 8.98%, and that for TT4 was 11.4%.

Sex determination

Sex determination by polymerase chain reaction and gel electrophoresis was performed for those individuals caught which had not been previously molecularly sexed at the Norwegian University of Science and Technology, Trondheim. The analysis was performed as described in Griffiths et al. (1998) and outlined in the Supplemental Data.

Data treatment and statistical analyses

From a total of 60 kittiwakes caught, 4 individuals were sampled in both years. All these individuals were females, and in the data treatment only one random (determined by a coin toss) observation for each of these individuals was included.

Furthermore, one male and one female kittiwake were excluded from the TT4 data sets because of too high variance between their duplicated plasma samples (>15%). Thus, the final sample size for male kittiwakes was 30 individuals for biometric variables and PFAS concentrations and 29 individuals for TT4 concentrations. For female kittiwakes the final sample size was 26 individuals for biometric variables and for TT3 and PFAS concentrations and 25 individuals for TT4 concentrations.

Statistical analyses were performed using R, Ver 3.6.1 (R Development Core Team 2019). The data (concentrations of PFAS and thyroid hormones) were resampled by bootstrapping (using the infer package, 5000 replicates). Some of the PFAS variables contained observations which were below the limit of detection (LOD; <0.05 ng/mL for all compounds) and were thus left‐censored. Specifically, branched perfluoroctane sulfonate anion (brPFOS), perfluorotetradecanoate (PFTeDA), perfluorohexanoate (PFHxA), perfluoroheptanoate (PFHpA), perfluorooctanoate (PFOA), perfluorohexane sulfonate (PFHxS), and perfluorodecane sulfonate (PFDS) contained observations below the LOD; PFHxA, PFHpA, PFOA, PFHxS, and PFDS contained >70% censored observations and were omitted from analyses. Because it is inappropriate to substitute nondetects with a value (Helsel 2006), the nonparametric Kaplan‐Meier technique was used on brPFOS and PFTeDA using the NADA package to obtain the descriptive statistics. Based on the bootstrapped data, mean values were calculated for the mean, the 95% confidence interval (CI) of the mean, and median.

The body condition of each individual bird was determined using the biometric variable (length of skull, wing, or tarsus) that was most correlated to body mass. Skull length was most strongly correlated to body mass (r = 0.50, p < 0.001). The standardized residual of linear regression plots of skull length as the independent variable against body mass as the dependent variable was determined to be the body condition.

The kittiwake data set was examined for differences between the years (2013 and 2014) by comparing the 95% CIs of the means, as described by Erickson and Rattner (2020; see specifically the section “Confidence intervals”). Differences in thyroid hormone and PFAS mean concentrations between male and female kittiwakes were also examined by comparing their 95% CIs. Potential associations between thyroid hormones and body condition were investigated using linear regression and considered meaningful if the 95% CI of the regression slope did not contain zero (Ismay and Kim 2020).

To examine if levels of thyroid hormones were associated with PFAS levels—thus indicating an effect of PFAS on thyroid hormones—a 2‐pronged approach was used: principal component analysis (PCA) and linear regression. Because of the issue of multicollinearity of the PFAS, PCA was performed to obtain principal components (PCs), which were subsequently regressed against the thyroid hormones. Because PCA is an unsupervised feature extraction method (Jolliffe 2002), the linear regression was performed with all PCs; and those which were significant were selected and are presented in the Supplemental Data. Afterward, linear regression of single PFAS was employed to further investigate which PFAS were correlated with thyroid hormones because effects of PFAS may depend on their carbon chain length (Costantini et al. 2019). To account for the censored PFAS—brPFOS and PFTeDA—linear regression with maximum likelihood estimation was employed, as described in Jaspers et al. (2013). Again, the correlations were considered meaningful if the 95% CI of the regression slope did not contain zero. The plotted linear regressions were bias‐corrected by subtracting the mean difference in “original” coefficients and bootstrapped coefficients from the “original” coefficients. Likewise, associations between PFAS and the body condition of the kittiwakes were also examined with linear regression models. R scripts can be found in the Supplemental Data.

RESULTS

Differences in PFAS and thyroid hormone concentrations between 2013 and 2014

All PFAS and thyroid hormone concentrations are given in nanograms per milliliter. Female kittiwakes sampled in 2014 had a higher mean TT3 concentration than those sampled in 2013 (95% CI 2013, 1.88–2.83; 2014, 3.01–3.95). There was no difference in mean concentrations of TT3 between 2013 and 2014 for males (95% CI 2013, 2.31–2.91; 2014, 2.60–3.36). There were no differences in the mean TT4 concentrations in male and female kittiwakes from 2013 (95% CI males, 29.97–43.88; females, 29.19–49.61) compared to 2014 (95% CI males, 32.11–40.49; females, 29.80–42.39). On average, perfluorononanoate (PFNA) levels were higher in both males and females sampled in 2013 (95% CI males, 1.44–2.75; females, 1.14–2.15) than 2014 (95% CI males, 0.866–1.19; females, 0.788–1.06). In 2013 mean concentrations of C12‐C14 perfluoroalkyl carboxylic acids (PFCAs) were higher in male and female kittiwakes than they were in 2014 (95% CI [2013/2014] males: perfluorododecanoate [PFDoDA], 0.833–1.32/1.87–2.34; perfluorotridecanoate [PFTrDA], 3.39–4.61/5.92–7.16; PFTeDA, 0.254–0.410/0.867–1.13; and females: PFDoDA, 0.537–0.970/1.35–1.68; PFTrDA, 2.01–3.25/4.20–5.01; PFTeDA, 0.077–0.160/0.529–0.692). The differences in concentrations are shown in Figure 1. Although there were some differences in mean concentrations between years, the PFAS' mechanism(s) of action should not differ. Thus, the data sets from the 2 yr were pooled to obtain a larger sample size and variation in PFAS and thyroid hormone concentrations.

Mean concentrations (with 95% confidence intervals) of per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male and female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014. F = female; M = male; brPFOS = branched PFOS anion; linPFOS = linear PFOS anion; PFDA = perfluorodecanoate; PFDoDA = perfluorododecanoate; PFNA = perfluorononanoate; PFOS = perfluorooctane sulfonate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; PFUnDA = perfluoroundecanoate.
Figure 1.

Mean concentrations (with 95% confidence intervals) of per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male and female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014. F = female; M = male; brPFOS = branched PFOS anion; linPFOS = linear PFOS anion; PFDA = perfluorodecanoate; PFDoDA = perfluorododecanoate; PFNA = perfluorononanoate; PFOS = perfluorooctane sulfonate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; PFUnDA = perfluoroundecanoate.

PFAS concentrations

The mean, 95% CI of the mean and median in PFAS concentrations after bootstrapping, as well as the range (of the original data set) in PFAS concentrations in kittiwakes are shown in Table 2. In both males and females linear (lin) PFOS, perfluoroundecanoate (PFUnDA), and PFTrDA were the PFAS present at highest concentrations. There is considerable overlap in the CIs of linPFOS (males, 6.13–8.65; females, 3.27–5.68) and PFUnDA (males, 6.10–7.82; females, 4.99–6.30), suggesting no difference in concentrations for both male and female kittiwakes. In males, both these chemicals (linPFOS and PFUnDA) were detected at higher levels than PFTrDA (95% CI 4.59–5.91). In female kittiwakes, however, only PFUnDA was detected at a higher concentration compared to PFTrDA (95% CI 3.14–4.16).

Table 2.

Concentrations of thyroid hormones (nanograms per milliliter of plasma) and per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male and female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014a

MalesFemales
MeanCI of meanMedianRangeMeanCI of meanMedianRange
TH (ng/mL)TT32.782.54–3.052.671.75–4.592.962.55–3.363.051.12–4.49
TT436.4132.49–40.4836.9918.61–70.0537.2631.73–43.1535.139.53–79.83
PFAS (ng/mL)brPFOS0.1810.116–0.2510.11nd–0.6690.1440.090–0.2150.083nd–0.164
linPFOS7.326.13–8.656.312.47–19.34.353.27–5.683.620.937–14.6
PFNA1.531.19–1.951.220.377–5.611.250.988–1.551.050.415–3.34
PFDA1.291.11–1.481.210.392–2.651.110.949–1.271.050.360–2.08
PFUnDA6.946.10–7.826.243.47–12.25.624.99–6.305.582.47–9.49
PFDoDA1.581.32–1.841.560.429–2.961.160.963–1.361.220.158–2.07
PFTrDA5.254.59–5.915.092.22–8.423.673.14–4.163.831.20–6.07
PFTeDA0.6590.521–0.8010.6390.084–1.640.380.277–0.4870.339nd–1.01
MalesFemales
MeanCI of meanMedianRangeMeanCI of meanMedianRange
TH (ng/mL)TT32.782.54–3.052.671.75–4.592.962.55–3.363.051.12–4.49
TT436.4132.49–40.4836.9918.61–70.0537.2631.73–43.1535.139.53–79.83
PFAS (ng/mL)brPFOS0.1810.116–0.2510.11nd–0.6690.1440.090–0.2150.083nd–0.164
linPFOS7.326.13–8.656.312.47–19.34.353.27–5.683.620.937–14.6
PFNA1.531.19–1.951.220.377–5.611.250.988–1.551.050.415–3.34
PFDA1.291.11–1.481.210.392–2.651.110.949–1.271.050.360–2.08
PFUnDA6.946.10–7.826.243.47–12.25.624.99–6.305.582.47–9.49
PFDoDA1.581.32–1.841.560.429–2.961.160.963–1.361.220.158–2.07
PFTrDA5.254.59–5.915.092.22–8.423.673.14–4.163.831.20–6.07
PFTeDA0.6590.521–0.8010.6390.084–1.640.380.277–0.4870.339nd–1.01
a

The data have been bootstrapped (5000 replicates), and the mean of the descriptive statistics is presented (mean, 95% confidence interval of the mean, median, and range, of the original data set): For males, branched perfluoroctane sulfonate anion (brPFOS) was left‐censored; and for females, brPFOS and PFTeDA were left‐censored, and the values presented are estimated based on the Kaplan‐Meier technique.

nd = nondetect; brPFOS = branched perfluorooctane sulfonate anion; linPFOS = linear perfluorooctane sulfonate anion; PFDA = perfluorodecanoate; PFDoDA = perfluorododecanoate; PFNA = perfluorononanoate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; PFUnDA = perfluoroundecanoate; T3 = triiodothyronine; T4 = thyroxine.

Table 2.

Concentrations of thyroid hormones (nanograms per milliliter of plasma) and per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male and female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard, in 2013 and 2014a

MalesFemales
MeanCI of meanMedianRangeMeanCI of meanMedianRange
TH (ng/mL)TT32.782.54–3.052.671.75–4.592.962.55–3.363.051.12–4.49
TT436.4132.49–40.4836.9918.61–70.0537.2631.73–43.1535.139.53–79.83
PFAS (ng/mL)brPFOS0.1810.116–0.2510.11nd–0.6690.1440.090–0.2150.083nd–0.164
linPFOS7.326.13–8.656.312.47–19.34.353.27–5.683.620.937–14.6
PFNA1.531.19–1.951.220.377–5.611.250.988–1.551.050.415–3.34
PFDA1.291.11–1.481.210.392–2.651.110.949–1.271.050.360–2.08
PFUnDA6.946.10–7.826.243.47–12.25.624.99–6.305.582.47–9.49
PFDoDA1.581.32–1.841.560.429–2.961.160.963–1.361.220.158–2.07
PFTrDA5.254.59–5.915.092.22–8.423.673.14–4.163.831.20–6.07
PFTeDA0.6590.521–0.8010.6390.084–1.640.380.277–0.4870.339nd–1.01
MalesFemales
MeanCI of meanMedianRangeMeanCI of meanMedianRange
TH (ng/mL)TT32.782.54–3.052.671.75–4.592.962.55–3.363.051.12–4.49
TT436.4132.49–40.4836.9918.61–70.0537.2631.73–43.1535.139.53–79.83
PFAS (ng/mL)brPFOS0.1810.116–0.2510.11nd–0.6690.1440.090–0.2150.083nd–0.164
linPFOS7.326.13–8.656.312.47–19.34.353.27–5.683.620.937–14.6
PFNA1.531.19–1.951.220.377–5.611.250.988–1.551.050.415–3.34
PFDA1.291.11–1.481.210.392–2.651.110.949–1.271.050.360–2.08
PFUnDA6.946.10–7.826.243.47–12.25.624.99–6.305.582.47–9.49
PFDoDA1.581.32–1.841.560.429–2.961.160.963–1.361.220.158–2.07
PFTrDA5.254.59–5.915.092.22–8.423.673.14–4.163.831.20–6.07
PFTeDA0.6590.521–0.8010.6390.084–1.640.380.277–0.4870.339nd–1.01
a

The data have been bootstrapped (5000 replicates), and the mean of the descriptive statistics is presented (mean, 95% confidence interval of the mean, median, and range, of the original data set): For males, branched perfluoroctane sulfonate anion (brPFOS) was left‐censored; and for females, brPFOS and PFTeDA were left‐censored, and the values presented are estimated based on the Kaplan‐Meier technique.

nd = nondetect; brPFOS = branched perfluorooctane sulfonate anion; linPFOS = linear perfluorooctane sulfonate anion; PFDA = perfluorodecanoate; PFDoDA = perfluorododecanoate; PFNA = perfluorononanoate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; PFUnDA = perfluoroundecanoate; T3 = triiodothyronine; T4 = thyroxine.

Male kittiwakes had higher mean concentrations of linPFOS (95% CI males, 6.13–8.65; females, 3.27–5.68), PFDoDA (95% CI males, 1.32–1.84; females, 0.963–1.36), PFTrDA (95% CI males, 4.59–5.91; females, 3.14–4.16), and PFTeDA (95% CI males, 0.521–0.801; females, 0.277–0.487) compared to females.

Thyroid hormone concentrations and body condition

The mean, 95% CI of the mean, and median of TT3 and TT4 concentrations in kittiwakes after bootstrapping are presented in Table 2, together with the range in thyroid hormone concentrations of the original data set. The 95% CIs for both mean TT3 (males, 2.54–3.05; females, 2.55–3.36) and mean TT4 (males, 32.49–40.48; females, 31.73–43.15) concentrations overlapped greatly for males and females, indicating no difference in thyroid hormone concentrations between the sexes.

For female kittiwakes, there were no meaningful relationships between the thyroid hormones and body condition (95% CI of the regression slope TT3, –0.12 to 0.58; TT4, –0.23 to 0.50). The same was the case for males: no meaningful relationships were found between thyroid hormones and body condition (95% CI TT3, –0.21 to 0.47; TT4, –0.15 to 0.52).

Associations between PFAS, thyroid hormones, and body condition

Regressing the PCs from the PCA (Supplemental Data, Figure S2) against TT4 (Supplemental Data, Figure S3) and TT3 (Supplemental Data, Figure S4) showed positive associations between PFAS and both thyroid hormones.

In male kittiwakes, only TT4 was found to be positively correlated with linPFOS (95% CI of the regression slope, 0.08–0.70) and perfluorodecanoate (PFDA; 95% CI 0.09–0.67; Figure 2).

Linear regression (±95% confidence interval [CI], shaded area) of total thyroxine (nanograms per milliliter of plasma) against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 29). The 95% CI for the correlation is shown for each subfigure. corr = correlation; PFDA = perfluorodecanoate; linPFOS = linear perfluorooctane sulfonate anion; TT4 = total thyroxine.
Figure 2.

Linear regression (±95% confidence interval [CI], shaded area) of total thyroxine (nanograms per milliliter of plasma) against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 29). The 95% CI for the correlation is shown for each subfigure. corr = correlation; PFDA = perfluorodecanoate; linPFOS = linear perfluorooctane sulfonate anion; TT4 = total thyroxine.

In female kittiwakes, TT3 was positively correlated with PFDoDA (95% CI 0.08–0.70), PFTrDA (95% CI 0.12–0.71), and PFTeDA (95% CI 0.08–0.75; Figure 3).

Linear regression (±95% confidence interval [CI], shaded area) of total triiodothyronine (nanograms per milliliter of plasma) against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 26). The 95% CI for the correlation is shown for each subfigure. The red line in (C) is the limit of detection (LOD, 0.05 ng/mL) for PFTeDA; the 5 points under the LOD have been assigned a value of one‐half the LOD for the purpose of the visual presentation only. corr = correlation; PFDoDA = perfluorododecanoate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; TT3 = total triiodothyronine.
Figure 3.

Linear regression (±95% confidence interval [CI], shaded area) of total triiodothyronine (nanograms per milliliter of plasma) against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 26). The 95% CI for the correlation is shown for each subfigure. The red line in (C) is the limit of detection (LOD, 0.05 ng/mL) for PFTeDA; the 5 points under the LOD have been assigned a value of one‐half the LOD for the purpose of the visual presentation only. corr = correlation; PFDoDA = perfluorododecanoate; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate; TT3 = total triiodothyronine.

For male kittiwakes PFTrDA and PFTeDA were positively correlated with body condition (95% CI 0.07–0.66 for both correlations), as shown in Figure 4. For female kittiwakes, body condition was negatively correlated with linPFOS and PFNA (95% CI –0.73 to –0.14 and –0.76 to –0.22, respectively), presented in Figure 5.

Linear regression (±95% confidence interval [CI], shaded area) of body condition against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 30). The 95% CI for the correlation is shown for each subfigure. corr = correlation; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate.
Figure 4.

Linear regression (±95% confidence interval [CI], shaded area) of body condition against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in male chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 30). The 95% CI for the correlation is shown for each subfigure. corr = correlation; PFTeDA = perfluorotetradecanoate; PFTrDA = perfluorotridecanoate.

Linear regression of body condition against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 26). The 95% confidence interval (CI) for the correlation is shown for each subfigure. corr = correlation; PFNA = perfluorononanoate; PFOS = perfluorooctane sulfonate.
Figure 5.

Linear regression of body condition against natural log–transformed per‐ and polyfluoroalkyl substances (nanograms per milliliter of plasma) in female chick‐rearing kittiwakes (Rissa tridactyla) caught in Kongsfjorden, Svalbard (2013–2014, n = 26). The 95% confidence interval (CI) for the correlation is shown for each subfigure. corr = correlation; PFNA = perfluorononanoate; PFOS = perfluorooctane sulfonate.

Supplemental Data, Table S3, shows the 95% CIs of all the linear regressions.

DISCUSSION

PFAS concentrations

Linear PFOS, PFUnDA, and PFTrDA were the dominant PFAS in both male and female kittiwakes in the present study. In another study of kittiwakes in the chick‐rearing part of their breeding cycle, Tartu et al. (2014) found PFTrDA to be the dominant PFAS detected in 20 kittiwakes caught in 2012 from the same colony (“Krykkjefjellet,” Kongsfjorden) as in the present study. Furthermore, the concentrations of PFTrDA in male and female kittiwakes reported by Tartu et al. (2014) were higher than those observed in the present study. Nevertheless, PFTrDA was the third most dominant substance in both males and females in the present study. For the other compounds, levels reported in the present study were comparable to the levels reported in kittiwakes by Tartu et al. (2014).

There have been indications that there are no meaningful differences in PFAS concentrations between the sexes in arctic wildlife (Verreault et al. 2005; Butt et al. 2007). This contradicts the results reported by Blévin et al. (2017a) and the results reported in the present study. Sex differences have also been reported for some PFAS in great skuas (Stercorarius skua; Leat et al. 2013), lesser black‐backed gulls (Larus fuscus; Bustnes et al. 2008), glaucous gulls (Larus hyperboreus; Melnes et al. 2017; Haarr et al. 2018), and European shags (Phalacrocorax aristotelis; Carravieri et al. 2020). There also appeared to be sex differences in PFAS concentrations for black guillemots (Cepphus grylle) and arctic skuas (Stercorarius parasiticus), although the differences were not statistically significant (Haarr et al. 2018). These studies generally found higher concentrations of PFAS in males compared to females, which is consistent with the present study. In several avian species, PFAS have been detected in eggs (Verreault et al. 2005; Miller et al. 2015). Thus, the sex differences observed in the present study are likely a result of the female bird being able to deposit some of her PFAS burden into the eggs, as suggested by Bustnes et al. (2008). Nevertheless, the difference in PFAS levels may also be caused by a difference in exposure if males and females target different prey species and/or if the sexes have differences in biotransformation capacities. To the best of our knowledge, there are no published results on whether male and female kittiwakes select different prey species. However, a personal communication from our team (Sabrina Tartu, Centre d'Etudes Biologiques de Chizé and Université de La Rochelle, Villiers‐en‐Bois, France) suggests that blood concentrations of the stable nitrogen isotope 15N relative to the isotope 14N (i.e., δ15N; used as a proxy of trophic level [Hobson 1999]) in chick‐rearing male kittiwakes are significantly higher than those measured in females (logistic regression, estimate 0.28‰, t = 3.12, p = 0.002). Although fluorotelomer alcohols can be biotransformed to PFCAs (Nilsson et al. 2013), there is currently a lack of information on the role of biotransformation processes of PFAS in birds and, in particular, whether there are sex differences. Thus, maternal deposition of PFAS into eggs in addition to sex‐specific differences in diet appear as reliable factors to explain the observed sex difference in the PFAS concentrations.

Newsted et al. (2005) calculated avian toxicity reference values for PFOS based on experiments with northern bobwhite quail (Colinus virginianus). The predicted‐no‐effect concentration (PNEC) of PFOS in serum of a trophic level 4 avian predator is 2.4 µg/mL for males and 0.15 µg/mL for females. The single highest ∑PFOS concentrations for male and female kittiwakes were 2 and 1 orders of magnitude lower than the reported PNECs, respectively (0.02 and 0.015 µg/mL). The observed concentrations in the present study are thus well below the PNECs calculated by Newsted et al. (2005), and it is unlikely that the concentrations of PFOS and other PFAS pose an overt, acute toxicological risk to the kittiwakes. Nonetheless, we must remain cautious because the concentrations experienced by the kittiwakes could be sufficient for endocrine disruption and because other factors not taken into account in laboratory experiments could make wild species more sensitive to PFAS toxicity.

Associations between PFAS and thyroid hormone concentrations

In the present study male kittiwakes had positive correlations between their concentrations of TT4 and concentrations of linPFOS and PFDA. In females, concentrations of TT3 were positively correlated to concentrations of PFDoDA, PFTrDA, and PFTeDA.

The literature on the effects of PFAS on the thyroid hormone system in free‐ranging seabirds is scarce. However, Nøst et al. (2012) reported positive associations between TT4 and perfluoroheptane sulfonate (PFHpS), brPFOS, linPFOS, and PFNA in kittiwake chicks from Svalbard. That particular study also investigated associations in northern fulmar (Fulmarus glacialis) chicks, and TT4 was found to have positive associations with PFHpS, brPFOS, and PFNA. No associations were found between TT3 and PFAS for either of the 2 species included in that study. A study on glaucous gulls from Kongsfjorden, Svalbard, found a positive association between PFOS and free T3 in females (Melnes et al. 2017). Another study of female glaucous gulls from Kongsfjorden found positive correlations between PFDA, PFUnDA, PFTrDA, and PFTeDA and TT3 and between PFDA and TT4 (Haugerud 2011). Although a study of adult kittiwakes sampled in 2012 did not find any significant relationships between PFAS and thyroid hormones, possibly because of lower sample sizes compared to the present study, they reported a positive association between some long‐chained PFAS and BMR (Blévin et al. 2017b). A positive association of BMR with T3 has been shown in kittiwakes (Welcker et al. 2013) and other birds (Chastel et al. 2003). Thus, there are several indications that PFAS are positively associated with concentrations of TT3 and TT4 in seabirds. The weak correlations found in the present study are hence strengthened and likely to be indicative of an association between PFAS and thyroid hormones in kittiwakes.

It is interesting that the particular individual PFAS found to be associated with thyroid hormones vary between different studies, even for the same species. One explanation for this might be that PFAS are not evenly distributed spatially and temporally (Rauert et al. 2018), and levels of other halogenated compounds that are thyroid‐disruptive, such as chlorinated and brominated compounds, may also vary across time and space. This could cause complicated interactive thyroid hormone‐related effects via similar or dissimilar mechanisms of actions of the various chemicals. Indeed, in female glaucous gulls, Melnes et al. (2017) reported a positive relationship between PFOS and T3 but negative relationships between chlorinated compounds and T3, indicating that that these 2 groups of contaminants could have interactive effects on thyroid hormones in a complex manner. Thus, the effect of the exposure chemicals on thyroid hormones would depend on the concentrations in birds and the relative distribution of the various compounds, which could vary between years.

The mechanism(s) of action of PFAS on the thyroid hormone system is, as of yet, largely unknown. However, one potential disruptive effect of PFAS is their competition with thyroid hormones for binding proteins: binding potencies of various PFAS in glaucous gulls have been investigated for both transthyretin and albumin (Ucán‐Marin et al. 2010; Mortensen et al. 2020), and a recent study has found that albumin is the main transport protein for PFOS, PFOA, PFHxS, PFNA, and PFDA in human plasma (Forsthuber et al. 2020). Albumin is also the most important transport protein for thyroid hormones in avian species (McNabb 2007). Most of the T4 and T3 in the plasma is thus bound to proteins, with much lower fractions of T4 and T3 being free in the plasma. Because both thyroid hormones and PFAS are associated with plasma proteins, it is possible that the positive relationship between TT4 and TT3 and the PFAS is caused by differing plasma protein concentrations in the birds. Thus, birds with low plasma protein concentrations would have low concentrations of both total thyroid hormones and PFAS, whereas birds with high plasma protein concentrations would have high thyroid hormone and PFAS concentrations. Information on the free fraction of thyroid hormones in the plasma would also be necessary in the present study to be able to better assess the importance of binding site competition between thyroid hormones and PFAS as an explanation for the increased thyroid hormone levels observed.

As for the implications of the increased thyroid hormone levels observed in the current work on kittiwake physiology, the use of plasma concentrations of thyroid hormones to diagnose hyper‐ and hypothyroidism in birds has been reported to be unreliable (Merryman and Buckles 1998). Moreover, although reference values of thyroid hormones exist for some species (Spano et al. 1987), there are considerable differences between the species, hindering any meaningful comparison. Thus, it is unknown whether the thyroid hormone concentrations measured in the kittiwakes in the present study are within the normal range of hormone concentrations or not. Nevertheless, a finely tuned regulation of the thyroid hormone system is crucial for free‐ranging species, and a potential disruption of thyroid hormones could lead to inappropriate regulation of the basal metabolism, leading to inappropriate storage of lipids, which could become highly problematic because females need to store lipids prior to egg formation. Furthermore, when considering the positive relationship reported between TT3 and BMR in adult kittiwakes (Welcker et al. 2013), it can be hypothesized that increased TT3 levels caused by high body burdens of certain PFAS could lead to an increased BMR, which could impair thermal acclimation or acclimatization in kittiwakes (Hulbert 2000; McKechnie 2008). Lastly, an increase in thyroid hormone concentrations could disrupt the adaptive decrease in BMR in kittiwakes during chick‐rearing (Bech et al. 2002).

Associations between PFAS and body condition

In male kittiwakes, there were positive correlations between body condition and 2 PFCAs: PFTrDA and PFTeDA. In females, however, body condition was negatively correlated to linPFOS and PFNA.

A previous study on kittiwakes found a positive correlation between PFNA and body condition in kittiwake males (Tartu et al. 2014), whereas PFNA was not correlated to body condition in males in the present study, although these variables were correlated in females. The concentrations of PFNA in these 2 studies are similar; 1.24 ± 0.55 ng/mL (mean ± standard deviation) in the study by Tartu et al. (2014) and 1.53 (95% CI 1.19–1.95) ng/mL in the present study. Furthermore, in contrast to the present study, Tartu et al. (2014) found no significant relationships of the investigated PFAS and body condition in female kittiwakes, whereas in the present study a negative correlation was found between body condition and linPFOS and PFNA in the females. It should also be noted that using a multivariate statistical modeling approach, Melnes et al. (2017) identified plasma levels of PFOA, PFHxS, and perfluorooctane sulfonamide as strong determinants of having a high body condition. These apparently contrasting sex‐ and species‐dependent responses in body condition to PFAS, as indicated by the correlations, highlight the need for more research on this topic before the biological effects, if any, of PFAS on the body condition and lipid storage in kittiwakes can be determined.

It is interesting that the relationship between PFAS and body condition is contrasted in male and female kittiwakes in the present study. Cheng et al. (2016) reported that male zebrafish (Danio rerio) had increased mRNA expression of 4 nuclear receptors involved in adipogenesis following PFOS exposure, whereas females only had a modest increase in expression of the pregnane X receptor. Furthermore, they noted a sex‐specific difference in concentration of serum lipids. Sex differences in energy homeostasis are also recognized in mammals (Mauvais‐Jarvis 2015; Link and Reue 2017), but the scientific literature on this topic in birds is scarce. Still, the contrasting associations between body condition and PFAS observed in the present study might be attributable to, for example, the role of sex hormones in adipogenesis (Palmer and Clegg 2015) and/or genetics (Link and Reue 2017).

Another possible explanation for the positive associations between body condition and the 2 PFCAs in male kittiwakes in the present study is that all 3 variables might be related to males feeding at a higher trophic level than females during chick‐rearing (S. Tartu, personal communication). Thus, it is possible that a high‐quality dietary intake from high–trophic level prey would lead to a higher body condition (MacArthur and Pianka 1966), accompanied by exposure to greater levels of PFAS because PFAS biomagnify along the food chain (Haukås et al. 2007). Yet, this does not explain the negative associations between body condition with linPFOS and PFNA in female kittiwakes in the present study.

It is, however, necessary to acknowledge the shortcomings of using body condition as a proxy for lipid storage in seabirds (Jacobs et al. 2012), and thus the results presented herein should be considered exploratory. Given the increasing scrutiny of environmental pollutants as so‐called obesogens (Heindel and Blumberg 2019), this is a topic which should be examined in greater detail in future studies.

CONCLUSION

In kittiwakes, PFAS were detected, half of which were at higher concentrations in males compared to females; and there were associations between some PFAS and circulating thyroid hormones. It should be noted that the positive correlations between PFAS and thyroid hormones in kittiwakes do not allow for conclusions on causal relationships. Future studies should include determination of the free fraction of thyroid hormones in the plasma. The mechanisms of how PFAS exert their effects on the thyroid hormone system are probably indirect, and future research efforts should be directed toward elucidating this. In male kittiwakes, there were positive associations between 2 PFCAs and body condition, but in female kittiwakes, the opposite was observed, with negative associations between body condition and both linPFOS and PFNA. Given the contrasting associations of PFAS on the kittiwakes' body condition, this is identified as an area where more research is needed before any meaningful conclusions can be drawn.

Supplemental Data

The Supplemental Data are available on the Wiley Online Library at https://doi.org/10.1002/etc.4978.

Acknowledgment

O. Chastel, F. Angelier, and S. Tartu were supported by Institut Polaire Français (IPEV project 330, to O. Chastel), Agence National de la Recherche (ANR project PolarTop, to O. Chastel), and INSU‐EC2CO 2014 (to O. Chastel). A.V. Ask received the Arctic Field Grant from the Research Council of Norway for the field season of 2014. At the Norwegian Institute for Air Research, we thank D. Herzke, A. Götsch, and L. Hanssen for their facilitation and help with the PFAS analyses. At Centre d'Etudes Biologiques de Chizé, we thank C. Parenteau for her help in the T3 and T4 assays. We thank S. Nilsen for her help with the kittiwake sexing. We also thank T. Cielsiski and P. Massicotte for helpful discussions on the manuscript. S. Dielschneider is gratefully acknowledged for his help with the language. Finally, we thank 2 anonymous reviewers for their helpful comments.

Disclaimer

The authors declare no conflicts of interest.

Data Availability Statement

The data are available from the corresponding author upon request ([email protected]).

REFERENCES

Anker‐Nilssen
 
T
,
Barrett
 
RT
,
Bustnes
 
JO
,
Bårdsen
 
B‐J
,
Christensen‐Dalsgaard
 
S
,
Dehnhard
 
N
,
Descamps
 
S
,
Erikstad
 
KE
,
Fauchald
 
P
,
Hanssen
 
SA
,
Langset
 
M
,
Lorentsen
 
S‐H
,
Moe
 
B
,
Reiertsen
 
TK
,
Strøm
 
H
,
Systad
 
GH
.  
2019
. Sjøfugl i Norge 2018. SEAPOP, Tromsø, Norway.

Bech
 
C
,
Langseth
 
I
,
Moe
 
B
,
Fyhn
 
M
,
Gabrielsen
 
GW
.  
2002
.
The energy economy of the arctic‐breeding kittiwake (Rissa tridactyla): A review
.
Comp Biochem Physiol A Mol Integr Physiol
 
133
:
765
770
.

Blévin
 
P
,
Angelier
 
F
,
Tartu
 
S
,
Bustamante
 
P
,
Herzke
 
D
,
Moe
 
B
,
Bech
 
C
,
Gabrielsen
 
GW
,
Bustnes
 
JO
,
Chastel
 
O
.  
2017a
.
Perfluorinated substances and telomeres in an Arctic seabird: Cross‐sectional and longitudinal approaches
.
Environ Pollut
 
230
:
360
367
.

Blévin
 
P
,
Tartu
 
S
,
Ellis
 
HI
,
Chastel
 
O
,
Bustamante
 
P
,
Parenteau
 
C
,
Herzke
 
D
,
Angelier
 
F
,
Gabrielsen
 
GW
.  
2017b
.
Contaminants and energy expenditure in an Arctic seabird: Organochlorine pesticides and perfluoroalkyl substances are associated with metabolic rate in a contrasted manner
.
Environ Res
 
157
:
118
126
.

Bourgeon
 
S
,
Riemer
 
AK
,
Tartu
 
S
,
Aars
 
J
,
Polder
 
A
,
Jenssen
 
BM
,
Routti
 
H
.  
2017
.
Potentiation of ecological factors on the disruption of thyroid hormones by organo‐halogenated contaminants in female polar bears (Ursus maritimus) from the Barents Sea
.
Environ Res
 
158
:
94
104
.

Braune
 
BM
,
Trudeau
 
S
,
Jeffrey
 
DA
,
Mallory
 
ML
.  
2011
.
Biomarker responses associated with halogenated organic contaminants in northern fulmars (Fulmarus glacialis) breeding in the Canadian Arctic
.
Environ Pollut
 
159
:
2891
2898
.

Bustnes
 
JO
,
Borgå
 
K
,
Erikstad
 
KE
,
Lorentsen
 
S‐H
,
Herzke
 
D
.  
2008
.
Perfluorinated, brominated, and chlorinated contaminants in a population of lesser black‐backed gulls (Larus fuscus)
.
Environ Toxicol Chem
 
27
:
1383
1392
.

Butt
 
CM
,
Mabury
 
SA
,
Muir
 
DCG
,
Braune
 
BM
.  
2007
.
Prevalence of long‐chained perfluorinated carboxylates in seabirds from the Canadian Arctic between 1975 and 2004
.
Environ Sci Technol
 
41
:
3521
3528
.

Carravieri
 
A
,
Burthe
 
SJ
,
de la Vega
 
C
,
Yonehara
 
Y
,
Daunt
 
F
,
Newell
 
MA
,
Jeffreys
 
RM
,
Lawlor
 
AJ
,
Hunt
 
A
,
Shore
 
RF
,
Pereira
 
MG
,
Green
 
JA
.  
2020
.
Interactions between environmental contaminants and gastrointestinal parasites: Novel insights from an integrative approach in a marine predator
.
Environ Sci Technol
 
54
:
8938
8948
.

Chastel
 
O
,
Lacroix
 
A
,
Kersten
 
M
.  
2003
.
Pre‐breeding energy requirements: Thyroid hormone, metabolism and the timing of reproduction in house sparrows Passer domesticus
.
J Avian Biol
 
34
:
298
306
.

Cheng
 
J
,
Lv
 
S
,
Nie
 
S
,
Liu
 
J
,
Tong
 
S
,
Kang
 
N
,
Xiao
 
Y
,
Dong
 
Q
,
Huang
 
C
,
Yang
 
D
.  
2016
.
Chronic perfluorooctane sulfonate (PFOS) exposure induces hepatic steatosis in zebrafish
.
Aquat Toxicol
 
176
:
45
52
.

Costantini
 
D
,
Blévin
 
P
,
Herzke
 
D
,
Moe
 
B
,
Gabrielsen
 
GW
,
Bustnes
 
JO
,
Chastel
 
O
.  
2019
.
Higher plasma oxidative damage and lower plasma antioxidant defences in an Arctic seabird exposed to longer perfluoroalkyl acids
.
Environ Res
 
168
:
278
285
.

Drent
 
RH
,
Daan
 
S
.  
1980
.
The prudent parent: Energetic adjustments in avian breeding
.
Ardea
 
55
:
225
252
.

Erickson
 
RA
,
Rattner
 
BA
.  
2020
.
Moving beyond p < 0.05 in ecotoxicology: A guide for practitioners
.
Environ Toxicol Chem
 
39
:
1657
1669
.

Forsthuber
 
M
,
Kaiser
 
AM
,
Granitzer
 
S
,
Hassl
 
I
,
Hengstschläger
 
M
,
Stangl
 
H
,
Gundacker
 
C
.  
2020
.
Albumin is the major carrier protein for PFOS, PFOA, PFHxS, PFNA and PFDA in human plasma
.
Environ Int
 
137
:
105324
.

Griffiths
 
R
,
Double
 
MC
,
Orr
 
K
,
Dawson
 
RJG
.  
1998
.
A DNA test to sex most birds
.
Mol Ecol
 
7
:
1071
1075
.

Haarr
 
A
,
Hylland
 
K
,
Eckbo
 
N
,
Gabrielsen
 
GW
,
Herzke
 
D
,
Bustnes
 
JO
,
Blévin
 
P
,
Chastel
 
O
,
Moe
 
B
,
Hanssen
 
SA
,
Sagerup
 
K
,
Borgå
 
K
.  
2018
.
DNA damage in Arctic seabirds: Baseline, sensitivity to a genotoxic stressor, and association with organohalogen contaminants
.
Environ Toxicol Chem
 
37
:
1084
1091
.

Hanssen
 
L
,
Dudarev
 
AA
,
Huber
 
S
,
Odland
 
,
Nieboer
 
E
,
Sandanger
 
TM
.  
2013
.
Partition of perfluoroalkyl substances (PFASs) in whole blood and plasma, assessed in maternal and umbilical cord samples from inhabitants of arctic Russia and Uzbekistan
.
Sci Total Environ
 
447
:
430
437
.

Haugerud
 
AJ
.  
2011
. Levels and effects of organohalogenated contaminants on thyroid hormone levels in glaucous gulls (Larus hyperboreus) from Kongsfjorden, Svalbard. Master's thesis. Norwegian University of Science and Technology, Trondheim, Norway.

Haukås
 
M
,
Berger
 
U
,
Hop
 
H
,
Gulliksen
 
B
,
Gabrielsen
 
GW
.  
2007
.
Bioaccumulation of per‐ and polyfluorinated alkyl substances (PFAS) in selected species from the Barents Sea food web
.
Environ Pollut
 
148
:
360
371
.

Heindel
 
JJ
,
Blumberg
 
B
.  
2019
.
Environmental obesogens: Mechanisms and controversies
.
Annu Rev Pharmacol Toxicol
 
59
:
89
106
.

Helsel
 
DR
.  
2006
.
Fabricating data: How substituting values for nondetects can ruin results, and what can be done about it
.
Chemosphere
 
65
:
2434
2439
.

Hobson
 
KA
.  
1999
.
Tracing origins and migration of wildlife using stable isotopes: A review
.
Oecologia
 
120
:
314
326
.

Hovden
 
TS
.  
2018
. Thyroid hormone status and thyroid gland histology in PFAS‐ and mercury contaminated glaucous gulls (Larus hyperboreus) from Svalbard. Master's thesis. Norwegian University of Science and Technology, Trondheim, Norway.

Hulbert
 
AJ
.  
2000
.
Thyroid hormones and their effects: A new perspective
.
Biol Rev
 
75
:
519
631
.

Ishihara
 
A
,
Nishiyama
 
N
,
Sugiyama
 
S
,
Yamauchi
 
K
.  
2003
.
The effect of endocrine disrupting chemicals on thyroid hormone binding to Japanese quail transthyretin and thyroid hormone receptor
.
Gen Comp Endocrinol
 
134
:
36
43
.

Ismay
 
C
,
Kim
 
AY
.  
2020
.
Statistical Inference via Data Science: A Modern Dive into R and the Tidyverse
.
CRC
,
Boca Raton, FL, USA
.

Jacobs
 
SR
,
Elliott
 
K
,
Guigueno
 
MF
,
Gaston
 
AJ
,
Redman
 
P
,
Speakman
 
JR
,
Weber
 
J‐M
.  
2012
.
Determining seabird body condition using nonlethal measures
.
Physiol Biochem Zool
 
85
:
85
95
.

Jaspers
 
VLB
,
Herzke
 
D
,
Eulaers
 
I
,
Gillespie
 
BW
,
Eens
 
M
.  
2013
.
Perfluoroalkyl substances in soft tissues and tail feathers of Belgian barn owls (Tyto alba) using statistical methods for left‐censored data to handle non‐detects
.
Environ Int
 
52
:
9
16
.

Jolliffe
 
IT
.  
2002
.
Principal Component Analysis
, 2nd ed.  
Springer Series in Statistics. Springer‐Verlag
,
New York, NY, USA
.

Kim
 
B
.  
2008
.
Thyroid hormone as a determinant of energy expenditure and the basal metabolic rate
.
Thyroid
 
18
:
141
144
.

Leat
 
EHK
,
Bourgeon
 
S
,
Eze
 
JI
,
Muir
 
DCG
,
Williamson
 
M
,
Bustnes
 
JO
,
Furness
 
RW
,
Borgå
 
K
.  
2013
.
Perfluoroalkyl substances in eggs and plasma of an avian top predator, great skua (Stercorarius skua), in the north Atlantic
.
Environ Toxicol Chem
 
32
:
569
576
.

Letcher
 
RJ
,
Bustnes
 
JO
,
Dietz
 
R
,
Jenssen
 
BM
,
Jørgensen
 
EH
,
Sonne
 
C
,
Verreault
 
J
,
Vijayan
 
MM
,
Gabrielsen
 
GW
.  
2010
.
Exposure and effects assessment of persistent organohalogen contaminants in arctic wildlife and fish
.
Sci Total Environ
 
408
:
2995
3043
.

Link
 
JC
,
Reue
 
K
.  
2017
.
Genetic basis for sex differences in obesity and lipid metabolism
.
Annu Rev Nutr
 
37
:
225
245
.

Liu
 
Y
,
Wang
 
J
,
Fang
 
X
,
Zhang
 
H
,
Dai
 
J
.  
2011
.
The thyroid‐disrupting effects of long‐term perfluorononanoate exposure on zebrafish (Danio rerio)
.
Ecotoxicology
 
20
:
47
55
.

MacArthur
 
RH
,
Pianka
 
ER
.  
1966
.
On optimal use of a patchy environment
.
Am Nat
 
100
:
603
609
.

Mauvais‐Jarvis
 
F
.  
2015
.
Sex differences in metabolic homeostasis, diabetes, and obesity
.
Biology of Sex Differences
 
6
:
14
.

McKechnie
 
AE
.  
2008
.
Phenotypic flexibility in basal metabolic rate and the changing view of avian physiological diversity: A review
.
J Comp Physiol B
 
178
:
235
247
.

McNabb
 
FMA
.  
2007
.
The hypothalamic–pituitary–thyroid (HPT) axis in birds and its role in bird development and reproduction
.
Crit Rev Toxicol
 
37
:
163
193
.

McNabb
 
FMA
,
Hughes
 
TE
.  
1983
.
The role of serum binding proteins in determining free thyroid hormone concentrations during development in quail
.
Endocrinology
 
113
:
957
963
.

Melnes
 
M
,
Gabrielsen
 
GW
,
Herzke
 
D
,
Sagerup
 
K
,
Jenssen
 
BM
.  
2017
.
Dissimilar effects of organohalogenated compounds on thyroid hormones in glaucous gulls
.
Environ Res
 
158
:
350
357
.

Merryman
 
JI
,
Buckles
 
EL
.  
1998
.
The avian thyroid gland
. Part two: A review of function and pathophysiology.
J Avian Med Surg
 
12
:
238
242
.

Miller
 
A
,
Elliott
 
JE
,
Elliott
 
KH
,
Lee
 
S
,
Cyr
 
F
.  
2015
.
Temporal trends of perfluoroalkyl substances (PFAS) in eggs of coastal and offshore birds: Increasing PFAS levels associated with offshore bird species breeding on the Pacific coast of Canada and wintering near Asia
.
Environ Toxicol Chem
 
34
:
1799
1808
.

Mortensen
 
Å‐K
,
Mæhre
 
S
,
Kristiansen
 
K
,
Heimstad
 
ES
,
Gabrielsen
 
GW
,
Jenssen
 
BM
,
Sylte
 
I
.  
2020
.
Homology modeling to screen for potential binding of contaminants to thyroid hormone receptor and transthyretin in glaucous gull (Larus hyperboreus) and herring gull (Larus argentatus)
.
Comput Toxicol
 
13
:
100120
.

Newsted
 
JL
,
Jones
 
PD
,
Coady
 
K
,
Giesy
 
JP
.  
2005
.
Avian toxicity reference values for perfluorooctane sulfonate
.
Environ Sci Technol
 
39
:
9357
9362
.

Nilsson
 
H
,
Kärrman
 
A
,
Rotander
 
A
,
van Bavel
 
B
,
Lindström
 
G
,
Westberg
 
H
.  
2013
.
Biotransformation of fluorotelomer compound to perfluorocarboxylates in humans
.
Environ Int
 
51
:
8
12
.

Nishiwaki‐Ohkawa
 
T
,
Yoshimura
 
T
.  
2016
.
Molecular basis for regulating seasonal reproduction in vertebrates
.
J Endocrinol
 
229
:
R117
R127
.

Nøst
 
TH
,
Helgason
 
LB
,
Harju
 
M
,
Heimstad
 
ES
,
Gabrielsen
 
GW
,
Jenssen
 
BM
.  
2012
.
Halogenated organic contaminants and their correlations with circulating thyroid hormones in developing Arctic seabirds
.
Sci Total Environ
 
414
:
248
256
.

Palmer
 
BF
,
Clegg
 
DJ
.  
2015
.
The sexual dimorphism of obesity
.
Mol Cell Endocrinol
 
402
:
113
119
.

R Development Core Team
.
2019
.
R: A Language and Environment for Statistical Computing
.
R Foundation for Statistical Computing
.
Vienna, Austria
.

Rauert
 
C
,
Shoieb
 
M
,
Schuster
 
JK
,
Eng
 
A
,
Harner
 
T
.  
2018
.
Atmospheric concentrations and trends of poly‐ and perfluoroalkyl substances (PFAS) and volatile methyl siloxanes (VMS) over 7 years of sampling in the Global Atmospheric Passive Sampling (GAPS) network
.
Environ Pollut
 
238
:
94
102
.

Ruuskanen
 
S
,
Hsu
 
B‐Y
,
Nord
 
A
.  
2021
.
Endocrinology of thermoregulation in birds in a changing climate
.
Mol Cell Endocrinol
 
519
:
111088
.

Shi
 
X
,
Liu
 
C
,
Wu
 
G
,
Zhou
 
B
.  
2009
.
Waterborne exposure to PFOS causes disruption of the hypothalamus–pituitary–thyroid axis in zebrafish larvae
.
Chemosphere
 
77
:
1010
1018
.

Sirsat
 
TS
,
Dzialowski
 
EM
.  
2020
.
Manipulating plasma thyroid hormone levels at hatching alters development of endothermy and ventilation in Pekin duck (Anas platyrhynchos domestica)
.
J Exp Biol
 
223
:
jeb237701
.

Spano
 
JS
,
Pedersoli
 
WM
,
Kemppainen
 
RJ
,
Krista
 
LM
,
Young
 
DW
.  
1987
.
Baseline hematologic, endocrine, and clinical chemistry values in ducks and roosters
.
Avian Dis
 
31
:
800
803
.

Tartu
 
S
,
Gabrielsen
 
GW
,
Blévin
 
P
,
Ellis
 
H
,
Bustnes
 
JO
,
Herzke
 
D
,
Chastel
 
O
.  
2014
.
Endocrine and fitness correlates of long‐chain perfluorinated carboxylates exposure in Arctic breeding black‐legged kittiwakes
.
Environ Sci Technol
 
48
:
13504
13510
.

Tartu
 
S
,
Lille‐Langøy
 
R
,
Størseth
 
TR
,
Bourgeon
 
S
,
Brunsvik
 
A
,
Aars
 
J
,
Goksøyr
 
A
,
Jenssen
 
BM
,
Polder
 
A
,
Thiemann
 
GW
,
Torget
 
V
,
Routti
 
H
.  
2017
.
Multiple‐stressor effects in an apex predator: Combined influence of pollutants and sea ice decline on lipid metabolism in polar bears
.
Sci Rep
 
7
:
16487
.

Ucán‐Marin
 
F
,
Arukwe
 
A
,
Mortensen
 
AS
,
Gabrielsen
 
GW
,
Letcher
 
RJ
.  
2010
.
Recombinant albumin and transthyretin transport proteins from two gull species and human: Chlorinated and brominated contaminant binding and thyroid hormones
.
Environ Sci Technol
 
44
:
497
504
.

Verreault
 
J
,
Houde
 
M
,
Gabrielsen
 
GW
,
Berger
 
U
,
Haukås
 
M
,
Letcher
 
RJ
,
Muir
 
DCG
.  
2005
.
Perfluorinated alkyl substances in plasma, liver, brain, and eggs of glaucous gulls (Larus hyperboreus) from the Norwegian Arctic
.
Environ Sci Technol
 
39
:
7439
7445
.

Welcker
 
J
,
Chastel
 
O
,
Gabrielsen
 
GW
,
Guillaumin
 
J
,
Kitaysky
 
AS
,
Speakman
 
JR
,
Tremblay
 
Y
,
Bech
 
C
.  
2013
.
Thyroid hormones correlate with basal metabolic rate but not field metabolic rate in a wild bird species
.
PLoS One
 
8
:
e56229
.

Young
 
CJ
,
Furdui
 
VI
,
Franklin
 
J
,
Koerner
 
RM
,
Muir
 
DCG
,
Mabury
 
SA
.  
2007
.
Perfluorinated acids in Arctic snow:  New evidence for atmospheric formation
.
Environ Sci Technol
 
41
:
3455
3461
.

Zhang
 
S
,
Guo
 
X
,
Lu
 
S
,
Sang
 
N
,
Li
 
G
,
Xie
 
P
,
Liu
 
C
,
Zhang
 
L
,
Xing
 
Y
.  
2018
.
Exposure to PFDoA causes disruption of the hypothalamus–pituitary–thyroid axis in zebrafish larvae
.
Environ Pollut
 
235
:
974
982
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The moral rights of the named author(s) have been asserted.