Associations of Dietary Fats with All-Cause Mortality and Cardiovascular Disease Mortality among Patients with Cardiometabolic Disease

Previous studies have shown distinct associations between specific dietary fats and mortality. However, evidence on specific dietary fats and mortality among patients with cardiometabolic disease (CMD) remains unclear. The aim of this study was to estimate the association between consumption of specific fatty acids and survival of patients with CMD and examine whether cardiometabolic biomarkers can mediate the above effects. The study included 8537 participants with CMD, from the Third National Health and Nutrition Examination Survey (NHANES III) and NHANES 1999–2014. Cox proportional hazards regression, restricted cubic spline regression, and isocaloric substitution models were used to estimate the associations of dietary fats with all-cause mortality and cardiovascular disease (CVD) mortality among participants with CMD. Mediation analysis was performed to assess the potential mediating roles of cardiometabolic biomarkers. During a median follow-up of 10.3 years (0–27.1 years), 3506 all-cause deaths and 882 CVD deaths occurred. The hazard ratios (HRs) of all-cause mortality among patients with CMD were 0.85 (95% confidence interval (CI), 95% CI, 0.73–0.99; p trend = 0.03) for ω-6 polyunsaturated fatty acids (ω-6 PUFA), 0.86 (95% CI, 0.75–1.00; p trend = 0.05) for linoleic acid (LA), and 0.86 (95% CI, 0.75–0.98; p trend = 0.03) for docosapentaenoic acid (DPA). Isocalorically replacing energy from SFA with PUFA and LA were associated with 8% and 4% lower all-cause mortality respectively. The HRs of CVD mortality among CMD patients comparing extreme tertiles of specific dietary fats were 0.60 (95% CI, 0.48–0.75; p trend = 0.002) for eicosapentaenoic acid (EPA), and 0.64 (95% CI, 0.48–0.85; p trend = 0.002) for DPA and above effects were mediated by levels of total cholesterol (TC), triglycerides (TG), low density lipoprotein cholesterol (LDL), and high density lipoprotein cholesterol (HDL). Restricted cubic splines showed significant negative nonlinear associations between above specific dietary fats and mortality. These results suggest that intakes of ω-6 PUFA, LA, and DPA or replacing SFA with PUFA or LA might be associated with lower all-cause mortality for patients with CMD. Consumption of EPA and DPA could potentially reduce cardiovascular death for patients with CMD, and their effects might be regulated by cardiometabolic biomarkers indirectly. More precise and representative studies are further needed to validate our findings.


Introduction
Cardiometabolic disease (CMD), mainly consisting of cardiovascular disease (CVD) and diabetes mellitus [1][2][3], are a leading contributor to death and disability in the world [4][5][6][7]. There were estimated 451 million adults living with diabetes worldwide in 2017 [8]. The prevalence of CVD increased from 271 million in 1990 to 523 million in 2019 [9]. Furthermore, CVD has been the leading cause of death globally. An estimated 17.9 million people die from CVD each year, accounting for 31% of all deaths worldwide [10]. Among various environmental exposures, diet is one of the most important modifiable risk factors for CMD [11], especially dietary fats [12].

Study Population
Detailed information of NHANES has been described elsewhere [19,20]. Briefly, NHANES is a major program of the National Center for Health Statistics (NCHS) conducted annually for the purpose of assessing the health and nutritional status of people in the United States by interviews and physical examinations. The complex and multi-stage survey examines about 5000 persons each year from various counties across the country, and they were assigned weights for their probability of being randomly selected into the study. To represent the U.S. population, their weighs were taken into account in our analysis. In our current study, we used data available on website from NHANES III and NHANES 1999-2014. Disease information was collected from self-reported questionnaires. CVD was defined if participants reported they had congestive heart failure, coronary heart disease, angina, heart attack, or stroke. Diabetes mellitus was defined if participants had self-reported diabetes or taking medication (including diabetes pills or insulin) to lower blood sugar at baseline. In this analysis, 11,249 patients with CMD were included (4428 with CVD, 4961 with diabetes only, and 1860 with CVD and diabetes simultaneously). Of the 11,124 patients with CMD, participants were excluded because: (a) they had extreme energy intake (total energy intake <800 kcal/day or >4200 kcal/day for men, and <500 kcal/day or >3500 kcal/day for women) (n = 395), (b) they had missing information on fatty acids consumption and covariates (n = 1252), or (c) they had no mortality follow-up (n = 1065). Thus, the remaining analysis consisted of 8537 individuals. The detailed inclusion and exclusion criteria have been presented in Figure 1. up (n = 1065). Thus, the remaining analysis consisted of 8537 individuals. The detailed inclusion and exclusion criteria have been presented in Figure 1.

