Associations between thyroid hormones and serum metabolite profiles in an euthyroid population.

The aim was to characterise associations between circulating thyroid hormones-free thyroxine (FT4) and thyrotropin (TSH)-and the metabolite profiles in serum samples from participants of the German population-based KORA F4 study. Analyses were based on the metabolite profile of 1463 euthyroid subjects. In serum samples, obtained after overnight fasting (≥8), 151 different metabolites were quantified in a targeted approach including amino acids, acylcarnitines (ACs), and phosphatidylcholines (PCs). Associations between metabolites and thyroid hormone concentrations were analysed using adjusted linear regression models. To draw conclusions on thyroid hormone related pathways, intra-class metabolite ratios were additionally explored. We discovered 154 significant associations (Bonferroni p < 1.75 × 10-04) between FT4 and various metabolites and metabolite ratios belonging to AC and PC groups. Significant associations with TSH were lacking. High FT4 levels were associated with increased concentrations of many ACs and various sums of ACs of different chain length, and the ratio of C2 by C0. The inverse associations observed between FT4 and many serum PCs reflected the general decrease in PC concentrations. Similar results were found in subgroup analyses, e.g., in weight-stable subjects or in obese subjects. Further, results were independent of different parameters for liver or kidney function, or inflammation, which supports the notion of an independent FT4 effect. In fasting euthyroid adults, higher serum FT4 levels are associated with increased serum AC concentrations and an increased ratio of C2 by C0 which is indicative of an overall enhanced fatty acyl mitochondrial transport and β-oxidation of fatty acids.

Introduction

Thyroid hormones such as thyroxine (T4) and triiodothyronine (T3) are responsible for the proper function of nearly all organ systems (Brix et al. 2011). They do not only modulate the metabolism of carbohydrates and lipids, influence protein synthesis, affect the cellular reaction to catecholamines in almost every tissue; they are also unique in their ability to influence all over energy expenditure, e.g., via the regulation of thermogenesis (Kim 2008; Panicker et al. 2010; Moreno et al. 2008). The release of T4 and T3 by the thyroid gland is controlled by thyroid stimulating hormone (thyrotropin; TSH) which is synthesized by the anterior pituitary gland in response to TSH-releasing hormone, which in turn is secreted by the hypothalamus. T4 is the quantitatively main product of the thyroid gland (Yen 2001). Almost all of T3 and T4 in the circulation are protein-bound (Boelaert and Franklyn 2005); only about 0.03 and 0.3 % of total serum T4 and T3, respectively, are free or unbound (Yen 2001). The free thyroid hormones (FT4 and FT3) actively enter target cells and trigger biological reactions (Yen 2001). Perturbations of thyroid hormone release and functioning belong to the most common endocrine disorders and are connected to major health problems (Brix et al. 2011). A prolonged imbalance of thyroid hormones leads to diseases, such as hyperthyroidism or hypothyroidism, and has been linked to an increased risk for cardiovascular morbidity and mortality (Brix et al. 2011; Boelaert and Franklyn 2005; Biondi 2012b). The circulating concentrations of TSH, FT4 and FT3 have a strong heritable component and are supposed to be controlled by multiple genes (Panicker et al. 2010). However, in the case of FT4 these genes are still unknown.

In the present work we investigate the association between FT4 concentrations and serum metabolite profiles. The concept of targeted metabolomics is the quantification of a defined set of metabolites in a body fluid, representing an image of the current metabolic state of the organism (Altmaier et al. 2009). It has previously been shown that this method has the power to identify perturbations of the body’s metabolic homeostasis and allows for the identification of and access to biomarkers of metabolic pathways that are impacted, for example, by diseases (Suhre et al. 2010; Altmaier et al. 2008; Vinayavekhin and Saghatelian 2009; Watson 2006; Dumas et al. 2007; Assfalg et al. 2008).

The focus of the present study was on euthyroid participants, i.e., participants with FT3, FT4, and TSH levels in the normal range, no history of known thyroid diseases, and no intake of thyroid-specific medication. We hypothesised that within the normal ranges of FT4 and TSH an effect on substrate oxidation exists, and that the activity of involved enzymes is reflected in serum metabolite concentrations. Therefore, we took a targeted quantitative metabolomics approach to identify FT4 and TSH triggered pathways and differences in varying population strata by, e.g., gender, obesity status, age, or physical activity.

Materials and methods

Ethics statement

The study was approved by the ethics committee of the Bavarian Medical Association. Written informed consent was obtained from each participant in accordance with institutional requirements and the Declaration of Helsinki Principles.

Study design and population

