- 1Hormone Laboratory, Department of Medical Biochemistry and Pharmacology, Haukeland University Hospital, Bergen, Norway
- 2Mohn Nutrition Research Laboratory, Department of Clinical Science, University of Bergen, Bergen, Norway
- 3Department of Clinical Science, University of Bergen, Bergen, Norway
- 4Department of Heart Disease, Haukeland University Hospital, Bergen, Norway
- 5Mohn Research Center for Diabetes Precision Medicine, Department of Clinical Science, University of Bergen, Bergen, Norway
Background: Clinical studies on effects of marine-derived omega-3 (n-3) polyunsaturated fatty acids (PUFAs), mainly eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), and the plant-derived omega-6 (n-6) PUFA linoleic acid (LA) on lipoprotein-lipid components and glucose-insulin homeostasis have shown conflicting results, which may partly be explained by differential responses in females and males. However, we have lacked data on sexual dimorphism in the response of cardiometabolic risk markers following increased consumption of n-3 or n-6 PUFAs.
Objective: To explore sex-specific responses after n-3 (EPA + DHA) or n-6 (LA) PUFA supplementation on circulating lipoprotein subfractions, standard lipids, apolipoproteins, fatty acids in red blood cell membranes, and markers of glycemic control/insulin sensitivity among people with abdominal obesity.
Methods: This was a randomized double-blind crossover study with two 7-week intervention periods separated by a 9-week washout phase. Females (n = 16) were supplemented with 3 g/d of EPA + DHA (fish oil) or 15 g/d of LA (safflower oil), while males (n = 23) received a dose of 4 g/d of EPA + DHA or 20 g/d of LA. In fasting blood samples, we measured lipoprotein particle subclasses, standard lipids, apolipoproteins, fatty acid profiles, and markers of glycemic control/insulin sensitivity.
Results: The between-sex difference in relative change scores was significant after n-3 for total high-density lipoproteins (females/males: −11%*/−3.3%, p = 0.036; *: significant within-sex change), high-density lipoprotein particle size (+2.1%*/−0.1%, p = 0.045), and arachidonic acid (−8.3%*/−12%*, p = 0.012), and after n-6 for total (+37%*/+2.1%, p = 0.041) and small very-low-density lipoproteins (+97%*/+14%, p = 0.021), and lipoprotein (a) (−16%*/+0.1%, p = 0.028). Circulating markers of glucose-insulin homeostasis differed significantly after n-3 for glucose (females/males: −2.1%/+3.9%*, p = 0.029), insulin (−31%*/+16%, p < 0.001), insulin C-peptide (−12%*/+13%*, p = 0.001), homeostasis model assessment of insulin resistance index 2 (−12%*/+14%*, p = 0.001) and insulin sensitivity index 2 (+14%*/−12%*, p = 0.001), and quantitative insulin sensitivity check index (+4.9%*/−3.4%*, p < 0.001).
Conclusion: We found sex-specific responses after high-dose n-3 (but not n-6) supplementation in circulating markers of glycemic control/insulin sensitivity, which improved in females but worsened in males. This may partly be related to the sex differences we observed in several components of the lipoprotein-lipid profile following the n-3 intervention.
Clinical trial registration: https://clinicaltrials.gov/, identifier [NCT02647333].
1. Introduction
The risk of cardiovascular disease (CVD), a leading cause of death worldwide (1), is lower in females than in males (2–4). Several lines of evidence indicate that this may partly be explained by sex-specific differences in the metabolism and levels of fatty acids (FAs), lipoproteins, lipids, and apolipoproteins (apo) due to the actions of sex hormones, genetic mechanisms, and other factors (4–9). Sex hormones may for example influence the metabolism of essential FAs and lead to sex-specific differences in levels of polyunsaturated fatty acids (PUFAs) in blood and cell membranes (6, 10, 11), and several genotype × sex interactions have been associated with the concentrations of total PUFAs and omega-3 (n-3) PUFAs in different blood lipid fractions (12, 13). Additionally, different measures of body adiposity have been related to blood levels of FAs, lipoprotein subfractions, and lipids in a sex-dependent manner (9, 14–17).
Elevated blood levels of triacylglycerols (TAGs) and TAG-rich lipoproteins (TRLs) or their remnants are strongly associated with abdominal obesity (18), atherosclerotic CVD (19–21), and all-cause mortality (19), and increasing cardiometabolic risk has been observed with increasing TAG levels within the widely accepted “normal” range (22, 23). The concentrations of different lipoprotein subclass particles, as well as their composition of lipids and proteins, including apoB and apoC-III, have been shown to be independent predictors of CVD comparable to or stronger than routinely measured levels of cholesterol in low-density lipoproteins (LDLs) and high-density lipoproteins (HDLs) (24–33). Premenopausal females are generally characterized by less atherogenic lipoprotein-lipid-apolipoprotein (LLA) patterns than similar-aged males (34–38). Cross-sectional analyses have shown less atherogenic lipoprotein subfraction profiles characterized by smaller-sized very-low-density lipoproteins (VLDLs), larger-sized LDLs, and/or larger-sized HDLs in females relative to males (16, 34–39). Studies have overall reported lower levels of large/medium VLDLs and higher concentrations of large/medium HDLs in females compared with males (16, 36–38, 40), as well as lower TAGs and higher HDL cholesterol (HDL-C) (14, 16, 34–38, 40–48), while the sex differences were more variable in other lipoprotein subfractions and lipids.
Lower CVD risk and cardiovascular mortality have been associated with higher intakes or blood levels of marine-derived n-3 PUFAs, mainly eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (49–56), and with higher intakes or circulating levels of plant-derived omega-6 (n-6) PUFAs, mainly linoleic acid (LA) (57–60). These inverse relationships have been partially explained by changes in blood fat fractions (58, 61, 62), primarily reduced TAG levels after n-3 intake (51, 63–67) and lower cholesterol levels after n-6 intake (62, 68–70). However, effects on CVD morbidity, mortality, and intermediate outcomes, including lipoproteins and lipids, have not been consistent across studies for n-3 (50, 53, 71–74), nor for n-6 (57, 59, 75–78), and it remains controversial whether dietary or supplemental PUFAs reduce cardiometabolic risk (70, 73, 79). Inconsistent results may partly be explained by differential responses to PUFA supplementation between sexes and other subgroups, and previous studies have suggested that n-3 and n-6 PUFAs modify CVD risk factors differently in females and males (80–82).
We previously reported that high-dose supplementation of high-quality oils with n-3 (EPA + DHA) or n-6 (LA) PUFAs was followed by reductions in primarily TAG-or cholesterol-related LLA components, respectively, and concluded that the responses after both interventions point to changes in the LLA profile that have been associated with reduced cardiometabolic risk, also among people with TAG or LDL cholesterol (LDL-C) levels within the normal range (83). However, to the best of our knowledge, no randomized controlled trials (RCTs) have previously explored sex-dependent lipoprotein subgroup responses after n-3 or n-6 PUFA interventions. Additionally, human studies have produced conflicting results regarding the effects of marine n-3 PUFAs on glucose responses and insulin sensitivity (84–89), which are important factors in the metabolism of LLA components (90–92), and findings for sex-dependent effects of n-3 PUFAs on measures of insulin resistance are inconsistent (86, 88). Notably, very few studies have investigated sexually dimorphic responses in glucose-insulin homeostasis after increased consumption of n-6 PUFAs (LA) (93), and none have directly compared females and males. Thus, we report here from a randomized crossover trial exploratory analyses of sex-specific responses in concentrations and sizes of all major lipoprotein subclass particles, determined by nuclear magnetic resonance (NMR) spectroscopy, and in levels of associated lipids and apolipoproteins, as well as in indices of glycemic control/insulin sensitivity, at baseline and after high-dose supplementation with n-3 or n-6 PUFAs in people with abdominal obesity. We demonstrate sex-dependent responses after n-3 (but not n-6) supplementation in glucose-insulin homeostasis that may be related to sex differences in several components of the lipoprotein-lipid profile following the n-3 intervention.
2. Materials and methods
This study reports exploratory analyses of sex-specific responses at baseline and after n-3 or n-6 PUFA interventions in an RCT with a crossover design previously described in detail (83). Briefly, the study was registered at ClinicalTrials.gov (NCT02647333) and conducted according to the guidelines in the Declaration of Helsinki. The trial procedures were reviewed and approved by the Regional Committee for Medical and Health Research Ethics (2014/2336/REK South-East). All participants provided written informed consent.
2.1. Participants
People with a sedentary lifestyle (<2 h/week of physical activity), 30–70 years of age, and abdominal obesity [waist circumference (WC) of ≥80 cm in females and ≥ 94 cm in males] but otherwise healthy, were screened for eligibility. Exclusion criteria were as follows: diagnosed diabetes, severe psychiatric illness, malabsorption disorders, or other chronic diseases; regular use of medications (except antibiotics, NSAIDs, antihistamines, diuretics, and hormone replacement therapy if prescribed by a physician); fasting serum levels of TAGs >5 mmol/L; cigarette smoking; alcohol abuse (>7 alcohol units/week for females; >14 alcohol units/week for males); previous coronary interventions; previous bariatric surgery; pregnancy or lactation; blood donation within 3 months before baseline; scheduled hospitalization during the study; and pacemaker or implantable cardioverter defibrillator (ICD). Participants taking dietary supplements, including fish oil, were asked to stop supplementation at least 3 months before pre-treatment baseline, except for those taking prescribed iron, calcium, or vitamin D for medical reasons. Changes in prescribed medication and supplementation during the study have been previously reported (83). Participants were instructed to maintain their usual lifestyle, dietary habits, and physical activity level during the study. The habitual food intake was obtained by dietary recordings from Monday to Sunday two times in the run-in period (baseline) and one time during both intervention periods (in week five). Participants were also questioned about their physical activity level at each baseline and follow-up visit, as previously described (83).
2.2. Study design
This randomized double-blind two-period crossover study was conducted at Haukeland University Hospital in Bergen, Norway, from May 2015 to March 2016. After a run-in period of 15 week without any dietary supplements, eligible participants were randomly assigned to one of two treatment sequences: supplementation with n-3 fatty acids (fish oil) in period one, followed by supplementation with n-6 fatty acids (safflower oil) in period two (sequence AB), or n-6 supplementation in period one, followed by n-3 supplementation in period two (sequence BA). The intervention periods lasted for 7 weeks and were separated by a 9-week washout phase. Participants were at pre-treatment baseline allocated to the sequences by stratified randomization using sex as a stratum. Both participants and study investigators were blinded to the intervention order until statistical analyses of all primary and most secondary outcomes were performed. Clinical measurements and tissue sampling were conducted 1 day before treatment period one started (baseline visit 1, B1; September 2015), the day after 7 weeks of the first intervention (follow-up visit 1, I1; November 2015), before treatment period two started (baseline visit 2, B2; January 2016), and after the second intervention (follow-up visit 2, I2; March 2016).
2.3. Interventions
The n-3 supplement was a hydrolyzed and re-esterified fish oil (TAG form) containing mainly EPA (C20:5n-3) and DHA (C22:6n-3), while the n-6 supplement was an organic, cold-pressed safflower oil containing mainly LA (C18:2n-6). Based on analyzed PUFA levels in these supplementation products, females should consume 5.6 mL/5.0 g (1.8 g EPA and 1.2 g DHA) and males 7.5 mL/6.7 g (2.4 g EPA and 1.6 g DHA) fish oil each day during the n-3 period to obtain a total dose with EPA + DHA of 3 and 4 g/d, respectively. In the n-6 period, the daily dose of safflower oil was for females 24.5 mL/22.6 g (15 g LA) and for males 32.7 mL/30.1 g (20 g LA). Different amounts of oils among females and males (a ratio of 3:4) allowed for similar intakes of EPA + DHA or LA per kg body weight. The total amount of supplement was divided into two daily dosages, one in the morning and one in the afternoon with meals, and the participants were instructed to take the supplement (using a syringe with the dose marked) in each period from the morning after the baseline visit to the evening before the follow-up visit. Both supplements were produced and sponsored by Pharmatech AS (Rolvsøy, Norway). The FA compositions of the oils were analyzed by the Lipid Research Group at the University of Bergen, while analyses of primary and secondary oxidation products and free fatty acids, as well as acid values, were performed by Multilab Østfold AS (Rolvsøy, Norway) and indicated high-quality products at delivery (83). Participants were instructed to store the bottles with oil supplements at 4°C in dark (i.e., in the refrigerator) during the entire intervention periods to reduce the rate of lipid oxidation. The trial was blinded for the participants and the study investigators by equal appearance of the study products and the bottles containing these oils. The taste was modified by adding 1.8% of a natural citrus aroma and 0.045% of steviol glycosides (E 960). Equal appearance was accomplished by using 0.85% calcium carbonate (E 170) and 0.2% turmeric (E 100). We have previously reported an analysis of the blinding, as well as compliance with the interventions and recorded adverse events (83).
2.4. Study visits
Prior to each study visit, participants were instructed to fast for at least 10–12 h and avoid alcohol for at least 48 h. Besides, they were asked to not exercise the day before, as this may also affect circulating levels of lipids and lipoproteins (94). In the morning between 07:30 and 11:30 a.m., venous blood samples were collected, and anthropometric variables were measured. The participants were also questioned about fasting status, alcohol intake, physical activity, stress level, use of medications and supplements, any health conditions, and compliance with the intervention since the last visit.
2.5. Anthropometrics
Body weight, height, and waist and hip circumferences were measured by standardized procedures, and body composition, i.e., body fat mass, body fat percent, visceral fat area, and fat-free mass, was further analyzed by using a segmental multifrequency bioelectrical impedance measurement system (InBody S10; InBody Co., Ltd., Seoul, South Korea) in a lying posture. Studies indicate that multifrequency bioelectrical impedance analysis is a simple and convenient method for estimating visceral fat area and significantly correlates with computed tomography scan measurements (95, 96), although its accuracy may vary depending on the population being measured (97). The means of two measurements were calculated for all the anthropometric measures.
2.6. Biochemical variables
2.6.1. Fatty acid profiles in RBCMs
All blood samples were aliquoted and stored at −80°C until the end of the study when samples from all visits were analyzed at the same time. Fasting levels of 47 FAs in red blood cell membranes (RBCMs) were measured by gas–liquid chromatography using the internal standard C21:0 as previously described (83, 98). The total levels of analyzed n-3 and n-6 PUFAs are reported as the n-6/n-3 ratio, while the n-3 index was calculated as the sum of EPA and DHA in RBCMs and expressed as a percentage by weight of the total fatty acid (TFA) content (g FA/100 g TFA; wt%). We also report the levels of other RBCM FAs measured as wt% to eliminate differences in RBCM amounts, which may vary between individual blood sampling, even in the same person. However, since there is no agreement on the choice of FA measure (absolute vs. relative) (99–101), data on the concentrations (μg/mL) of RBCM FAs are also reported.
2.6.2. Lipoprotein particle subclasses
Concentrations and sizes of lipoprotein particle subclasses, and a lipoprotein-based insulin resistance index (LP-IR), were determined in EDTA-plasma samples by LipoScience (now LabCorp, Inc., Raleigh, NC, United States) using an automated NMR spectroscopy assay according to the LipoProfile-3 algorithm (102, 103). The following nine subfractions were measured by NMR (estimated ranges of particle diameter): large VLDLs (including chylomicrons if present; >60 nm), medium-sized (medium) VLDLs (42–60 nm), small VLDLs (29–42 nm), intermediate-density lipoproteins (IDLs; 23–29 nm), large LDLs (20.5–23.0 nm), small LDLs (18–20.5 nm), large HDLs (9.4–14 nm), medium HDLs (8.2–9.4 nm), and small HDLs (7.3–8.2 nm). Total VLDL and HDL particle concentrations were calculated as the sum of small, medium, and large subclass particle concentrations, and total LDL particle concentration is the sum of small LDL, large LDL, and IDL concentrations.
2.6.3. Lipids
We measured fasting serum TAGs, total cholesterol (TC), LDL-C, and HDL-C at the Department of Medical Biochemistry and Pharmacology (DMBP), Haukeland University Hospital, according to standardized procedures. Serum phospholipids, free cholesterol, and non-esterified fatty acids (NEFAs) were analyzed on a Hitachi 917 Chemistry Analyzer (Boehringer Mannheim GmbH, Mannheim, Germany) using kits from DiaSys Diagnostic Systems GmbH (Holzheim, Germany). Non-HDL cholesterol (non-HDL-C) was calculated by subtracting HDL-C from TC, while TRL cholesterol (TRL-C) was derived as non-HDL-C minus LDL-C.
2.6.4. Lipoprotein (a) and apolipoproteins
Serum concentrations of lipoprotein (a) [Lp(a)], apoB, and apoA-I were analyzed at DMBP. Serum levels of apoA-II, apoC-II, apoC-III, and apoE were measured using the MILLIPLEX MAP Human Apolipoprotein Magnetic Bead Panel—Cardiovascular Disease Multiplex Assay (APOMAG-62 K; Merck Millipore, Billerica, MA, United States) and detected by the Bio-Plex 200 System (Bio-Rad, Hercules, CA, United States).
2.6.5. Indices of glycemic control and insulin sensitivity
Fasting serum glucose, insulin, and insulin C-peptide (INCP), as well as glycated hemoglobin (HbA1c) in whole blood, were analyzed at DMBP. The homeostasis model assessment (HOMA) is reported as an index of insulin sensitivity, and we used the HOMA2 calculator developed by the University of Oxford1 to estimate insulin resistance (homeostasis model assessment of insulin resistance index 2, HOMA2-IR), insulin sensitivity (HOMA2-%S), and beta-cell function (HOMA2-%B) based on the updated computerized model (104). Fasting serum levels of glucose and INCP were used as inputs in these calculations. Another measure of insulin sensitivity, the quantitative insulin sensitivity check index (QUICKI), was calculated from the measurements of fasting glucose and insulin concentrations (105): 1/(log[fasting glucose (mg/dL)] + log[fasting insulin (μU/mL)]). We also report the revised version of QUICKI (rQUICKI), which incorporates the fasting level of NEFAs (106): 1/(log[fasting glucose (mg/dL)] + log[fasting insulin (μU/mL)] + log[fasting NEFAs (mmol/L)]). LP-IR, a marker of insulin resistance ranging from 0 (least) to 100 (most insulin resistant), was calculated by LipoScience as a weighted combination of lipoprotein subclass (large VLDLs, small LDLs, and large HDLs) and size (VLDL, LDL, and HDL sizes) parameters most closely associated with HOMA-IR (107).
2.6.6. Other biomarkers
Alanine aminotransferase, albumin, alkaline phosphatase, aspartate aminotransferase, bile acids, bilirubin, creatine kinase, γ-glutamyl transpeptidase, lactate dehydrogenase, estrogen (17β-estradiol), testosterone, sex hormone-binding globulin (SHBG), thyroid stimulating hormone, free thyroxine, and insulin-like growth factor-1 were analyzed at DMBP. The free androgen index, often used as a surrogate marker of free testosterone, was calculated as (testosterone/SHBG) × 100. Serum 25-OH vitamin D3 was analyzed by liquid chromatography–tandem mass spectrometry at BEVITAL AS, Bergen.2 The ketone bodies 3-hydroxybutyrate and acetoacetate were measured by gas chromatography-tandem mass spectrometry at BEVITAL AS.
2.7. Statistical analyses
2.7.1. Sex-specific responses in relative change scores (primary analysis)
The reported outcomes in this study are between-sex differences within treatments for relative change scores in fasting blood levels of LLA components: NMR-measured lipoprotein subclass particle concentrations and mean sizes, Lp(a), blood lipids (TAGs, NEFAs, phospholipids, TC, free cholesterol, LDL-C, and HDL-C), cholesterol levels calculated from these measures (non-HDL-C and TRL-C), apolipoproteins (apoB, apoA-I, apoA-II, apoC-II, apoC-III, and apoE), and ratios between these variables, as well as RBCM FAs. We also report between-treatment differences within females and males. The results presented here are derived from an intention-to-treat analysis including all randomized participants (16 females and 23 males), and the sample size calculation for the primary outcomes was previously reported (83).
Data are presented as raw unadjusted means (SDs), geometric means (1 SD ranges), or mean score differences (±absolute/relative effect estimates [95% CIs]). The geometric SD ranges used in descriptive statistics were calculated by dividing and multiplying the geometric means with the geometric SD factors to obtain the lower and upper limits, respectively (108). The distribution of data points from different measurements is shown by violin and error bar plots in Supplementary Figures S1–S4. All inferential tests were two-tailed with a nominal alpha level of 0.05. Adjustments for multiplicity were not performed due to the exploratory nature of the analyses and because a general adjustment method for mixed modeling of repeated measurements has not yet been developed due to difficulties associated with the correlation structure, which has to be taken into account (109). The statistical analyses were conducted with R v4.1.2,3 data transformation and exploration were done by using the tidyverse packages,4 and plots were made by the ggplot2 package v3.3.0 and GraphPad Prism v.8.0.1 (GraphPad Software, San Diego, CA, United States).
Study outcomes were analyzed by period-and baseline-adjusted constrained linear mixed-effects models (cLMMs), a constrained longitudinal data analysis technique (110–114), with “subjects” as the random factor by using the lme function in the nlme package v3.1–157. We followed a design-driven approach and defined a fixed effects structure including “time,” “sex,” “period,” “treatment × time,” “time × sex,” and “treatment × time × sex,” chose a random effects structure with random intercepts and slopes for “time,” and used the correlation structure compound symmetry, which is appropriate for a crossover design (115). In case of heterogeneity, the model included a variance structure allowing for different variances per stratum of “treatment,” “time,” and/or “sex.” In the between-sex and between-treatment comparisons, the male sex and the n-6 intervention were defined as the reference groups, respectively.
By excluding the main term “treatment” and thereby constraining the baseline values to be equal across treatment groups, a reasonable assumption in RCTs, the cLMM inherently adjusts for baseline differences (111, 113, 114), i.e., period-specific baselines (psb; within-subject effect) in a crossover design. We also included subject-averaged baselines (sab; between-subject effect) as a fixed covariable to control for cross-level bias in this design (115). To obtain both baseline-adjusted within-treatment changes and between-treatment differences in change scores for each stratum of sex, we replaced the two “treatment × time” interaction terms in the cLMM with a time-varying covariable that describes the treatment effect. This computational trick was done by adding the new variable to the data in a long format and setting this covariable equal to the “time” variable for one treatment group and to zero (i.e., baseline) at all time points for the other (reference) group. When specified correctly, the intercept from the model output will be equal regardless of the chosen reference treatment group. Since the cLMM assumes that the baseline means are equal across arms, the new time-varying covariable essentially estimates the mean difference between treatments at the follow-up time, which is analytically equivalent to an ANCOVA approach (113). However, in contrast to this approach, the cLMM also provides baseline-adjusted within-treatment change scores.
In a separate model, we additionally adjusted for age and covariables that can differ during the trial (period level factors) or between sexes. Total energy intake normalized to body weight was included because of different amounts of oils between the intervention periods. Vitamin D3 was controlled for due to possible effects on blood lipids, such as LDL-C (116), and its often-varying tissue concentrations between seasons. Additionally, we adjusted for alcohol intake normalized to body weight (although participants were instructed to avoid alcohol for at least 48 h before study visits) and physical activity level (although we recruited people with a sedentary lifestyle), which impact LLA components (94, 117–120). To assess the effect of menopausal status, we fitted a separate model including a binary variable for menopausal status among females (0: premenopausal; 1: postmenopausal) and age categories among males (0: <55 years; 1: ≥55 years). Females were classified as postmenopausal based on the following criteria: follicle-stimulating hormone (FSH) ≥30 IU/L or age ≥ 55 years.
Values were transformed by the natural logarithm before analyses of responses in relative terms (108). Relative within-treatment changes from baseline to follow-up for each stratum of sex and between-sex differences in change scores within treatments are reported in the main text and tables as percentages calculated from the regression coefficients (i.e., the average of log-ratios) by the formula 100 × (expestimate − 1). In the figures, however, we show results in relative terms as sympercents (s%), of reasons previously explained (83). This relative measure is calculated as the difference between the natural logs of two numbers multiplied by 100, i.e.,100 × ln(a) – 100 × ln(b), making it straightforward to present and interpret without back transformation (121). As part of the model validation procedure, the Shapiro–Wilk test for normality, the D’Agostino test for skewness, and graphical tools (boxplots, quantile-quantile plots, and histograms) were used to assess the distribution of standardized residuals [see Supplementary material (Supplementary Text, Materials and Methods), for further details about the model validation procedure].
2.7.2. Sex-specific responses in relative follow-up scores (sensitivity analysis)
When no outcome data are missing, the test for treatment difference over time in the cLMM is essentially equivalent to a test for treatment difference in an ANCOVA mixed model (AMM) (110, 111, 113, 114), which is often regarded as the most robust and powerful method (122–124). However, the cLMM is at least as efficient and powerful as the AMM (111–114). Yet, the results of cLMM slightly differ from the results of AMM because of the random part of the models (111, 114). Additionally, when outcome data are missing, the two analyses are based on different populations because participants with missing baseline or follow-up measurements are deleted in the AMM, while the cLMM uses all available data (111–114). We, therefore, conducted a sensitivity analysis of follow-up scores adjusted for the main effects of period, period-specific baselines, and subject-averaged baselines. This AMM included the fixed terms “treatment,” “sex,” “period,” “psb,” “sab,” and “treatment × sex,” a random effects structure with random intercepts-only, the correlation structure compound symmetry, and a data-driven variance structure.
2.7.3. Missing data handling
Linear mixed-effects modeling has been shown to efficiently deal with data sets containing missing outcome values and to be the optimal estimator in trials with repeated outcome measurements comprising a large portion of missing data (125–127). Therefore, we did not pre-specify another strategy for dealing with potential intermittent missing repeated outcome measurements or missing data resulting from dropouts, such as multiple imputation before mixed modeling, as this has shown to add no obvious gain compared to a standard mixed model approach without imputed values (126, 127).
2.7.4. Linear regression-determined associations
Sex-specific associations between continuous variables were analyzed with linear regression models by using the lm function in the R stats package v4.1.2. Before the analyses, data were transformed by the natural logarithm and standardized by the scale function in the R base package v4.1.2. The models included LLA components/secondary outcomes as univariate outcomes and anthropometric measures/biochemical variables as univariable predictors of interest. In one model, we adjusted for age. Another model was additionally adjusted for body mass index (BMI) and WC. We performed split sample analyses (sex-stratified multivariable models) and interaction analyses [augmented product term models (128)] to estimate within-sex associations and between-sex differences, respectively. From these models, we obtained two different effect sizes: (1) standardized regression coefficients (95% CIs), and (2) partial Cohen’s f2 by using the partial_f2 function in the R package sensemakr v0.1.4. Cohen’s f2 ≥ 0.02 (2%), f2 ≥ 0.15 (15%), and f2 ≥ 0.35 (35%) represent a small effect/weak association, medium effect/moderate association, and large effect/strong association, respectively (129–131).
3. Results
3.1. Study participants
A flow diagram of the 39 randomized individuals (16 females and 23 males) included in the final intention-to-treat analysis is reported in Supplementary Figure S5 [CONSORT 2019 format (132)]. Pre-treatment clinical characteristics of the subjects randomly assigned to intervention sequence AB or BA are shown by sequence and period for each stratum of sex in Supplementary Table S1 (geometric means for biochemical variables) and Supplementary Table S2 (arithmetic means for biochemical variables). The study participants were at baseline middle-aged [females: 56.3 (SD 9.3, range 38–69) years; males: 55.2 (SD 9.5, range 37–69) years], overweight [BMI: 28.5 (SD 4.5) in females, 29.8 (SD 3.8) in males], had abdominal obesity [WC: 100 (SD 11.2) cm in females, 107 (SD 8.7) cm in males; visceral fat area: 134 (SD 51) cm2 in females, 183 (SD 61) cm2 in males], and were normolipidemic (n = 7; all males) or hyperlipidemic (n = 32; 16 males) (Supplementary Figure S6) (83). Among the females, we identified 10 (63%) to be postmenopausal at B1. None of the females with FSH levels <30 IU/L were 55 years or older, and none with FSH levels ≥30 IU/L were younger than 55 years. 12 (52%) of the males were 55 years or older.
The pre-treatment TAG levels were measured to be <1.13 mmol/L in 9 females and 8 males, ≥1.7 mmol/L in 2 females and 9 males, and ≥ 2.3 mmol/L in 1 female and 2 males. All females and males had at B1 an n-3 index in RBCMs higher than 4 wt%. According to proposed n-3 index risk zones for coronary heart disease (133), 4 females and 11 males were at intermediate (4–8 wt%) risk, while 12 females and 12 males were at low (>8 wt%) risk. In 3 females and 1 male, the n-3 index was >10 wt%. The mean B1 levels of LA in RBCMs were similar to the wt% of EPA + DHA (Supplementary Table S1) and lower compared with circulating levels of LA measured as wt% in erythrocyte phospholipids in seven prospective observational studies (59). The relative content of LA in RBCMs was measured to be >6 wt% in all participants, >8 wt% in 15 females and 18 males, and > 10 wt% in 1 female and 2 males. Thus, our sample was characterized by a relatively high EPA + DHA/LA ratio in blood in both females and males.
3.2. Dietary and supplemental oil intakes
The cross-sectional analyses of dietary data at the pre-treatment baseline showed no significant sex differences in absolute intakes of energy and macronutrients per unit of body weight or fat-free mass, except for a significantly higher fiber intake per unit of fat-free mass among females compared with males (Supplementary Table S3). During the n-3 intervention, the calculated average intake of supplemental EPA and DHA was in total 3.68 (SD 0.78) g/d in females and 4.04 (0.68) g/d in males (females vs. males [95% CI]: −0.35 g/d [−0.69, −0.013], p = 0.042; from linear regression). The average intake per 100 kg body weight was 4.63 (1.15) g/d in females and 4.21 (0.91) g/d in males (0.40 g/d [−0.08, 0.87], p = 0.101). When normalized to 100 kg fat-free mass, the average intake was 7.70 (1.81) g/d in females and 6.08 (1.15) g/d in males (1.48 g/d [0.82, 2.15], p < 0.001). During the n-6 period, the calculated average intake of supplemental LA was 13.7 (2.53) g/d in females and 17.7 (2.25) g/d in males (females vs. males [95% CI]: −3.92 g/d [−5.03, −2.81], p < 0.001). The average intake per 100 kg body weight was 17.2 (3.46) g/d in females and 18.4 (3.06) g/d in males (−1.24 g/d [−2.77, 0.28], p = 0.110). When normalized to 100 kg fat-free mass, the average intake was 28.1 (4.63) g/d in females and 26.7 (3.70) g/d in males (1.45 g/d [−0.58, 3.47], p = 0.160).
We observed non-significant reductions in energy intake from the habitual diet after both interventions (Table 1). However, the caloric contribution from the supplemental safflower oil resulted in significantly higher total energy intake after the n-6 intervention in both females and males. The significantly lower amount of supplemental PUFA oil during the n-3 periods compared with the n-6 periods, which was recorded to be on average 16 g/d (~144 kcal/d) less in females (n-3 vs. n-6 [95% CI]: −16.2 g/d [−18.4, −14.1], p < 0.001) and 22 g/d (~198 kcal/d) less in males (n-3 vs. n-6 [95% CI]: −22.1 g/d [−24.0, −20.3], p < 0.001), explains most of the approximately 170 less kcal that were consumed in total during the n-3 periods versus the n-6 periods in females (n-3 vs. n-6 [95% CI]: −170 kcal [−307, −32.9], p = 0.016) and males (−169 kcal [−287, −49.8], p = 0.006). However, and most importantly for our analyses of sex-dependent responses, changes from baseline in total energy intake per 100 kg body weight did not differ significantly between sexes after the n-3 (females vs. males [95% CI]: −31.1 kcal [−272, 210], p = 0.799) and n-6 (5.37 kcal [−234, 245], p = 0.965) interventions (Table 1). This was also the case when we normalized total energy intake to 100 kg fat-free mass.
Table 1. Recorded daily intakes of energy and macronutrients at baseline and during the intervention periods.1
The energy percentage of fat increased significantly with the n-6 supplementation in both females and males compared with the n-3 intervention, and a significant decrease in protein intake was recorded after n-6 versus n-3 (Table 1). The changes in these nutrients, as well as carbohydrates and fiber, did not differ significantly between sexes after n-3 or n-6, nor when the intakes were normalized to 100 kg body weight, while we observed a significant between-sex difference in change scores for energy percentage of alcohol after n-3, also when normalized to body weight. Notably, the intakes of saturated, monounsaturated, and trans FAs did not change significantly with any treatment within both sexes and were not significantly different between sexes within treatments. Importantly, the recorded model-adjusted changes in total, n-3, and n-6 PUFA intakes were not significantly different between the sexes (Table 1). Additionally, the Likert scale data on physical activity level showed no significant differences between sexes at B1 (females vs. males [95% CI]: 3.49 [−2.05, 9.04], p = 0.224), and no significant sex differences in change scores after n-3 (3.99 [−0.62, 8.60], p = 0.089) or n-6 (0.40 [−4.22, 5.03], p = 0.863) supplementation.
3.3. FA profiles in RBCMs
We measured a comprehensive panel of 47 FAs in RBCMs because cardiometabolic effects of supplemental PUFAs are largely mediated through changes in the FA composition of different tissues, and circulating FA levels can reflect such changes (133, 134), as well as adherence to the interventions (99). The cross-sectional analyses of pre-treatment data showed that the contents of FAs in RBCMs differed significantly between sexes only for γ-linolenic acid (Table 2 and Supplementary Table S4). In the period-and baseline-adjusted mixed modeling, we found that the relative changes from baseline to follow-up differed significantly between sexes after the n-3 intervention in the n-6/n-3 ratio and wt% of arachidonic acid (AA) (more reduced in males than in females), while there were no significant sex-specific responses after n-6 supplementation for any of the variables (Table 3 and Supplementary Table S5). The relative content (wt%) of supplemental fatty acids in RBCMs increased less, although statistically non-significant, among females than males for EPA and DHA after the n-3 intervention, and for LA after n-6 supplementation. Furthermore, our analyses demonstrated significant between-treatment differences within both females and males for all the reported FAs, except α-linolenic acid. Notably, the observed changes in RBCM fatty acid profiles indicated good compliance with the interventions for both sexes. In the sensitivity analyses of follow-up scores adjusted for the main effects of period and baselines, we found after the n-3 intervention significant between-sex differences in the RBCM n-6/n-3 ratio (females vs. males [95% CI]: 7.83% [0.63, 15.5], p = 0.033) and wt% of EPA (−11.6% [−21.0, −1.11], p = 0.032), docosapentaenoic acid (−4.07% [−6.92, −1.14], p = 0.008), and AA (4.12% [1.21, 7.11], p = 0.007), but not DHA (−2.49% [−5.86, 0.99], p = 0.153), while no significant sex differences were observed after n-6 supplementation (Supplementary Table S6).
Table 3. Sex-specific responses in relative changes for RBCM fatty acid levels after 7 weeks of supplementation with n-3 or n-6 PUFAs.1
3.4. Lipoprotein particle subclasses
We found in the cross-sectional analyses of sex-specific differences at the pre-treatment baseline that several LLA components differed significantly between females and males (Figure 1 and Supplementary Tables S7, S8). An overall assessment of these variables indicates that females, compared with males, had a less atherogenic profile based on most of the LLA components with significant sex differences (more large LDLs and HDLs, fewer small LDLs and HDLs, more total HDLs, larger mean HDL size, higher HDL-C, lower TRL-C, more apoA-I, lower apoB/apoA-I ratio) but not all (more TC, more apoC-III, lower apoC-II/apoC-III ratio).
Figure 1. Comparison of pre-treatment characteristics. The dot plots show relative between-sex differences (females vs. males) at the pre-treatment baseline in lipoprotein subclass particle concentrations and mean sizes (A), lipids and lipoproteins (B), and indices of glycemic control/insulin sensitivity and anthropometric measures (C). Data were analyzed with GLS models using the gls function in the R package nlme v3.1–157. In case of heterogeneity, the model included a variance structure allowing for different variances per stratum of “sex.” Before the analysis, values were transformed by the natural logarithm and multiplied by 100 to show differences as additive, symmetric percentages (sympercents; see main text). Error bars represent 95% CIs. Number of participants included: n = 39 (16 females, 23 males). Apo, apolipoprotein; FC, free cholesterol; GLS, generalized least squares; HbA1c, glycated hemoglobin; HDLs, high-density lipoprotein particles; HDL-C, HDL cholesterol; HOMA2-IR, homeostasis model assessment of insulin resistance index 2 (computer model); HOMA2-%B, homeostasis model assessment of 𝛽-cell function index 2 (computer model); HOMA2-%S, homeostasis model assessment of insulin sensitivity index 2 (computer model); IDLs, intermediate-density lipoprotein particles; INCP, insulin C-peptide; LDLs, low-density lipoprotein particles; LDL-C, LDL cholesterol; Lp(a), lipoprotein (a); LP-IR, lipoprotein-based insulin resistance index; PLs, phospholipids; s%, sympercent; TC, total cholesterol; TRL-C, TAG-rich lipoprotein cholesterol; VLDLs, vey-low-density lipoprotein particles. Total VLDLs and large VLDLs also include chylomicrons if present.
The period- and baseline-adjusted mixed modeling showed that the relative changes from baseline to follow-up differed significantly between sexes after the n-3 intervention in total HDLs (females vs. males: −11.3% [−3.64 μmol/L] vs. −3.29% [−0.96 μmol/L], p = 0.036) and mean HDL size (+2.09% [+0.21 nm] vs. −0.093% [−0.006 nm], p = 0.045), and after n-6 supplementation in total VLDLs (+37.0% [+10.1 nmol/L] vs. +2.13% [+0.11 nmol/L], p = 0.041) and small VLDLs (+97.1% [+10.1 nmol/L] vs. +14.2% [+2.26 nmol/L], p = 0.021) (Figure 2 and Supplementary Tables S9, S10). Large but non-significant sex differences after mean changes in opposite directions for females and males were observed after n-3 supplementation in small LDLs (−13.7% [−46.1 nmol/L] vs. +16.3% [+94.3 nmol/L], p = 0.140) and medium HDLs (−31.2% [−2.99 μmol/L] vs. +7.96 [−0.34 μmol/L], p = 0.097), and after the n-6 intervention in large VLDLs (+34.9% [−0.16 nmol/L] vs. −7.17% [−1.25 nmol/L], p = 0.254). Smaller non-significant sex differences after mean changes in opposite directions were observed for mean VLDL size after n-3 (+1.69% [+0.55 nm] vs. −3.83% [−1.97 nm], p = 0.323) and for small LDLs after n-6 (−2.59% [+35.0 nmol/L] vs. +7.92% [+43.4 nmol/L], p = 0.611). In the analyses of between-treatment differences within sexes, we found that females and males differed for medium HDLs, LDL size, and HDL size, all of them demonstrating a significant difference between interventions only in females. Among these variables, only medium HDLs showed a between-sex difference after n-3 that was significantly different from the between-sex difference after n-6 (n-3 vs. n-6 in females [95% CI]: −44.9 [−62.9, −18.3], p = 0.003; n-3 vs. n-6 in males [95% CI]: 4.07 [−25.7, 45.7], p = 0.815; females vs. males after n-3 vs. n-6: −47.1 [−68.5, −11.1], p = 0.017).
Figure 2. Differences between females and males in change scores for lipoprotein subfractions after 7 weeks of supplementation with n-3 or n-6 PUFAs. The bar plots show relative changes from baseline to follow-up and between-sex differences in change scores of particle concentrations and mean sizes for VLDLs (A-E), LDLs (F-J), and HDLs (K-O). VLDLs and LDLs were measured in nmol/L, HDLs in μmol/L, and particle sizes in nm. Data were analyzed with period- and baseline-adjusted cLMMs (see main text). Before the analysis, values were transformed by the natural logarithm and multiplied by 100 to show within-sex changes and between-sex differences as additive, symmetric percentages (sympercents; see main text). Error bars represent 95% CIs. Relative changes from baseline to follow-up within females and males are shown in orange and yellow bars, respectively. Between-sex differences in relative change scores are shown in gray bars. Number of participants included: n = 39 (16 females, 23 males) at baseline in P1; n = 38 (16 females, 22 males) at follow-up in P1; n = 38 at baseline in P2 (16 females, 22 males); n = 38 at follow-up in P2 (16 females, 22 males). cLMMs, constrained linear mixed-effects models; HDL, high-density lipoprotein particles; IDLs, intermediate-density lipoprotein particles; LDLs, low-density lipoprotein particles; n-3, omega-3 PUFAs; n-6, omega-6 PUFAs; PUFAs, polyunsaturated fatty acids; P1, the first intervention period; P2, the second intervention period; s%, sympercent; VLDLs, very-low-density lipoprotein particles. Total VLDLs and large VLDLs also include chylomicrons if present.
The sensitivity analysis of follow-up scores adjusted for the main effects of period and baselines showed a significant between-sex difference in mean HDL size (females vs. males [95% CI]: 2.50% [0.57, 4.47], p = 0.012) after the n-3 intervention, but not in total HDLs (−5.66% [−11.6, 0.72], p = 0.079) after n-3 and not in total VLDLs (4.76% [−15.3, 29.5], p = 0.659) or small VLDLs (33.8% [−7.21, 93.0], p = 0.115) after n-6 (Supplementary Table S6). In contrast to the primary analysis, we found here a significant sex-specific difference in small LDLs after n-3 (−42.9% [−61.0, −16.5], p = 0.005). Moreover, large but non-significant sex differences were found for medium HDLs after n-3 (−39.2% [−64.7, 4.70], p = 0.072) and for small LDLs (−27.5% [−50.5, 6.05], p = 0.095) after n-6.
3.5. Lp(a) and standard lipids
The relative changes from baseline to follow-up differed significantly between sexes after the n-6 intervention in Lp(a) (females vs. males: −15.5% [−23.1 mg/L] vs. +0.52% [+8.60 mg/L], p = 0.028), while there were no significant sex-specific responses after n-3 supplementation (Figure 3 and Supplementary Tables S9, S10). However, a large but non-significant sex difference was observed in NEFAs (+6.43% [+0.030 mmol/L] vs. −26.7% [−0.10 mmol/L], p = 0.065) after mean changes in opposite directions for females and males following the n-3 intervention. Smaller non-significant sex differences after mean changes in opposite directions were observed after n-6 for TAGs (+5.92% [+0.026 mmol/L] vs. −2.64% [−0.080 mmol/L], p = 0.380) and NEFAs (+0.87% [+0.033 mmol/L] vs. −16.6% [−0.063 mmol/L], p = 0.301). The analyses of between-treatment differences within sexes demonstrated that females and males differed for Lp(a), TC, LDL-C, and non-HDL-C, all of which showed a significant difference between interventions only in females. Among these variables, only Lp(a) showed a between-sex difference after n-3 that was significantly different from the between-sex difference after n-6 (n-3 vs. n-6 in females [95% CI]: 21.9 [6.25, 39.8], p = 0.005; n-3 vs. n-6 in males [95% CI]: −0.37 [−11.4, 12.0], p = 0.950; females vs. males after n-3 vs. n-6: 22.3 [2.15, 46.5], p = 0.029). The sensitivity analysis showed significant between-sex differences in NEFAs (females vs. males [95% CI]: 51.3% [10.1, 108], p = 0.012) following the n-3 intervention and in Lp(a) (−11.6% [−21.8, −0.23], p = 0.046) after n-6 supplementation (Supplementary Table S6). Moderate to large but non-significant sex differences were found for NEFAs (25.3% [−9.02, 72.5], p = 0.162) and TRL-C (−15.4% [−33.9, 8.27], p = 0.177) after n-6. In the exploratory analyses of linear regression-determined associations, we found after the n-3 intervention that relative change scores in NEFA and TAG levels were inversely (non-significantly) associated and stronger in males (standardized regression coefficient [95% CI]: −0.38 [−0.81, 0.058], p = 0.086, f2 = 16.4%) than females (−0.098 [−0.67, 0.47], p = 0.717, f2 = 0.98%), while positive (non-significant) relationships were observed after n-6 in both females (0.052 [−0.52, 0.62], p = 0.850, f2 = 0.27%) and males (0.36 [−0.070, 0.80], p = 0.096, f2 = 15.3%).
Figure 3. Differences between females and males in change scores for lipoprotein (a), blood lipids, and apolipoproteins after 7 weeks of supplementation with n-3 or n-6 PUFAs. The bar plots show relative changes from baseline to follow-up and between-sex differences in change scores for lipoprotein (a) (A), blood lipids (B–L), and apolipoproteins (M–T). Blood lipids were measured in mmol/L and apolipoproteins in g/L. Apo, apolipoprotein; HDL-C, HDL cholesterol; LDL-C, LDL cholesterol; Lp(a), lipoprotein (a); non-HDL-C, non-HDL cholesterol; TC, total cholesterol; TRL-C, TAG-rich lipoprotein cholesterol. See further explanation and abbreviations in Figure 2.
3.6. Apolipoproteins
We found no significant between-sex differences for any of the apolipoproteins in relative or absolute change scores (Figure 3 and Supplementary Tables S9, S10), or in the sensitivity analysis of follow-up scores (Supplementary Table S6). Notably, we observed that females and males had changes in opposite directions for apoA-I (females vs. males: −4.25% [−0.059 g/L] vs. +1.16% [+0.010 g/L], p = 0.052) after the n-3 intervention (Supplementary Tables S9, S10). Additionally, the analyses of between-treatment differences within sexes demonstrated that only apoA-I differed for females (significant) and males (non-significant), although the between-sex differences after n-3 and n-6 were not significantly different (n-3 vs. n-6 in females [95% CI]: −5.66 [−9.88, −1.24], p = 0.013; n-3 vs. n-6 in males [95% CI]: −1.63 [−5.39, 2.29], p = 0.406; females vs. males after n-3 vs. n-6: −4.09 [−9.69, 1.85], p = 0.171).
Taken together, this exploratory analysis of sex-specific responses in change scores of LLA components after supplementation with n-3 or n-6 PUFAs showed in females compared with males relatively fewer total HDLs and larger mean HDL particle size after the n-3 intervention, and relatively more of total and small VLDLs after n-6 supplementation, as well as more AA in RBCMs after n-3 and fewer Lp(a) after n-6, while we found no significant sex differences for the apolipoproteins. The sensitivity analyses of follow-up scores supported this for mean HDL particle size, AA, Lp(a), and apolipoproteins, but not for total HDLs, total VLDLs, or small VLDLs. Additionally, after n-3 supplementation, large non-significant between-sex differences in relative change scores were observed for small LDLs (non-significantly decreased in females and significantly increased in males) and NEFAs (non-significantly increased in females and significantly decreased in males), while the analyses of follow-up scores showed significantly fewer small LDLs and more NEFAs in females than in males following this intervention. We also found in the sensitivity analyses significantly less EPA in RBCMs from females compared with males after n-3 supplementation. Notably, the analyses of between-treatment differences in change scores within sexes demonstrated for several LLA components significant differences between the n-3 and n-6 interventions only in females. Finally, covariable-adjusted effect estimates additionally controlled for covariables that can differ during the trial (period level factors) or between sexes, are reported in Supplementary material (Supplementary text, Results).
3.7. Glycemic control and insulin sensitivity
We measured several indices of glycemic control and insulin sensitivity because differences in glucose-insulin homeostasis have shown to affect the metabolism and circulating levels of LLA components (90–92). The cross-sectional analyses of sex-specific differences at B1 showed that females compared with males had lower fasting glucose levels and higher insulin sensitivity (Figure 1 and Supplementary Tables S11, S12). Relative changes from baseline to follow-up differed significantly between sexes after the n-3 intervention for glucose (females vs. males: −2.07% [−0.12 mmol/L] vs. +3.90% [+0.21 mmol/L], p = 0.029), insulin (−30.5% [−2.48 mU/L] vs. +16.2% [+1.78 mU/L], p < 0.001), INCP (−11.8% [−0.077 nmol/L] vs. +12.5% [+0.084 nmol/L], p = 0.001), HOMA2-IR (−12.4% [−0.18] vs. +13.9% [+0.21], p = 0.001), HOMA2-%S (+14.2% [+13.2] vs. −12.2% [−8.60], p = 0.001), and QUICKI (+4.92% [+0.007] vs. −3.42% [−0.005], p < 0.001), while there were no significant sex-specific responses after n-6 supplementation (Figure 4 and Supplementary Tables S13, S14). The sensitivity analysis supported these findings (Supplementary Table S6). Thus, the results indicated that males were significantly more insulin resistant at baseline than females and that the circulating markers of glucose and insulin responses improved in females but worsened in males after n-3 supplementation. Notably, LP-IR, another measure of insulin resistance calculated from lipoprotein parameters associated with HOMA-IR (107), decreased after n-3 in both females (non-significantly) and males (significantly) and showed no significant between-sex difference. Besides, rQUICKI demonstrated no significant sex-dependent effect. Finally, in the analyses of between-treatment differences within sexes, we found that females and males differed for insulin, INCP, HOMA2-IR, HOMA2-%S, and QUICKI, of which all showed a significant difference between interventions only in females, and that the between-sex differences after n-3 were significantly different from the between-sex differences after n-6.
Figure 4. Differences between females and males in change scores for circulating markers of glycemic control/insulin sensitivity after 7 weeks of supplementation with n-3 or n-6 PUFAs. The bar plots show relative changes from baseline to follow-up and between-sex differences in change scores for glycemic control (A,B) and insulin sensitivity (C–J). Glucose was measured in mmol/L, HbA1c in %, insulin in mU/L, and INCP in nmol/L. HbA1c, glycated hemoglobin; HOMA2-IR, homeostasis model assessment of insulin resistance index 2 (computer model); HOMA2-%B, homeostasis model assessment of 𝛽-cell function index 2; HOMA2-%S, homeostasis model assessment of insulin sensitivity index 2; INCP, insulin C-peptide; LP-IR, lipoprotein-based insulin resistance index; QUICKI, quantitative insulin sensitivity check index; rQUICKI, revised QUICKI. See further explanation and abbreviations in Figure 2.
3.8. Anthropometrics
Body weight and composition were measured because the sexual dimorphism in anthropometric measures typically observed between females and males has been related to differences in glucose-insulin homeostasis and LLA components (9, 16, 18), and because we wanted to investigate if different anthropometric traits would be affected differently between sexes after high-dose PUFA supplementation. The cross-sectional analyses of pre-treatment data showed that females compared with males had higher body fat percent but less abdominal adiposity (Supplementary Tables S11, S12). No significant between-sex differences in relative or absolute change scores were observed (Supplementary Tables S13, S14), while the sensitivity analysis of follow-up scores showed significant sex differences in body weight (females vs. males [95% CI]: −1.28% [−2.25, −0.29], p = 0.013), BMI (−1.02% [−1.94, −0.093], p = 0.032), and body fat mass (−3.68% [−6.38, −0.91], p = 0.011) following the n-3 intervention (Supplementary Table S6). Notably, all anthropometric measures except fat-free mass changed after n-3 supplementation in opposite directions for females (decreased) and males (increased), while this was observed only for visceral fat area after the n-6 treatment (Supplementary Tables S13, S14). Additionally, small but significant increases in body weight and BMI were found among females and males after n-6 supplementation, probably related to the significantly higher total energy intake caused by this intervention, although the exploratory analyses of linear regression-determined associations did not show any significant relationships between relative changes in total energy intake and body weight for females (standardized regression coefficient [95% CI]: 0.28 [−0.15, 0.71], p = 0.182, f2 = 15.3%) or males (0.066 [−0.41, 0.54], p = 0.772, f2 = 0.48%).
3.9. Markers of liver function
The liver is a major regulator of glucose and lipid metabolism (135, 136), and circulating levels of liver enzymes may be related to hepatic lipid content (9, 137), which has been associated with hepatic production and secretion of TRLs (138–143). A standard liver panel was also part of the monitoring of potential adverse effects of the high-dose PUFA supplementation, which has been reported previously (83). All circulating markers of liver function (alanine aminotransferase, albumin, alkaline phosphatase, aspartate aminotransferase, bile acids, bilirubin, creatine kinase, γ-glutamyl transpeptidase, and lactate dehydrogenase) were at B1 found to be lower in females compared with males (Supplementary Tables S11, S12), indicating less hepatic lipid content in females compared with males (9, 137), possibly related to higher hepatic insulin sensitivity (9, 137, 144). We observed no significant between-sex differences in relative or absolute change scores in any of these measures after the PUFA interventions (Supplementary Tables S13, S14), while the sensitivity analysis showed a significant sex difference in γ-glutamyl transpeptidase (females vs. males [95% CI]: −15.8% [−28.6, −0.73], p = 0.041) after n-3 supplementation (Supplementary Table S6).
3.10. Additional hormones
Several hormones, including the sex steroids estrogen (17β-estradiol) and testosterone, were measured because they are known to affect the metabolism of n-3 and n-6 PUFAs (6), as well as the metabolism of LLA components (5, 8, 145). The secosteroid vitamin D3 was also measured due to possible effects on blood lipids, such as LDL-C (116). Our cross-sectional analyses of sex-specific differences at B1 showed, as expected, significantly lower testosterone levels in females than in males, higher SHBG, and lower free androgen index (Supplementary Tables S11, S12). Estrogen levels did not differ significantly, but it should be noted that both pre- and postmenopausal females were included in the study. Additionally, we found significantly lower free thyroxine levels in females than males, but no significant sex differences in thyroid stimulating hormone, insulin-like growth factor-1, or vitamin D3. None of the measured hormones showed significant between-sex differences in relative or absolute change scores (Supplementary Tables S13, S14), or in follow-up scores (Supplementary Table S6), after n-3 or n-6 supplementation. In the exploratory analyses of linear regression-determined associations, we found among males after the n-3 intervention a strong inverse relationship between relative changes of testosterone levels and the n-3 index (as well as in EPA and DHA) adjusted for age, BMI, WC, and SHBG (standardized regression coefficient [95% CI]: −0.62 [−1.06, −0.17], p = 0.009, f2 = 54.7%), while this was not evident among females (−0.022 [−0.90, 0.85], p = 0.955, f2 = 0.040%). We did not find any significant associations in relative change scores for testosterone versus LA after n-6, for estrogen versus n-3 index after n-3, or for estrogen versus LA after n-6.
3.11. Ketone bodies
Circulating levels of ketone bodies were measured as markers of hepatic beta-oxidation of FAs (146, 147), which has been shown to be affected by n-3 PUFAs (89). The fasting serum levels of 3-hydroxybutyrate and acetoacetate did not differ significantly between females and males at the pre-treatment baseline (Supplementary Tables S11, S12). These ketone bodies changed non-significantly after n-3 supplementation in opposite directions for females (increased) and males (decreased) and showed large but non-significant sex differences in relative change scores for 3-hydroxybutyrate (females vs. males: +13.1% [+16.5 μmol/L] vs. −25.8% [−27.8 μmol/L], p = 0.154) and acetoacetate (+17.1% [+10.3 μmol/L] vs. −11.2% [−7.69 μmol/L], p = 0.254) (Supplementary Tables S13, S14), and in follow-up scores (Supplementary Table S6). After the n-6 intervention, these ketones were non-significantly lowered in both females and males and demonstrated no significant sex-specific responses. Finally, the exploratory analyses of linear regression-determined relationships showed after adjustment for age, BMI, and WC strong positive associations between relative changes in 3-hydroxybutyrate and NEFA levels after the n-3 intervention for females (standardized regression coefficient [95% CI]: 0.85 [0.54, 1.16], p < 0.001, f2 = 322%) and males (0.75 [0.40, 1.11], p < 0.001, f2 = 118%), and after n-6 supplementation for females (0.64 [0.22, 1.06], p = 0.006, f2 = 103%) and males (0.56 [0.16, 0.96], p = 0.009, f2 = 50.8%). Similar associations were found for acetoacetate.
4. Discussion
Although biological sex is an important modifier of cardiometabolic risk factors (7, 9), including body adiposity, glucose-insulin homeostasis, blood fat fractions, and LLA profiles (5, 8–11, 18, 35, 36, 38, 91, 144, 145, 148–154), and several of these factors have shown sex-specific responses to dietary changes (5, 9, 45, 155), we have lacked analyses of sexual dimorphism in cardiometabolic risk following increased consumption of n-3 or n-6 PUFAs. We here report sex-dependent responses after PUFA interventions in LLA components, glycemic control, and insulin sensitivity based on data that have directly compared females and males using n-3 (EPA + DHA) or n-6 (LA) PUFA supplementation with high-quality oils and similar doses per kg body weight across sexes. In our analyses, we observed sex-dependent responses (females vs. males) following n-3 supplementation for small LDLs (fewer), total HDLs (fewer), mean HDL particle size (larger), and NEFAs (more), as well as EPA (less) and AA (more) in RBCMs, and after n-6 supplementation for small VLDLs (more), total VLDLs (more), and Lp(a) (fewer). Importantly, we found after the n-3 (but not n-6) intervention significant sex differences in circulating markers of glycemic control/insulin sensitivity, which improved in females but worsened in males.
The significant changes in markers of insulin sensitivity among females after n-3 supplementation are in line with a systematic review and meta-analysis of 31 RCTs, showing that interventions with n-3 PUFAs (α-linolenic acid, EPA, and/or DHA) for at least 6 weeks significantly improve indices of insulin resistance in females (86). Also, in two cohorts of people with overweight/obesity, the n-3 status (EPA + DHA) in erythrocytes or plasma phospholipids was inversely associated with type 2 diabetes in females (156, 157). In contrast to our findings, the meta-analysis and cohort studies did not demonstrate significant effects or relationships in males, and a previous trial in males with overweight showed that a higher n-3 index was associated with increased insulin sensitivity (84). Notably, none of the included studies in the meta-analysis directly compared females and males. In a more recent 12-week parallel-arm study specifically investigating sex-dependent effects of a DHA-enriched fish oil compared with corn oil, fasting insulin levels and HOMA-IR were significantly reduced in adults with abdominal obesity, but no significant sex differences were evident (88). In our trial, the significant increase in fasting glucose levels observed among males after the n-3 intervention agrees with a study in healthy males that found after 8 weeks of n-3 supplementation significantly elevated fasting and postprandial plasma glucose concentrations but no significant changes in insulin levels (89), while the DHA-enriched fish oil study showed no increase in fasting glucose levels and no sex-dependent effects (88). Although it is questionable, within the context of short-term increased n-3 PUFA intake, to extrapolate our observed changes in glucose and insulin markers to sex-specific effects on CVD endpoints, previous research has shown that fasting glucose, fasting insulin, and HOMA-IR were positively associated with incident CVD (158–161). A prospective cohort study in males found an odds ratio of 1.6 (95% CI: 1.1–2.3) for ischemic heart disease with each increase of 1 SD in fasting insulin concentration after adjustment for plasma TAGs, apoB, LDL-C, and HDL-C (158), while a meta-analysis of cohort studies and nested case–control studies showed that the relative risk of CVD per 1 SD increase in HOMA-IR was 1.37 (1.05, 1.80) for females and 1.41 (1.12, 1.77) for males (159).
Possible mechanisms explaining how n–3 PUFAs may affect insulin sensitivity have been previously described (82, 86, 162, 163), but reasons behind the sex differences we observed in fasting glucose levels and insulin sensitivity after the n-3 intervention are not clear (9, 82, 86, 89). Available data indicate that the sexual dimorphism may be caused or mediated, at least partly, by changes in steroid receptor signaling pathways (8, 164, 165). Previous n-3 PUFA intervention studies have shown inconsistent results regarding responses in sex hormones, including testosterone levels in males (166, 167). The n-3 intervention in the current study was not followed by significant changes or sex differences in testosterone levels or free androgen index, but we found among males a strong inverse association between relative changes in the n-3 index and testosterone levels adjusted for age, BMI, WC, and SHBG. In contrast, a DHA-enriched fish oil intervention (860 mg DHA and 120 mg EPA/d) for 12 weeks showed in males that change in testosterone was positively correlated with change in DHA in erythrocyte membranes after controlling for age and BMI (167).
Regardless of sex hormone levels, higher intakes of EPA and DHA may decrease the expression of the androgen receptor (AR) (168, 169) or modulate the actions of AR (170). This nuclear transcription factor plays a major role in glucose-insulin homeostasis (164, 171), potentially leading to sexually dimorphic effects on glucose, lipid, and energy metabolism in females and males (164, 165, 172). Additionally, the activity of 5-alpha-reductase, which catalyzes the conversion of testosterone to the more biologically active form dihydrotestosterone (DHT), may also be hampered by n-3 PUFAs and EPA’s 15-lipoxygenase metabolite (173). Furthermore, AR-induced activity has shown to be inhibited by small heterodimer partner (SHP) (174, 175), and animal data indicate that fish oil feeding increases the hepatic expression of SHP (176). This nuclear receptor regulates genes involved in, e.g., glucose and lipid metabolism (177–179) and may have sex-dependent and tissue-specific functions leading to disparate effects on TAG metabolism (180). Notably, the bidirectional modulation of glucose, lipid, and energy metabolism by testosterone, DHT, and AR-induced activity in females and males are among the most important sexually dimorphic aspects of metabolic regulation (164, 172, 181). While a reduction in testosterone, DHT, and/or AR signaling has been associated with abdominal obesity, hyperglycemia, insulin resistance, leptin resistance, metabolic syndrome, and type 2 diabetes in males, it has been related to an improved metabolic phenotype in females (164, 172, 181–184).
We observed similar reductions in TRLs, TAGs, and TRL-C in both sexes after the n-3 intervention, but the mechanisms may have been different in females and males because of the sexual dimorphism in glucose-insulin homeostasis, possibly affecting the balance between hepatic secretion, peripheral lipolysis, and vascular clearance of TRLs. Several lines of evidence have shown that NEFAs contribute the largest fraction to VLDL-TAG production in both the fasted and fed states in different metabolic conditions (185–187), yet indirectly via a cytosolic TAG pool (188–190). It has also been proposed that the TAG-lowering effect of n-3 PUFAs is best explained by an effect on the NEFA pool (187), which, at least in obesity, may partly occur via suppression of adipose tissue inflammation and thereby reducing intracellular lipolysis in adipocytes mediated by adipose triacylglycerol lipase and hormone-sensitive lipase. Recent experimental evidence from an animal model of aging, where rats received a fish oil intervention, point to potential mechanisms involved in n-3 PUFA-induced changes in adipose tissue function, inflammation, and lipolysis, and circulating lipid levels (191). This study showed that the expression of several proteins involved in lipid metabolism, notably ER lipid raft-associated protein 1 (Erlin1) and fatty acid binding protein 1 (Fabp1), were significantly changed (down-regulated Erlin1 and Fabp1) after the fish oil-derived n-3 PUFA intervention. However, the spatiotemporal sequence of events from hepatic NEFA influx to VLDL efflux is complex and incompletely understood (92, 138, 190, 192–195), and not all studies showing lowered TAGs concurrently found a significant reduction in NEFAs (89). While NEFAs in the present study increased non-significantly in females and decreased significantly in males after n-3 supplementation, TAG levels were significantly reduced in both sexes and to a similar degree, in contrast to a previous study of healthy adults showing a greater TAG-lowering effect of EPA + DHA in males compared with females (45).
In the present study, elevated insulin levels in males after the n-3 intervention probably reduced the NEFA flux from adipose tissues to the liver (196–198), as insulin suppresses the activity of adipose triacylglycerol lipase and hormone-sensitive lipase in insulin-sensitive states. By blocking lipolysis in adipocytes, as well as restraining intrahepatic ketogenesis and enhancing peripheral ketone body clearance, insulin also lowers circulating ketone body concentrations (146, 147, 199). This may explain the large (non-significant) sex differences in fasting serum levels of ketone bodies (increased in females and decreased in males), possibly reflecting sexual dimorphism in hepatic beta-oxidation of FAs and ketogenesis after n-3 supplementation. Additionally, we observed after the n-3 intervention that measures of body adiposity changed (non-significantly) in opposite directions among females (decreased) and males (increased), which may be related to sexual dimorphism in glucose and insulin responses.
A widely held hypothesis, which was initially based on short-term experiments with acute insulin infusions (194, 196, 200, 201), is that insulin inhibits hepatic production and secretion of larger, TAG-rich VLDLs (VLDL1) (141, 189, 202), and that this effect is partly lost if the liver becomes insulin resistant, leading to increased hepatic VLDL1-TAG output (139–141, 202–204), despite efficient suppression of serum NEFAs (205). However, others have questioned that hypertriacylglycerolemia is due to hepatic insulin resistance (196). Studies have demonstrated that higher insulin signaling in the long term (chronic hyperinsulinemia) does not inhibit but rather stimulates hepatic VLDL1-TAG synthesis and secretion, even when NEFA levels are lowered (138, 192–194, 196). The stimulation of de novo lipogenesis, decreased FA oxidation, and increased TAG esterification associated with chronic hyperinsulinemia presumably override the effects of reduced NEFA flux to the liver (194). Other lines of evidence point to the presence of differential insulin resistance between tissues or metabolic pathways within the liver, in which insulin fails to suppress gluconeogenesis but continues to activate lipogenesis, thereby increasing the hepatic output of both glucose and TAGs (206), although recent animal and human data have questioned the importance of selective hepatic insulin resistance (207, 208). Furthermore, elevated glucose levels have been shown to decrease FA oxidation and increase hepatic TAG synthesis and secretion (194, 209–212), of which a large proportion may originate from other sources than plasma NEFAs in hyperglycemic hyperinsulinemic individuals (194, 209). In an analysis of factors that predicted overproduction of TRLs, only liver fat and plasma glucose were identified as significant predictors in the multiple regression modeling (140). Importantly, recent experiments with primary hepatocytes, involving diverse in vivo and in vitro assays, have shown that insulin enhances in the fed state VLDL1-TAG production and secretion via the phosphatidic acid (PA) pathway (92, 213). This response is downregulated when fasting causes insulin levels to drop, thereby preventing the release of TAGs into an already NEFA-rich bloodstream in the fasted state (195). These experiments showed that higher insulin signaling elevated dramatically in lipid droplets the content of PA (generated by phospholipase-D1 after its activation by ADP-ribosylation factor 1), and that PA signaling induced transport of lipid droplets to smooth endoplasmic reticulum inside hepatocytes, where lipid droplets were catabolized to produce lipoproteins (92, 213).
Based on these lines of evidence, and other supporting studies (189, 190, 214, 215), it is tempting to speculate that the higher fasting glucose and insulin levels observed in the current study among males after n-3 supplementation, possibly in part caused by modulated AR expression or actions, (1) reduced the NEFA flux from adipose tissues, but (2) enhanced the output of hepatic VLDL1-TAG, possibly recruited at least in part from an intracellular phospholipid pool via the PA pathway (190, 214), such as PA-rich activated LDs (92, 213), and (3) stimulated lipolytic conversion of TRLs and subsequent FA uptake and accumulation of TAGs in peripheral tissues by glucose-, insulin-, and eventually n-3 PUFA-stimulated actions of lipoprotein lipase (LPL) (198, 210, 216–221). These effects would ultimately lead to reduced circulating TRL, TAG, and NEFA levels, as well as ketone bodies, and increased adiposity. In females, the TAG-lowering effects of n-3 PUFAs may be partly explained by reduced hepatic secretion of TRLs involving mechanisms normally attributed to EPA and DHA, i.e., less available substrates for TAG synthesis, decreased activity of TAG-synthesizing enzymes, and/or stimulated degradation of apoB-100 (49, 187, 222–227). It is known that females, compared with males, secrete more TAG-rich VLDLs from the liver (228). This is matched with accelerated VLDL1-TAG clearance rates, collectively contributing to lower plasma VLDL1-TAG levels in obesity in females (228). If n-3 PUFAs increase the hepatic expression of SHP (176), which inhibits estrogen receptor-induced activity (178, 180, 229), then this mechanism may also contribute to lower TAG levels, as estrogen has been shown to increase hepatic production of VLDL1-TAG (230).
The sex-specific responses in LDL and HDL measures after the n-3 intervention may also be related to the sexual dimorphism in glucose-insulin homeostasis and associated changes in VLDL metabolism. Higher circulating TRL levels increase the activity of cholesteryl ester transfer protein (CETP) (202, 231–234), which exchanges TAGs from TRLs for cholesteryl esters from LDLs and HDLs (231–234). LPL promotes the exchange of lipids between lipoproteins (218, 235), and the compositional changes make LDL and HDL particles better substrates for hepatic lipase (HL) and to a lesser degree LPL (232–234, 236). HL activity is normally higher in males than females (43, 237, 238) and is elevated in abdominal obesity and hyperglycemic and insulin-resistant states (234, 237–240). Increased hepatic production of TRLs and subsequent CETP-mediated lipid exchange and lipolytic actions of LPL and HL typically result in fewer larger cholesterol-enriched LDLs and HDLs, more numerous smaller TAG-enriched LDLs and HDLs, and smaller mean LDL and HDL particle sizes (141, 202, 232, 233, 236, 241–243). The effects of n-3 PUFAs on these pathways may be modulated by glucose and insulin responses and thus explain some of the sex differences observed in the present study, of which the sexual dimorphism in small, dense LDLs may translate into a lower CVD risk in females compared with males (244–248). Finally, while multiple lines of evidence point to several specific mechanisms possibly explaining the sex-dependent responses we observed after n-3 supplementation, we are not aware of experimental data that could indicate reasons behind the sex differences in VLDL-related measures or Lp(a) after the n-6 intervention.
Strengths and limitations of our crossover trial have been discussed previously (83). Other strengths related specifically to the present work include the supplementation of similar PUFA doses per kg body weight across sexes. Moreover, we used cLMMs to obtain baseline-adjusted effect estimates, and we controlled for cross-level bias in the mixed modeling. Among additional limitations, we were not able to recruit the same number of females and males, and the underrepresentation of females in this study may have limited our ability to detect some clinically significant sex differences. Furthermore, our sample included both pre- and postmenopausal females and we lacked data on the time of last menstruation. However, controlling for menopausal status in the mixed modeling had a negligible impact on effect and precision estimates. It should also be recognized that the present study was exploratory since sex-specific responses were not defined as a primary outcome. Finally, while our results show that sex-dependent responses to PUFA supplementation occur even in people with a relatively high pre-treatment n-3 index and low proportion of LA in RBCMs (83), which is typically observed in Nordic countries (59, 249), the findings may not be generalizable to people with a lower EPA + DHA/LA ratio in blood or with other anthropometric traits.
5. Conclusion
Our findings demonstrate sexually dimorphic responses after n-3 (but not n-6) supplementation in circulating markers of glycemic control/insulin sensitivity, which improved in females but worsened in males. These differential changes in glucose-insulin homeostasis may partly be related to the sex differences we observed in several components of the lipoprotein-lipid profile following the n-3 intervention. After n-6 supplementation, Lp(a) decreased significantly in females and increased non-significantly in males. These insights motivate greater attention to individual differences in PUFA-mediated regulation of glucose and lipid metabolism, at least in part dependent on differences in body composition and sex-specific hormonal signaling.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors upon request, without undue reservation.
Ethics statement
This study involved human participants and was reviewed and approved by The Regional Committee for Medical and Health Research Ethics (2014/2336/REK South-East). The patients/ participants provided their written informed consent to participate in this study.
Author contributions
JL-B, EG, BB, JS, JN, RB, ER, and ON: conceptualization. JL-B, EG, PB, BB, JS, JN, RB, ER, and ON: methodology. JL-B: software, formal analysis, and writing—original draft preparation. JL-B, EG, PB, JN, ER, GM, SD, and ON: validation. JL-B, EG, PB, BB, ES, and ER: investigation. JN, ER, and ON: resources. JL-B and EG: data curation and visualization. JL-B, EG, PB, BB, ES, JN, GM, SD, and ON: writing—review and editing. ER, GM, SD, and ON: supervision. JL-B, EG, BB, ES, JS, JN, RB, ER, and ON: project administration. ER and ON: funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the Western Norway Regional Health Authority (Bergen), Trond Mohn Foundation (Bergen), and the University of Bergen.
Acknowledgments
We thank all participants of the study for their dedication to the project. We thank the Norwegian company Pharmatech AS (NO-1664 Rolvsøy) for providing the supplementation products used in the study. We also thank our project nurse, Hilde Netland, for her contributions during the planning and implementation of the project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Glossary
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2023.1020678/full#supplementary-material
Footnotes
References
1. Roth, GA, Abate, D, Abate, KH, Abay, SM, Abbafati, C, Abbasi, N, et al. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the global burden of disease study 2017. Lancet. (2018) 392:1736–88. doi: 10.1016/S0140-6736(18)32203-7
2. Vaidya, D, Becker, DM, Bittner, V, Mathias, RA, and Ouyang, P. Ageing, menopause, and ischaemic heart disease mortality in England, Wales, and the United States: modelling study of national mortality data. BMJ. (2011) 343:d5170. doi: 10.1136/bmj.d5170
3. Walli-Attaei, M, Joseph, P, Rosengren, A, Chow, CK, Rangarajan, S, Lear, SA, et al. Variations between women and men in risk factors, treatments, cardiovascular disease incidence, and death in 27 high-income, middle-income, and low-income countries (pure): a prospective cohort study. Lancet. (2020) 396:97–109. doi: 10.1016/s0140-6736(20)30543-2
4. Virani, SS, Alonso, A, Aparicio, HJ, Benjamin, EJ, Bittencourt, MS, Callaway, CW, et al. Heart disease and stroke statistics – 2021 update: a report from the American Heart Association. Circulation. (2021) 143:e254–743. doi: 10.1161/CIR.0000000000000950
5. Knopp, RH, Paramsothy, P, Retzlaff, BM, Fish, B, Walden, C, Dowdy, A, et al. Gender differences in lipoprotein metabolism and dietary response: basis in hormonal differences and implications for cardiovascular disease. Curr Atheroscler Rep. (2005) 7:472–9. doi: 10.1007/s11883-005-0065-6
6. Decsi, T, and Kennedy, K. Sex-specific differences in essential fatty acid metabolism. Am J Clin Nutr. (2011) 94:1914S–9S. doi: 10.3945/ajcn.110.000893
7. Regitz-Zagrosek, V, and Kararigas, G. Mechanistic pathways of sex differences in cardiovascular disease. Physiol Rev. (2017) 97:1–37. doi: 10.1152/physrev.00021.2015
8. Palmisano, BT, Zhu, L, Eckel, RH, and Stafford, JM. Sex differences in lipid and lipoprotein metabolism. Mol Metab. (2018) 15:45–55. doi: 10.1016/j.molmet.2018.05.008
9. Goossens, GH, Jocken, JWE, and Blaak, EE. Sexual dimorphism in Cardiometabolic health: the role of adipose tissue. Muscle and Liver Nat Rev Endocrinol. (2021) 17:47–66. doi: 10.1038/s41574-020-00431-8
10. Childs, CE, Romeu-Nadal, M, Burdge, GC, and Calder, PC. Gender differences in the N-3 fatty acid content of tissues. Proc Nutr Soc. (2008) 67:19–27. doi: 10.1017/S0029665108005983
11. Lohner, S, Fekete, K, Marosvölgyi, T, and Decsi, T. Gender differences in the long-chain polyunsaturated fatty acid status: systematic review of 51 publications. Ann Nutr Metab. (2013) 62:98–112. doi: 10.1159/000345599
12. Fisk, HL, Irvine, M, Miles, EA, Lietz, G, Mathers, JC, Packard, CJ, et al. Association of oily fish intake, sex, age, BMI and APOE genotype with plasma long-chain N-3 fatty acid composition. Br J Nutr. (2018) 120:23–32. doi: 10.1017/S000711451800106X
13. Coltell, O, Sorlí, JV, Asensio, EM, Barragán, R, González, JI, Giménez-Alba, IM, et al. Genome-wide association study for serum omega-3 and omega-6 polyunsaturated fatty acids: exploratory analysis of the sex-specific effects and dietary modulation in Mediterranean subjects with metabolic syndrome. Nutrients. (2020) 12:310. doi: 10.3390/nu12020310
14. Lyu, L-C, Yeh, C-Y, Lichtenstein, AH, Li, Z, Ordovas, JM, and Schaefer, EJ. Association of sex, adiposity, and diet with HDL subclasses in middle-aged Chinese. Am J Clin Nutr. (2001) 74:64–71. doi: 10.1093/ajcn/74.1.64
15. Howe, PRC, Buckley, JD, Murphy, KJ, Pettman, T, Milte, C, and Coates, AM. Relationship between erythrocyte omega-3 content and obesity is gender dependent. Nutrients. (2014) 6:1850–60. doi: 10.3390/nu6051850
16. Hatch-Stein, JA, Kelly, A, Gidding, SS, Zemel, BS, and Magge, SN. Sex differences in the associations of visceral adiposity, homeostatic model assessment of insulin resistance, and body mass index with lipoprotein subclass analysis in obese adolescents. J Clin Lipidol. (2016) 10:757–66. doi: 10.1016/j.jacl.2016.02.002
17. Nakanishi, S, Nagano, C, Miyahara, M, and Sawano, F. Sex differences in the association between the eicosapentaenoic acid/arachidonic acid ratio and the visceral fat area among patients with type 2 diabetes. Intern Med. (2016) 55:1269–74. doi: 10.2169/internalmedicine.55.5624
18. Tchernof, A, and Després, J-P. Pathophysiology of human visceral obesity: an update. Physiol Rev. (2013) 93:359–404. doi: 10.1152/physrev.00033.2011
19. Nordestgaard, BG. Triglyceride-rich lipoproteins and atherosclerotic cardiovascular disease: new insights from epidemiology, genetics, and biology. Circ Res. (2016) 118:547–63. doi: 10.1161/CIRCRESAHA.115.306249
20. Sandesara, PB, Virani, SS, Fazio, S, and Shapiro, MD. The forgotten lipids: triglycerides, remnant cholesterol, and atherosclerotic cardiovascular disease risk. Endocr Rev. (2018) 40:537–57. doi: 10.1210/er.2018-00184
21. Hussain, A, Ballantyne, CM, Saeed, A, and Virani, SS. Triglycerides and ASCVD risk reduction: recent insights and future directions. Curr Atheroscler Rep. (2020) 22:25. doi: 10.1007/s11883-020-00846-8
22. Miller, M, Stone, NJ, Ballantyne, C, Bittner, V, Criqui, MH, Ginsberg, HN, et al. Triglycerides and cardiovascular disease: a scientific statement from the American Heart Association. Circulation. (2011) 123:2292–333. doi: 10.1161/CIR.0b013e3182160726
23. Aberra, T, Peterson, ED, Pagidipati, NJ, Mulder, H, Wojdyla, DM, Philip, S, et al. The association between triglycerides and incident cardiovascular disease: what is “optimal”? J Clin Lipidol. (2020) 14:438–47.e3. doi: 10.1016/j.jacl.2020.04.009
24. Freedman, DS, Otvos, JD, Jeyarajah, EJ, Barboriak, JJ, Anderson, AJ, and Walker, JA. Relation of lipoprotein subclasses as measured by proton nuclear magnetic resonance spectroscopy to coronary artery disease. Arterioscler Thromb Vasc Biol. (1998) 18:1046–53. doi: 10.1161/01.atv.18.7.1046
25. Blake, GJ, Otvos, JD, Rifai, N, and Ridker, PM. Low-density lipoprotein particle concentration and size as determined by nuclear magnetic resonance spectroscopy as predictors of cardiovascular disease in women. Circulation. (2002) 106:1930–7. doi: 10.1161/01.cir.0000033222.75187.b9
26. Colhoun, HM, Otvos, JD, Rubens, MB, Taskinen, MR, Underwood, SR, and Fuller, JH. Lipoprotein subclasses and particle sizes and their relationship with coronary artery calcification in men and women with and without type 1 diabetes. Diabetes. (2002) 51:1949–56. doi: 10.2337/diabetes.51.6.1949
27. Alaupovic, P. The concept of apolipoprotein-defined lipoprotein families and its clinical significance. Curr Atheroscler Rep. (2003) 5:459–67. doi: 10.1007/s11883-003-0036-8
28. Cromwell, WC, Otvos, JD, Keyes, MJ, Pencina, MJ, Sullivan, L, Vasan, RS, et al. LDL particle number and risk of future cardiovascular disease in the framingham offspring study: implications for LDL management. J Clin Lipidol. (2007) 1:583–92. doi: 10.1016/j.jacl.2007.10.001
29. Mora, S, Otvos, JD, Rifai, N, Rosenson, RS, Buring, JE, and Ridker, PM. Lipoprotein particle profiles by nuclear magnetic resonance compared with standard lipids and apolipoproteins in predicting incident cardiovascular disease in women. Circulation. (2009) 119:931–9. doi: 10.1161/circulationaha.108.816181
30. Musunuru, K, Orho-Melander, M, Caulfield, MP, Li, S, Salameh, WA, Reitz, RE, et al. Ion mobility analysis of lipoprotein subfractions identifies three independent axes of cardiovascular risk. Arterioscler Thromb Vasc Biol. (2009) 29:1975–80. doi: 10.1161/atvbaha.109.190405
31. Krauss, RM. Lipoprotein subfractions and cardiovascular disease risk. Curr Opin Lipidol. (2010) 21:305–11. doi: 10.1097/MOL.0b013e32833b7756
32. Arsenault, BJ, Boekholdt, SM, and Kastelein, JJP. Lipid parameters for measuring risk of cardiovascular disease. Nat Rev Cardiol. (2011) 8:197–206. doi: 10.1038/nrcardio.2010.223
33. Wyler von Ballmoos, MC, Haring, B, and Sacks, FM. The risk of cardiovascular events with increased apolipoprotein CIII: a systematic review and meta-analysis. J Clin Lipidol. (2015) 9:498–510. doi: 10.1016/j.jacl.2015.05.002
34. Li, Z, McNamara, JR, Fruchart, JC, Luc, G, Bard, JM, Ordovas, JM, et al. Effects of gender and menopausal status on plasma lipoprotein subspecies and particle sizes. J Lipid Res. (1996) 37:1886–96. doi: 10.1016/s0022-2275(20)37553-2
35. Pascot, A, Lemieux, I, Bergeron, J, Tremblay, A, Nadeau, A, Prud'homme, D, et al. HDL particle size: a marker of the gender difference in the metabolic risk profile. Atherosclerosis. (2002) 160:399–406. doi: 10.1016/s0021-9150(01)00579-2
36. Freedman, DS, Otvos, JD, Jeyarajah, EJ, Shalaurova, I, Cupples, LA, Parise, H, et al. Sex and age differences in lipoprotein subclasses measured by nuclear magnetic resonance spectroscopy: the framingham study. Clin Chem. (2004) 50:1189–200. doi: 10.1373/clinchem.2004.032763
37. Vaidya, D, Dobs, A, Gapstur, SM, Golden, SH, Hankinson, A, Liu, K, et al. The Association of endogenous sex hormones with lipoprotein subfraction profile in the multi-ethnic study of atherosclerosis. Metabolism. (2008) 57:782–90. doi: 10.1016/j.metabol.2008.01.019
38. Lu, Y, Zhou, S, Dreyer, RP, Caulfield, M, Spatz, ES, Geda, M, et al. Sex differences in lipid profiles and treatment utilization among young adults with acute myocardial infarction: results from the VIRGO study. Am Heart J. (2017) 183:74–84. doi: 10.1016/j.ahj.2016.09.012
39. Sancho-Rodríguez, N, Avilés-Plaza, FV, Granero-Fernández, E, Hernández-Martínez, AM, Albaladejo-Otón, MD, Martínez-Mernández, P, et al. Observational study of lipid profile and LDL particle size in patients with metabolic syndrome. Lipids Health Dis. (2011) 10:162. doi: 10.1186/1476-511X-10-162
40. Hansson, P, Holven, KB, Øyri, LKL, Brekke, HK, Gjevestad, GO, Thoresen, M, et al. Sex differences in postprandial responses to different dairy products on lipoprotein subclasses: a randomised controlled cross-over trial. Br J Nutr. (2019) 122:780–9. doi: 10.1017/S0007114519001429
41. McNamara, JR, Campos, H, Ordovas, JM, Peterson, J, Wilson, PW, and Schaefer, EJ. Effect of gender, age, and lipid status on low density lipoprotein subfraction distribution: results from the framingham offspring study. Arteriosclerosis. (1987) 7:483–90. doi: 10.1161/01.atv.7.5.483
42. Williams, PT, Krauss, RM, Vranizan, KM, Stefanick, ML, Wood, PD, and Lindgren, FT. Associations of lipoproteins and apolipoproteins with gradient gel electrophoresis estimates of high density lipoprotein subfractions in men and women. Arterioscler Thromb. (1992) 12:332–40. doi: 10.1161/01.atv.12.3.332
43. Tan, CE, Foster, L, Caslake, MJ, Bedford, D, Watson, TD, McConnell, M, et al. Relations between plasma lipids and postheparin plasma lipases and VLDL and LDL subfraction patterns in normolipemic men and women. Arterioscler Thromb Vasc Biol. (1995) 15:1839–48. doi: 10.1161/01.atv.15.11.1839
44. Tilly-Kiesi, M, Lichtenstein, AH, Joven, J, Vilella, E, Cheung, MC, Carrasco, WV, et al. Impact of gender on the metabolism of apolipoprotein A-I in HDL subclasses LpAI and LpAI:AII in older subjects. Arterioscler Thromb Vasc Biol. (1997) 17:3513–8. doi: 10.1161/01.atv.17.12.3513
45. Caslake, MJ, Miles, EA, Kofler, BM, Lietz, G, Curtis, P, Armah, CK, et al. Effect of sex and genotype on cardiovascular biomarker response to fish oils: the FINGEN study. Am J Clin Nutr. (2008) 88:618–29. doi: 10.1093/ajcn/88.3.618
46. Matthan, NR, Jalbert, SM, Barrett, PHR, Dolnikowski, GG, Schaefer, EJ, and Lichtenstein, AH. Gender-specific differences in the kinetics of nonfasting TRL, IDL, and LDL apolipoprotein B-100 in men and premenopausal women. Arterioscler Thromb Vasc Biol. (2008) 28:1838–43. doi: 10.1161/ATVBAHA.108.163931
47. Anagnostis, P, Stevenson, JC, Crook, D, Johnston, DG, and Godsland, IF. Effects of menopause, gender and age on lipids and high-density lipoprotein cholesterol subfractions. Maturitas. (2015) 81:62–8. doi: 10.1016/j.maturitas.2015.02.262
48. Klingel, SL, Roke, K, Hidalgo, B, Aslibekyan, S, Straka, RJ, An, P, et al. Sex differences in blood HDL-C, the total cholesterol/HDL-C ratio, and palmitoleic acid are not associated with variants in common candidate genes. Lipids. (2017) 52:969–80. doi: 10.1007/s11745-017-4307-5
49. Mozaffarian, D, and Wu, JHY. Omega-3 fatty acids and cardiovascular disease: effects on risk factors, molecular pathways, and clinical events. J Am Coll Cardiol. (2011) 58:2047–67. doi: 10.1016/j.jacc.2011.06.063
50. Casula, M, Soranna, D, Catapano, AL, and Corrao, G. Long-term effect of high dose omega-3 fatty acid supplementation for secondary prevention of cardiovascular outcomes: a meta-analysis of randomized, double blind, placebo controlled trials. Atheroscler Suppl. (2013) 14:243–51. doi: 10.1016/S1567-5688(13)70005-9
51. AbuMweis, S, Jew, S, Tayyem, R, and Agraib, L. Eicosapentaenoic acid and docosahexaenoic acid containing supplements modulate risk factors for cardiovascular disease: a meta-analysis of randomised placebo-control human clinical trials. J Hum Nutr Diet. (2018) 31:67–84. doi: 10.1111/jhn.12493
52. Bhatt, DL, Steg, PG, Miller, M, Brinton, EA, Jacobson, TA, Ketchum, SB, et al. Cardiovascular risk reduction with icosapent ethyl for hypertriglyceridemia. N Engl J Med. (2019) 380:11–22. doi: 10.1056/NEJMoa1812792
53. Hu, Y, Hu, FB, and Manson, JE. Marine omega-3 supplementation and cardiovascular disease: an updated meta-analysis of 13 randomized controlled trials involving 127 477 participants. J Am Heart Assoc. (2019) 8:e013543. doi: 10.1161/JAHA.119.013543
54. Innes, JK, and Calder, PC. Marine omega-3 (N-3) fatty acids for cardiovascular health: an update for 2020. Int J Mol Sci. (2020) 21:1362. doi: 10.3390/ijms21041362
55. Harris, WS, Tintle, NL, Imamura, F, Qian, F, Korat, AVA, Marklund, M, et al. Blood N-3 fatty acid levels and Total and cause-specific mortality from 17 prospective studies. Nat Commun. (2021) 12:2329. doi: 10.1038/s41467-021-22370-2
56. Alfaddagh, A, Kapoor, K, Dardari, ZA, Bhatt, DL, Budoff, MJ, Nasir, K, et al. Omega-3 fatty acids, subclinical atherosclerosis, and cardiovascular events: implications for primary prevention. Atherosclerosis. (2022) 353:11–9. doi: 10.1016/j.atherosclerosis.2022.06.1018
57. Farvid, MS, Ding, M, Pan, A, Sun, Q, Chiuve, SE, Steffen, LM, et al. Dietary linoleic acid and risk of coronary heart disease: a systematic review and meta-analysis of prospective cohort studies. Circulation. (2014) 130:1568–78. doi: 10.1161/circulationaha.114.010236
58. Sacks, FM, Lichtenstein, AH, Wu, JHY, Appel, LJ, Creager, MA, Kris-Etherton, PM, et al. Dietary fats and cardiovascular disease: a presidential advisory from the American Heart Association. Circulation. (2017) 136:e1–e23. doi: 10.1161/CIR.0000000000000510
59. Marklund, M, Wu, JHY, Imamura, F, Gobbo, LCD, Fretts, A, Goede, J, et al. Biomarkers of dietary omega-6 fatty acids and incident cardiovascular disease and mortality: an individual-level pooled analysis of 30 cohort studies. Circulation. (2019) 139:2422–36. doi: 10.1161/CIRCULATIONAHA.118.038908
60. Marangoni, F, Agostoni, C, Borghi, C, Catapano, AL, Cena, H, Ghiselli, A, et al. Dietary linoleic acid and human health: focus on cardiovascular and cardiometabolic effects. Atherosclerosis. (2020) 292:90–8. doi: 10.1016/j.atherosclerosis.2019.11.018
61. Lewington, S, Whitlock, G, Clarke, R, Sherliker, P, Emberson, J, Halsey, J, et al. Blood cholesterol and vascular mortality by age, sex, and blood pressure: a meta-analysis of individual data from 61 prospective studies with 55 000 vascular deaths. Lancet. (2007) 370:1829–39. doi: 10.1016/S0140-6736(07)61778-4
62. Froyen, E, and Burns-Whitmore, B. The effects of linoleic acid consumption on lipid risk markers for cardiovascular disease in healthy individuals: a review of human intervention trials. Nutrients. (2020) 12:2329. doi: 10.3390/nu12082329
63. Harris, WS. N-3 fatty acids and serum lipoproteins: human studies. Am J Clin Nutr. (1997) 65:1645S–54S. doi: 10.1093/ajcn/65.5.1645S
64. Chan, DC, Watts, GF, Mori, TA, Barrett, PHR, Redgrave, TG, and Beilin, LJ. Randomized controlled trial of the effect of N–3 fatty acid supplementation on the metabolism of apolipoprotein B-100 and chylomicron remnants in men with visceral obesity. Am J Clin Nutr. (2003) 77:300–7. doi: 10.1093/ajcn/77.2.300
65. Skulas-Ray, AC, Kris-Etherton, PM, Harris, WS, Vanden Heuvel, JP, Wagner, PR, and West, SG. Dose-response effects of omega-3 fatty acids on triglycerides, inflammation, and endothelial function in healthy persons with moderate hypertriglyceridemia. Am J Clin Nutr. (2011) 93:243–52. doi: 10.3945/ajcn.110.003871
66. Wei, MY, and Jacobson, TA. Effects of eicosapentaenoic acid versus docosahexaenoic acid on serum lipids: a systematic review and meta-analysis. Curr Atheroscler Rep. (2011) 13:474–83. doi: 10.1007/s11883-011-0210-3
67. Leslie, MA, Cohen, DJA, Liddle, DM, Robinson, LE, and Ma, DWL. A review of the effect of omega-3 polyunsaturated fatty acids on blood triacylglycerol levels in normolipidemic and borderline hyperlipidemic individuals. Lipids Health Dis. (2015) 14:53. doi: 10.1186/s12944-015-0049-7
68. Kris-Etherton, PM, and Yu, S. Individual fatty acid effects on plasma lipids and lipoproteins: human studies. Am J Clin Nutr. (1997) 65:1628S–44S. doi: 10.1093/ajcn/65.5.1628S
69. Mensink, RP, Zock, PL, Kester, AD, and Katan, MB. Effects of dietary fatty acids and carbohydrates on the ratio of serum total to HDL cholesterol and on serum lipids and apolipoproteins: a meta-analysis of 60 controlled trials. Am J Clin Nutr. (2003) 77:1146–55. doi: 10.1093/ajcn/77.5.1146
70. Maki, KC, Eren, F, Cassens, ME, Dicklin, MR, and Davidson, MH. Ω-6 polyunsaturated fatty acids and cardiometabolic health: current evidence, controversies, and research gaps. Adv Nutr. (2018) 9:688–700. doi: 10.1093/advances/nmy038
71. Aung, T, Halsey, J, Kromhout, D, Gerstein, HC, Marchioli, R, Tavazzi, L, et al. Associations of omega-3 fatty acid supplement use with cardiovascular disease risks: meta-analysis of 10 trials involving 77 917 individuals. JAMA Cardiol. (2018) 3:225–34. doi: 10.1001/jamacardio.2017.5205
72. Nicholls, SJ, Lincoff, AM, Garcia, M, Bash, D, Ballantyne, CM, Barter, PJ, et al. Effect of high-dose omega-3 fatty acids vs corn oil on major adverse cardiovascular events in patients at high cardiovascular risk: the strength randomized clinical trial. JAMA. (2020) 324:2268–80. doi: 10.1001/jama.2020.22258
73. Sharma, G, Martin, SS, and Blumenthal, RS. Effects of omega-3 fatty acids on major adverse cardiovascular events: what matters most: the drug, the dose, or the placebo? JAMA. (2020) 324:2262–4. doi: 10.1001/jama.2020.22387
74. Moosavi, D, Vuckovic, I, Kunz, HE, and Lanza, IR. A randomized trial of omega-3 fatty acid supplementation and circulating lipoprotein subclasses in healthy older adults. J Nutr. (2022) 152:1675–89. doi: 10.1093/jn/nxac084
75. Ramsden, CE, Hibbeln, JR, Majchrzak, SF, and Davis, JM. N-6 fatty acid-specific and mixed polyunsaturate dietary interventions have different effects on CHD risk: a meta-analysis of randomised controlled trials. Br J Nutr. (2010) 104:1586–600. doi: 10.1017/S0007114510004010
76. Hamley, S. The effect of replacing saturated fat with mostly N-6 polyunsaturated fat on coronary heart disease: a meta-analysis of randomised controlled trials. Nutr J. (2017) 16:30. doi: 10.1186/s12937-017-0254-5
77. Hooper, L, Al-Khudairy, L, Abdelhamid, AS, Rees, K, Brainard, JS, Brown, TJ, et al. Omega-6 fats for the primary and secondary prevention of cardiovascular disease. Cochrane Database Syst Rev. (2018) 11:CD011094. doi: 10.1002/14651858.CD011094.pub4
78. Abdelhamid, AS, Brown, TJ, Brainard, JS, Biswas, P, Thorpe, GC, Moore, HJ, et al. Omega-3 fatty acids for the primary and secondary prevention of cardiovascular disease. Cochrane Database Syst Rev. (2020) 3:CD003177. doi: 10.1002/14651858.CD003177.pub5
79. Wang, DD, and Hu, FB. Dietary fat and risk of cardiovascular disease: recent controversies and advances. Annu Rev Nutr. (2017) 37:423–46. doi: 10.1146/annurev-nutr-071816-064614
80. Marklund, M, Leander, K, Vikström, M, Laguzzi, F, Gigante, B, Sjögren, P, et al. Polyunsaturated fat intake estimated by circulating biomarkers and risk of cardiovascular disease and all-cause mortality in a population-based cohort of 60-year-old men and women. Circulation. (2015) 132:586–94. doi: 10.1161/CIRCULATIONAHA.115.015607
81. Innes, JK, and Calder, PC. The differential effects of eicosapentaenoic acid and docosahexaenoic acid on cardiometabolic risk factors: a systematic review. Int J Mol Sci. (2018) 19:532. doi: 10.3390/ijms19020532
82. Thota, RN, Ferguson, JJA, Abbott, KA, Dias, CB, and Garg, ML. Science behind the cardio-metabolic benefits of omega-3 polyunsaturated fatty acids: biochemical effects vs. clinical outcomes. Food Funct. (2018) 9:3576–96. doi: 10.1039/c8fo00348c
83. Grytten, E, Laupsa-Borge, J, Bohov, P, Bjørndal, B, Strand, E, Skorve, J, et al. Changes in lipoprotein particle subclasses, standard lipids, and apolipoproteins after supplementation with N-3 or N-6 PUFAs in abdominal obesity: a randomized double-blind crossover study. Clin Nutr. (2021) 40:2556–75. doi: 10.1016/j.clnu.2021.03.040
84. Albert, BB, Derraik, JGB, Brennan, CM, Biggs, JB, Smith, GC, Garg, ML, et al. Higher omega-3 index is associated with increased insulin sensitivity and more favourable metabolic profile in middle-aged overweight men. Sci Rep. (2014) 4:6697. doi: 10.1038/srep06697
85. Flachs, P, Rossmeisl, M, and Kopecky, J. The effect of N-3 fatty acids on glucose homeostasis and insulin sensitivity. Physiol Res. (2014) 63:S93–S118. doi: 10.33549/physiolres.932715
86. Abbott, KA, Burrows, TL, Thota, RN, Acharya, S, and Garg, ML. Do Ω-3 Pufas affect insulin resistance in a sex-specific manner? A systematic review and meta-analysis of randomized controlled trials. Am J Clin Nutr. (2016) 104:1470–84. doi: 10.3945/ajcn.116.138172
87. Lalia, AZ, and Lanza, IR. Insulin-sensitizing effects of omega-3 fatty acids: lost in translation? Nutrients. (2016) 8:329. doi: 10.3390/nu8060329
88. Abbott, KA, Burrows, TL, Acharya, S, Thota, RN, and Garg, ML. Dha-enriched Fish oil reduces insulin resistance in overweight and obese adults. Prostaglandins Leukot Essent Fatty Acids. (2020) 159:102154. doi: 10.1016/j.plefa.2020.102154
89. Green, CJ, Pramfalk, C, Charlton, CA, Gunn, PJ, Cornfield, T, Pavlides, M, et al. Hepatic de novo lipogenesis is suppressed and fat oxidation is increased by omega-3 fatty acids at the expense of glucose metabolism. BMJ Open Diabetes Res Care. (2020) 8:e000871. doi: 10.1136/bmjdrc-2019-000871
90. Hodson, L, and Fielding, BA. Trafficking and partitioning of fatty acids: the transition from fasted to fed state. Clin Lipidol. (2010) 5:131–44. doi: 10.2217/clp.09.72
91. Varlamov, O, Bethea, CL, and Roberts, CT. Sex-specific differences in lipid and glucose metabolism. Front Endocrinol. (2015) 5:241. doi: 10.3389/fendo.2014.00241
92. Kumar, M, Ojha, S, Rai, P, Joshi, A, Kamat, SS, and Mallik, R. Insulin activates intracellular transport of lipid droplets to release triglycerides from the liver. J Cell Biol. (2019) 218:3697–713. doi: 10.1083/jcb.201903102
93. Raeisi-Dehkordi, H, Amiri, M, Zimorovat, A, Moghtaderi, F, Zarei, S, Forbes, SC, et al. Canola oil compared with sesame and sesame-canola oil on glycaemic control and liver function in patients with type 2 diabetes: a three-way randomized triple-blind cross-over trial. Diabetes Metab Res Rev. (2021) 37:e3399. doi: 10.1002/dmrr.3399
94. Wang, Y, and Xu, D. Effects of aerobic exercise on lipids and lipoproteins. Lipids Health Dis. (2017) 16:132. doi: 10.1186/s12944-017-0515-5
95. Ogawa, H, Fujitani, K, Tsujinaka, T, Imanishi, K, Shirakata, H, Kantani, A, et al. Inbody 720 as a new method of evaluating visceral obesity. Hepatogastroenterology. (2011) 58:42–4.
96. Nagai, M, Komiya, H, Mori, Y, Ohta, T, Kasahara, Y, and Ikeda, Y. Estimating visceral fat area by multifrequency bioelectrical impedance. Diabetes Care. (2010) 33:1077–9. doi: 10.2337/dc09-1099
97. Lee, LC, Hsu, PS, Hsieh, KC, Chen, YY, Chu, LP, Lu, HK, et al. Standing 8-electrode bioelectrical impedance analysis as an alternative method to estimate visceral fat area and body fat mass in athletes. Int J Gen Med. (2021) 14:539–48. doi: 10.2147/ijgm.S281418
98. Strand, E, Bjorndal, B, Nygard, O, Burri, L, Berge, C, Bohov, P, et al. Long-term treatment with the pan-PPAR agonist tetradecylthioacetic acid or fish oil is associated with increased cardiac content of N-3 fatty acids in rat. Lipids Health Dis. (2012) 11:82. doi: 10.1186/1476-511X-11-82
99. Hodson, L, Skeaff, CM, and Fielding, BA. Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res. (2008) 47:348–80. doi: 10.1016/j.plipres.2008.03.003
100. Mocking, RJT, Assies, J, Lok, A, Ruhé, HG, Koeter, MWJ, Visser, I, et al. Statistical methodological issues in handling of fatty acid data: percentage or concentration, imputation and indices. Lipids. (2012) 47:541–7. doi: 10.1007/s11745-012-3665-2
101. Miura, K, Hughes, MCB, Ungerer, JPJ, Smith, DD, and Green, AC. Absolute versus relative measures of plasma fatty acids and health outcomes: example of phospholipid omega-3 and omega-6 fatty acids and all-cause mortality in women. Eur J Nutr. (2018) 57:713–22. doi: 10.1007/s00394-016-1358-y
102. Otvos, JD, Jeyarajah, EJ, and Cromwell, WC. Measurement issues related to lipoprotein heterogeneity. Am J Cardiol. (2002) 90:22–9. doi: 10.1016/s0002-9149(02)02632-2
103. Jeyarajah, EJ, Cromwell, WC, and Otvos, JD. Lipoprotein particle analysis by nuclear magnetic resonance spectroscopy. Clin Lab Med. (2006) 26:847–70. doi: 10.1016/j.cll.2006.07.006
104. Levy, JC, Matthews, DR, and Hermans, MP. Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care. (1998) 21:2191–2. doi: 10.2337/diacare.21.12.2191
105. Katz, A, Nambi, SS, Mather, K, Baron, AD, Follmann, DA, Sullivan, G, et al. Quantitative insulin sensitivity check index: a simple, accurate method for assessing insulin sensitivity in humans. J Clin Endocrinol Metab. (2000) 85:2402–10. doi: 10.1210/jcem.85.7.6661
106. Perseghin, G, Caumo, A, Caloni, M, Testolin, G, and Luzi, L. Incorporation of the fasting plasma FFA concentration into QUICKI improves its association with insulin sensitivity in nonobese individuals. J Clin Endocrinol Metab. (2001) 86:4776–81. doi: 10.1210/jcem.86.10.7902
107. Shalaurova, I, Connelly, MA, Garvey, WT, and Otvos, JD. Lipoprotein insulin resistance index: a lipoprotein particle-derived measure of insulin resistance. Metab Syndr Relat Disord. (2014) 12:422–9. doi: 10.1089/met.2014.0050
108. Limpert, E, and Stahel, WA. Problems with using the normal distribution – and ways to improve quality and efficiency of data analysis. PLoS One. (2011) 6:e21403. doi: 10.1371/journal.pone.0021403
109. Bender, R, and Lange, S. Adjusting for multiple testing—when and how? J Clin Epidemiol. (2001) 54:343–9. doi: 10.1016/s0895-4356(00)00314-0
110. Liang, K-Y, and Zeger, SL. Longitudinal data analysis of continuous and discrete responses for pre-post designs. Sankhyā. (2000) 62:134–48.
111. Liu, GF, Lu, K, Mogg, R, Mallick, M, and Mehrotra, DV. Should baseline be a covariate or dependent variable in analyses of change from baseline in clinical trials? Stat Med. (2009) 28:2509–30. doi: 10.1002/sim.3639
112. Lu, K. On efficiency of constrained longitudinal data analysis versus longitudinal analysis of covariance. Biometrics. (2010) 66:891–6. doi: 10.1111/j.1541-0420.2009.01332.x
113. Coffman, CJ, Edelman, D, and Woolson, RF. To condition or not condition? Analysing ‘change’ in longitudinal randomised controlled trials. BMJ Open. (2016) 6:e013096. doi: 10.1136/bmjopen-2016-013096
114. Twisk, J, Bosman, L, Hoekstra, T, Rijnhart, J, Welten, M, and Heymans, M. Different ways to estimate treatment effects in randomised controlled trials. Contemp Clin Trials Commun. (2018) 10:80–5. doi: 10.1016/j.conctc.2018.03.008
115. Kenward, MG, and Roger, JH. The use of baseline covariates in crossover studies. Biostatistics. (2010) 11:1–17. doi: 10.1093/biostatistics/kxp046
116. Wang, H, Xia, N, Yang, Y, and Peng, D-Q. Influence of vitamin D supplementation on plasma lipid profiles: a meta-analysis of randomized controlled trials. Lipids Health Dis. (2012) 11:42. doi: 10.1186/1476-511X-11-42
117. Hannuksela, ML, Liisanantti, MK, and Savolainen, MJ. Effect of alcohol on lipids and lipoproteins in relation to atherosclerosis. Crit Rev Clin Lab Sci. (2002) 39:225–83. doi: 10.1080/10408360290795529
118. Kraus, WE, Houmard, JA, Duscha, BD, Knetzger, KJ, Wharton, MB, McCartney, JS, et al. Effects of the amount and intensity of exercise on plasma lipoproteins. N Engl J Med. (2002) 347:1483–92. doi: 10.1056/NEJMoa020194
119. Brien, SE, Ronksley, PE, Turner, BJ, Mukamal, KJ, and Ghali, WA. Effect of alcohol consumption on biological markers associated with risk of coronary heart disease: systematic review and meta-analysis of interventional studies. BMJ. (2011) 342:d636. doi: 10.1136/bmj.d636
120. Vu, KN, Ballantyne, CM, Hoogeveen, RC, Nambi, V, Volcik, KA, Boerwinkle, E, et al. Causal role of alcohol consumption in an improved lipid profile: the atherosclerosis risk in communities (ARIC) study. PLoS One. (2016) 11:e0148765. doi: 10.1371/journal.pone.0148765
121. Cole, TJ. Sympercents: symmetric percentage differences on the 100 loge scale simplify the presentation of log transformed data. Stat Med. (2000) 19:3109–25. doi: 10.1002/1097-0258(20001130)19:22<3109::AID-SIM558>3.0.CO;2-F
122. Vickers, AJ, and Altman, DG. Statistics notes: analysing controlled trials with baseline and follow up measurements. BMJ. (2001) 323:1123–4. doi: 10.1136/bmj.323.7321.1123
123. Vickers, AJ. Parametric versus non-parametric statistics in the analysis of randomized trials with non-normally distributed data. BMC Med Res Methodol. (2005) 5:35. doi: 10.1186/1471-2288-5-35
124. Senn, S. Change from baseline and analysis of covariance revisited. Stat Med. (2006) 25:4334–44. doi: 10.1002/sim.2682
125. Beunckens, C, Molenberghs, G, and Kenward, MG. Direct likelihood analysis versus simple forms of imputation for missing data in randomized clinical trials. Clin Trials. (2005) 2:379–86. doi: 10.1191/1740774505cn119oa
126. Peters, SAE, Bots, ML, den Ruijter, HM, Palmer, MK, Grobbee, DE, Crouse, JR, et al. Multiple imputation of missing repeated outcome measurements did not add to linear mixed-effects models. J Clin Epidemiol. (2012) 65:686–95. doi: 10.1016/j.jclinepi.2011.11.012
127. Twisk, J, de Boer, M, de Vente, W, and Heymans, M. Multiple imputation of missing values was not necessary before performing a longitudinal mixed-model analysis. J Clin Epidemiol. (2013) 66:1022–8. doi: 10.1016/j.jclinepi.2013.03.017
128. Buckley, JP, Doherty, BT, Keil, AP, and Engel, SM. Statistical approaches for estimating sex-specific effects in endocrine disruptors research. Environ Health Perspect. (2017) 125:067013. doi: 10.1289/EHP334
129. Cohen, J. Statistical power analysis for the behavioral sciences. Hillsdale, NJ: Lawrence Erlbaum Associates, Inc. (1988).
130. Selya, A, Rose, J, Dierker, L, Hedeker, D, and Mermelstein, R. A practical guide to calculating Cohen’s F2, a measure of local effect size, from proc mixed. Front Psychol. (2012) 3:111. doi: 10.3389/fpsyg.2012.00111
131. Lorah, J. Effect size measures for multilevel models: definition, interpretation, and TIMSS example. Large Scale Assess Educ. (2018) 6:8. doi: 10.1186/s40536-018-0061-2
132. Dwan, K, Li, T, Altman, DG, and Elbourne, D. Consort 2010 statement: extension to randomised crossover trials. BMJ. (2019) 366:l4378. doi: 10.1136/bmj.l4378
133. Harris, WS. The omega-3 index as a risk factor for coronary heart disease. Am J Clin Nutr. (2008) 87:1997S–2002S. doi: 10.1093/ajcn/87.6.1997S
134. Fenton, JI, Gurzell, EA, Davidson, EA, and Harris, WS. Red blood cell Pufas reflect the phospholipid Pufa composition of major organs. Prostaglandins Leukot Essent Fatty Acids. (2016) 112:12–23. doi: 10.1016/j.plefa.2016.06.004
135. Postic, C, Dentin, R, and Girard, J. Role of the liver in the control of carbohydrate and lipid homeostasis. Diabetes Metab. (2004) 30:398–408. doi: 10.1016/s1262-3636(07)70133-7
136. Han, H-S, Kang, G, Kim, JS, Choi, BH, and Koo, S-H. Regulation of glucose metabolism from a liver-centric perspective. Exp Mol Med. (2016) 48:e218. doi: 10.1038/emm.2015.122
137. Kotronen, A, and Yki-Jarvinen, H. Fatty liver: a novel component of the metabolic syndrome. Arterioscler Thromb Vasc Biol. (2008) 28:27–38. doi: 10.1161/ATVBAHA.107.147538
138. Adeli, K, Taghibiglou, C, Van Iderstine, SC, and Lewis, GF. Mechanisms of hepatic very low-density lipoprotein overproduction in insulin resistance. Trends Cardiovasc Med. (2001) 11:170–6. doi: 10.1016/s1050-1738(01)00084-6
139. Avramoglu, RK, Basciano, H, and Adeli, K. Lipid and lipoprotein dysregulation in insulin resistant states. Clin Chim Acta. (2006) 368:1–19. doi: 10.1016/j.cca.2005.12.026
140. Adiels, M, Taskinen, MR, Packard, C, Caslake, MJ, Soro-Paavonen, A, Westerbacka, J, et al. Overproduction of large VLDL particles is driven by increased liver fat content in man. Diabetologia. (2006) 49:755–65. doi: 10.1007/s00125-005-0125-z
141. Adiels, M, Olofsson, S-O, Taskinen, M-R, and Borén, J. Overproduction of very low–density lipoproteins is the hallmark of the dyslipidemia in the metabolic syndrome. Arterioscler Thromb Vasc Biol. (2008) 28:1225–36. doi: 10.1161/atvbaha.107.160192
142. Fabbrini, E, Magkos, F, Mohammed, BS, Pietka, T, Abumrad, NA, Patterson, BW, et al. Intrahepatic fat, not visceral fat, is linked with metabolic complications of obesity. Proc Natl Acad Sci U S A. (2009) 106:15430–5. doi: 10.1073/pnas.0904944106
143. Klop, B, Elte, J, and Cabezas, M. Dyslipidemia in obesity: mechanisms and potential targets. Nutrients. (2013) 5:1218–40. doi: 10.3390/nu5041218
144. Ter Horst, KW, Gilijamse, PW, de Weijer, BA, Kilicarslan, M, Ackermans, MT, Nederveen, AJ, et al. Sexual dimorphism in hepatic, adipose tissue, and peripheral tissue insulin sensitivity in obese humans. Front Endocrinol. (2015) 6:182. doi: 10.3389/fendo.2015.00182
145. Robinson, GA, Peng, J, Peckham, H, Radziszewska, A, Butler, G, Pineda-Torra, I, et al. Sex hormones drive changes in lipoprotein metabolism. iScience. (2021) 24:103257. doi: 10.1016/j.isci.2021.103257
146. Grabacka, M, Pierzchalska, M, Dean, M, and Reiss, K. Regulation of ketone body metabolism and the role of Pparα. Int J Mol Sci. (2016) 17:2093. doi: 10.3390/ijms17122093
147. Puchalska, P, and Crawford, PA. Multi-dimensional roles of ketone bodies in fuel metabolism, signaling, and therapeutics. Cell Metab. (2017) 25:262–84. doi: 10.1016/j.cmet.2016.12.022
148. Lovejoy, JC, and Sainsbury, A, Group tSCW. Sex differences in obesity and the regulation of energy homeostasis. Obes Rev. (2009) 10:154–67. doi: 10.1111/j.1467-789X.2008.00529.x
149. Sugiyama, MG, and Agellon, LB. Sex differences in lipid metabolism and metabolic disease risk. Biochem Cell Biol. (2012) 90:124–41. doi: 10.1139/o11-067 %M 22221155
150. Fuente-Martín, E, Argente-Arizón, P, Ros, P, Argente, J, and Chowen, JA. Sex differences in adipose tissue: it is not only a question of quantity and distribution. Adipocytes. (2013) 2:128–34. doi: 10.4161/adip.24075
151. Fried, SK, Lee, M-J, and Karastergiou, K. Shaping fat distribution: new insights into the molecular determinants of depot- and sex-dependent adipose biology. Obesity. (2015) 23:1345–52. doi: 10.1002/oby.21133
152. Palmer, BF, and Clegg, DJ. The sexual dimorphism of obesity. Mol Cell Endocrinol. (2015) 402:113–9. doi: 10.1016/j.mce.2014.11.029
153. Lin, C, Rajalahti, T, Mjøs, SA, and Kvalheim, OM. Predictive associations between serum fatty acids and lipoproteins in healthy non-obese Norwegians: implications for cardiovascular health. Metabolomics. (2016) 12:6. doi: 10.1007/s11306-015-0886-4
154. Chella Krishnan, K, Mehrabian, M, and Lusis, AJ. Sex differences in metabolism and cardiometabolic disorders. Curr Opin Lipidol. (2018) 29:404–10. doi: 10.1097/mol.0000000000000536
155. Tierney, AC, McMonagle, J, Shaw, DI, Gulseth, HL, Helal, O, Saris, WHM, et al. Effects of dietary fat modification on insulin sensitivity and on other risk factors of the metabolic syndrome – lipgene: a European randomized dietary intervention study. Int J Obes. (2011) 35:800–9. doi: 10.1038/ijo.2010.209
156. Abbott, KA, Veysey, M, Lucock, M, Niblett, S, King, K, Burrows, T, et al. Sex-dependent association between erythrocyte N-3 PUFA and type 2 diabetes in older overweight people. Br J Nutr. (2016) 115:1379–86. doi: 10.1017/S0007114516000258
157. Abbott, KA, Burrows, TL, Thota, RN, Alex, A, Acharya, S, Attia, J, et al. Association between plasma phospholipid omega-3 polyunsaturated fatty acids and type 2 diabetes is sex dependent: the Hunter Community study. Clin Nutr. (2020) 39:1059–66. doi: 10.1016/j.clnu.2019.04.007
158. Després, J-P, Lamarche, B, Mauriège, P, Cantin, B, Dagenais, GR, Moorjani, S, et al. Hyperinsulinemia as an independent risk factor for ischemic heart disease. N Engl J Med. (1996) 334:952–8. doi: 10.1056/NEJM199604113341504
159. Gast, KB, Tjeerdema, N, Stijnen, T, Smit, JWA, and Dekkers, OM. Insulin resistance and risk of incident cardiovascular events in adults without diabetes: meta-analysis. PLoS One. (2012) 7:e52036. doi: 10.1371/journal.pone.0052036
160. Xun, P, Wu, Y, He, Q, and He, K. Fasting insulin concentrations and incidence of hypertension, stroke, and coronary heart disease: a meta-analysis of prospective cohort studies. Am J Clin Nutr. (2013) 98:1543–54. doi: 10.3945/ajcn.113.065565
161. Adeva-Andany, MM, Martínez-Rodríguez, J, González-Lucán, M, Fernández-Fernández, C, and Castro-Quintela, E. Insulin resistance is a cardiovascular risk factor in humans. Diabetes Metab Syndr. (2019) 13:1449–55. doi: 10.1016/j.dsx.2019.02.023
162. Titos, E, and Clària, J. Omega-3-derived mediators counteract obesity-induced adipose tissue inflammation. Prostaglandins Other Lipid Mediat. (2013) 107:77–84. doi: 10.1016/j.prostaglandins.2013.05.003
163. Coelho, OGL, da Silva, BP, Rocha, DMUP, Lopes, LL, and Alfenas, RCG. Polyunsaturated fatty acids and type 2 diabetes: impact on the glycemic control mechanism. Crit Rev Food Sci Nutr. (2017) 57:3614–9. doi: 10.1080/10408398.2015.1130016
164. Navarro, G, Allard, C, Xu, W, and Mauvais-Jarvis, F. The role of androgens in metabolism, obesity, and diabetes in males and females. Obesity. (2015) 23:713–9. doi: 10.1002/oby.21033
165. Ortiz-Huidobro, RI, Velasco, M, Larqué, C, Escalona, R, and Hiriart, M. Molecular insulin actions are sexually dimorphic in lipid metabolism. Front Endocrinol. (2021) 12:690484. doi: 10.3389/fendo.2021.690484
166. Giltay, EJ, Geleijnse, JM, Heijboer, AC, de Goede, J, Oude Griep, LM, Blankenstein, MA, et al. No effects of N−3 fatty acid supplementation on serum total testosterone levels in older men: the alpha omega trial. Int J Androl. (2012) 35:680–7. doi: 10.1111/j.1365-2605.2012.01255.x
167. Abbott, K, Burrows, TL, Acharya, S, Thota, RN, and Garg, ML. Dietary supplementation with docosahexaenoic acid rich fish oil increases circulating levels of testosterone in overweight and obese men. Prostaglandins Leukot Essent Fatty Acids. (2020) 163:102204. doi: 10.1016/j.plefa.2020.102204
168. Friedrichs, W, Ruparel, SB, Marciniak, RA, and deGraffenried, L. Omega-3 fatty acid inhibition of prostate cancer progression to hormone independence is associated with suppression of mTOR signaling and androgen receptor expression. Nutr Cancer. (2011) 63:771–7. doi: 10.1080/01635581.2011.570892
169. Hu, Z, Qi, H, Zhang, R, Zhang, K, Shi, Z, Chang, Y, et al. Docosahexaenoic acid inhibits the growth of hormone-dependent prostate cancer cells by promoting the degradation of the androgen receptor. Mol Med Rep. (2015) 12:3769–74. doi: 10.3892/mmr.2015.3813
170. Chung, BH, Mitchell, SH, Zhang, J-S, and Young, CYF. Effects of docosahexaenoic acid and eicosapentaenoic acid on androgen-mediated cell growth and gene expression in LNCaP prostate Cancer cells. Carcinogenesis. (2001) 22:1201–6. doi: 10.1093/carcin/22.8.1201
171. Xu, W, Morford, J, and Mauvais-Jarvis, F. Emerging role of testosterone in pancreatic Β cell function and insulin secretion. J Endocrinol. (2019) 240:R97–R105. doi: 10.1530/joe-18-0573
172. Morford, J, and Mauvais-Jarvis, F. Sex differences in the effects of androgens acting in the central nervous system on metabolism. Dialogues Clin Neurosci. (2016) 18:415–24. doi: 10.31887/DCNS.2016.18.4/fmauvais
173. Pham, H, and Ziboh, VA. 5α-reductase-catalyzed conversion of testosterone to dihydrotestosterone is increased in prostatic adenocarcinoma cells: suppression by 15-lipoxygenase metabolites of gamma-linolenic and eicosapentaenoic acids. J Steroid Biochem Mol Biol. (2002) 82:393–400. doi: 10.1016/s0960-0760(02)00217-0
174. Gobinet, J, Auzou, G, Nicolas, J-C, Sultan, C, and Jalaguier, S. Characterization of the interaction between androgen receptor and a new transcriptional inhibitor, SHP. Biochemistry. (2001) 40:15369–77. doi: 10.1021/bi011384o
175. Xiao, J, Gong, A-Y, Eischeid, AN, Chen, D, Deng, C, Young, CYF, et al. Mir-141 modulates androgen receptor transcriptional activity in human prostate cancer cells through targeting the small heterodimer partner protein. Prostate. (2012) 72:1514–22. doi: 10.1002/pros.22501
176. Miyata, M, Shinno, K, Kinoshita, T, Kinoshita, Y, and Sugiura, Y. Fish oil feeding reverses hepatomegaly and disrupted hepatic function due to the lack of FXR signaling. J Toxicol Sci. (2017) 42:671–81. doi: 10.2131/jts.42.671
177. Chanda, D, Park, J-H, and Choi, H-S. Molecular basis of endocrine regulation by orphan nuclear receptor small heterodimer partner. Endocr J. (2008) 55:253–68. doi: 10.1507/endocrj.K07E-103
178. Zhang, Y, Hagedorn, CH, and Wang, L. Role of nuclear receptor SHP in metabolism and cancer. Biochim Biophys Acta. (2011) 1812:893–908. doi: 10.1016/j.bbadis.2010.10.006
179. Garruti, G, Wang, HH, Bonfrate, L, de Bari, O, Wang, DQH, and Portincasa, P. A pleiotropic role for the orphan nuclear receptor small heterodimer partner in lipid homeostasis and metabolic pathways. J Lipids. (2012) 2012:304292. doi: 10.1155/2012/304292
180. Palmisano, BT, Zhu, L, Litts, B, Burman, A, Yu, S, Neuman, JC, et al. Hepatocyte small heterodimer partner mediates sex-specific effects on triglyceride metabolism via androgen receptor in male mice. Meta. (2021) 11:330. doi: 10.3390/metabo11050330
181. Morford, JJ, Wu, S, and Mauvais-Jarvis, F. The impact of androgen actions in neurons on metabolic health and disease. Mol Cell Endocrinol. (2018) 465:92–102. doi: 10.1016/j.mce.2017.09.001
182. Ding, EL, Song, Y, Malik, VS, and Liu, S. Sex differences of endogenous sex hormones and risk of type 2 diabetes: a systematic review and meta-analysis. JAMA. (2006) 295:1288–99. doi: 10.1001/jama.295.11.1288
183. Yu, I-C, Lin, H-Y, Sparks, JD, Yeh, S, and Chang, C. Androgen receptor roles in insulin resistance and obesity in males: the linkage of androgen-deprivation therapy to metabolic syndrome. Diabetes. (2014) 63:3180–8. doi: 10.2337/db13-1505
184. Schiffer, L, Kempegowda, P, Arlt, W, and O’Reilly, MW. Mechanisms in endocrinology: the sexually dimorphic role of androgens in human metabolic disease. Eur J Endocrinol. (2017) 177:R125–43. doi: 10.1530/eje-17-0124
185. Barrows, BR, and Parks, EJ. Contributions of different fatty acid sources to very low-density lipoprotein-triacylglycerol in the fasted and fed states. J Clin Endocrinol Metab. (2006) 91:1446–52. doi: 10.1210/jc.2005-1709
186. Vedala, A, Wang, W, Neese, RA, Christiansen, MP, and Hellerstein, MK. Delayed secretory pathway contributions to VLDL-triglycerides from plasma NEFA, diet, and de novo lipogenesis in humans. J Lipid Res. (2006) 47:2562–74. doi: 10.1194/jlr.M600200-JLR200
187. Shearer, GC, Savinova, OV, and Harris, WS. Fish oil: how does it reduce plasma triglycerides? Biochim Biophys Acta. (2012) 1821:843–51. doi: 10.1016/j.bbalip.2011.10.011
188. Gibbons, GF, Islam, K, and Pease, RJ. Mobilisation of triacylglycerol stores. Biochim Biophys Acta. (2000) 1483:37–57. doi: 10.1016/s1388-1981(99)00182-1
189. Gibbons, GF, Brown, AM, Wiggins, D, and Pease, R. The roles of insulin and fatty acids in the regulation of hepatic very-low-density lipoprotein assembly. J R Soc Med. (2002) 95:23–32.
190. Gibbons, GF, Wiggins, D, Brown, A-M, and Hebbachi, A-M. Synthesis and function of hepatic very-low-density lipoprotein. Biochem Soc Trans. (2004) 32:59–64. doi: 10.1042/bst0320059
191. Pan, D, Yang, L, Yang, X, Xu, D, Wang, S, Gao, H, et al. Potential nutritional strategies to prevent and reverse sarcopenia in aging process: role of fish oil-derived Ω-3 polyunsaturated fatty acids, wheat oligopeptide and their combined intervention. J Adv Res. (2023). doi: 10.1016/j.jare.2023.04.005
192. Steiner, G, and Lewis, G. Hyperinsulinemia and triglyceride-rich lipoproteins. Diabetes. (1996) 45:S24–6. doi: 10.2337/diab.45.3.s24
193. Lewis, GF. Fatty acid regulation of very low density lipoprotein production. Curr Opin Lipidol. (1997) 8:146–53. doi: 10.1097/00041433-199706000-00004
194. Julius, U. Influence of plasma free fatty acids on lipoprotein synthesis and diabetic dyslipidemia. Exp Clin Endocrinol Diabetes. (2003) 111:246–50. doi: 10.1055/s-2003-41284
195. Schulze, RJ, and McNiven, MA. Fasting inhibits the recruitment of Kinesin-1 to lipid droplets and stalls hepatic triglyceride secretion. Hepatology. (2019) 69:444–6. doi: 10.1002/hep.30104
196. Reaven, GM. Compensatory hyperinsulinemia and the development of an Atherogenic lipoprotein profile: the Price paid to maintain glucose homeostasis in insulin-resistant individuals. Endocrinol Metab Clin N Am. (2005) 34:49–62. doi: 10.1016/j.ecl.2004.12.001
197. Schreiber, R, Xie, H, and Schweiger, M. Of mice and men: the physiological role of adipose triglyceride lipase (ATGL). Biochim Biophys Acta Mol Cell Biol Lipids. (2019) 1864:880–99. doi: 10.1016/j.bbalip.2018.10.008
198. Carpentier, AC. 100th anniversary of the discovery of insulin perspective: insulin and adipose tissue fatty acid metabolism. Am J Physiol Endocrinol Metab. (2021) 320:E653–70. doi: 10.1152/ajpendo.00620.2020
199. Keller, U, Lustenberger, M, Müller-Brand, J, Gerber, PPG, and Stauffacher, W. Human ketone body production and utilization studied using tracer techniques: regulation by free fatty acids, insulin, catecholamines, and thyroid hormones. Diabetes Metab Rev. (1989) 5:285–98. doi: 10.1002/dmr.5610050306
200. Lewis, GF, Uffelman, KD, Szeto, LW, and Steiner, G. Effects of acute hyperinsulinemia on VLDL triglyceride and VLDL apoB production in normal weight and obese individuals. Diabetes. (1993) 42:833–42. doi: 10.2337/diab.42.6.833
201. Malmström, R, Packard, CJ, Watson, TDG, Rannikko, S, Caslake, M, Bedford, D, et al. Metabolic basis of hypotriglyceridemic effects of insulin in normal men. Arterioscler Thromb Vasc Biol. (1997) 17:1454–64. doi: 10.1161/01.atv.17.7.1454
202. Packard, CJ, Boren, J, and Taskinen, M-R. Causes and consequences of hypertriglyceridemia. Front Endocrinol. (2020) 11:252. doi: 10.3389/fendo.2020.00252
203. Taskinen, M-R, and Borén, J. New insights into the pathophysiology of dyslipidemia in type 2 diabetes. Atherosclerosis. (2015) 239:483–95. doi: 10.1016/j.atherosclerosis.2015.01.039
204. Hirano, T. Pathophysiology of diabetic dyslipidemia. J Atheroscler Thromb. (2018) 25:771–82. doi: 10.5551/jat.RV17023
205. Malmström, R, Packard, CJ, Caslake, M, Bedford, D, Stewart, P, Yki-Järvinen, H, et al. Defective regulation of triglyceride metabolism by insulin in the liver in NIDDM. Diabetologia. (1997) 40:454–62. doi: 10.1007/s001250050700
206. Brown, MS, and Goldstein, JL. Selective versus total insulin resistance: a pathogenic paradox. Cell Metab. (2008) 7:95–6. doi: 10.1016/j.cmet.2007.12.009
207. Vatner, DF, Majumdar, SK, Kumashiro, N, Petersen, MC, Rahimi, Y, Gattu, AK, et al. Insulin-independent regulation of hepatic triglyceride synthesis by fatty acids. Proc Natl Acad Sci U S A. (2015) 112:1143–8. doi: 10.1073/pnas.1423952112
208. Ter Horst, KW, Vatner, DF, Zhang, D, Cline, GW, Ackermans, MT, Nederveen, AJ, et al. Hepatic insulin resistance is not pathway selective in humans with nonalcoholic fatty liver disease. Diabetes Care. (2021) 44:489–98. doi: 10.2337/dc20-1644
209. Sidossis, LS, Mittendorfer, B, Walser, E, Chinkes, D, and Wolfe, RR. Hyperglycemia-induced inhibition of splanchnic fatty acid oxidation increases hepatic triacylglycerol secretion. Am J Phys. (1998) 275:E798–805. doi: 10.1152/ajpendo.1998.275.5.E798
210. Brown, A-M, Wiggins, D, and Gibbons, GF. Glucose phosphorylation is essential for the turnover of neutral lipid and the second stage assembly of triacylglycerol-rich apoB-containing lipoproteins in primary hepatocyte cultures. Arterioscler Thromb Vasc Biol. (1999) 19:321–9. doi: 10.1161/01.atv.19.2.321
211. Adiels, M, Borén, J, Caslake, MJ, Stewart, P, Soro, A, Westerbacka, J, et al. Overproduction of VLDL1 driven by hyperglycemia is a dominant feature of diabetic dyslipidemia. Arterioscler Thromb Vasc Biol. (2005) 25:1697–703. doi: 10.1161/01.ATV.0000172689.53992.25
212. Zammit, VA. Hepatic triacylglycerol synthesis and secretion: DGAT2 as the link between glycaemia and triglyceridaemia. Biochem J. (2013) 451:1–12. doi: 10.1042/bj20121689
213. Rai, P, Kumar, M, Sharma, G, Barak, P, Das, S, Kamat, SS, et al. Kinesin-dependent mechanism for controlling triglyceride secretion from the liver. Proc Natl Acad Sci U S A. (2017) 114:12958–63. doi: 10.1073/pnas.1713292114
214. Wiggins, D, and Gibbons, GF. Origin of hepatic very-low-density lipoprotein triacylglycerol: the contribution of cellular phospholipid. Biochem J. (1996) 320:673–9. doi: 10.1042/bj3200673
215. Asp, L, Claesson, C, Borén, J, and Olofsson, S-O. ADP-ribosylation factor 1 and its activation of phospholipase D are important for the assembly of very low density lipoproteins. J Biol Chem. (2000) 275:26285–92. doi: 10.1074/jbc.M003520200
216. Semenkovich, CF, Wims, M, Noe, L, Etienne, J, and Chan, L. Insulin regulation of lipoprotein lipase activity in 3t3-L1 adipocytes is mediated at posttranscriptional and posttranslational levels. J Biol Chem. (1989) 264:9030–8. doi: 10.1016/S0021-9258(18)81898-1
217. Fried, SK, Russell, CD, Grauso, NL, and Brolin, RE. Lipoprotein lipase regulation by insulin and glucocorticoid in subcutaneous and omental adipose tissues of obese women and men. J Clin Invest. (1993) 92:2191–8. doi: 10.1172/JCI116821
218. Wang, H, and Eckel, RH. Lipoprotein lipase: from gene to obesity. Am J Physiol Endocrinol Metab. (2009) 297:E271–88. doi: 10.1152/ajpendo.90920.2008
219. Li, Y, He, P-P, Zhang, D-W, Zheng, X-L, Cayabyab, FS, Yin, W-D, et al. Lipoprotein lipase: from gene to atherosclerosis. Atherosclerosis. (2014) 237:597–608. doi: 10.1016/j.atherosclerosis.2014.10.016
220. Khan, S, Minihane, A-M, Talmud, PJ, Wright, JW, Murphy, MC, Williams, CM, et al. Dietary long-chain N-3 PUFAs increase LPL gene expression in adipose tissue of subjects with an Atherogenic lipoprotein phenotype. J Lipid Res. (2002) 43:979–85. doi: 10.1016/S0022-2275(20)30473-9
221. Park, Y, and Harris, WS. Omega-3 fatty acid supplementation accelerates chylomicron triglyceride clearance. J Lipid Res. (2003) 44:455–63. doi: 10.1194/jlr.M200282-JLR200
222. Davidson, MH. Mechanisms for the hypotriglyceridemic effect of marine omega-3 fatty acids. Am J Cardiol. (2006) 98:27–33. doi: 10.1016/j.amjcard.2005.12.024
223. Harris, WS, and Bulchandani, D. Why do omega-3 fatty acids lower serum triglycerides? Curr Opin Lipidol. (2006) 17:387–93. doi: 10.1097/01.mol.0000236363.63840.16
224. Bays, HE, Tighe, AP, Sadovsky, R, and Davidson, MH. Prescription omega-3 fatty acids and their lipid effects: physiologic mechanisms of action and clinical implications. Expert Rev Cardiovasc Ther. (2008) 6:391–409. doi: 10.1586/14779072.6.3.391
225. Harris, WS, Miller, M, Tighe, AP, Davidson, MH, and Schaefer, EJ. Omega-3 fatty acids and coronary heart disease risk: clinical and mechanistic perspectives. Atherosclerosis. (2008) 197:12–24. doi: 10.1016/j.atherosclerosis.2007.11.008
226. Jacobson, TA. Role of N-3 fatty acids in the treatment of hypertriglyceridemia and cardiovascular disease. Am J Clin Nutr. (2008) 87:1981S–90S. doi: 10.1093/ajcn/87.6.1981S
227. Oscarsson, J, and Hurt-Camejo, E. Omega-3 fatty acids eicosapentaenoic acid and docosahexaenoic acid and their mechanisms of action on apolipoprotein B-containing lipoproteins in humans: a review. Lipids Health Dis. (2017) 16:149. doi: 10.1186/s12944-017-0541-3
228. Palmisano, BT, Zhu, L, and Stafford, JM. Role of estrogens in the regulation of liver lipid metabolism In: F Mauvais-Jarvis, editor. Sex and gender factors affecting metabolic homeostasis, diabetes and obesity. Cham: Springer International Publishing (2017). 227–56.
229. Johansson, L, Thomsen, JS, Damdimopoulos, AE, Spyrou, G, Gustafsson, J-Å, and Treuter, E. The orphan nuclear receptor SHP inhibits agonist-dependent transcriptional activity of estrogen receptors ERalpha and ERbeta. J Biol Chem. (1999) 274:345–53. doi: 10.1074/jbc.274.1.345
230. Walsh, BW, Schiff, I, Rosner, B, Greenberg, L, Ravnikar, V, and Sacks, FM. Effects of postmenopausal estrogen replacement on the concentrations and metabolism of plasma lipoproteins. N Engl J Med. (1991) 325:1196–204. doi: 10.1056/NEJM199110243251702
231. Mann, CJ, Yen, FT, Grant, AM, and Bihain, BE. Mechanism of plasma cholesteryl Ester transfer in hypertriglyceridemia. J Clin Invest. (1991) 88:2059–66. doi: 10.1172/JCI115535
232. Packard, CJ, and Shepherd, J. Lipoprotein heterogeneity and apolipoprotein B metabolism. Arterioscler Thromb Vasc Biol. (1997) 17:3542–56. doi: 10.1161/01.atv.17.12.3542
233. Packard, CJ. Triacylglycerol-rich lipoproteins and the generation of small, dense low-density lipoprotein. Biochem Soc Trans. (2003) 31:1066–9. doi: 10.1042/bst0311066
234. Rashid, S, Watanabe, T, Sakaue, T, and Lewis, GF. Mechanisms of HDL lowering in insulin resistant, hypertriglyceridemic states: the combined effect of HDL triglyceride enrichment and elevated hepatic lipase activity. Clin Biochem. (2003) 36:421–9. doi: 10.1016/s0009-9120(03)00078-x
235. Murdoch, SJ, and Breckenridge, WC. Influence of lipoprotein lipase and hepatic lipase on the transformation of VLDL and HDL during lipolysis of VLDL. Atherosclerosis. (1995) 118:193–212. doi: 10.1016/0021-9150(95)05606-8
236. Berneis, KK, and Krauss, RM. Metabolic origins and clinical significance of LDL heterogeneity. J Lipid Res. (2002) 43:1363–79. doi: 10.1194/jlr.R200004-JLR200
237. Carr, MC, Hokanson, JE, Zambon, A, Deeb, SS, Barrett, PHR, Purnell, JQ, et al. The contribution of intraabdominal fat to gender differences in hepatic lipase activity and low/high density lipoprotein heterogeneity. J Clin Endocrinol Metab. (2001) 86:2831–7. doi: 10.1210/jcem.86.6.7586
238. Deeb, SS, Zambon, A, Carr, MC, Ayyobi, AF, and Brunzell, JD. Hepatic lipase and dyslipidemia: interactions among genetic variants, obesity, gender, and diet. J Lipid Res. (2003) 44:1279–86. doi: 10.1194/jlr.R200017-JLR200
239. Després, JP, Ferland, M, Moorjani, S, Nadeau, A, Tremblay, A, Lupien, PJ, et al. Role of hepatic-triglyceride lipase activity in the association between intra-abdominal fat and plasma HDL cholesterol in obese women. Arteriosclerosis. (1989) 9:485–92. doi: 10.1161/01.ATV.9.4.485
240. van Deursen, D, Jansen, H, and Verhoeven, AJM. Glucose increases hepatic lipase expression in HepG2 liver cells through upregulation of upstream stimulatory factors 1 and 2. Diabetologia. (2008) 51:2078. doi: 10.1007/s00125-008-1125-6
241. Taskinen, MR. Diabetic dyslipidaemia: from basic research to clinical practice. Diabetologia. (2003) 46:733–49. doi: 10.1007/s00125-003-1111-y
242. Lewis, GF, and Rader, DJ. New insights into the regulation of HDL metabolism and reverse cholesterol transport. Circ Res. (2005) 96:1221–32. doi: 10.1161/01.RES.0000170946.56981.5c
243. Vergès, B. New insight into the pathophysiology of lipid abnormalities in type 2 diabetes. Diabetes Metab. (2005) 31:429–39. doi: 10.1016/s1262-3636(07)70213-6
244. St-Pierre, AC, Cantin, B, Dagenais, GR, Mauriège, P, Bernard, P-M, Després, J-P, et al. Low-density lipoprotein subfractions and the long-term risk of ischemic heart disease in men: 13-year follow-up data from the Québec cardiovascular study. Arterioscler Thromb Vasc Biol. (2005) 25:553–9. doi: 10.1161/01.ATV.0000154144.73236.f4
245. Rizzo, M, Pernice, V, Frasheri, A, Di Lorenzo, G, Rini, GB, Spinas, GA, et al. Small, dense low-density lipoproteins (LDL) are predictors of cardio- and cerebro-vascular events in subjects with the metabolic syndrome. Clin Endocrinol. (2009) 70:870–5. doi: 10.1111/j.1365-2265.2008.03407.x
246. Williams, PT, Zhao, X-Q, Marcovina, SM, Brown, BG, and Krauss, RM. Levels of cholesterol in small LDL particles predict atherosclerosis progression and incident CHD in the HDL-atherosclerosis treatment study (HATS). PLoS One. (2013) 8:e56782. doi: 10.1371/journal.pone.0056782
247. Hoogeveen, RC, Gaubatz, JW, Sun, W, Dodge, RC, Crosby, JR, Jiang, J, et al. Small dense low-density lipoprotein-cholesterol concentrations predict risk for coronary heart disease: the atherosclerosis risk in communities (ARIC) study. Arterioscler Thromb Vasc Biol. (2014) 34:1069–77. doi: 10.1161/atvbaha.114.303284
248. Shiffman, D, Louie, JZ, Caulfield, MP, Nilsson, PM, Devlin, JJ, and Melander, O. Ldl subfractions are associated with incident cardiovascular disease in the Malmö prevention project study. Atherosclerosis. (2017) 263:287–92. doi: 10.1016/j.atherosclerosis.2017.07.003
Keywords: sex-specific responses, omega-3 PUFAs, lipoprotein subfractions, blood lipids, glucose-insulin homeostasis, abdominal obesity, crossover design, constrained linear mixed-effects modeling
Citation: Laupsa-Borge J, Grytten E, Bohov P, Bjørndal B, Strand E, Skorve J, Nordrehaug JE, Berge RK, Rostrup E, Mellgren G, Dankel SN and Nygård OK (2023) Sex-specific responses in glucose-insulin homeostasis and lipoprotein-lipid components after high-dose supplementation with marine n-3 PUFAs in abdominal obesity: a randomized double-blind crossover study. Front. Nutr. 10:1020678. doi: 10.3389/fnut.2023.1020678
Edited by:
Leslie Shen, Texas Tech University Health Sciences Center, United StatesReviewed by:
Salvatore Vaccaro, IRCCS Local Health Authority of Reggio Emilia, ItalyDa Pan, Southeast University, China
Copyright © 2023 Laupsa-Borge, Grytten, Bohov, Bjørndal, Strand, Skorve, Nordrehaug, Berge, Rostrup, Mellgren, Dankel and Nygård. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Johnny Laupsa-Borge, am9obm55LmxhdXBzYS1ib3JnZUB1aWIubm8=
†PRESENT ADDRESSES: Bodil Bjørndal, Department of Sports, Physical activity and Food, Western Norway University of Applied Sciences, Bergen, Norway; Elin Strand, Department of Immunology and Transfusion Medicine, Haukeland University Hospital, Bergen, Norway; Espen Rostrup, VIS Vestlandets Innovasjonsselskap AS, Department of Technology Transfer, Bergen, Norway