Dietary Assessment and Other Covariates
Our dietary data were collected from dietary recall from NHANES III and NHANES 1999-2014. In NHANES III and NHANES 1999-2002, one piece of 24 h diet data were collected to estimate intakes of food during the 24-h period prior to the interview (midnight to midnight), and the method has been validated as a useful method for dietary assessment [21]. After 2002, the second 24 h diet data were collected by telephone 3 to 10 days later. An average of two pieces of 24 h dietary intake data during NHANES 2003-2014 was used to represent usual dietary intake status of participants, and it is a more accurate estimation to use a minimum of two nonconsecutive days of dietary intake data

Dietary Assessment and Other Covariates
Our dietary data were collected from dietary recall from NHANES III and NHANES 1999-2014. In NHANES III and NHANES 1999-2002, one piece of 24 h diet data were collected to estimate intakes of food during the 24-h period prior to the interview (midnight to midnight), and the method has been validated as a useful method for dietary assessment [21]. After 2002, the second 24 h diet data were collected by telephone 3 to 10 days later. An average of two pieces of 24 h dietary intake data during NHANES 2003-2014 was used to represent usual dietary intake status of participants, and it is a more accurate estimation to use a minimum of two nonconsecutive days of dietary intake data [22]. Based on obtained dietary consuming information, intakes of total energy and nutrients were calculated based on the United States Department of Agriculture (USDA) Survey Nutrient Database [23,24]. Other covariates were obtained by questionnaires, examinations, and laboratory tests, including gender, age, race (non-Hispanic white, non-Hispanic black, Mexican American, and others), education level (below high school, accomplished high school and higher than high school), family income-poverty ratio (tertiles of PIR, <1.3, 1.3-2.4, ≥2.4), Body Mass Index (BMI, <18.5, 18.5-30, ≥30), use of medication to lower blood sugar, self-reported health status, baseline history of hypertension, hyperlipidemia and cancer, total cholesterol (TC), triglycerides (TG), high density lipoprotein (HDL), and low density lipoprotein (LDL). Never smoking was defined in the questionnaire as smoking fewer than 100 cigarettes in life. Participants who had smoked over 100 cigarettes during their lifetime and reported smoking at the time of the interview were considered as current smokers. Former smokers were defined as those who smoked over 100 cigarettes and had quit smoking. Alcohol drinking status was categorized as none drinking (0 g/day), moderate drinking (0.1 to 27.9 g/day for men and 0.1 to 13.9 g/day for women), and heavy drinking (28 g/day or more for men and 14 g/day or more for women). Physical activity was classified as inactive, insufficiently active, and active according to tertiles of metabolic equivalent of task (MET).

Ascertainment of Death
Mortality status was determined by probabilistic matching to the National Death Index (NDI) through 31 December 2015 based on a unique study identifier [25]. Briefly, the linkage between the survey data and the NDI combined both deterministic and probabilistic approaches. The deterministic approach is based on primary identifiers of Social Security number 9 (SSN9) or Social Security number 4 (SSN4) (depending on the survey year or cycle of the survey), first name, middle initial, last name, or father's surname, race, sex, and so on. The probabilistic approach performs weighting and link adjudication following the Fellegi-Sunter method, which is the foundational methodology used for record linkage, showing a high accuracy by a prior validation study [26]. The underlying causes-of-death were classified based on the international statistical classification of diseases, 10th revision (ICD-10). In addition, data on mortality caused by stroke (I60-I69) were available only until 31 December 2011. Therefore, the main outcomes of our study were all-cause mortality and CVD (not including stroke) mortality.

