Associations of Whole Blood n-3 and n-6 Polyunsaturated Fatty Acids with Blood Pressure in Children and Adolescents – Results from the IDEFICS/I.Family Cohort

Background Polyunsaturated n-3 and n-6 polyunsaturated fatty acids (PUFA) are precursors of biologically active metabolites that affect blood pressure (BP) regulation. This study investigated the association of n-3 and n-6 PUFA and BP in children and adolescents. Methods In a subsample of 1267 children aged 2–9 years at baseline of the European IDEFICS (Identification and prevention of dietary- and lifestyle-induced health effects in children and infants) cohort whole blood fatty acids were measured by a validated gas chromatographic method. Systolic and diastolic BP was measured at baseline and after two and six years. Mixed-effects models were used to assess the associations between fatty acids at baseline and BP z-scores over time adjusting for relevant covariables. Models were further estimated stratified by sex and weight status. Results The baseline level of arachidonic acid was positively associated with subsequent systolic BP (β = 0.08, P = 0.002) and diastolic BP (β = 0.07, P<0.001). In thin/normal weight children, baseline alpha-linolenic (β = -1.13, P = 0.003) and eicosapentaenoic acid (β = -0.85, P = 0.003) levels were inversely related to baseline and also to subsequent systolic BP and alpha-linolenic acid to subsequent diastolic BP. In overweight/obese children, baseline eicosapentaenoic acid level was positively associated with baseline diastolic BP (β = 0.54, P = 0.005). Conclusions Low blood arachidonic acid levels in the whole sample and high n-3 PUFA levels in thin/normal weight children are associated with lower and therefore healthier BP. The beneficial effects of high n-3 PUFA on BP were not observed in overweight/obese children, suggesting that they may have been overlaid by the unfavorable effects of excess weight.


Methods
In a subsample of 1267 children aged 2-9 years at baseline of the European IDEFICS (Identification and prevention of dietary-and lifestyle-induced health effects in children and infants) cohort whole blood fatty acids were measured by a validated gas chromatographic method. Systolic and diastolic BP was measured at baseline and after two and six years. Mixed-effects models were used to assess the associations between fatty acids at baseline and BP z-scores over time adjusting for relevant covariables. Models were further estimated stratified by sex and weight status.

Introduction
Hypertension is a major public health issue in industrialized countries. Given the high prevalence of overweight and obesity in all age groups, already children are affected by elevated blood pressure (BP) [1][2][3][4]. Further, elevated childhood BP has been shown to predict elevated BP in adolescence and adulthood [5][6][7] and other adulthood cardiovascular diseases [8][9][10]. Serum and erythrocyte fatty acids (FA) have been suggested to be associated with BP in adults [11][12][13] although not all studies confirmed this association [14].
In particular, polyunsaturated fatty acids (PUFA) and their metabolites, the long-chain PUFA (LC PUFA) such as eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) of the n-3 series and arachidonic acid (ARA, 20:4n-6) of the n-6 series have been linked to BP [11,12,14]. As precursors for the production of prostaglandins and thromboxanes with effects on vasodilation and platelet aggregation LC PUFA can affect BP [15]. Additionally, ARA and EPA are substrates for the cytochrome P450 (CYP450) catalyzed biosynthesis of metabolites of which the ARA metabolite 20-hydroxyeicosatetraenoic acid  has been shown to act as a vasoconstrictor. In contrast, epoxides of ARA like epoxyeicosatrienoic acids (EETs) and of n-3 LC PUFA like epoxyeicosatetraenoic acids seem to exert BP lowering effects, e.g. by increasing nitric oxide production and vasodilation [16,17]. Beneficial effects of n-3 LC PUFA on BP have also been attributed to their role as competitors of n-6 LC PUFA in the biosynthesis of eicosanoids and lipid mediators including those catalyzed by CYP450 [18]. Additionally, 20-HETE and EETs are involved in renal tubular and vascular function that may affect BP [19][20][21].
Intervention studies with fish oil supplements or n-3 rich fatty fish showed a small blood pressure lowering effect [22]. Also a meta-analysis of randomized controlled trials with EPA and DHA supplementation concluded that n-3 LC PUFA lowers systolic BP (SBP) and in high doses also diastolic BP (DBP) [23]. However, data from observational studies and from the limited number of intervention studies in children are less consistent [18]. Three studies reported positive associations between childhood n-3 LC PUFA and BP. In Danish children, cross-sectional data indicated whole blood EPA to be positively associated with DBP in boys but not in girls [24] while in Danish adolescents DHA was positively associated with SBP [25]. In a Finnish cohort, the sum of serum cholesterylester n-3 PUFA in childhood was positively associated with BP in adulthood after 27 years in males but not in females [26]. Against the background of the presumed biological mechanisms and results of intervention studies in adults, the reported positive associations of n-3 PUFA and BP were unexpected. Therefore, this study investigates the cross-sectional and prospective associations between whole blood n-3 and n-6 PUFA and BP in a large cohort of European children.