The present analyses are based on subsamples of the metabolically characterised subjects who participated in the KORA F4 study, the seven year follow-up of the KORA S4 study. KORA (Cooperative Health Research in the Region of Augsburg) is a research platform performing population-based surveys and subsequent follow-ups in the region of Augsburg in Southern Germany and has been described in detail elsewhere (Holle et al. 2005). The metabolic profiles of 3044 participants in KORA F4, aging 31–82 years, were determined. For the present analyses only euthyroid participants were included. A person was defined as euthyroid if their FT3, FT4 and TSH hormone levels were in the ranges of 3.3–6.1 pmol/L, 9.8–18.8 pmol/L and 0.3–4.0 mIU/L according to the assay manuals and Brabant et al. (2006) as well as Surks et al. (2004), respectively, their thyroid peroxidase antibody level did not exceed 200 IU/mL (Demers and Spencer 2003; Knudsen et al. 1999) and there was no known history of thyroid disease or thyroid medication use in accordance with the anatomical therapeutic chemical code for thyroid therapy (H03). Further, subjects with a known history of myocardial infarction, stroke, diabetes or cancer were excluded as well as subjects taking ACE inhibitors or anti-lipidemic drugs, and women being pregnant at the time of blood collection. The selection process of the study participants is displayed in Fig. 1. Participants were excluded for several reasons if applicable. For a subanalysis, data on the participants’ body composition was used. However, the bioelectrical impedance analysis measurements were only conducted in the KORA S4 study which was seven years prior to the KORA F4 examinations. To account for this, only weight-stable participants of the KORA F4 population were included in this subanalysis. Persons were defined as weight-stable if their weight gain or loss did not exceed more than 0.5 % per year between measurements in KORA S4 and F4 (Nimptsch et al. 2010). This resulted in a sample size of n = 1463 for KORA F4 and of n = 616 for the weight-stable subgroup in KORA F4.

Fig. 1

Flowchart of the selection process of the study participants for the present analyses. (Participants were excluded for multiple reasons.)

Blood sampling

Blood samples were drawn between 8:00 and 10:30 a.m. after a period of overnight fasting (minimum of 8 h). S-Monovette® serum gel tubes (SARSTEDT AG & Co., Nümbrecht, Germany) were used. After blood was drawn, tubes were inverted twice and then allowed to rest for 30 min at room temperature (18–25 °C) to achieve total coagulation. Tubes were centrifuged for 10 min at 2,660 g at 15 °C and the obtained serum was divided into aliquots, chilled for a maximum of 6 h at 4 °C before being stored at −80 °C until analysis.

Blood analyses

Concentrations of TSH, FT3, and FT4 were analysed by immunochemiluminescent procedures (Dimension Vista System, Siemens, Germany). The functional sensitivity of the TSH assay was 0.005 mIU/L; the TSH, FT3, and FT4 working ranges were 0.005–100 mIU/L, 0.77–46.1 pmol/L, and 1.3–103 pmol/L, respectively. The (low or high) inter-assay coefficients of variations were 2.0 or 2.2 % for TSH, and 1.7 or 4.1 % for FT4. All assays were performed according to the manufacturer’s recommendations. However, the applied analytic method for FT3 did not fulfil the strict quality control criteria in the end (Thienpont et al. 2010). Thus, data on FT3 was only used for the definition of euthyroid subjects but not for analyses. The thyroid peroxidase antibodies were determined by an enzyme immunoassay (VARELISA, Elias Medizintechnik GmbH, Freiburg, Germany) with a functional sensitivity of 1 IU/mL.

The applied analyses methods to determine glutamate pyruvate transaminase, glutamic oxaloacetic transaminase, gamma-glutamyl transpeptidase, C-reactive protein (Ruckert et al. 2011), cystatin C (Goek et al. 2012), and creatinine (Seissler et al. 2012), concentrations have been described earlier. White blood cell count was determined from fresh venous EDTA blood using impedance measurements (Coulter® LH 780, Beckman Coulter Inc., USA).

Serum metabolite quantification

Metabolic characterisation of the KORA F4 serum samples was done in three batches of approximately 1000 samples at three different time points with a recalibration of the equipment in between.

The targeted metabolomics approach for KORA F4 was based on ESI–MS/MS measurements by AbsoluteIDQ™ p150 kit (BIOCRATES, Life Sciences AG, Innsbruck, Austria). The assay allows the simultaneous quantification of 163 metabolites out of 10 μl serum. The AbsoluteIDQ™ p150 kit has been described in detail previously (Illig et al. 2010; Römisch-Margl et al. 2011) Sample handling was performed by a Hamilton Micro Lab Star robot (Hamilton Bonaduz AG, Bonaduz, Switzerland) and a Ultravap nitrogen evaporator (Porvair, Sciences, Leatherhead, United Kingdom). Mass spectrometric (MS) analyses were done on a 4000 QTRAP mass spectrometer (AB Sciex) coupled to Promincence HPLC (Shimadzu) apparatus equipped with a SIL-20AC autosampler and controlled by the software Analyst 1.4. Data evaluation for quantification of metabolite concentrations and quality assessment was performed with the MetIQ™ software package, which is an integral part of the AbsoluteIDQ™ kit. Internal standards serve as reference for the calculation of metabolite concentrations. The methods of the AbsoluteIDQ™ p150 kits has been proven to be in conformance with the FDA guideline “Guidance for Industry—Bioanalytical Method Validation, May 2001” (U.S. Department of Health and Human Services et al. May 2001), which implies proof of reproducibility within a given error range. Measurements were performed as described in the manufacturer manuals. Concentrations are reported in μmol/L.

Metabolite spectrum

In total, 163 different metabolites were detected, including 14 amino acids (13 proteinogenic + ornithine), hexose (sum of hexoses—about 90–95 % glucose), free carnitine and 40 acylcarnitines (Cx:y), 15 sphingomyelins (SMx:y), 77 phosphatidylcholines (PCs, aa = diacyl, ae = acyl–alkyl) and 15 lyso-phosphatidylcholines (lysoPCs). The lipid side chain composition is abbreviated as Cx:y, with x denoting the number of carbons in the side chain and y denoting the number of double-bonds. Full biochemical names are provided in the Supplementary Table S1.

Quality control of the metabolomics dataset