Statistical Analyses
All analyses accounted for the complex survey design, including sample weights, clustering and stratification [19]. In our analyses, the amount of macronutrients intake was represented as percent of total energy (%E) by the nutrient density method, which was computed by dividing energy from the nutrient by the total daily energy intake (%E = ((fats (g) × 4)/total energy intake (kcal)) × 100) [27]. We categorized participants into tertiles of intake levels. Cox proportional hazards regression was applied to investigate the associations of dietary fats with all-cause and CVD mortality among patients with CMD. Person-years of follow-up were calculated from baseline to the earliest of time of death, loss to or unavailability for follow-up, or the end of follow-up, whichever came first. Hazard ratios (HRs) and the corresponding 95% confidence intervals (CIs) were calculated by using the lowest tertile group as reference. Tests for trends were performed by assigning median values to corresponding categories of intakes as continuous variables, which were adjusted for all potential confounding variables. Additionally, restricted cubic splines with three knots were used to detect dose-response relationships between dietary fats and mortality. Furthermore, we evaluated the effect of replacing SFA with other types of fats by establishing isocaloric substitution models that simultaneously included total energy intake and percentages of energy from proteins, carbohydrates, and all fats except SFAs. By leaving SFAs out of the model, regression coefficients of the other fats could be interpreted as estimated effects of replacing SFAs with specific fats while leaving other fatty acids unchanged. Mediation analysis was performed to assess the mediating roles of cardiometabolic biomarkers (including total cholesterol (TC), total triglyceride (TG), low density lipoprotein cholesterol (LDL), and high density lipoprotein cholesterol (HDL)) and Nutrients 2022, 14, 3608 5 of 16 beneficial fatty acids in CVD mortality among participants with CMD. Furthermore, we also conducted additional analysis to explore the association of specific dietary fats and all-cause mortality and CVD mortality in the general population.
Sensitivity analyses were conducted to test the robustness of our results, including evaluating the associations of dietary fats with all-cause and CVD mortality among CMD patients after excluding patients who died within the first two years of follow-up and further adjustment of the use of medication (including diabetes pills and insulin) from NHANES III and NHANES 1999-2014. All statistical analyses were conducted using SAS software version 9.4 (SAS Institute, Cary, NC, USA) and R software. Mediation analysis was generated by SAS Mediation Macro. A two-sided p-value < 0.05 was considered statistically significant.

Population Characteristics
During a median follow-up length of 10.3 years (0-27.1 years), 3506 (41.07%) deaths occurred, including 882 deaths from CVD. Table 1 shows baseline characteristics of patients with CMD in NHANES III and NHANES 1999-2014. Among 8537 patients with CMD, participants with higher SFA, MUFA, and PUFA intakes were more likely to be non-Hispanic white, richer, and fatter. However, they were less likely to drink. The participants with higher consumption of SFA and MUFA were more likely to be male, younger, current smokers, and with good self-reported health status. For participants consuming more PUFA, they tended to be higher educated, and with more prevalent hypertension and hyperlipidemia. Intakes of SFA, MUFA, and PUFA were positively associated with total energy intake but inversely related to carbohydrate intake. In addition, those with higher PUFA intake seemed to have more protein consumption (Table 1). For specific PUFA, the baseline characteristics of patients with CMD according to the tertiles of energy percentage from ω-6 PUFA and ω-3 PUFA are presented in the Table S1. The participants with more ω-6 PUFA consumption had the similar baseline characteristics with those having more consumption of PUFA. Those with higher ω-3 PUFA consumption were more likely to be higher educated, richer, and with good self-reported health status. They also tended to consume more total energy but less carbohydrates (Table S1).  Tertile of PUFA Intake   T1  T2  T3  T1  T2  T3  T1  T2  T3 BMI group (%) <18. 5