Study group
In the IDEFICS (Identification and prevention of dietary-and lifestyle-induced health effects in children and infants) baseline survey (T0) in 2007/2008, a population-based sample of 16 228 children aged 2 to 9.9 years from eight European countries (Belgium, Cyprus, Estonia, Germany, Hungary, Italy, Spain, Sweden) was examined. Follow-up examinations were conducted two (T1) and six (T3, I.Family study) years later; the study design has been described in detail elsewhere [27,28]. In brief, the IDEFICS/I.Family study comprises one of the largest prospective European child cohorts. The main focus of the IDEFICS study was to investigate diet-and lifestyle related disorders in young children with a focus on overweight and obesity. A further follow-up of the health status including 7105 children of the original IDEFICS cohort was performed in the I.Family study to explore the determinants of dietary and other lifestyle behaviors in these children and adolescents as they have grown older. A large set of anthropometric and clinical examinations was conducted and biosamples were taken as well as several questionnaires and interviews were collected. Lifestyle behaviors like dietary intake, physical activity and screen time of young children up to the age of 11 years were reported by parents whereas this information was self-reported by adolescents. Parents reported their own educational level and lifestyle behaviors like smoking in self-completion questionnaires. In a face-to-face interview the health conditions of the participating child and the family history of diseases were collected by trained interviewers.
Since in IDEFICS community-and setting-oriented intervention programs were implemented, in each country two mutually comparable areas have been allocated to either serve as intervention or control region [29].
Parents gave written informed consent for themselves and for their participating children. Children aged 12 years and over also gave written informed consent while small children gave oral consent prior to the start of the examinations. The participating centers in each country obtained ethical approval by the competent Institutional Review Boards prior to starting the study. These were the following ethics committees: This study is based on a subsample of IDEFICS T0 participants with an oversampling of overweight and obese children [27]. FA profiles were analyzed in 2600 children of which 10 children were excluded because of blood drawing in a non-fasting state. Furthermore, only children with at least one follow-up measurement were included such that additional 1297 children had to be excluded due to missing BP, missing follow-up or missing covariate data. Because the remaining sample of children from Belgium and Cyprus consisted of only 15 and 11 children, respectively, these were not considered in the final analysis resulting in an analysis sample of 1267 children (Fig 1). The age range of the children was 2-9 years at T0, 4-11 years at T1 and 9-16 years at T3, respectively.

