Research Article  Open Access
Structural Equation Modeling for Analyzing Erythrocyte Fatty Acids in Framingham
Abstract
Research has shown that several types of erythrocyte fatty acids (i.e., omega3, omega6, and trans) are associated with risk for cardiovascular diseases. However, there are complex metabolic and dietary relations among fatty acids, which induce correlations that are typically ignored when using them as risk predictors. A latent variable approach could summarize these complex relations into a few latent variable scores for use in statistical models. Twentytwo red blood cell (RBC) fatty acids were measured in Framingham (N = 3196). The correlation matrix of the fatty acids was modeled using structural equation modeling; the model was tested for goodnessoffit and gender invariance. Thirteen fatty acids were summarized by three latent variables, and gender invariance was rejected so separate models were developed for men and women. A score was developed for the polyunsaturated fatty acid (PUFA) latent variable, which explained about 30% of the variance in the data. The PUFA score included loadings in opposing directions among three omega3 and three omega6 fatty acids, and incorporated the biosynthetic and dietary relations among them. Whether the PUFA factor score can improve the performance of risk prediction in cardiovascular diseases remains to be tested.
1. Introduction
Higher blood levels of the essential omega3 polyunsaturated fatty acids (PUFA) are associated with reduced risk for sudden cardiac death [1, 2] and allcause mortality [2, 3]. There is also evidence that the essential omega6 PUFA intakes and blood levels are inversely associated with risk for coronary heart disease [4]. Other fatty acids, such as the trans fatty acids found in partially hydrogenated vegetable oils, are believed to increase risk for cardiovascular disease [5]. Hence, the study of these fatty acids is of vital importance.
PUFA are “essential” since they cannot be produced in vivo and must be consumed. Foods are composed of multiple fatty acids, and dietary habits manifest themselves as correlated fatty acid levels in the blood. Once consumed, simpler PUFA species can be acted upon by enzymes that convert them into more complex PUFA which have a wide variety of metabolic functions. Desaturase enzymes insert double bonds (points of “desaturation”) into fatty acid molecules, and elongase enzymes are needed to increase the carbon chain length [6]. Importantly, competition among PUFA species exists for these enzymes such that different ratios in the diet can affect overall PUFA patterns (Figure 1) [7]. While elongase enzymes are readily available, the desaturase enzymes are rate limiting, and thus their levels may impact the amount of the 20 and 22carbon fatty acids present in the system [8]. Additionally, omega3 fatty acids are the preferential substrates over omega6 for these desaturase enzymes [6]. These biochemical and dietary relations induce a correlation structure in the blood fatty acids.
Fatty acids have been reported as weight%, mol%, or concentration (by volume or cell count). Since there are no laboratory standards in the USA to uniformly report fatty acid data, these multiple presentations exist. The main debate is relative versus absolute amounts, whose quantities become more divergent as they increase [9]. Chow argues in favor of absolute concentrations, citing the obvious drawback of relative amounts being the imposed linear constraint (i.e, summation to 100%) [10]. However, Crowe prefers relative amounts since absolute concentrations rise and fall with total cholesterol, which is made up of lipoproteins composed of fatty acids [11]. Fatty acid nomenclature is as follows: C#:#n# = the number of carbon (C) atoms in the molecule, the number of double bonds, and the omega family (n), whether 3, 6, 7, or 9. The latter indicates on which carbon the final double bond resides. In Bradbury et al. concentrations and mol% are compared for C14:0 and C18:2n6 in plasma cholesterol ester and phospholipids [12]. The study results show that C18:2n6 is significantly directly correlated with total cholesterol when represented as a concentration () but significantly inversely related as mol%. The paper lists several references that support the cholesterol lowering effect of C18:2n6 and concludes that the metabolic pathways are influenced by the percentage of total fatty acids and not by concentration. The other advantage of weight% representation is that RBC fatty acids have a strong correlation with myocardium tissue fatty acids () [13] and dietary intake [14]. This preferred technique of using relative weight% of total fatty acids also induces a correlation structure in the data.
Structural equation modeling (SEM) is well suited to incorporate the metabolic, dietary, and measurement correlations observed in fatty acid data. Only in the last decade has SEM been applied to fatty acids [15–18]. SEM allows complex highdimensional relations to be simplified into a few latent variable scores, which can be evaluated as novel risk markers. Sexspecific risk prediction models have been implemented for coronary heart disease [19] to account for gender differences in the amount of risk attributable to cholesterol and blood pressure levels. Similarly the observed fatty acids may relate to the underlying latent constructs differentially for men and women, and this potential gender invariance needs to be evaluated. A recent taxonomy has been developed to specifically test measurement invariance over multiple groups in SEM by comparing models with different constraints applied to the correlation structure [20, 21].
The objective of the present study was to reduce the dimensionality of the complex fatty acid correlation structure by incorporating dietary intake patterns and biosynthesis processes as constraints in a structural equation model. This technique was applied to the Framingham Offspring/Omni RBC samples, and differences in fatty acid means, loadings, residuals, and latent variable covariance structures between men and women were tested.
2. Materials and Methods
2.1. Materials
The Framingham Heart Study (FHS) was established in 1948 to research the factors that contribute to cardiovascular disease. Its study design and methods are described at http://www.nhlbi.nih.gov/about/framingham. In 1971, the children (and their spouses) of the original FHS were recruited; they constitute the Framingham Offspring cohort [22]. In 1994, to better reflect the changing demographics of the area, recruitment began for Framingham residents aged 40–74 who described themselves as members of a minority group, that is, Omni cohort [23]. The Offspring and Omni cohorts were scheduled together for comprehensive examinations every 4–8 years. These included anthropometric measurements, biochemical assessment for CVD risk factors, medical history, and physical examination by a study physician. RBC samples taken from Offspring Examination 8 and Omni Exam 3 (2005–2007) were collected and subsequently 22 fatty acids were analyzed using gas chromatography (GC), and their content was expressed as a weight% of total fatty acids [24]. These data are publicly available as part of the National Heart Lung and Blood Institute (NHLBI) SNP Health Association Resource (SHARe) project (release date: March 26, 2013, dataset name: l_rbcfa_2008_m_0420s, dataset accession: pht002568). Written informed consent was provided by all participants, and the Institutional Review Board at the Boston University Medical Center approved the study protocol.
The study participants had a mean age (SD) of 66 (9) years, 55% were female, and 91% were white (see Table 1in Supplementary Material available online at http://dx.doi.org/10.1155/2014/160520). The prevalence of chronic disease was diabetes (14%), heart disease (10%), and congestive heart failure (2%). The participants were taking hypertension medications (49%), lipid pharmacotherapy (43%), aspirin 3+ per week (43%), and fish oil supplements (10%).
2.2. Methods
2.2.1. Model Notation
SEM is a twopart modeling process. The first part defines a measurement model, which specifies the relations between the fatty acids and the latent variables, given by () [25] where is a vector of observed fatty acids measured on subject , is a vector of fatty acid means, is a matrix of unknown loading parameters, is a vector of latent variable scores for subject , and is vector of normal random errors with covariance matrix that is independent of . The second part defines a path model for the latent variables , which allows regressing one latent variable on the set of other latent variables , given by , where is a vector of latent variable means, is a vector of unknown regression parameters, and is vector of normal random errors with covariance matrix .
2.2.2. Comparing SEM to Other Multivariate Techniques
The proportion of variance explained in the th fatty acid by the latent variables is defined as the th communality. The variance in each fatty acid is equal to its communality plus its unique variance; that is, . Principal components analysis (PCA) ignores the specific variance (measurement error) and uses the identity matrix for . This results in factoring the total variance instead of the common variance; the latter is the proportion of variance shared by the fatty acids. If communalities < 1 are used with principal components method of decomposing the observed correlation matrix using eigenvalues and eigenvectors, then the method is principal factoring. This is akin to using the reduced sample correlation matrix where the main diagonals are less than one. Likewise, an exploratory factor analysis imposes no structure and also assumes independence in the residuals matrix, but the common variance is extracted and the factors can be correlated through oblique rotations. Moving to confirmatory factor analysis requires restrictions on the model parameter, which allows testing the model goodnessoffit. However, only structural equation modeling (SEM) allows regression paths among the observed and unobserved variables, also a sparse correlation matrix for the residuals can be specified. Therefore, SEM allows the most model flexibility for implementing the dietary intake patterns and metabolic processes among the fatty acids.
2.2.3. Data Preparation
SEM requires multivariate normality for maximum likelihood (ML) estimation. We started by assessing univariate normality which is implied by multivariate normality. To examine univariate normality, skewness and kurtosis measures were calculated for each fatty acid. If a fatty acid had an absolute kurtosis >10, considered problematic [26], or skewness >3, then a natural logarithm transformation was employed. Next, the null hypothesis of multivariate normality was tested using the SAS macro %MULTNORM which calculates the squared Mahalanobis distances (); for large samples, is distributed as a [27]. When multivariate normality is not tenable, robust ML estimation should be implemented [28].
The scales of the RBC fatty acids differed by two orders of magnitude, on average Palmitic acid (C16:0) accounted for 20% of total fatty acid abundance, whereas alphalinolenic acid (C18:3n3) accounted for only 0.2%. However, there are meaningful fatty acids even at small relative weight%. For example, the mean levels of C20:5n3 and C16:1t are 0.7% and 0.2%, but these are widely studied biomarkers of fish and dairy intake, respectively. Therefore, all fatty acids were standardized in order to have similar effect sizes, which prevented the large abundance fatty acids from dominating the variance extraction. The measure of sampling adequacy developed by Kaiser [29] MSA was used to identify fatty acids that were not sufficiently related to the core latent structure. Individual fatty acids with a MSA were considered “unacceptable” [30] and dropped from analysis.
2.2.4. Gender Invariance Testing
The primary motivation of the study was to reduce the dimensions of the fatty acid correlation matrix and to develop latent variable scores for each subject using the regression scoring method [31]. An assumption when using latent variable scores is that the indicators (i.e, fatty acids) have the same relations with the underlying latent variables between groups of interest, in this case gender. To this end, the correlation matrix was partitioned as , and it along with the mean profile of the observed fatty acids was tested for gender invariance using the likelihood ratio test (LRT) by comparing models with different imposed constraints. When using robust maximum likelihood the LRT has been modified by the deviance scaling [32]. The specific hypotheses are as follows.(1)The fatty acid means are equal between genders, .(2)The loadings matrix is equal between genders, .(3)The fatty acid variances are equal between genders, .(4)The latent variable covariances are equal between genders, .(5)The fatty acid covariances are equal between genders, .
The number of underlying dimensions was examined using exploratory factory analysis where eigenvalues >1 were retained. Then a SEM was built using the same number of latent variables, and it was used to test if the fatty acid means were equal between men and women (H_{1}). This was done by testing the model between a model with intercepts, loadings, and unique variances freely estimated for men and women (model M0) versus one with a single set of intercepts imposed for both genders (model M1). To test for equality in the loading matrix between genders (H_{2}), a model with equal loading constraints, but freely estimated intercepts and unique variances for each gender (model M2), was compared to model M0. Hypothesis H_{3} was tested for gender differences in fatty acid residual variances by comparing a model with freely estimated intercepts and loadings for each gender, but with constrained variances (Model M3) versus model M0. To test the latent variable covariance structure , model M2 was used for comparison with additional constraints placed on the six latent covariances to be equal between genders (Model M4). Lastly the fatty acid covariance matrix was tested by constraining loadings, latent covariances, and unique variances to be equal for men and women (Model M5). Each model used the direct Quartimin oblique rotation (available in Mplus and SAS).
2.2.5. Comparing Model Fits
To evaluate the fit of the SEM there are several indices, and the following is the minimal set established by current practice: (1) model chisquare, (2) SteigerLind root mean square error of approximation (RMSEA), (3) Bentler comparative fix index (CFI), and (4) standardized root mean square residual (SRMSR) [26]. The RMSEA indicates the discrepancy in model fit per degree of freedom as defined by [33]. The RMSEA follows a noncentral distribution, which allows reversing the role of the null hypothesis to testing a poorly fitting model and then a larger sample provides evidence of good fit. RMSEA is not used to test for perfect fit but to test the alternative hypothesis of “close fit” or “reasonable fit” [26]. Bentler’s CFI and the TuckerLewis Index (TLI) are relative fit indexes; these measures should be >0.90 [34], and CFI differences of 0.01 between models are considered relevant [35]. The absolute model fit was assessed by calculating the SRMSR between the fatty acids’ observed correlations and the correlations predicted by the latent variables; these residuals should be less than 0.10 for a good fitting model [26]. The Schwarz Bayesian Criterion [36], which includes a larger penalty for lack of parsimony than Akaike Information Criteria [37], was also reported. Analyses were performed using SAS software (version 9.2; SAS Institute Inc., Cary, NC) and Mplus (version 6.12; Muthen & Muthen, Los Angeles, CA).
3. Results
3.1. Exploratory Factor Analysis
Table 1 indicates gender differences in mean levels, in 12 out of 22 RBC fatty acids. The greatest relative differences were higher levels of C16:1 and C18:3n3 in females. The largest absolute differences were that females had about 0.3 percentage point higher and lower levels of C18:2n6 and C22:4n6 than males, respectively. Skewness and kurtosis were calculated for the individual fatty acids, and the following had distributions with an absolute kurtosis index >10 and/or a skew index >3, that is, C20:1, C18:3n3, C20:5n3, and C18:3n6, which became approximately Gaussian using a natural logarithm transformation. However, about 60% of C18:3n6 measurements were <0.1%, which is considered as the reliable detection limit for the GC method, and it appears that the log transformation simply produced normally distributed noise. Therefore, C18:3n6 was excluded from latent variable analysis. Even though univariate normality was reasonable for the individual fatty acids, multivariate normality was rejected () so robust ML method was implemented in Mplus [28].
 
Twosample test, the critical level alpha was set to 0.05/22 = 0.0023 for statistical significance using Bonferroni correction (shown in bold). Use the following mean (SD) of the logtransformed values when standardizing data as explained in the discussion section for C18:3n3 −6.38 (0.41) and C20:5n3 −5.04 (0.48). 
Fatty acid concentrations were standardized to produce a correlation matrix (Supplemental Table 2), and the MSA was calculated for each fatty acid and overall (the latter was initially 0.20). The fatty acid with the lowest MSA value was dropped from analysis until all fatty acid MSA values were >0.60. This resulted in the following fatty acids being sequentially excluded from the latent variable analysis: C18:1, C20:1, C18:2n6, C20:2n6, C20:3n6, C24:0, and C24:1. After excluding these fatty acids the overall MSA for the correlation matrix increased to 0.75. Additionally, C22:5n3 needed to be removed from the correlation matrix because it was causing the explained variability in C20:5n3 to be greater than 100%, that is, Heywood condition [38]; hence, 13 fatty acids remained. Afterwards three dimensions were identified for men and women with eigenvalues greater than one.
3.2. Confirmatory Factor Analysis
Model M0 allowed intercepts, loadings, and unique variances to be freely estimated for men and women. The absolute fit was good with a SRMSR of 0.035, and the fit was much better than a model with zero correlations since (Table 2). However, the fit measures which adjust for parsimony, that is, RMSEA and TLI, were not near acceptable ranges. In model M1 when the 13 fatty acid means were held constant between gender , the SBC increased by over 400 and the hypothesis was rejected using the chisquared difference testing between nested models , . As pointed out earlier since robust ML method was used due to lack of multivariate normality, the chisquared test statistic was modified [32]. Specifically, to compare model M1 nested in M0, the scaling correction factor was computed as . Then (Table 2). All other model fit measures deteriorated as well. These results importantly show that the multivariate fatty acid mean profile was not the same between genders for these fatty acids.
 
The total degrees of freedom (DF) = 2 * (91 fatty acid variances/covariances + 13 fatty acid means) = 208 parameters in all models; hence, the number of estimated parameters equals DF; SRMSR = standardized root mean square residual; SBC = Schwarz Bayesian criteria; RMSEA = root mean square error of approximation; CFI = Bentler Comparative Fit Index; TLI = TuckerLewis Index. 
Next the fatty acid correlation structure was tested for gender invariance in multiple steps. When comparing model M2 to M0 we concluded that the factor loadings were different between men and women = 92, (Table 2). To test for gender differences in the fatty acids’ variances, models M3 and M0 were compared. The fatty acids’ variances were not different between genders, = 20.7, respectively, . Next the latent variable covariance structure was tested , by comparing model M2 with M4 and found to be different between men and women (). Likewise the overall covariance structure was different ().
3.3. Structural Equation Modeling Constraints
The above results suggested differential fatty acid functioning for men and women, so separate models were developed by gender that allows the standardized latent variables scores to be compared between men and women. Model M3 had a good fit compared to M0 shown by chisquared difference testing, SRMSR, and CFI; however, the parsimony measures suggest there were still too many parameters. Model M3 is shown for men and women in Supplemental Tables 3 and 4, respectively; the latent variables were named for the fatty acids with the strongest correlations as PUFA, SATURATED, and TRANS FACTORS. When examining the loading estimates between men and women they were quite similar; there were only 3 parameters that differ by >0.10 which included the saturated fatty acids C14:0 and C18:0. These two fatty acids have slightly stronger correlations with the underlying latent variables in women than men. To further reduce the model complexity, constraints were placed on the loading matrix. A threshold of 0.15 was chosen, and parameters were constrained to zero if they had loadings below this threshold.
Correlations among dietary intakes of fatty acids were determined for the subset of 2332 participants with valid food frequency questionnaires [39] (Supplemental Table 5). Dietary intake (g/d) was available for 11 out of 13 RBC fatty acids included in the latent variable model, and C22:4n6 and C22:5n6 were not calculated from the diet. Since RBC fatty acids were correlated with corresponding dietary intakes of fatty acids, covariances were added to the fatty acids residual matrix, , to account for foods being composed of many different fatty acids. Being able to specify which residual covariances to include is a feature unique to structural equation modeling and cannot be accomplished in the context of confirmatory factor analysis. There were 14 strong dietary correlations () that were added to the model. The correlation between C18:1t and C18:2t was extremely high () and caused model convergence issues; therefore, it was subsequently removed.
The biosynthesis process is well known for omega3 and omega6 fatty acids [7]. Delta6 and delta5 desaturase activity is needed to convert C18:3n3 into C20:5n3 (Figure 1). Then delta6 desaturase (D6D) is required again to further convert C20:5n3 into C22:6n3. The amount of D6D available for the second conversion to synthesize C22:6n3 may be limited by a function of what is initially consumed to support converting C18:3n3 [40]. Therefore, the amount and variability of both C20:5n3 and C22:6n3 depend (to some extent) on the intake of the parent n3 fatty acid C18:3n3. In the omega6 fatty acid family, D6D is needed for C22:4n6 to synthesize into C22:5n6. These biochemical steps introduce structural elements into the SEM model, so these three covariances were added to fatty acid residuals matrix . However, omega3 and omega6 fatty acids also compete for the desaturase enzymes, and the omega3 fatty acids are the preferential substrates [6]. So with higher levels of C20:5n3 (whether by biosynthesis or fish oil consumption), D5D activity is inhibited (feedback inhibition, whereby the enzyme senses when enough product has been made and then shuts down). This slows the synthesis of C20:4n6 from C20:3n6. Likewise C22:6n3 and C22:5n6 compete for D6D. These two additional fatty acid covariances were added to the SEM residual matrix as well.
3.4. Final Structural Equation Model
After the above constraints were placed on the model, the resulting SEM (Model M6) had a significantly better fit than the unrestricted model by gender (M0), = 506 (). Additionally model M6 was the only model to have a “reasonable” fit with RMSEA <0.08. Also it was the only model to have CFI and TLI > 0.90. Model M6 had a total of estimated parameters, including the following gender specific 2 (23 loadings, 9 means, and 3 latent variable covariances) = 70 and gender invariant (13 residual variances and 13 dietaryrelated and 5 desaturaserelated residual covariances) . The loadings and latent variable correlations (Supplemental Table 6) and fatty acid residual matrix (Supplemental Table 7) are shown for model M6. One loading and three residual covariances were <0.05; these were set to zero for a more parsimonious model (M7). The nested fit between the reduced SEM model M7 was similar to model M6, = 4.3 () and all the parsimony fit measures (i.e, SBC, RMSEA, and TLI) were improved.
There were four fatty acids with gender mean differences which were not significantly different than zero (i.e, C22:6n3, C20:4n6, C16:0, and C18:1t). The greatest mean differences (all 0.30 to 0.40 SD) between genders were found in two PUFA [ln(C18:3n3) and C22:4n6], one saturated fatty acid [C14:0] and one monounsaturated fatty acid [C16:1]. The loadings, latent variable covariance structure, mean profile, and fatty acid residual matrix for model M7 are shown in Tables 3 and 4. The Mplus code for model M7 is given in Algorithm 1.
 
Direction of signs is arbitrary; the signs for the PUFA loadings and factor correlations have been reversed in this table to make more n3 fatty acid positively associated with the PUFA FACTOR. 
 
indicates dietary intakerelated covariances; indicates biosynthesisrelated covariances; fatty acid residual variances are shown on the main diagonal. 

4. Discussion
Although 22 individual fatty acids were measured during the GC process, 9 were removed from the latent variable analysis because they were not related to the core structure as explained above. However, the individual fatty acids that were removed may still have clinical utility as individual predictor variables. For example, C18:2n6 (Linoleic acid) was fairly independent of the PUFA, SATURATED, and TRANS latent variable scores, and all correlations were around 0.10 (Supplemental Table 8). Since C18:2n6 has been reported inversely related with heart disease [4], it could still have clinical utility as in independent predictor variable in combination with these newly defined latent variables. C22:5n3 (DPA) is intermediate of C20:5n3 (EPA) and C22:6n3 (DHA) in the biosynthesis process (Figure 1) and was the only excluded fatty acid that had a correlation >0.10 with the PUFA FACTOR (). However, DPA has less biological activity than the other marine fish oils [41], so it is not anticipated that DPA would be useful as an individual predictor variable of clinical outcomes. Afterwards, the remaining 13 fatty acids were found to be represented by three dimensions, which constituted nearly 70% of the total fatty acid abundance. In the SEM model, the residual variances were equal between genders. Likewise the structural dietary correlations and biosynthesis processes, which were accounted for with residual correlations, did not vary between men and women. The final SEM model fit the data well by all measures.
The PUFA FACTOR includes the following fatty acids: ln(C18:3n3), ln(C20:5n3), C22:6n3, C20:4n6, C22:4n6, and C22:5n6 (Figure 2); all of which are found in the PUFA biosynthesis processes shown in Figure 1. The loading directions of the fatty acids included in the PUFA FACTOR are also supported by many of the competing metrics being used in fatty acid research. The omega3 index is implemented in clinical laboratory testing and is defined as RBC C20:5n3 + C22:6n3 [42]. The omega3 index was an independent risk factor for allcause mortality in a study of stable coronary heart disease patients, with higher levels indicating reduced risk [3]. Lower amounts of the omega3 index were associated with depression in a casecontrol study of adolescents [43]. In the PUFA FACTOR, for both men and women, these two fatty acids operate in the same direction with similar magnitudes, which supports their summation as a biomarker (although C20:5n3 has been log transformed in the PUFA FACTOR).
(a)
(b)
The n6/n3 ratio [44], n6 HUFA/total HUFA ratio [45], and C20:4n6/C20:5n3 ratio [46] are all metrics that seek to combine individual fatty acids into more powerful predictors of risk. Although the goal is reasonable, these approaches are criticized as being imprecise and impractical [46]. All of these ratios may be flawed in that the same ratio can be obtained by increasing the numerator or decreasing the denominator, when these fatty acids do not have the same physiological properties. An improvement to these ratios may be the PUFA FACTOR. It is a more nuanced metric since it does not simply add up the masses of different PUFA families and create a ratio; it takes into account the relative strengths of relationship among these linearly “opposing” and interrelated fatty acids and reduces this nexus into a single number. The algorithm for scoring these latent variables from raw fatty acid data is given in the following.
Algorithm (algorithm for scoring latent variables)
Step 1. Measure fatty acids as a % of total fatty acids.
Step 2. For C18:3n3 and C20:5n3 transform the values using natural logarithm.
Step 3. Standardize all fatty acids using corresponding overall means and standard deviations from Table 1 into a row vector (as shown below, the headings indicate the required order).
Step 4. Calculate the latent variable scores for subject using the appropriate male or female matrices as where is the standardized fatty acid mean column vector given by gender in Table 3. is the latent variable covariance matrix given by gender in Table 3. is the loading matrix given by gender in Table 3. is the fatty acid residual covariance matrix given in Table 4. means to transpose the matrix. −1 means to take the inverse of the matrix.
Example. Measure RBC fatty acid percent weight composition using gas chromatography as detailed in Harris et al. [41] or similar, then log transform C18:3n3 and C20:5n3, and then standardize all raw data by (value − mean)/SD from Table 1 to produce a row vector:
Lastly the latent variable scores are derived by using, with the appropriate vector and matrices used for men or women given in Tables 3 and 4. If the blood sample was from a man or woman, the respective latent variables scores would be
Another approach has been to construct “desaturase ratios” which are based on the known biosynthetic relationships among PUFA [7]. Since it is far too invasive (requiring liver biopsy) to measure the activity of these enzymes directly, they have been estimated empirically by dividing RBC levels of product fatty acids by levels of precursor fatty acids. Thus the delta6 desaturase (D6D) activity can be estimated by the ratio of 20:3n6/C18:2n6 and the delta5 desaturase (D5D) activity by the ratio of C20:4n6/C20:3n6 (Figure 1). Interestingly, both of these desaturase ratios have been associated with risk for the development of Type 2 diabetes mellitus in a recent metareview [47]. The ultimate clinical utility of the PUFA FACTOR (versus desaturase or other fatty acid ratios) will be determined in future studies by comparing these metrics as predictors of disease outcomes for mortality, CHD events, development of type 2 diabetes or dementia, and so forth.
5. Conclusion
The PUFA FACTOR has much supporting evidence based on fatty acid metabolism and dietary patterns. It was also the first dimension extracted from the data, due to explaining the most variability (about 30% of total in men and women) for these 13 fatty acids. In a previous study these same Framingham subjects were included in a heritability analysis, and it was found that about 25% and 40% of the variance in two of the fatty acids included in the PUFA factor (i.e, EPA and DHA) was due to genetic and environment, respectively [41]. The PUFA factor can also be seen as a unifying theme among the various n3 and n6 metrics typically used in fatty acid research. Since n3 and n6 fatty acids have been implicated in cardiovascular diseases [4, 47], cognitive function [48], brain magnetic resonance imaging [49, 50], depression [43], mortality [1–3], and cellular aging [51] it is reasonable to expect the PUFA FACTOR to have clinical utility for predicting these outcomes. In contrast, the SATURATED and TRANS FACTORS had several cross loadings between them and even include some PUFA. Thus, their interpretations are unclear, which will likely limit their usefulness.
The strengths of this study include a wellcharacterized structural equation model applied to RBC fatty acid data which incorporates elements of both fatty acid metabolism and dietary intake patterns in defining the model. Additionally the correlation structure of the SEM was decomposed, and the separate components were tested for gender invariance. Another benefit was the use of a large, extensively studied cohort with enrichment for minorities (Framingham). Limitations include that the RBC measurements were from a particular GC method, and since national standards have not been established for measuring fatty acids the sensitivity of these results to other GC methods is unknown. This study measured erythrocytes; other blood fractions or sample types (e.g., whole blood, plasma, and plasma phospholipids) have different rank orders of fatty acid abundances and these may require unique structural equation models. The fit of this SEM needs be tested in independent samples to determine its generalizability beyond the Framingham Study.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge Ramachandran S. Vasan and Sander J. Robins at the Framingham Heart Study for their assistance in developing the overall research plan. The authors are supported in part by the National Heart Lung and Blood Institute (NHLBI; R01HL089590) and by contract N01HC25195, the Framingham Heart Study (NHLBI) and Boston University School of Medicine.
Supplementary Materials
The supplementary material includes tables of subject’s characteristics, RBC fatty acid and dietary intake correlation matrices, factor correlations with excluded RBC fatty acids, and interim model results.
References
 C. M. Albert, H. Campos, M. J. Stampfer et al., “Blood levels of longchain n3 fatty acids and the risk of sudden death,” The New England Journal of Medicine, vol. 346, no. 15, pp. 1113–1118, 2002. View at: Publisher Site  Google Scholar
 C. Wang, W. S. Harris, M. Chung et al., “n3 Fatty acids from fish or fishoil supplements, but not αlinolenic acid, benefit cardiovascular disease outcomes in primary and secondaryprevention studies: a systematic review,” American Journal of Clinical Nutrition, vol. 84, no. 1, pp. 5–17, 2006. View at: Google Scholar
 J. V. Pottala, S. Garg, B. E. Cohen, M. A. Whooley, and W. S. Harris, “Blood eicosapentaenoic and docosahexaenoic acids predict allcause mortality in patients with stable coronary heart disease: the heart and soul study,” Circulation: Cardiovascular Quality and Outcomes, vol. 3, no. 4, pp. 406–412, 2010. View at: Publisher Site  Google Scholar
 W. S. Harris, D. Mozaffarian, E. Rimm et al., “Omega6 fatty acids and risk for cardiovascular disease: a science advisory from the American Heart Association nutrition subcommittee of the council on nutrition, physical activity, and metabolism; council on cardiovascular nursing; and council on epidemiology and prevention,” Circulation, vol. 119, no. 6, pp. 902–907, 2009. View at: Publisher Site  Google Scholar
 N. T. Bendsen, S. Stender, P. B. Szecsi et al., “Effect of industrially produced trans fat on markers of systemic inflammation: evidence from a randomized trial in women,” Journal of Lipid Research, vol. 52, no. 10, pp. 1821–1828, 2011. View at: Publisher Site  Google Scholar
 J. B. Barham, M. B. Edens, A. N. Fonteh, M. M. Johnson, L. Easter, and F. H. Chilton, “Addition of eicosapentaenoic acid to γlinolenic acidsupplemented diets prevents serum arachidonic acid accumulation in humans,” Journal of Nutrition, vol. 130, no. 8, pp. 1925–1931, 2000. View at: Google Scholar
 J. R. Marszalek and H. F. Lodish, “Docosahexaenoic acid, fatty acidinteracting proteins, and neuronal function: breastmilk and fish are good for you,” Annual Review of Cell and Developmental Biology, vol. 21, pp. 633–657, 2005. View at: Publisher Site  Google Scholar
 D. Mozaffarian, A. Ascherio, F. B. Hu et al., “Interplay between different polyunsaturated fatty acids and risk of coronary heart disease in men,” Circulation, vol. 111, no. 2, pp. 157–164, 2005. View at: Publisher Site  Google Scholar
 R. J. T. Mocking, J. Assies, A. Lok et al., “Statistical methodological issues in handling of fatty acid data: percentage or concentration, imputation and indices,” Lipids, vol. 47, no. 5, pp. 541–547, 2012. View at: Publisher Site  Google Scholar
 C. K. Chow, “Fatty acid composition of plasma phospholipids and risk of prostate cancer,” The American Journal of Clinical Nutrition, vol. 89, no. 6, pp. 1946–1947, 2009. View at: Google Scholar
 F. L. Crowe, “Reply to CK Chow,” American Journal of Clinical Nutrition, vol. 89, no. 6, pp. 1946–1947, 2009. View at: Publisher Site  Google Scholar
 K. E. Bradbury, C. Murray Skeaff, T. J. Green, A. R. Gray, and F. L. Crowe, “The serum fatty acids myristic acid and linoleic acid are better predictors of serum cholesterol concentrations when measured as molecular percentages rather than as absolute concentrations,” American Journal of Clinical Nutrition, vol. 91, no. 2, pp. 398–405, 2010. View at: Publisher Site  Google Scholar
 H. Harris, “Omega3 fatty acids in cardiac biopsies from heart transplantation patients—correlation with erythrocytes and response to supplementation,” Circulation, vol. 110, no. 12, pp. 1645–1649, 2004. View at: Google Scholar
 L. Hodson, C. M. Skeaff, and B. A. Fielding, “Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake,” Progress in Lipid Research, vol. 47, no. 5, pp. 348–380, 2008. View at: Publisher Site  Google Scholar
 M. A. Beydoun, J. S. Kaufman, J. Ibrahim, J. A. Satia, and G. Heiss, “Measurement error adjustment in essential fatty acid intake from a food frequency questionnaire: alternative approaches and methods,” BMC Medical Research Methodology, vol. 7, article 41, 2007. View at: Publisher Site  Google Scholar
 M. Tournoud, R. Ecochard, J. Iwaz, J.P. Steghens, G. Bellon, and I. Durieu, “Structural equations to model relationships between pulmonary function, fatty acids and oxidation in cystic fibrosis,” Scandinavian Journal of Clinical and Laboratory Investigation, vol. 69, no. 1, pp. 36–44, 2009. View at: Publisher Site  Google Scholar
 J. C. N. Chan, P. C. Y. Tong, and J. A. J. H. Critchley, “The insulin resistance syndrome: mechanisms of clustering of cardiovascular risk,” Seminars in Vascular Medicine, vol. 2, no. 1, pp. 45–57, 2002. View at: Google Scholar
 J. Bradbury, L. Brooks, and S. P. Myers, “Are the adaptogenic effects of omega 3 fatty acids mediated via inhibition of proinflammatory cytokines?” EvidenceBased Complementary and Alternative Medicine, vol. 2012, Article ID 209197, 14 pages, 2012. View at: Publisher Site  Google Scholar
 P. W. F. Wilson, R. B. D'Agostino, D. Levy, A. M. Belanger, H. Silbershatz, and W. B. Kannel, “Prediction of coronary heart disease using risk factor categories,” Circulation, vol. 97, no. 18, pp. 1837–1847, 1998. View at: Google Scholar
 H. W. Marsh, B. Muthén, T. Asparouhov et al., “Exploratory structural equation modeling, integrating CFA and EFA: application to students' evaluations of university teaching,” Structural Equation Modeling, vol. 16, no. 3, pp. 439–476, 2009. View at: Publisher Site  Google Scholar
 H. W. Marsh, O. Lüdtke, B. Muthén et al., “A new look at the big five factor structure through exploratory structural equation modeling,” Psychological Assessment, vol. 22, no. 3, pp. 471–491, 2010. View at: Publisher Site  Google Scholar
 W. B. Kannel, M. Feinleib, and P. M. McNamara, “An investigation of coronary heart disease in families. The Framingham offspring study,” American Journal of Epidemiology, vol. 110, no. 3, pp. 281–290, 1979. View at: Google Scholar
 S. F. Quan, B. V. Howard, C. Iber et al., “The Sleep Heart Health Study: design, rationale, and methods,” Sleep, vol. 20, no. 12, pp. 1077–1085, 1997. View at: Google Scholar
 W. S. Harris, J. V. Pottala, R. S. Vasan et al., “Changes in erythrocyte membrane trans and marine fatty acids between 1999 and 2006 in older Americans,” Journal of Nutrition, vol. 142, no. 7, pp. 1297–1303, 2012. View at: Google Scholar
 A. Skrondal and S. RabeHesketh, Generalized Latent Variable Modeling: Multilevel, Longitudinal, and Structural Equation Modeling, Chapman & Hall, Boca Raton, Fla, USA, 2004.
 R. B. Kline, Principles and Practice of Structural Equation Modeling, The Guilford Press, New York, NY, USA, 2nd edition, 2005.
 K. V. Mardia, “Applications of some measures of multivariate skewness and kurtosis in testing normality and robustness studies,” Sankhya: The Indian Journal of Statistics B, vol. 36, pp. 115–128, 1974. View at: Google Scholar
 L. K. Muthen and B. O. Muthen, Mplus User's Guide, Muthen & Muthen, Los Angeles, Calif, USA, 6th edition, 2010.
 H. F. Kaiser, “An index of factorial simplicity,” Psychometrika, vol. 39, no. 1, pp. 31–36, 1974. View at: Publisher Site  Google Scholar
 M. A. Pett, N. R. Lackey, and J. J. Sullivan, Making Sense of Factor Analysis: The Use of Factor Analysis for Instrument Development in Health Care Research, Sage Publications, Thousand Oaks, Calif, USA, 2003.
 R. Johnson and D. Wichern, Applied Multivariate Statistical Analysis, Pearson Education, Upper Saddle River, NJ, USA, 6th edition, 2007.
 A. Satorra and P. M. Bentler, “A scaled difference chisquare test statistic for moment structure analysis,” Psychometrika, vol. 66, no. 4, pp. 507–514, 2001. View at: Google Scholar
 J. H. Steiger and J. M. Lind, “Statistically based tests for the number of common factors,” in Proceedings of the Annual Meeting Psychometric Society, Iowa City, IA, Iowa, USA, 1980. View at: Google Scholar
 P. M. Bentler, “Comparative fit indexes in structural models,” Psychological Bulletin, vol. 107, no. 2, pp. 238–246, 1990. View at: Google Scholar
 G. W. Cheung and R. B. Rensvold, “The effects of model parsimony and sampling error on the fit of structural equation models,” Organizational Research Methods, vol. 4, no. 3, pp. 236–264, 2001. View at: Google Scholar
 G. Schwarz, “Estimating the dimension of a model,” Annals of Statistics, vol. 6, pp. 461–464, 1978. View at: Google Scholar
 H. Akaike, “Information theory and an extension of maximum likelihood principal,” in Proceedings of the 2nd International Symposium of Information Theory and Control, Akademia Kiado, Budapest, Hungary, 1973. View at: Google Scholar
 F. Chen, K. A. Bollen, P. Paxton, P. J. Curran, and J. B. Kirby, “Improper solutions in structural equation models: causes, consequences, and strategies,” Sociological Methods and Research, vol. 29, no. 4, pp. 468–508, 2001. View at: Google Scholar
 M. E. Rumawas, J. T. Dwyer, N. M. Mckeown, J. B. Meigs, G. Rogers, and P. F. Jacques, “The development of the mediterraneanstyle dietary pattern score and its application to the american diet in the framingham offspring cohort,” Journal of Nutrition, vol. 139, no. 6, pp. 1150–1156, 2009. View at: Publisher Site  Google Scholar
 R. Portolesi, B. C. Powell, and R. A. Gibson, “Competition between 24:5n3 and ALA for Δ6 desaturase may limit the accumulation of DHA in HepG2 cell membranes,” Journal of Lipid Research, vol. 48, no. 7, pp. 1592–1598, 2007. View at: Publisher Site  Google Scholar
 W. S. Harris, J. V. Pottala, S. M. Lacey et al., “Clinical correlates and heritability of erythrocyte eicosapentaenoic and docosahexaenoic acid content in the Framingham Heart Study,” Atherosclerosis, vol. 225, no. 2, pp. 425–431, 2012. View at: Google Scholar
 W. S. Harris and C. Von Schacky, “The Omega3 Index: a new risk factor for death from coronary heart disease?” Preventive Medicine, vol. 39, no. 1, pp. 212–220, 2004. View at: Publisher Site  Google Scholar
 J. V. Pottala, J. A. Talley, S. W. Churchill, D. A. Lynch, C. von Schacky, and W. S. Harris, “Red blood cell fatty acids are associated with depression in a casecontrol study of adolescents,” Prostaglandins Leukotrienes and Essential Fatty Acids, vol. 86, no. 45, pp. 161–165, 2012. View at: Publisher Site  Google Scholar
 A. P. Simopoulos, “The importance of the omega6/omega3 fatty acid ratio in cardiovascular disease and other chronic diseases,” Experimental Biology and Medicine, vol. 233, no. 6, pp. 674–688, 2008. View at: Publisher Site  Google Scholar
 W. E. M. Lands, “Diets could prevent many diseases,” Lipids, vol. 38, no. 4, pp. 317–321, 2003. View at: Publisher Site  Google Scholar
 W. S. Harris, “The omega6/omega3 ratio and cardiovascular disease risk: uses and abuses,” Current Atherosclerosis Reports, vol. 8, no. 6, pp. 453–459, 2006. View at: Publisher Site  Google Scholar
 J. Kröger and M. B. Schulze, “Recent insights into the relation of Δ5 desaturase and Δ6 desaturase activity to the development of type 2 diabetes,” Current Opinion in Lipidology, vol. 23, no. 1, pp. 4–10, 2012. View at: Publisher Site  Google Scholar
 J. G. Robinson, N. Ijioma, and W. Harris, “Omega3 fatty acids and cognitive function in women,” Women's Health (London, England), vol. 6, no. 1, pp. 119–134, 2010. View at: Publisher Site  Google Scholar
 J. V. Pottala, K. Yaffe, J. G. Robinson et al., “Higher RBC EPA + DHA corresponds with larger total brain and hippocampal volumes: WHIMSMRI Study,” Neurology, vol. 82, no. 5, pp. 435–442, 2014. View at: Google Scholar
 Z. S. Tan, W. S. Harris, A. S. Beiser et al., “Red blood cell omega3 fatty acid levels and markers of accelerated brain aging,” Neurology, vol. 78, no. 9, pp. 658–664, 2012. View at: Publisher Site  Google Scholar
 R. FarzanehFar, J. Lin, E. S. Epel, W. S. Harris, E. H. Blackburn, and M. A. Whooley, “Association of marine omega3 fatty acid levels with telomeric aging in patients with coronary heart disease,” JAMA—Journal of the American Medical Association, vol. 303, no. 3, pp. 250–257, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 James V. Pottala et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.