All-Cause Mortality
Among patients with CMD, total SFA, MUFA, and PUFA were not associated with allcause mortality, but ω-6 PUFA, and the majority ω-6 PUFA, LA were inversely associated with all-cause mortality in multivariable-adjusted models, with the HRs of 0.85 (95% CI: 0.73-0.99; p trend = 0.03) and 0.86 (95% CI: 0.75-1.00; p trend = 0.03) for highest vs. lowest tertile, respectively. Similarly, DPA were inversely associated with all-cause mortality with the HRs of 0.86 (95% CI: 0.75-0.98; p trend = 0.03) for highest vs. lowest tertile, respectively ( Table 2). No association of residual fatty acids components with all-cause mortality was observed ( Table 2). In the substitution model, a 5% replacement of energy from SFA with PUFA, and a 2% replacement of energy from SFA with LA were associated with an 8% lower all-cause mortality and a 4% lower all-cause mortality, respectively ( Figure 2). Figure 3 showed a nonlinear association between ω-6 PUFA, LA, DPA, and all-cause mortality (all p nonlinear < 0.05). The HRs for all-cause mortality decreased as the intake of above dietary fats increased and reached a plateau when intakes were approximately 7.7%E, 7.5%E, and 0.02%E for ω-6 PUFA, LA, and DPA, respectively. This suggests that once dietary fats reach a certain level, the all-cause mortality hardly decreases with further increases in ω-6 PUFA, LA, and DPA. For the general population, the intakes of ω-6 PUFA, ω-3 PUFA, marine ω-3 PUFA, and its components (DHA, EPA, DPA) were inversely associated with all-cause mortality. (Table S2).    The multivariable-adjusted model was adjusted for gender, age, race, educational levels, PIR, BMI, smoking status, alcohol drinking status, physical activity status, self-reported health status, baseline history of hypertension, hyperlipidemia and cancer, total energy intake, energy intake derived from carbohydrate, protein intake (in tertiles), cholesterol levels (in tertiles), and energy intake derived from remaining fatty acids (SFA, MUFA, and PUFA) where appropriate. Abbreviation and acronyms: CMD, cardiometabolic disease; CVD, cardiovascular disease; ω-6 PUFA, ω-6 polyunsaturated fatty acids; LA, linoleic acid; EPA, eicosapentaenoic acid; DPA, docosapentaenoic acid.

CVD Mortality
Among patients with CMD, total SFA, MUFA, and PUFA were not associated with CVD mortality. For specific PUFAs, marine ω-3 PUFA components-EPA and DPA were strongly associated with reduced CVD mortality, with HRs for highest vs. lowest tertile of 0.60 (95% CI, 0.48-0.75; p trend = 0.002) and 0.64 (95% CI, 0.48-0.85; p trend = 0.002), respectively (Table 3). No association of residual PUFA components with CVD mortality was observed ( Table 3). The nonlinear dose-response relationships of EPA and DPA with CVD mortality among patients with CMD are shown in Figure 3. A U-shaped curve was observed between EPA and CVD mortality. The CVD mortality decreased dramatically along with the increased intakes of EPA at a low level and then kept rising at a relatively higher level. The protective effect for CVD mortality approached to saturation when EPA accounted for 0.08%E. For DPA in CVD mortality, we can see the same trend as in allcause mortality. Mediating effect analysis indicated that the health effect of EPA and DPA in CVD mortality could be regulated by cardiometabolic biomarkers through indirect effects with different mediating proportions. When higher intake of EPA was compared with lower intake of EPA, the proportion mediated by cardiometabolic biomarkers was 5.33% for TC, and 6.18% for TG, 1.84% for HDL, and 3.94% for LDL (Table 4). Similarly, when higher intake of DPA was compared with lower intake of DPA, the proportion mediated by cardiometabolic biomarkers was 9.55% for TC, and 6.53% for TG, 2.74% for HDL, and 11.12% for LDL (Table 4). For the general population, the intakes of ω-3 PUFA, were estimated by multivariable adjusted Cox regression models based on restricted cubic splines. Solid splines represented estimated hazard ratios, and cloud areas represented the corresponding 95% CI. The multivariable-adjusted model was adjusted for gender, age, race, educational levels, PIR, BMI, smoking status, alcohol drinking status, physical activity status, self-reported health status, baseline history of hypertension, hyperlipidemia and cancer, total energy intake, energy intake derived from carbohydrate, protein intake (in tertiles), cholesterol levels (in tertiles), and energy intake derived from remaining fatty acids (SFA, MUFA, and PUFA) where appropriate. Abbreviation and acronyms: CMD, cardiometabolic disease; CVD, cardiovascular disease; ω-6 PUFA, ω-6 polyunsaturated fatty acids; LA, linoleic acid; EPA, eicosapentaenoic acid; DPA, docosapentaenoic acid.