Exposures
Blood was obtained by collecting a drop of blood from a fingertip or by venipuncture which was applied directly to a test strip prepared with butylated hydroxytoluene (BHT). FA in whole blood were separated and determined by gas-liquid chromatography after direct derivatisation to their methyl esters without prior extraction of total lipids from the samples in a central laboratory as previously described [30,31]. This simple and rapid analysis method was validated by several laboratories [31][32][33]. Whole blood includes FA from all lipid classes and is representative for the total FA pool [32]. In the present analysis we used FA values of the n-3 PUFA alpha-linolenic acid (ALA, 18:3n-3), eicosapentaenoic acid and docosahexaenoic acid and the n-6 PUFA linoleic acid (LA, 18:2n-6), gamma-linolenic acid (GLA, 18:3n-6), dihomo-gammalinolenic acid (DGLA, 20:3n-6) and arachidonic acid. FA are expressed as weight percentage of all FA detected (% wt/wt).
Outcomes SBP and DBP was measured at the right arm with a cuff of appropriate size for arm circumference after at least 5 minutes of rest in seated position according to a standardized protocol with an automated oscillometric device (Welch Allyn, Inc., 4200B-E2, Skaneateles Falls, NY, USA) [34]. Two recordings were taken, with a 2-min interval between each, plus a third measurement in case of a >5% difference between the first two readings. The mean value of the two measurements with the smallest difference was used as outcome measure. If the first and second measurement deviated by >5% and no third measurement value was obtained BP was set to missing. For SBP, DBP and BP classification to normal BP, pre-hypertension and hypertension, age, sex and height specific reference values from the US Department of Health and Human Services, National Institutes of Health (NIH, 2005) [35] were applied. SBP and DBP were further transformed to age-, sex-and height-specific z-scores according to the NIH references and were used as outcome variables in the subsequent analysis.

Anthropometric data
Height of the children was measured to the nearest 0.

Confounders
Potential confounders such as age (continuous), sex, country (6 categories: Estonia, Germany, Hungary, Italy, Spain, Sweden), maximum ISCED (International Standard Classification of Education [36]) level of parents (three categories: low level (ISCED 0,1,2); medium (ISCED 3,4); high (ISCED 5,6)) [36], parental smoking (5 categories, only assessed at T3: current daily smoker; current occasional smoker; ex-smoker; ex-smoker with missing pack years; never smoker), family history of hypertension (yes vs. no; as reported for biological parents/siblings at any of the three time points), BMI z-scores calculated based on Cole & Lobstein (2012) [37], height, birth weight and pubertal stage (pubertal, pre-pubertal, no information available; defined based on voice change in boys and first menstrual period in girls) were considered. To account for individual variations in the whole blood FA content depending on the FA distribution in lipid classes of plasma, lipoproteins and red blood cells of whole blood [32], we included weight percentage of total saturated (SFA) and total monounsaturated FA (MUFA) of total FA as covariates. A binary variable indicating control versus intervention regions was used to adjust for potential differences resulting from the intervention [29]. As indicator of dietary salt intake we used the weekly frequency of processed food consumption which includes Hamburger, hot dog, kebab, wrap, falafel, salty snacks such as savory pastries and fritters as well as cold cuts and ready to cook meat products. Television and computer time calculated as average screen time per week was also considered as a potential confounder.
Physical activity (PA) was measured objectively in a subset of children using the uniaxial Actigraph accelerometer (Actigraph MTI, model GT1M, Manufacturing Technology Inc., Fort Walton Beach, FL, USA) and the ActiTrainer which was the same technology as the Actigraph. We used time spent in moderate-to-vigorous physical activity (MVPA) according to the cutpoints of Evenson et al. (2008) [38] including children with at least 3 accelerometer measurement days and at least 360 minutes of valid time per day as indicator of PA in sensitivity analyses. measurements data as they do not require individuals to be measured at the same ages and as subjects with a varying number of repeated measurements can be considered under a missing at random assumption [39]. They further allow for correlations among repeated measurements taken from the same subject. In the present analysis, all children with a baseline and at least one follow-up BP measurement were included (N = 1267; 1267 baseline measurements; 1069 T1 measurements and 667 T3 measurements). A variable indicating time since baseline was constructed to assess how associations change over time (takes the values 0 for the baseline survey, approx. 2 years for the first follow-up and approx. 6 years for the second follow-up). The model was built as follows: where N denotes the number of study subjects; j = 0,1,2 denotes the survey wave), β 0 denotes the intercept, u i,0 a random subject-specific intercept, β 1 indicates the yearly change in average BP levels since baseline, u i,1 is a random subject-specific slope allowing BP to change differently in individuals over time, β 2 is the effect of the T0 exposure on the outcome at T0, β 3 describes the change per year in the effect of the FA exposure at baseline on BP over time and finally β covars was used to indicate all effect estimates of the covariates included in the different models (e.g. if age and sex were included as covariates, β covars would indicate the two effect estimates β age and β sex for age and sex). ε i,j is the error term for individual i at time j. An unstructured covariance matrix was selected for the random effects, i.e. each variance/covariance could take the value that the data demand.
For all exposures and the two outcomes (SBP z-score, DBP z-score), a basic model adjusting only for age (continuous), sex, height, country and intervention vs. control region (dummy variable) was estimated in a first step. In a second step, the following covariates were added to the basic model: maximum ISCED level of parents, parental smoking, family history of hypertension, BMI z-score, pubertal status, birth weight, weight percentage of the sum of SFA and of MUFA of total FA, weekly frequency of processed food consumption and weekly hours of screen time (fully adjusted model). Due to the significant interaction between EPA and weight status, models were further estimated stratified by sex and stratified by BMI category (thin/normal weight vs. overweight/obese; categorization according to Cole & Lobstein (2012) [37]).
In the subgroup with available accelerometer data (N = 1035), the fully adjusted model was additionally adjusted for minutes in MVPA.
All analyses were performed using SAS 1 statistical software version 9.3 (SAS Institute, Inc., Cary, NC). All models were run using SAS Proc Mixed. P values were not adjusted for multiple comparisons because of the explorative nature of this study. A P value of 0.01 was used as the level of statistical significance to account at least partially for multiple testing.