The quality control of the metabolomics dataset of KORA F4 was done in a two-step procedure. First, the quality of all metabolites was assessed using a reference blood which was measured five times on ten plates. With this data, a coefficient of variation was calculated for every metabolite and plate. All metabolites having a mean coefficient of variation over all ten plates bigger than 25 % were removed from the dataset (11 in total). One further metabolite was excluded as the number of missing values exceeded 5 %. In the second step, the dataset was controlled for outliers. A subject’s metabolite concentration was defined as an outlier if the concentration was greater or less than the mean ± five standard deviations of the particular metabolite over the whole population. All subjects having more than three independent outlying metabolite concentrations were excluded from the dataset. An outlier was defined as independent if the correlation with all other outliers was less than 70 %. If there were three or less independent outliers, only the data points were removed. All missing values (<0.09 %) were imputed with the R package “mice” (van Buuren and Groothuis-Oudshoorn 2011) which uses a linear regression approach (predictive mean matching; vector size = 151). This left us with 151 metabolites.

Anthropometry, physical activity and body composition assessment

Height and weight were measured to the nearest 0.1 cm and 0.1 kg, respectively. The body mass index (BMI) was calculated as weight in kg/(height in m)2 and according to the world health organisation (WHO) obesity was defined as having a BMI ≥ 30 kg/m2 (WHO 2000). Waist circumference was measured to the closest 0.1 cm at the smallest position between the lower rip and the upper margin of the iliac crest. Hip size was determined exactly to 0.1 cm as the widest circumference measured between the upper margin of the iliac crest and the crotch. Physical activity was assessed on a four-level graded scale by the amount of regular leisure time exercise per week during summer and wintertime. Based on those assessments, Meisinger et al. (2007) defined four levels of physical activity: (i) ‘No sports activities in leisure time’—almost no sports activity or no activity in summer and in winter; (ii) ‘Low level of sports activities in leisure time’—irregular exercise of 1 h per week at least in summer or winter; (iii) ‘Moderate level of sports activity in leisure time’—regular exercise of 1 h per week at least in summer or winter; (iv) ‘High level of sports activities in leisure time’—more than 2 h per week of regular exercise in summer and winter. Based on this variable, the physical activity variable which was used in the present analyses was created representing two levels: ‘physically inactive’—category (i) or (ii), and ‘physically active’—category (iii) or (iv). For the assessment of body composition in KORA S4, two bioelectrical impedance analysis measurements of resistance (R), reactance (Xc) and the phase angle (α) were taken between the dominant hand wrist and dorsum and the dominant foot angle and dorsum (placement of the electrodes) by means of a body impedance analyser (BIA 2000-S; Data Input GmbH, Frankfurt, Germany) while subjects were spreading their arms and legs and lying in a relaxed and supine position on a nonconductive surface with 50 kHz. Fat free mass was calculated by means of Kyle’s equation (Kyle et al. 2001). In accordance to the BMI calculation, a fat free mass index (FFMI) (fat free mass in kg/(height in m)2) was determined.

Statistical analyses

The descriptive data is presented as mean and standard deviation for the continuous variables and as absolute quantities and percentages for the discrete parameters. Differences in hormone concentrations across different categories were tested by means of pairwise t-tests. Besides the absolute metabolite concentrations, all pairs of intra-metabolite class ratios (n = 2188) were part of the metabolomics dataset. Associations between metabolite concentrations or metabolite ratios and FT4 or TSH hormone levels were assessed by means of linear models. First, metabolite variables were standardised with an inverse log-normal transformation to allow for comparison of the estimates derived from the linear regression models. Then, linear regression models were applied with the metabolite variables as dependent variables and the hormone concentrations as independent variable. The linear models were adjusted for age, sex, BMI, and batch. To control for multiple testing and with regard to the dependencies among the various metabolites and ratios, a p value of 1.75 × 10−4 was considered statistically significant at α = 5 %. For an association with a metabolite ratio to be regarded as statistically significant, a further criterion (p-gain) had to be fulfilled next to a significant p value. The p-gain is defined as the fold decrease in the p value of the association for the pair of metabolites compared to the lowest of the two p values for the single metabolites (Suhre et al. 2010). Thus, an association between a metabolite ratio (M1/M2) and FT4 or TSH is considered to be significant, if the p value of this association is significant and the p-gain exceeds a certain value. To calculate the p-gain we first have to determine the minimum p value of the associations between metabolite M1 and FT4 or TSH as well as metabolite M2 and FT4 or TSH. The p-gain is the quotient of this minimum to the p value of the association between the metabolite ratio and FT4 or TSH. This p-gain had to exceed 2.86 × 10+03. This cut-off is seen as Bonferroni-related correction in order to identify metabolite concentration pairs for which the strength of association improves considerably by using ratios (Suhre et al. 2010; Jourdan et al. 2012).

Subanalyses

In a subanalysis, the effect of FFMI on the association between FT4 and the metabolomics data was examined. Therefore, the weight-stable set was used and the applied linear models were adjusted for FFMI instead of BMI. Stratified analyses of the entire sample were also conducted by sex, obesity status (BMI ≥ 30 kg/m2 versus BMI < 30 kg/m2), age groups (≤60 years versus >60 years), and physical activity as potential confounder or effect modifiers. The equality of the strata-specific regression coefficients was tested by means of an approximately normally distributed test-statistic as described by Paternoster et al. (Paternoster et al. 1998). Further, the influence of serum insulin, different liver, kidney, and inflammation parameters on the identified associations between the metabolites and FT4 were analysed by including an interaction term between FT4 and those parameters in the described main model.