CVD Mortality
Among patients with CMD, total SFA, MUFA, and PUFA were not associated with CVD mortality. For specific PUFAs, marine ω-3 PUFA components-EPA and DPA were strongly associated with reduced CVD mortality, with HRs for highest vs. lowest tertile of 0.60 (95% CI, 0.48-0.75; p trend = 0.002) and 0.64 (95% CI, 0.48-0.85; p trend = 0.002), respectively (Table 3). No association of residual PUFA components with CVD mortality was observed ( Table 3). The nonlinear dose-response relationships of EPA and DPA with CVD mortality among patients with CMD are shown in Figure 3. A U-shaped curve was observed between EPA and CVD mortality. The CVD mortality decreased dramatically along with the increased intakes of EPA at a low level and then kept rising at a relatively higher level. The protective effect for CVD mortality approached to saturation when EPA accounted for 0.08%E. For DPA in CVD mortality, we can see the same trend as in all-cause mortality. Mediating effect analysis indicated that the health effect of EPA and DPA in CVD mortality could be regulated by cardiometabolic biomarkers through indirect effects with different mediating proportions. When higher intake of EPA was compared with lower intake of EPA, the proportion mediated by cardiometabolic biomarkers was 5.33% for TC, and 6.18% for TG, 1.84% for HDL, and 3.94% for LDL (Table 4). Similarly, when higher intake of DPA was compared with lower intake of DPA, the proportion mediated by cardiometabolic biomarkers was 9.55% for TC, and 6.53% for TG, 2.74% for HDL, and 11.12% for LDL (Table 4). For the general population, the intakes of ω-3 PUFA, marine ω-3 PUFA, and its components (DHA, EPA, DPA) were inversely associated with CVD mortality (Table S2).   a Values were presented as HR (95% CI) of total effect, direct effect, indirect effect, and the mediating proportion of cardiometabolic biomarkers between beneficial fatty acids and CVD mortality using a mediating analysis which was adjusted for the same potential confounders as the Cox proportional hazards regression.

Sensitivity Analyses
We observed similar relationships of dietary fatty acids with CVD mortality among two sensitivity analyses (Tables S3 and S4). For all-cause mortality, we can only see the protective effect of DPA, not including ω-6 PUFA and LA after excluding patients who died within the first two years of follow-up (Table S3). After further adjustment of use of medication (including diabetes pills or insulin), the protective effect of DPA disappeared but similar associations of ω-6 PUFA and LA with all-cause mortality still maintained (Table S4).