Results
The baseline characteristics of the study population are presented in Table 1. Of the girls, 20.7% and of the boys 20.2% were classified as pre-hypertensive or hypertensive. BP z-scores were higher in girls than in boys. Weight percentages of PUFA were similar among BP categories. 42% of the children were overweight or obese with slightly more boys than girls being classified as obese. BMI z-scores were higher whereas birth weight was slightly lower in children with normal BP compared to hypertensive children. Among the 220 children with a family history of hypertension, 27.3% were pre-hypertensive or hypertensive while this was true for only 19.0% of the 1047 children without family history of hypertension. 21.5% of children whose parents were current smokers were pre-hypertensive or hypertensive whereas this applied to 18.0% of children with never or ex-smoker parents. Table 2 shows the associations of the different PUFA with SBP and DBP z-scores at baseline, after two years and after six years of follow-up. In the fully adjusted model, baseline ARA was positively associated with subsequent SBP and DBP indicating that 1 unit increase of ARA was associated with a 0.08 units higher SBP and a 0.07 units higher DBP z-score after 6 years.
Results stratified by BMI category revealed that associations differ by baseline BMI ( Table 3). The observed positive association between ARA and BP pointed to the same direction in all subgroups. In thin/normal weight children, inverse associations of baseline ALA and EPA with SBP z-score at baseline and after 2 years of follow up and of baseline ALA with subsequent DBP z-score were observed. Children classified as overweight/obese showed a positive association of baseline EPA with DBP z-score at baseline. For instance, one unit increase of Table 2. Associations of polyunsaturated fatty acids with systolic and diastolic blood pressure z-scores at baseline, after 2 years and after 6 years of follow-up estimated based on fully adjusted mixed-effects models 1  baseline EPA was associated with 0.85 units lower SBP z-score at baseline in thin/normal weight children but with 0.54 higher DBP z-score at baseline in overweight/obese children. Table 4 shows the associations stratified by sex which confirms the positive association of baseline ARA with subsequent DBP in boys though not reaching significance in girls (p = 0.010). In boys, baseline GLA was positively associated with subsequent DBP indicating that 1 unit increase in GLA was associated with a 0.78 unit higher DBP z-score.
In the subgroup analyses with additional adjustment for MVPA the overall trends of associations remained unaltered (data not shown).