All statistical analyses were performed with the R software version 2.12.0 [R Development Core Team, 2010, http://www.r-project.org].

Results

Table 1 summarises the characteristics of all study participants and the weight-stable subset. On average, men and women were 51 ± 12 years old and had a BMI of 27 ± 4 and 26 ± 5 kg/m², respectively; 19–20 % of men and women were classified as obese. Mean serum TSH and FT4 concentrations were 1.40 ± 0.65 mIU/L and 13.45 ± 1.67 pmol/L in men and 1.50 ± 0.71 mIU/L and 13.55 ± 1.62 pmol/L in women. Comparing serum FT4 and TSH hormone concentrations across different categories of age, sex, BMI, and physical activity for all KORA F4 participants and the two subsets, weight-stable and obese participants, we hardly detected any differences (Table 2). Applying ANOVA and pairwise t-tests, mean FT4 values of all KORA F4 participants for the age strata 61–70 and 71–80 years were significantly different to the other strata. Also the mean FT4 value in the BMI category <18.5 consisting of 4 subjects was significantly higher as compared to the other BMI categories. No differences among categories were found for the TSH concentrations (Table 2).

Table 1

Characteristics (Mean ± SD; % or n absolute) of male and female subjects in the KORA F4 sample and the weight-stable subset

		All (n = 1463)		Weight-stablea (n = 616)
		Men (n = 792)	Women (n = 671)	Men (n = 343)	Women (n = 273)
Parameter	Unit	Mean ± SD	Mean ± SD	Mean ± SD	Mean ± SD
Age	years	50.94 ± 12.08	50.98 ± 12.39	53.34 ± 12.06	54.03 ± 12.7
Weight	kg	85.31 ± 13.02	69.76 ± 13.13	83.29 ± 11.50	67.80 ± 12.56
Height	cm	177.11 ± 6.85	162.93 ± 6.68	176.66 ± 7.02	162.48 ± 6.44
Body-mass-index	kg/m2	27.19 ± 3.76	26.32 ± 4.96	26.68 ± 3.31	25.7 ± 4.63
Waist	cm	97.19 ± 11.35	85.23 ± 12.24	95.89 ± 10.09	83.94 ± 12.17
Hip	cm	105.61 ± 8.06	104.86 ± 9.92	104.31 ± 6.26	103.46 ± 9.43
Waist-to-hip ratio		0.92 ± 0.06	0.81 ± 0.07	0.92 ± 0.06	0.81 ± 0.07
FT4	pmol/L	13.45 ± 1.67	13.55 ± 1.62	13.43 ± 1.62	13.57 ± 1.61
TSH	mIU/L	1.40 ± 0.65	1.50 ± 0.71	1.37 ± 0.63	1.47 ± 0.69
Fat free massb	kg	60.85 ± 6.12	43.28 ± 5.21	60.83 ± 6.10	43.13 ± 5.19
Body fat massb	%	26.29 ± 4.99	35.62 ± 5.72	26.41 ± 4.94	35.54 ± 5.62
FFMIb,c	kg/m2	19.55 ± 1.64	16.39 ± 1.84	19.61 ± 1.58	16.40 ± 1.78
BFMIb,d	kg/m2	7.17 ± 2.24	9.44 ± 3.24	7.21 ± 2.08	9.39 ± 3.14
Insulin	μIU/ml	6.87 ± 12.04	5.63 ± 9.58	6.21 ± 9.07	5.22 ± 11.22
Parameter	Category	% (n absolute)	% (n absolute)	% (n absolute)	% (n absolute)
Body-mass-index (kg/m²)	<18.5	0.12 (1)	0.45 (3)	0.29 (1)	0.37 (1)
	18.5 to <25	28.28 (224)	47.84 (321)	30.61 (105)	51.65 (141)
	25 to <30	52.53 (416)	31.74 (213)	51.60 (177)	32.23 (88)
	≥30	19.07 (151)	19.82 (133)	17.49 (60)	15.75 (43)
Physical activitye	Active	57.95 (459)	68.41 (375)	60.35 (207)	60.07 (164)

aWeight-stable subjects of KORA F4, whose weight gain or loss did not exceed more than 0.5 % per year since their body weight was measured in KORA S4 (on average seven years prior F4 examination)

bParameters derived from the BIA-measurements which were only conducted in KORA S4

cFat free mass index

dBody fat mass index

ePhysical Activity (active: >1 h sport per week on a regular basis; inactive: <1 h sport per week)

Table 2

Mean (±SD) of FT4 hormone concentration (pmol/l) and TSH hormone concentration (mIU/L) across different strata of age, BMI (kg/m2), or physical activity for the KORA F4 participants (n = 1463), the weight-stable subset (n = 616) and the obese participants (BMI ≥ 30 kg/m2; n = 284)

		KORA F4 (n = 1463)				Weight-stablea (n = 616)				Obese subset (n = 284)
		FT4		TSH		FT4		TSH		FT4		TSH
Parameter	Category	Mean ± SD		Mean ± SD		Mean ± SD		Mean ± SD		Mean ± SD		Mean ± SD
Age (years)	31–40	13.34 ± 1.59		1.54 ± 0.67		13.43 ± 1.59		1.52 ± 0.54		13.00 ± 1.54		1.44 ± 0.66
	41–50	13.23 ± 1.61		1.47 ± 0.67		13.10 ± 1.65		1.46 ± 0.68		13.35 ± 1.65		1.54 ± 0.58
	51–60	13.41 ± 1.54		1.40 ± 0.70		13.34 ± 1.36		1.42 ± 0.71		13.34 ± 1.45		1.48 ± 0.61
	61–70	13.88 ± 1.68	*	1.41 ± 0.64		13.80 ± 1.60	*of 2	1.34 ± 0.62		13.93 ± 1.71		1.58 ± 0.68
	71–80	14.45 ± 1.72	*	1.33 ± 0.70		14.43 ± 1.70	*of 1-3	1.25 ± 0.66		14.49 ± 1.85	*of 1-3	1.48 ± 0.80
	≥81	14.22 ± 2.05		0.91 ± 0.59		14.22 ± 2.05		0.91 ± 0.59		15.45 ± 3.18		0.43 ± 0.01
BMI	<18.5	16.65 ± 1.04	*	1.23 ± 0.59		16.90 ± 1.41	*	1.23 ± 0.89
(kg/m2)	18.5 to <25	13.50 ± 1.62		1.46 ± 0.73		13.49 ± 1.66		1.41 ± 0.71
	25 to <30	13.44 ± 1.63		1.41 ± 0.65		13.49 ± 1.50		1.36 ± 0.59
	≥30	13.55 ± 1.69		1.50 ± 0.65		13.46 ± 1.73		1.54 ± 0.69
Physical activity	Active	13.45 ± 1.57	#	1.44 ± 0.69		13.55 ± 1.55	#	1.38 ± 0.65		13.51 ± 1.55		1.49 ± 0.61
	Inactive	13.55 ± 1.73		1.45 ± 0.66		13.41 ± 1.71		1.46 ± 0.67		13.59 ± 1.80		1.51 ± 0.68

aWeight-stable subjects of KORA F4, whose weight gain or loss did not exceed more than 0.5 % per year since their body weight was measured in KORA S4 (on average seven years prior F4 examination)

bPhysical Activity (active: >1 h sport per week on a regular basis; inactive: less than 1 h sport per week)

*pairwise t test was significant

#levene test resulted in rejection of equal variances

After correction for multiple testing we obtained 154 statistically significant associations between various serum metabolite concentrations, metabolite ratios, and FT4. Significant associations between metabolites and TSH concentrations were entirely lacking in the main analysis as well as in the subanalyses. Selected statistically significant associations of FT4 with different metabolite concentrations for the KORA F4 population are shown in Table 3. Different parameters such as the β-estimate derived from the adjusted linear model are given, as well as the adjusted p value, the adjusted R2 and the p-gain. The full list of all obtained statistically significant associations can be found in Supplementary Table S2. The results of the subanalyses are given in the Supplementary Tables S3-S7.

Table 3

Selected metabolic traits significantly associated with FT4 in a linear regression model adjusted for age, sex, BMI, and batch in the KORA F4 sample

Trait	Mean (μmol/l)	SD	Beta	Adj.P valuea	R2 adjb	P-Gainc
Acylcarnitines
Σ AC	10.81	3.24	0.11	1.23 × 10−10	0.14
Σ long AC	0.40	0.09	0.10	1.25 × 10−10	0.21
Σ even AC	9.97	3.16	0.11	1.29 × 10−10	0.14
Σ poly AC	8.32	2.85	0.11	2.41 × 10−10	0.13
Σ mono AC	1.73	0.46	0.10	2.12 × 10−08	0.16
Σ med. AC	1.12	0.36	0.09	3.00 × 10−06	0.14
C18:1	0.13	0.04	0.12	7.04 × 10−14	0.18
C16:1	0.04	0.01	0.11	6.15 × 10−11	0.17
C14:2	0.03	0.01	0.11	3.17 × 10−10	0.14
C2	7.85	2.77	0.11	4.46 × 10−10	0.13
C18:2	0.05	0.01	0.09	6.44 × 10−07	0.18
C12:1	0.14	0.05	0.09	1.32 × 10−06	0.11
C12	0.13	0.05	0.09	1.81 × 10−06	0.12
C10	0.34	0.16	0.08	3.66 × 10−05	0.09
C10:1	0.16	0.06	0.07	2.37 × 10−04	0.13
C14:1	0.14	0.05	0.06	1.69 × 10−02	0.25
C8	0.21	0.09	0.06	2.37 × 10−02	0.10
(C2 + C3)/C0	0.24	0.08	0.11	1.72 × 10−10	0.14
C2/C0	0.23	0.08	0.11	3.12 × 10−10	0.15
Σ PCs	1984.61	355.11	−0.11	3.55 × 10−11	0.14
Σ mono PCs	348.02	86.20	−0.11	8.18 × 10−11	0.07
Σ poly PCs	1592.39	279.37	−0.11	9.87 × 10−10	0.15
Σ sat. PCs	44.21	7.90	−0.09	2.61 × 10−06	0.12
Σ short PCs	759.09	161.58	−0.11	3.37 × 10−11	0.25
Σ PC aa	1786.89	333.22	−0.12	8.01 × 10−12	0.14
Σ short PC aa	720.81	156.95	−0.11	2.82 × 10−11	0.25
Σ mono PC aa	323.42	83.04	−0.11	9.75 × 10−11	0.07
Σ poly PC aa	1435.77	261.53	−0.11	2.28 × 10−10	0.15
Σ sat. PC aa	27.71	5.60	−0.09	9.30 × 10−06	0.11
Σ long PC aa	47.11	12.34	−0.06	1.08 × 10−02	0.07
PC aa C34:4	2.26	0.83	−0.14	2.02 × 10−16	0.07
PC aa C36:1	53.43	13.58	−0.13	7.41 × 10−15	0.10
PC aa C36:6/PC aa C42:0	1.95	0.80	−0.14	7.46 × 10−16	0.06	7.90 × 10+05
PC aa C36:5/PC aa C42:0	51.17	24.93	−0.13	2.97 × 10−14	0.10	5.84 × 10+04
PC aa C36:0/PC aa C38:0	0.83	0.13	−0.11	3.23 × 10−09	0.04	1.34 × 10+08
Σ sat. PC ae	16.50	2.81	−0.08	7.06 × 10−05	0.12
Σ mono PC ae	24.60	4.98	−0.07	6.06 × 10−04	0.17
Σ short PC ae	38.28	8.23	−0.06	1.03 × 10−02	0.19
PC ae C40:1	1.69	0.40	−0.11	2.11 × 10−10	0.08
PC ae C42:2	0.68	0.15	−0.09	4.58 × 10−07	0.12
PC ae C42:2/PC ae C42:5	0.29	0.06	−0.15	9.27 × 10−20	0.07	4.94 × 10+12
PC ae C40:1/PC ae C40:5	0.48	0.10	−0.14	2.92 × 10−17	0.09	7.22 × 10+06
PC ae C38:2/PC ae C42:5	0.92	0.22	−0.14	2.97 × 10−16	0.08	9.92 × 10+06
Σ unsat. lysoPCs	63.67	17.36	−0.08	2.92 × 10−05	0.24
Σ lysoPCs	191.45	40.77	−0.08	6.32 × 10−05	0.19
Σ sat. lysoPCs	127.78	26.46	−0.07	2.02 × 10−03	0.14
lysoPC a C14:0	3.23	0.84	−0.09	1.37 × 10−06	0.05
lysoPC a C18:1	19.77	5.79	−0.08	3.40 × 10−05	0.22

The full results are given in the Supplementary Table S2

AC acylcarnitines, even even numbered, long long chain, med. medium chain, short short chain, sat. saturated, unsat. unsaturated

aFor multiple testing adjusted p value

bAdjusted R2 of the linear model

cp-gain, fold decrease in the p value association for the pair of metabolites, compared to the lowest of two p values for the single metabolites

FT4 and serum acylcarnitines

For a series of different serum acylcarnitines, statistically significant associations with serum FT4 concentrations were detected, and all of them showed a positive direction (of the β-estimate). This could be indicative of an overall elevated transport and subsequent β-oxidation of fatty acids with higher serum FT4 concentration. For acetylcarnitine (C2) as the quantitatively most important serum acylcarnitine the adjusted p value was 4.46 × 10−10. Similarly strong associations were detected for C16:1, C14:2, C18:2, the sums of all acylcarnitines, long-chain acylcarnitines, even-numbered acylcarnitines, and acylcarnitines carrying monounsaturated or polyunsaturated fatty acids. Special ratios, such as C2 by C0 or the sum of C2 and C3 by C0 were suggested as valid measures of β-oxidation of fatty acids (Pande 1975; Roschinger et al. 2000) and they were both positively associated with FT4 hormone concentrations in our study (Table 3, Fig. 2a).

Fig. 2

Boxplot of serum C2/C0 (acetylcarnitine/free carnitine) concentrations (a) and the sum of PC aa concentrations (b) by quintiles of FT4 (pmol/L), for the KORA F4 population

FT4 and serum phosphatidylcholines

The group of phosphatidylcholines consists of different PC diacyl (aa), PC acyl-alkyl (ae), and lysoPC acyl (a) compounds. Numerous associations between FT4 and PCs or PC ratios were observed. Generally, an increase in FT4 levels was associated with lower concentrations of PCs. This was reflected by the sum of all PCs, and the different sums of PCs containing saturated, monounsaturated, polyunsaturated, or short-chain fatty acids (Fig. 2b). These results are in line with the associations found for the sums of PC aa, PC ae, or lysoPC with FT4 concentrations. Furthermore, single PCs such as PC aa C34:4, PC aa 36:1, PC ae C40:1, PC ae C42:2, lysoPC a C14:0, or lysoPC a C18:1 were all inversely associated with increasing serum FT4 concentrations. Regarding metabolite ratios, we received some indication for an increase in saturation and chain length for PC aa by PC aa ratios, such as PC aa C36:6 by PC aa C42:0, PC aa C36:5 by PC aa C42:0 or PC aa 36:0 by PC aa C38:0 with elevated FT4 levels. For PC ae by PC ae ratios, an increase in desaturation and chain-length with higher FT4 concentrations was observed, e.g., PC ae C42:2 by PC ae C42.5, PC ae C40:1 by PC ae C40:5 or PC ae C38:2 by PC ae C42:5.

Subanalyses

Weight-stable subsample

Restricting our sample to the weight-stable subset and adjusting for FFMI instead of BMI resulted in a substantial loss of significant associations with FT4 as compared to the results in the entire sample (Supplementary Table S3). Although β-estimates are similar in both analyses, statistical power was lower due to the smaller size of the weight-stable group. Still, long-chain acylcarnitines, the sums of long-chain acylcarnitines and of monounsaturated acylcarnitines were significantly positively associated with FT4 levels in this subgroup. The significant associations between single PC concentrations and FT4 were similar to the results of the entire sample. However, significant associations with PC ratios are largely missing, indicating loss of statistical power as a consequence of the reduced sample size. Adjusting the linear models for BMI instead of FFMI in the weight-stable subgroup resulted in almost the same significant associations with FT4, suggesting no distinct effect of the more precise adjustment for lean body mass.

Males versus females

The acylcarnitine results for men were very much comparable to the results in the whole sample while the results for the female participants varied slightly (Supplementary Table S4). Here, the positive associations with the sums of acylcarnitines were still present, but regarding the single concentrations, only C18:1, C14:2, C18:2, C2, and C16:1 are positively associated with elevated FT4 hormone concentrations. Testing for differences of the estimated coefficients between men and women revealed no statistical significance. The associations between PCs and FT4 levels in men and women were in line with the results for the entire study population.

Obese versus non-obese

The associations with FT4 levels found for the whole population and the non-obese participants were similar for both, the acylcarnitine and the PC concentrations (see Supplementary Table S5). In the obese group, many significant associations identified in the entire group were lacking which is likely due to sample size reasons. This is supported by the results of the β-estimates which were largely comparable in size and direction to those obtained in the entire group or the non-obese participants, although not reaching statistical significance (data not shown). Still, significant associations were identified for long-chain acylcarnitines with FT4 in the group of obese participants.

Physical activity

The acylcarnitine results for physically inactive participants paralleled the results in the whole sample, showing an enhanced oxidation activity with no preference of acylcarnitine type. The results for the active participants varied slightly (Supplementary Table S6). Here, the positive associations with most sums of acylcarnitines were also present, but regarding individual compounds, only C18:1, C14:2, C16:1 and C2 were still positively associated with elevated FT4 hormone concentrations. Testing for differences of the estimated coefficients between groups revealed no statistical significance. The associations between PCs and FT4 levels in active and inactive participants were in line with the results for the entire study population. For the inactive participants, however, most associations regarding metabolite ratios were missing.

Age groups (≤60 years versus >60 years)

The results for the ≤60 age group (n = 1132) are very similar to those of the entire study population concerning significance and estimate size. In the >60 age group, associations with most sums of acylcarnitines and sums of PCs were still present. However on the individual compound level and concerning PC ratios several significant associations with FT4 present in the entire population and the younger age group were lacking. This could be due to sample size reasons (>60 years; n = 331). However, testing for differences of the estimated coefficients between age groups revealed no statistical significance. Results are shown in Supplementary Table S7.

Influence of insulin, liver, kidney function, or inflammation parameters

Almost all identified associations between FT4 and different serum metabolites were independent of serum insulin, and parameters of liver and kidney function, or inflammatory biomarkers such as glutamate pyruvate transaminase, glutamic oxaloacetic transaminase, gamma-glutamyl transpeptidase, cystatin C, creatinine, C-reactive protein, or white blood cell count. Only the associations between FT4 and C18:1, PC aa C38:5, PC aa C40:4, PC aa C40:5, PC aa C42:6, PC ae C40:0, PC ae C42:1 and lysoPC a C16:1 were influenced by gamma-glutamyl transpeptidase, and the association between FT4 and lysoPC a C18:2 was affected by C-reactive protein.

Discussion

In this euthyroid population, after overnight fasting, strong associations between serum FT4 and serum metabolite concentrations were found. With higher FT4 levels, the concentration of acylcarnitines as well as the ratio of C2 by C0 increased, indicating an overall increase of β-oxidation of fatty acids. In addition, the inverse associations between FT4 and PC metabolites are reflective of a uniform decline in serum PC concentrations. Figure 3 illustrates a schematic description of cellular fatty acid transport, β-oxidation, and how these are linked to our found associations between FT4, acylcarnitines, and PCs. These results were supported by the associations found in stratified analyses and strengthened by their independency of serum insulin and different liver, kidney function, and inflammation parameters. In euthyroid adults, TSH concentrations were not associated with serum metabolites. Also, no effect of FT4 on serum amino acids concentrations was found.

Fig. 3

Schematic description of the cellular fatty acid transport and β-oxidation as the metabolic basis for the observed associations between FT4 and serum acylcarnitines and phosphatidylcholines (PCs)

Knowledge of how thyroid hormones exert their effects has further increased since the successful molecular cloning of nuclear thyroid hormone receptors in the mid 90’s (Flamant et al. 2007; Ribeiro et al. 1998; Wagner et al. 1995; Yen 2001). Most of the physiological effects of thyroid hormones are exerted through thyroid hormone receptor complexes and the response elements of regulatory regions of target genes (Biondi 2012a; Brix et al. 2011). Nevertheless, not all thyroid hormone effects are attributable to nuclear mediated pathways (Moreno et al. 2008) and an increasing number of non-genetic effects such as the transportation of ions, e.g., calcium, sodium, and potassium across plasma membranes or glucose and amino acids have been described lately (Biondi et al. 2002; Fazio et al. 2004). From a clinical point of view, it is clear that minor changes in thyroid hormone concentrations cause distinct metabolic effects (Kim 2008). In a recent microarray study by Clement et al. (2002) on the transcriptome of the vastus lateralis, a daily dosage of 75 μg of T4 was administered to healthy male volunteers with the result of 22 genes involved in mitochondrial energy metabolism being up-regulated. At the cellular level, T4 is converted to T3, and T3 is considered the active hormone. It is evident that intracellular T3 concentrations do not necessarily correspond to blood T3 concentrations. Although serum FT3 concentrations were measured in our study as well, the applied analytic method did not fulfil the strict quality control criteria in the end (Thienpont et al. 2010). Thus, FT3 data could not be used in this project. However, our results and the findings of the microarray study argue for an independent effect of FT4 on the mitochondrial level. There is no indication of an association between serum TSH concentration and serum metabolites. It is serum FT3 and FT4 that reveal biological effects at the cellular level. In an euthyroid situation, the serum TSH concentration was not strongly linked to the concentration of FT4 (correlation coefficient r = −0.15) in our study.

In a human experimental setting, it has recently been described that all plasma acylcarnitines (except for C3, C4, and C5) showed a strong increase in times of elevated β-oxidation, e.g., during fasting and or exercise (Krug et al. 2012). Concomitantly, plasma free carnitine (C0) decreases. In this project blood samples from all participants were drawn in fasting state, i.e., after overnight fast (minimum of 8 h). Thus, all study participants were in a catabolic metabolic state which is characterised by low plasma insulin concentrations, increased plasma concentrations of non-esterified fatty acids and ketone bodies, and an enhanced gluconeogenesis. Under this fasting condition, our results strongly indicate elevated β-oxidation with increasing plasma FT4 concentrations. Results from single acylcarnitine concentrations, sums of different acylcarnitine groups, as well as the ratio of C2 by C0 (Pande 1975; Roschinger et al. 2000) support this conclusion. The observed decreasing serum PC concentrations with higher FT4 levels are probably the consequence of this enhanced lipid oxidation, and therefore underline our hypothesis.

Krug et al. (2012) further demonstrated that plasma concentrations of specific amino acids, especially branched-chain amino acids, increased during fasting status, indicating enhanced gluconeogenesis. Although, a strong relationship between branched-chain amino acid serum concentrations and fat free mass in fasted state has been described for the present population (Jourdan et al. 2012), we did not observe significant associations between serum amino acid concentrations and FT4 levels in this project. Thus, we conclude that FT4 does not affect gluconeogenesis in euthyroid and fasting subjects. This finding affirms the notion that we do not simply see the effects of fasting status in our study but a specific effect attributable to the thyroid hormone FT4 which is characterised by enhanced transport of fatty acids in form of carnitine acylesters to the mitochondrion and subsequent β-oxidation. To the best of our knowledge no study directly comparable to ours has been performed so far. Most recent studies using transcriptomics or proteomics techniques investigated effects of T3 administration. In these studies effects of T3 on proteins of gluconeogenic pathways have been described; a review is given by Silvestri et al. (2011).

For treatment of obesity, T3 and T4 in varying doses and duration have been administered to euthyroid obese subjects during caloric deprivation to enhance weight loss and a systematic review on this issue is given by Kaptein et al. (2009). The authors summarised that available studies are insufficient to draw conclusions on beneficial or adverse effects of T3 or T4 administration; however, because of possible side-effect they do not advise on treatment of obesity with T4 or T3 (Kaptein et al. 2009). The results of the subanalyses presented here did not suggest specific differences by fat free mass (analysis in weight-stable subjects) or between obese and non-obese subjects. Also, subanalyses by sex, age, and physical activity provided largely similar results across subgroups. Thus, we conclude that the positive association between FT4 and β-oxidation of fatty acids after overnight fast is not dependent on lean body mass, body fat mass, age, sex, or habitual physical activity. In other words, we gained no indication for obesity being a consequence of altered fatty acid oxidation in fasting status as related to serum FT4 concentrations in euthyroid subjects. Further, no impact of serum insulin, or biomarkers of liver or kidney function or inflammation on the associations between FT4 and AC were observed. This supports the stability of our results and argues for an independent effect of FT4. However, this lack of impact could also be due to the characteristics of our study population. The selection of our participants was accompanied by an exclusion of elderly participants which might represent a risk group, e.g., for inflammation.

Strength and limitations

The human serum metabolome is currently characterised in many studies using different analytic approaches (Psychogios et al. 2011), including the method applied here. To the best of our knowledge, the association of FT4, TSH, and serum metabolites was not explored before in an euthyroid population. The participants of our study represent a random sample of the underlying population. In such a population with a heterogeneous metabolic makeup influenced by different environmental factors, genetic predispositions, and lifestyles, the interpretation of metabolic profiling is quite ambitious (Suhre et al. 2010). The metabolite concentrations which are reported in the present paper were identified by means of the AbsoluteIDQ™ p150 kit. This technology allows the automated quantification of many metabolites with a focus on amino acids, acylcarnitines, and PCs for many samples at a time (Altmaier et al. 2011; Wenk 2005). Thus, the kit is well suitable to assess a suggested impact of thyroid hormones on energy metabolism, especially lipid metabolism. The BIA measurements did only take place in KORA S4. For the analyses of associations with FT4 and the metabolite concentrations in the linear models with the FFMI instead of the BMI as adjustment variable only the results of the BIA measurements at baseline examination (KORA S4) were available. Assuming that the participants’ body composition is fairly stable in weight-stable subjects, those analyses were restricted to the weight-stable subjects in order to account for this.

Conclusion

Our results indicate that in a euthyroid population after overnight fast, serum FT4 concentrations are strongly linked to serum acylcarnitines and PCs, indicating enhanced transport of fatty acids to the mitochondrion and subsequent β-oxidation. These findings were stable as they could be largely reproduced in different subsets of the population.

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 168 kb)