Discussion
In this study, we investigated the associations between intakes of dietary fatty acids with all-cause mortality and CVD mortality among representative US patients with CMD from NHANES III and NHANES 1994-2014. The ω-6 PUFA, ω-6 PUFA-LA, and marine ω-3 PUFA-DPA were associated with reduced all-cause mortality risk. The components of marine ω-3 PUFA (EPA and DPA) were strongly associated with a reduced risk of CVD mortality among patients with CMD. The findings underscore the importance of maintaining PUFA consumption, especially marine ω-3 PUFA, to reduce the risk of CVD mortality in CMD patients.
To our knowledge, this is the first prospective analysis that investigated the associations between specific dietary fats and all-cause mortality and CVD mortality among people with CMD. In our study, ω-6 PUFA (accounting for 88.72% of total PUFA) and LA (accounting for 98.86% of ω-6 PUFA) were related with reduced all-cause mortality among CMD patients. Substitution analyses also revealed the protective effect of PUFA and LA in all-cause mortality. The dose-response relationship clearly sees the declining trend with increasing intake of ω-6 PUFA and LA. Although similar research about effects of fatty acids and mortality is scarce in the CMD population, studies on CVD or diabetes mellitus could provide related information. A review from randomized controlled trials summarized that the replacement of SFA with ω-6 PUFA had potential benefits on CVD incidence and mortality [28]. A previous study by Nagai et al. indicated that lower circulating ω-6 PUFA levels were significantly associated with all-cause death in patients with acute decompensated heart failure [29]. A meta-analysis by Wu et al. concluded that LA had long-term benefits for the prevention of type 2 diabetes (T2D) [30]. Furthermore, two prospective, longitudinal cohort studies also found inverse associations of the intake of LA with all-cause mortality among patients with T2D [17]. However, a randomized controlled trial concluded that LA led to increased mortality among patients with coronary heart disease [31]. The observation might be confounded by high trans-fat contents of margarines, which were used as the food sources of LA in this intervention [32]. It reminds us that more accurate and effective intervention studies are needed to validate our findings. It is worth noting that, when we excluded patients who died within the first two years of follow-up to reduce selecting bias, we found that the protective effects of ω-6 PUFA and LA disappeared. It warns us to draw conclusions about the above associations cautiously.
In the current study, we found that dietary DPA could lead to decreased all-cause mortality among CMD patients. However, after further adjustment of the use of medication (including diabetes pills or insulin), the protective effect of DPA disappeared. We also observed reduced CVD mortality among CMD patients with higher intake of EPA and DPA (components of marine ω-3 PUFA). Numerous studies have recognized the potential cardiovascular benefits of marine ω-3 PUFA, including DHA, EPA, and DPA [33,34]. Our findings are consistent with a meta-analysis based on 13 random clinical trials which concluded that marine ω-3 fatty acids were associated with reduced CVD death, with the HR of 0.92 (95% CI: 0.88-0.97) compared with placebo [35]. Moreover, a multicenter prospective cohort also revealed that the circulating levels of ω-3 PUFA (including DHA, EPA, and DPA) were inversely associated with all-cause mortality, especially the CHD death in older adults [36]. Two cohort studies in the United States concluded that marine ω-3 PUFA could reduce CVD mortality among patients with type 2 diabetes [17]. A randomized controlled trial also found in patients with a history of coronary artery disease who were given EPA treatment, and major coronary events were reduced by 19% in the EPA supplement group vs. in the control group [37]. A total of 67 prospective studies that is comprised of 310,955 participants identified that EPA and DPA could lead to a lower risk of T2D [38]. For EPA, the mechanisms of the cardiac-protect effects may include decreasing plasma triglyceride, resting heart rate and blood pressure, inhibiting platelet activity, and stabilizing/regressing of coronary plaque [35,[39][40][41][42][43]. Compared with EPA, research concerning the protective effects of DPA in CVD development is limited. A study combining 19 cohort studies reported that DPA could reduce the risk of coronary heart disease in the general population [44]. The potential mechanisms involved include anti-inflammation, inhibition of cytokine synthesis, anti-thrombosis, and improvement of plasma lipids [36,[45][46][47][48]. However, we did not observe a significant association between ALA and CVD mortality among patients with CMD, which accounted for 92.48% of total ω-3 PUFA. Epidemiological studies concerning cardiovascular benefits of ALA remain mixed [49][50][51], indicating that more precise and representative studies are needed. In mediating analysis, the protective effect of EPA and DPA in CVD mortality can work through TC, TG, HDL, and LDL indirectly. Many studies have identified that ω-3 PUFA could influence the concentrations of TC, TG, LDL, and HDL [52][53][54]. However, how the potential influence of these bioactive lipids is regulated by marine ω-3 PUFA is still unclear.
Our study has several strengths. First, the prospective study design, well-controlled and rigorous protocol, large sample size, long follow-up period, and relatively accurate estimates ensured the validity of the current study. Second, besides categorizing intakes of dietary fatty acids into tertiles, we also regarded the dietary fatty acids as continuous variables to estimate the dose-response relationships of fatty acids with mortality. Third, besides ω-6 PUFA and ω-3 PUFA, we further estimated the protective effects of their specific components. However, there are still some limitations. First, dietary data were obtained by recalling, so inaccuracies may exist. However, well trained workers, an abundant dietary database from USDA, and the reliable method and system ensured that the measures of diets were as accurate as possible. Second, dietary behaviors may be altered during follow up. Although day-to-day variations could not be estimated, we have used the average intake level to represent their usual intake status to minimize the measurement errors. Third, even though we adjusted many confounders in regression models, residual confounding effects from unmeasured or unknown factors could not be excluded entirely such as the quality of diabetes management or insulin resistance. More precise and representative studies are further needed. Fourth, mortality status was based on the deterministic and probabilistic matching to the NDI, which might result in misclassification. However, the linkage algorithm can produce high-quality matches with a low degree of linkage error [26]. Lastly, the study was an observational study, and we could not verify causality.

Conclusions
These results suggest that intakes of ω-6 PUFA, LA, and DPA or replacing SFA with PUFA or LA might be associated with lower all-cause mortality for patients with CMD. Consumption of EPA and DPA could potentially reduce cardiovascular death for patients with CMD, and their effects might be regulated by cardiometabolic biomarkers indirectly. More precise and representative studies are further needed to validate our findings.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/nu14173608/s1, Table S1: Baseline characteristics of patients with cardiometabolic diseases in NHANES III and NHANES 1999-2014 according to the tertiles of energy percentage from ω-6 PUFA and ω-3 PUFA; Table S2. Associations between specific fatty acids and all-cause mortality and CVD mortality in the general population (n = 43,450); Table S3: Associations between specific fatty acids and all-cause mortality and CVD mortality after excluding patients who died within the first two years of follow-up among patients with CMD (n = 7492); Table S4. Associations between specific fatty acids and all-cause mortality and CVD mortality after further adjustment of use of medication (including diabetes pills or insulin) (n = 8537).