Discussion
In our large sample of European children, we observed that high blood levels of ARA are positively associated with subsequent SBP and DBP z-scores indicating a small detrimental effect of high ARA on BP. ALA and EPA showed a beneficial effect on BP in the subgroup of thin/  normal weight children in contrast to the results of previous observational studies investigating childhood blood n-3 PUFA (EPA and/or DHA) and BP [24][25][26]40]. However, in overweight/ obese children an unfavorable effect of EPA was observed. Thus, this study shows that already in children the weight status seems to modify the inverse association of n-3 FA with BP.

n-6 PUFA
The adverse effect of high ARA on BP observed in our study confirms previous results in adults [41,42] which to our knowledge have not been shown in children before. In contrast, in a Finnish cohort childhood serum cholesterylester levels of ARA were inversely associated with the BP after 27 years of follow-up but fat intake in Finland has strongly changed over time [26] and is probably different from our study population. Additionally, other factors may play a more important role after this long time period. ARA metabolites like prostaglandins and 20-HETE have vasoconstrictive, and thus blood pressure increasing, effects [15][16][17]. Additionally, ARA is a reliable indicator for animal-derived foods because the key desaturase steps in the conversion of LA to ARA, delta-6 desaturase (D6D) and delta-5 desaturase Table 4. Associations of fatty acids with systolic and diastolic blood pressure z-scores at baseline, after two years and after six years of followup stratified by sex 1 . Fatty Acids and Blood Pressure in Children (D5D), are inhibited by high concentrations of the substrate [15]. Thus, high blood ARA reflects a higher consumption of animal-derived foods with high intake of animal fats. Such a dietary pattern has also been shown to be associated with elevated BP because of the corresponding high sodium intake [43]. We observed the association of baseline ARA with BP after 6 years of follow up. This may indicate that the adverse effects of the dietary pattern may become stronger with increasing age as dietary patterns have been shown to track from childhood to adolescence [44].

Fatty acid
As observed in boys in our study, a positive association of GLA with BP was also found in Chinese and Western Alaskan Native adults [13,45]. GLA is generated from LA by D6D and can-after being elongated to DGLA-further be desaturated by D5D to ARA. As D6D activity is higher and D5D activity is lower in male subjects [46,47], the positive association in boys may be explained by sex differences in desaturase activity.

n-3 PUFA
Our results indicate a protective role of ALA and EPA against hypertension in thin/normal weight children whereas in overweight/obese children high EPA seemed to affect DBP unfavourably. Beneficial effects of n-3 PUFA on BP have been shown in several intervention studies investigating the effect of supplemental FA intake [48][49][50] and in observational studies investigating dietary intake [8, 23,[51][52][53]. While observational studies investigating the association of n-3 PUFA blood levels and BP showed beneficial effects in adults [11,42], results in children and adolescents were not consistent and indicated unfavorable effects [24,25,40]. In a Danish cross-sectional study including 8-11 year old children, whole-blood EPA was associated with a 2.7 mmHg higher DBP per weight% EPA in boys but not in girls [24]. Also in a large Finnish study an unexpected positive association of childhood cholesterylester EPA with BP after 27 years was observed indicating a stronger association in men than in women [26]. We did not observe sex-specific but a weight status-specific disparity in the n-3 PUFA-BP association in our study indicating a protective role of n-3 PUFA only in thin/normal weight children. Thus, the unfavorable effects of a high BMI and an altered n-3 PUFA status may have overlaid the protective effects of n-3 PUFA on BP. In line with our results, in Alaskan Native adults an inverse association of a marine food biomarker in blood and BP was observed in non-obese but not in obese subjects. The interaction of weight status and the biomarker was only significant for DBP but not for SBP [54]. In our analysis group, a high proportion of children (42%) were overweight or obese. Already in children, BMI has a strong independent blood pressure increasing effect [55,56], and thus in our study sample obese children were more often hypertensive than thin or normal weight children. Additionally, overweight and obese children enter puberty usually earlier than thin or normal weight children [57,58]. As BP increases with puberty [59,60] and early puberty has shown to be a risk factor for later high BP [57], this may have also contributed to higher BP in overweight/obese children. Also, inflammatory cytokines and adipocytokines from adipose tissue in overweight and obese subjects were shown to increase blood pressure [61,62] such that obesity-related effects on BP may have overlaid the BP lowering effects of n-3 PUFA observed in thin/normal weight children.
Data on the influence of BMI on blood n-3 LC PUFA concentrations are inconsistent [40,63] whereas lower estimated D5D activity has consistently been observed in obese subjects compared to normal weight subjects indicating a less effective last desaturation step to ARA and EPA, respectively [63,64]. Therefore, obese children may synthesize less EPA from its precursor. However, in our study, EPA percentages were slightly higher in overweight/obese children compared to thin/normal weight children (data not shown). This may have resulted from a higher total energy intake and a dietary pattern with high consumption of animal-derived foods including fish with high n-3 LC PUFA as also ARA percentages were higher in the upper BMI categories reflecting higher consumption of animal-derived fats as discussed before and as shown for overweight/obese subjects in previous studies [65][66][67].
Higher dietary EPA intake may be required to achieve BP lowering effects, particularly regarding DBP, as shown in a meta-analysis on the effects of fish oil or n-3 LC PUFA supplements on BP [23]. As children and adolescents generally do not meet the recommended dietary intake of n-3 LC PUFA [68], an increase of dietary EPA may beneficially influence BP.
According to a recent meta-analysis of prospective cohort studies in adult populations higher circulating DHA levels may prevent from elevated blood pressure [69]. In agreement with this, our data point to an inverse association of DHA with SBP in thin/normal weight children at baseline though the effect estimate did not reach statistical significance. In contrast, two smaller cross-sectional studies showed that whole-blood DHA was positively associated with DBP in 8-11 year old boys [40] and that erythrocyte DHA was positively related with SBP in 17 year old adolescents [25].
In our analysis, we used "time since baseline" (instead of age) as the main time axis in order to assess how associations between the FA exposures and BP change over time. However, when using age as the main time axis instead, we found that the association between ALA and SBP decreases with age whereas associations between ARA and DHA with BP increase with age.

Limitations and strengths
A limitation of our study is that the repeated blood pressure measurements were taken only on a single occasion at each survey wave whereas confirmation of BP values on repeated visits or even ambulatory 24-hour monitoring is recommended for the diagnosis of hypertension. Therefore intra-individual variability which is well documented in children could not be considered and may have led to an overestimation of BP values [35,70].
We adjusted only partially for multiple testing by choosing a significance level of 0.01 but no formal multiple testing correction was performed.
As detailed dietary intake data was only available in a smaller subgroup and would not have reliably reflected dietary salt intake, we used the intake frequency of processed foods as a proxy for sodium intake because processed foods have shown to be the main source of salt [71].
The large analysis group of >1200 children from 6 European countries as well as the longitudinal design and detailed phenotyping of the children are important strengths of our study. In all countries standardized protocols for data assessment, measurements and biosampling were applied at all survey waves and adherence to protocols was verified by a central quality control.

Conclusions
Our results suggest that ARA adversely affects BP already in children and may contribute to the development of hypertension whereas n-3 PUFA exert BP lowering effects in thin and normal weight children. In overweight and obese children the detrimental effects of increased weight status may have overlaid the protective effects of n-3 fatty acids such that an unfavorable relation of EPA and DBP became apparent in this subgroup. Higher EPA, and especially higher DHA, combined with a healthier food pattern might be required to lower BP in both, thin/normal weight and overweight/obese children. Health interventions need to consider effect modification by weight status and should focus particularly on overweight and obese children who are more susceptible to elevated BP despite high n-3 PUFA levels.