Sara Beigrezaei1,2
Sara Jambarsang3Sayyed Saeid Khayyatzadeh1,2Masoud Mirzaei4Amir Houshang Mehrparvar5Amin Salehi-Abargouei1,2,4*†- 1Nutrition and Food Security Research Center, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
- 2Department of Nutrition, Faculty of Health, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
- 3Center for Healthcare Data Modeling, Department of Biostatistics and Epidemiology, School of Public Health, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
- 4Yazd Cardiovascular Research Center, Non-communicable Disease Research Institute, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
- 5Industrial Diseases Research Center, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
Limited data exist on the advantage of data reduction hybrid methodologies for evaluating the relationship between dietary patterns (DPs) and chronic diseases and they have led to inconsistent results. This study aimed to investigate the association between DPs extracted using principal component analysis (PCA), partial least-squares (PLS), and reduced-rank regressions (RRRs) in identifying DPs associated with hypertension (HTN) risk. The current study was conducted in the context of two cohort studies accomplished in Iran. DPs were generated using PCA, PLS, and RRR methods. Log-binomial logistic regression test was used to assess the association between DPs and the risk of HTN. From a total of 12,403 included participants aged 20–70 years, 507 incident cases of confirmed HTN were identified. The PCA-DP2 was associated with HTN in the fully adjusted model (T3 vs. T1: RR: 0.737, 95% CI: 0.57–0.93, Ptrend = 0.013). The PLS-DP2 and HTN risk were inversely associated in the multivariate model (T3 vs. T1: RR: 0.704, 95% CI: 0.54–0.91, Ptrend = 0.013). The RRR-DP2 was associated with an increased risk of HTN (T3 vs. T1: RR: 1.412, 95% CI: 1.11–1.80, Ptrend = 0.007). Our findings suggest that the RRR method reveals stronger results in association with HTN risk. However, further investigations are required to confirm the association between DPs derived by PLS and RRR methods by incorporating biomarkers related to HTN as the response variables.
Introduction
Hypertension (HTN), or elevated blood pressure, is a leading risk factor for cardiovascular diseases (CVDs) and premature death (1). It is estimated that 1.28 billion adults aged 30–79 years have HTN, and its prevalence is increasing worldwide (2). Evidence has well established that HTN is associated with genetic and environmental factors such as obesity, tobacco use, sedentary lifestyle, and diet (3). A healthful lifestyle is a key strategy for the prevention of HTN, and diet is known as a modifiable lifestyle component with the strongest impact on blood pressure (4).
The majority of previous studies that have investigated the relationship between diet and HTN focused on the consumption of a single nutrient, food item, or food group (5, 6); however, more recently, nutritional epidemiologic investigations of chronic diseases have gone beyond these approaches, focusing on evaluating the effect of the overall diet (7).
A dietary patterns (DPs) approach, as a complementary approach, can provide helpful findings, translating into dietary guidelines without fully understanding the mechanisms underlying the observed relationships (7). Most previous studies have investigated DPs generated by principal component analysis (PCA) (8). Recently, two alternative reduction techniques have been proposed for DPs research: reduced rank regression (RRRs) and partial least-square analysis (PLS) (9). The RRR creates linear combinations of predictor variables with as much variation as possible in a combined set of response variables including biomarkers or macro/micronutrients (9). The PLS is a compromise between the PCA and RRR methods which explain variation in predictor variables and selected intermediate response variables (9). The main advantage of these methods is that the derived DPs are correlated with biological pathways related to the etiology of diseases (9).
Few studies have assessed the relationship between DPs and HTN via the RRR method (10, 11). In a cross-sectional study by da Silva et al. (10) was found that healthier DP from the RRR method was directly related to HTN in women 40–60 years. There is no previous study investigating HTN development risk and the PLS method, specifically. However, Naja et al. (12) reported that the western and traditional Lebanese DPs were associated with higher HTN risk in Lebanese men. Previous studies suggested that using the RRR method resulted in a more significant link with the HTN compared to PCA and PLS methods (13, 14), whereas, PCA-derived patterns were more associated with the real dietary habits in their study population (9).
Therefore, the main objective of the present study was to assess whether two hybrid methods with intermediate variables for HTN generate different DPs (RRR and PLS) compared to a common method (PCA).
Materials and methods
Study design and study population
The Yazd Health Study (YaHS) was conducted in September 2014 in Yazd greater area located in central Iran, including 9,962 adults aged 20–70 years were participated in the enrollment phase. The assessment of dietary intake was separately collected in Taghzieh Mardom Yazd (TAMYZ) study by using a validated semi-quantitative food frequency questionnaire (FFQ) (15). The Shahedieh cohort study is a part of a large Persian multicentral cohort study (Persian cohort) performed on 180,000 participants in 18 various geographical areas of Iran (16). The Shahedieh study was established in 2014 and 9,977 adults aged 35–70 years were enrolled to the study at baseline. To collect the dietary intake of participants, a semi-quantitative FFQ was completed. Data on demographic characteristics, medical history, smoking status, and physical activity was collected in both cohort studies. The study protocol for YaHS-TAMYZ (17) and Shahedieh (16) studies are completely described elsewhere.
The exclusion criteria were: participants who reported incomplete dietary intakes data or an implausible total energy intake < 800 kcal/d or ≥ 6,000 kcal/d (YaHS-TAMYZ study, n = 1,692, Shahedieh study, n = 1,708), those who had a previous diagnosis of HTN (YaHS-TAMYZ study, n = 1,773, Shahedieh study, n = 2,054), those who had not provided data on National Identifier (YaHS-TAMYZ study, n = 287, Shahedieh study, n = 0), had no follow-up data (YaHS-TAMYZ study, n = 0, Shahedieh study, n = 0), and people who died (YaHS-TAMYZ study, n = 50, Shahedieh study, n = 28) were excluded, which left 12,403 participants (YaHS-TAMYZ study: 6,216, Shahedieh study: 6,187) for current analyses.
All participants filled an informed consent before enrolling the cohort studies. Both studies were approved by the research Council of Shahid Sadoughi University of Medical Sciences. The present study was ethically approved by Shahid Sadoughi University’s Ethics Committee (approval code: IR.SSU.SPH.REC.1399.197).
Dietary assessment and food groups
Dietary intake in the YaHS-TAMYZ study was examined by a validated, multiple-choice semi-quantitative FFQ (15). The food item list consisted of 178 food or food groups representing commonly consumed foods in Yazd. Participants were asked by trained interviewers to report the frequency of food item intake during the past year by answering 10-multiple-choice frequency responses ranging from “never or less than once a month” to “10 or more times per day.” Moreover, the FFQ had five choices for portion size for estimating the intake of each food item (18). In Shahedieh study, dietary intake information was collected by a semi-quantitative open-ended FFQ based on 134-items through face-to-face interviews with the participants. Participants of the Shahedieh study were asked to report how often on average over the past year they consumed a typical portion size of each food item with multiple possible responses on a “daily,” “weekly,” or “monthly” basis. The frequency and portion size reported for each food item were converted to grams per day using household measures (19). The United States Department of Agriculture (USDA) food composition database was used to estimate daily intake of energy and nutrient for each participant (20). Food items were merged into 32 food groups based on food items similarity in their nutrient profiles and are presented in Table 1.
Table 1. Food grouping used in the dietary pattern (DP) analyses in the Yazd Health Study-Taghzieh Mardom Yazd (YaHS-TAMYZ) and Shahedieh cohort studies.
Assessment of other covariates
To collect data on age, sex, educational status, marital status, smoking status, and the history of chronic diseases, a similar demographic questionnaire in both YaHS-TAMYZ and Shahedieh cohort studies was used.
Five categories (20–30, 30–40, 40–50, 50–60, and ≥ 60 years) was considered for age variable. Educational status of participants was classified into four levels (Uneducated, Elementary or guidance school, High school diploma, and B.Sc. or higher academic degree). Smoker participants were defined as current smokers, while non-smoker participants categorized as former smokers and never smokers. Marital status was defined into three categories (Single, Married, and Widowed or divorced).
The height and body weight of the study participants were measured in both YaHS-TAMYS and Shahedieh studies. In the YaHS-TAMYS study, body weight (kg) was measured by using an Omron BF511 portable digital scale (Omron Inc., Nagoya, Japan) with the accuracy of 0.1 kg, while standing on the middle of the scale, without assistance and with minimum clothing and height (cm) was measured in a standing position using a tape measure on a straight wall to the nearest centimeter. In the Shahedieh study, body weight (kg) and height (cm) were measured by using the National Institute of Health protocols by trained staffs. Body weight was measured while the participants were with minimum clothing and without shoes by using a digital scale (SECA, model 755, Germany). Height was measured by using a measure tape attached to a flat wall with the accuracy of 0.5 cm. Body mass index (BMI) was calculated as weight (kg)/height squared (m2).
In the YaHS-TAMYZ cohort study, the short version of the International Physical Activity Questionnaire (IPAQ) was used to measure physical activity level (21). In the Shahedieh cohort study, participants were asked about their usual physical activity levels in the last year and in case they had seasonal jobs (22). Physical activity was expressed as metabolic equivalent hours per week (MET-h/wk).
Clinical assessment
Follow-up of study participants and case confirmation
In the YaHS-TAMYZ cohort study, information on death events was collected using data from population-based registries and linked outcome information from the aggregated hospital information system (Samanah Electronici PArvandeh Salamat-SEPAS) which covers 100% of public hospitals and the majority of private hospitals in Yazd province. The data was obtained from SEPAS by using the National Identifier number of each participant to link data.
During follow-up time in the Shahedieh cohort study, participants received annual phone calls and follow-up questionnaires were completed in terms of the death events. Moreover, a verbal autopsy form validated in the Iranian population was completed during the death events. Two trained internists examined the medical documents to determine the cause of death. In case of inconsistency, a third internist performed a final assessment of the documents to reach a final decision.
In both studies, new cases of HTN were defined as individuals reporting a physician-based diagnosis of HTN and antihypertensive medication use in the annual follow-up phone call who did not have HTN at baseline.
Statistical analysis
Dietary patterns analysis
To identify DPs out of 32 food groups, three complementary data reduction techniques, including PCA, PLS, and RRR were used. An explained variance calculates the proportion of variation in a DP which can be attributable to the food groups or responses. In the PCA method, the DPs explain as much variation as possible of the food groups. To obtain the variance of the response variables explained by the DPs derived from the PCA method, linear regression analysis of the DPs scores, and response variables was used. The RRR method identifies linear functions of predictors (food groups) which explain as much intermediate responses variation as possible using a covariance matrix of predictors and responses in calculating the DPs scores. The PLS method combines PCA and RRR methods and calculates DPs scores considering both the predictor and response matrices; therefore, the explained variance of both food groups and response variables is expected to be between the PCA and RRR methods.
The number of DPs initially generated by PCA is constrained by the number of food groups used (23); nevertheless, we retained only three DPs from PCA for subsequent analyses based on the scree plot, an eigenvalue (>1), and the interpretability of the principal DPs (24). To obtain orthogonal DPs and increase the interpretability of principle DPs, Varimax rotation was applied. Kaiser–Meyer–Olkin (KMO) test was performed for checking sample adequacy.
Based on previous studies, dietary intake of sodium, potassium, and saturated fatty acids were selected as the response variables for PLS and RRR.
The SAS procedure PLS was used to conduct PLS and RRR analysis. The number of DPs derived by PLS and RRR is restricted by the number of intermediate response variables used; thus, three DPs were specified in each method. For both PLS and RRR, we calculated the continuous DPs scores (the linear functions of food groups) in the subsequent analyses and interpretations.
Descriptive analyses and modeling
Quantitative and qualitative variables were compared between participants with and without HTN by using independent sample t-test and chi-square tests, respectively. For further analyses, DP scores were categorized into tertiles. To assess the association between DPs derived by PCA, PLS, and RRR analyses and the risk of HTN incidence, binomial logistic regression was used. In all analyses, in addition to the crude model, the first model was adjusted for age, sex, and energy intake; the second model was additionally controlled for education, marital status, smoking status, and physical activity; and in the third model, BMI was additionally adjusted.
All statistical analyses were performed with SAS Version 8.02 (SAS Institute, Cary, NC, USA) software. P values less than 0.05 were considered statistically significant.
Comparison of dietary pattern methods
A comparison between PCA, PLS, and RRR methods was conducted according to the relative factor loading within each DPs and its association with HTN risk. In addition, these three methods were examined based on the magnitude of variation of each method which explained the food groups and responses.
Results
Of the 12,403 participants included in the study analyses, about 54% were men (Table 2). After 6 years of follow-up for the YaHS-TAMYZ study and 4 years for the Shahedieh study, 507 participants were diagnosed with HTN. In addition, Table 2 represents the baseline characteristics of the study population. In addition, age, sex, and energy-adjusted dietary macronutrients and micronutrients intake of the participants according to tertiles of DPs derived by three methods are presented in the Supplementary Tables 1–3.
Table 2. Baseline characteristics of participants with and without hypertension (HTN)a,b.
Dietary patterns derived by PCA, PLS, and RRR methods
Dietary patterns factor loadings of the 32 food groups for the PCA, PLS, and RRR are presented in Table 3. Using the PCA method, three DPs were generated characterized by: PCA-DP1: high intake of processed meats, organ meats, fish, margarine, fruit juice, pizza, snacks, sweet dessert, and soft drinks, and low intake of whole grains; PCA-DP2: high intake of red meats, dairy products, fruits, tomatoes, potatoes, refined grains, vegetable oils, and pickles; PCA-DP3: high intake of tea, nuts, mayonnaise, snacks, sweet dessert, soft drinks, and sugars.
Table 3. Factor loadings of food groups in dietary patterns identified using principal component analysis (PCA), partial least-squares (PLS), and reduced-rank regression (RRR) methods.
The percentage of explained variance by the DPs is presented in Table 4. The three PCA-derived DPs accounted for 11.165, 7.024, and 5.251%, respectively of the variance in food groups and 10.524, 25.338, and 4.086% of the variance in the response variables, respectively.
Table 4. Explained variation in food groups and responses using principal component analysis (PCA), partial least-squares (PLS), and reduced-rank regression (RRR).
The first PLS DP, labeled “PLS-DP1,” was characterized by a high intake of organ meats, dairy products, fruits, fruit juice, other vegetables, dried fruit, nuts, and sweet dessert. A high intake of tea, tomatoes, potatoes, refined grains, whole grains, and vegetable oils and a low intake of sweet dessert, processed meats, organ meats, and margarine characterized the second PLS DP, “PLS-DP2.” The third PLS DP, identified as “PLS-DP3” had high intakes of refined grains, whole grains, hydrogenated fats, and sugars and a low intake of fruits, fruit juice, other vegetables, and dried fruit. The three PLS DPs explained 39.863, 6.645, and 7.324% of the variance in the response variables, respectively.
Regarding PLS results, all DPs accounted for 21.064% of the variation in food groups (8.818, 8.440, and 3.805%, respectively).
The first DP derived by the RRR method (RRR-DP1) was characterized by high intakes of dairy products, fruits, other vegetables, mayonnaise, sweet dessert, and low intake of pizza, and snacks explained 46.840% of the variation in the response variables. The second RRR-derived DP labeled “RRR-DP2” was characterized by a high intake of processed meats, red meats, margarine, hydrogenated fats, and sweet dessert, and a low intake of fruits, tomatoes, other vegetables, whole grains, and mayonnaise, which explained 9.147% of the variance in the responses. The third DP (RRR-DP3) was characterized by a high intake of whole grains, nuts, and sugars and a low intake of red meats, organ meats, margarine, dairy products, fruits, and fruit juice, which explained 2.718% of the variance in the responses.
Findings from the RRR method indicated that three DPs explained 15.212% of the total variation in food groups (6.798, 4.074, and 4.339%, respectively).
Association between dietary patterns and HTN
Table 5 represents the risk of HTN for each tertile of the DP scores. Highest adherence to the PCA-DP1 and the PCA-DP3 in comparison with lowest adherence were associated with a decreased risk of HTN only in the crude model (RR: 0.788, 95% CI: 0.63–0.98, Ptrend = 0.013 and RR: 0.728, 95% CI: 0.57–0.93, Ptrend = 0.013, respectively). PCA-DP2 was significantly associated with HTN in the crude and all adjusted models (crude modelT3vs.T1: RR = 0.733, 95% CI: 0.57–0.94, Ptrend = 0.011; Model IT3vs.T1: RR = 0.770, 95% CI: 0.60–0.98, = 0.034; Model IIT3vs.T1: RR = 0.756, 95% CI: 0.59–0.96, Ptrend = 0.024, and Model IIIT3vs.T1: RR = 0.737, 95% CI: 0.57–0.93, Ptrend = 0.013).
Table 5. Risk ratios and 95% confidence intervals for hypertension (HTN) and dietary patterns (DPs) derived using principal component analysis (PCA), partial least-squares (PLS), and reduced-rank regression (RRR).
The RRs for the HTN in the crude model was 0.689 (95% CI: 0.52–0.90, Ptrend = 0.006) and 0.734 (95% CI: 0.57–0.95, Ptrend = 0.029) for the participants in the highest tertile compared to those in the lowest tertile of adherence to the PLS-DP1 and the PLS-DP3, respectively; however, these associations were not observed in all adjusted models. The association the PLS-DP2 was inversely significant in all adjusted models with HTN risk (Model IT3vs.T1: RR = 0.784, 95% CI: 0.61–0.99, Ptrend = 0.056; Model IIT3vs.T1: RR = 0.738, 95% CI: 0.56–0.95, Ptrend = 0.031, and Model IIIT3vs.T1: RR = 0.704, 95% CI: 0.54–0.91, Ptrend = 0.013).
Compared with the participants in the lowest tertile, those in the third tertile of RRR-DP2 adherence had 1.412 times (95% CI: 1.11–1.80, Ptrend = 0.007) higher risks of HTN after adjustment for all confounders. In addition, this association was also significant in Model IT3vs.T1 (RR = 1.299, 95% CI: 1.03–1.64, Ptrend = 0.033) and Model IIT3vs.T1 (RR = 1.348, 95% CI: 1.05–1.72, Ptrend = 0.020). Highest adherence to RRR-DP1 in comparison to lowest adherence was associated with a reduced risk of HTN only in crude modelT3vs.T1 (RR = 0.692, 95% CI: 0.52–0.92, Ptrend = 0.019).
Discussion
The aim of our study was to investigate the potential association between DPs and the risk of HTN incidence using three different methods to generate DPs. The PCA-DP2 and the PLS-DP2 were significantly associated with HTN, however, compared to the participants in the lowest tertile, those in the third tertile of the RRR-DP2 adherence had 1.414 times higher odds of HTN incidence. Though all three DPs have an association with HTN risk, the RRR method showed a stronger statistical relationship.
Our findings showed that the highest adherence to the RRR-DP2 in comparison to the lowest adherence was associated with an increased risk of HTN. The RRR-DP2 was characterized by high intakes of processed and red meats, margarine, hydrogenated fats, and sweet dessert and a low intake of fruits, tomatoes, other vegetables, whole grains, and mayonnaise. In line with our findings, Du and Cao (25) reported that RRR derived-DP associated with C-reactive protein (CRP), which was high in red meat, snacks, and nuts, and low in grains, poultry, and fish led to a higher risk of high-normal and high levels of systolic blood pressure. In contrast, da Silva et al. (10) found that the RRR-derived DP associated with sodium, potassium, and saturated fat is characterized by a high intake of healthy food items such as banana, papaya, apple, orange, pumpkin, green leafy vegetables, other vegetables, and natural fruit juice in contrast with a low intake of unhealthy foods such as processed meat presented a direct association with HTN women between 40 and 60 years.
The mechanisms of red meat intake in increasing blood pressure are complex and attributed to multiple factors. The high content of saturated fatty acids and cholesterol in red meat has been directly associated with HTN in previous studies (26, 27). It is demonstrated that dietary cholesterol is directly associated with endothelial dysfunction and arterial stiffness (28). Evidence from animal studies has shown that increased oxidative stress in cholesterol-fed rabbits was strongly attributed to endothelial dysfunction (29). In addition, some evidence has demonstrated that Maillard reaction products during cooking or processing of meat, such as advanced glycosylation end-products and heterocyclic amines, may induce mediators of inflammation and oxidation pathways and impair insulin activity which can contribute to elevated blood pressure (23, 30). A direct association between dietary advanced glycosylation endproducts (AGEs) and CRP, tumor necrosis factor-a (TNF-a), vascular cell adhesion molecule 1, and 8-isoprostane (24, 31) is proposed. Vascular damage caused by AGEs can lead to increased vascular and myocardial stiffening, oxidative stress, and inflammation (32). Moreover, HTN is related to the low circulating level of the soluble receptor for AGE (23). The metabolism of L-carnitine, as abundant trimethylamine in red meat, by gut microbiota produces trimethylamine-N-oxide and led to atherosclerosis in mice (33). Though, processed meats are rich in sodium and nitrite contents, which are used as food preservatives (34); it is established that sodium intake has the main effect on increasing blood pressure (35).
Several prospective cohort studies have shown that long-term consumption of all types of fruits was relative to a reduced risk of HTN in adults. Tsubota-Utsugi et al. (36) found that subjects in the highest quintile in comparison to those in the lowest quintile of fruit consumption had 60% lower odds of HTN. There are several mechanisms of fruits and vegetables related to reducing blood pressure. The high content of polyphenols like flavonoids and anthocyanins in some fruits and vegetables can reduce the risk of HTN (37, 38). In addition, it has been revealed that grape polyphenols may contribute to vasorelaxation and decrease blood pressure (39).
The findings of a meta-analysis by Schwingshackl et al. (40) showed a significant inverse relationship between whole grain consumption and HTN in the linear and non-linear dose-response analysis. Sawicki et al. (41) reported that whole grain intake was inversely associated with systolic blood pressure, but not refined grains intake. In addition, a meta-analysis of prospective studies revealed an inverse association between whole grain intake and HTN risk, although no significant link was observed between refined grains intake and the risk of HTN (40). According to the previous evidence, the components of whole grains, including phytochemicals, fiber, and minerals like potassium, magnesium, zinc, and selenium can decrease blood pressure (42). Moreover, the high content of fiber in whole grains can prevent the development of HTN through the alteration in gut microbiota (43).
In the present study, high dairy product intake is one of the features of the second PLS-derived DP associated with a lower risk of HTN. This finding supports the previous evidence on the inverse relationship between dairy intake and HTN risk. Several mechanisms may involve in contributing to the beneficial effect of dairy products on blood pressure. Dairies are rich in calcium, magnesium, potassium, and phosphorus and all of them are contributed to reduce blood pressure through vascular resistance regulation and promotion of vasodilation (44). Potential evidence proposed that bioactive peptides from dairy protein may inhibit the effect of angiotensin I–converting enzyme on constriction of blood vessels (45).
In the present study, the PLS-DP2 which was high in tea, tomatoes, potatoes, refined grains, whole grains, and vegetable oils, but low in processed meats, organ meats, margarine, and sweet dessert was inversely related to a lower risk of HTN. The food groups loaded in this DP are similar to the RRR-DP2, and both of them are associated with our selected response variables including sodium, potassium, and saturated fatty acids. One of the specific food items of the PLS-DP2 related to blood pressure is tea consumption. Evidence has shown that habitual tea drinkers had a 14% lower risk of HTN incidence in comparison with non-habitual tea drinkers (46). It is demonstrated that tea had anti-inflammatory and endothelial function improvement effects through a wide range of antioxidants and polyphenols like flavonoids (47).
The DP derived by the PCA method (PCA-DP2) associated with a lower risk of HTN was a mix of healthy and unhealthy food groups of usual DPs, including a high intake of red meats, dairy products, fruits, and tomatoes, potatoes, refined grains, vegetable oils, and pickles. In a similar way, in a cross-sectional study of India, the “fruits-vegetables-sweets-snacks” pattern had a positive link with diastolic blood pressure (48). Given the interpretation of the DPs, the PCA-derived DPs reflected a wide extent of the dietary behavior of the studied population (9). The behavior-related DPs derived by the PCA method cannot predict the risk of diseases (9).
There are some previous studies that have found that RRR-derived DPs had a stronger link with the disease compared to those derived by PCA or PLS; it clearly could be explained because DPs from PLS and RRR methods are generated via disease-related response variables. Naja et al. (12) reported that a western DP from all PCA, PLS, and RRR methods was associated with an increased risk of HTN in the Lebanese population. However, a cross-sectional study by DiBello et al. (49) showed that PCA generated more myocardial infarction–associated DPs than RRR that used nutrient intake as response variables.
There are several limitations to the present study. In our study, the dietary intakes were evaluated using the FFQ which is prone to measurement error. Moreover, data on dietary intake was collected only at the baseline of both cohort studies and we could not warrant that participants have not changed their diets during the follow-up time. The follow-up time was 6 years for the YaHS-TAMYZ study and 4 years for the Shahedieh study; although the follow-up time is considered in the analyses, its difference between the two studies might be considered as a limitation. Moreover, both cohort studies had less than a 6-year follow-up time; hence, we might not reveal the long-term effects of DPs on HTN risk. Additionally, the limited number of incident HTN cases could affect the observed associations in the present study. In addition, using different evaluation instruments in the two cohorts could cause bias in the analyses. Furthermore, the IPAQ is only validated on individuals up to 65 years of age, however, the age range of the present study is higher and physical activity data collected on participants over 65 are not valid. We tried to adjust a wide variety of confounders, however, the role of residual confounding factors cannot be completely ruled out in observational studies.
Altogether, DP derived by PCA associated with HTN risk was more similar to common dietary behavior of the studied population and DPs derived by PLS and RRR methods reflected food groups in association with response variables. Though all three DPs have an association with HTN risk, the RRR method showed a stronger statistical relationship. However, further investigations are required to confirm the association between DPs derived by the PLS and the RRR methods by incorporating biomarkers related to HTN as the response variables.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Shahid Sadoughi University of Medical Sciences. The patients/participants provided their written informed consent to participate in this study.
Author contributions
AS-A and SB conceived and designed the study. SJ, SB, and AS-A conducted and interpreted the statistical analyses. SB wrote the first draft of the manuscript. AS-A, SJ, SK, MM, and AM critically reviewed the manuscript. All authors read and approved the final version of the manuscript.
Funding
This present study was funded by the Shahid Sadoughi University of Medical Sciences.
Acknowledgments
We thank all participants in the YaHS-TAMYZ and Shahedieh cohort studies, and authorities of Shahid Sadoughi University of Medical Sciences for their cooperation.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.1076723/full#supplementary-material
References
1. GBD 2017 Risk Factor Collaborators. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990-2017: a systematic analysis for the global burden of disease study 2017. Lancet. (2018) 392:1923–94.
3. Whelton P, Einhorn P, Muntner P, Appel L, Cushman W, Diez Roux A, et al. Research needs to improve hypertension treatment and control in African Americans. Hypertension. (2016) 68:1066–72.
5. Grillo A, Salvi L, Coruzzi P, Salvi P, Parati G. Sodium intake and hypertension. Nutrients. (2019) 11:1970.
6. Zhang Y, Zhang D. Red meat, poultry, and egg consumption with the risk of hypertension: a meta-analysis of prospective cohort studies. J Hum Hypertens. (2018) 32:507–17. doi: 10.1038/s41371-018-0068-8
7. Hu F. Dietary pattern analysis: a new direction in nutritional epidemiology. Curr Opin Lipidol. (2002) 13:3–9.
8. Tatoli R, Lampignano L, Bortone I, Donghia R, Castellana F, Zupo R, et al. Dietary patterns associated with diabetes in an older population from southern italy using an unsupervised learning approach. Sensors. (2022) 22:2193. doi: 10.3390/s22062193
9. Hoffmann K, Boeing H, Boffetta P, Nagel G, Orfanos P, Ferrari P, et al. Comparison of two statistical approaches to predict all-cause mortality by dietary patterns in German elderly subjects. Br J Nutr. (2005) 93:709–16. doi: 10.1079/bjn20051399
10. da Silva B, Neutzling M, Camey S, Olinto M. Dietary patterns and hypertension: a population-based study with women from Southern Brazil. Cad Saude Publica. (2014) 30:961–71.
11. Sauvageot N, Leite S, Alkerwi A, Sisanni L, Zannad F, Saverio S, et al. Association of empirically derived dietary patterns with cardiovascular risk factors: a comparison of PCA and RRR methods. PLoS One. (2016) 11:e0161298. doi: 10.1371/journal.pone.0161298
12. Naja F, Itani L, Hwalla N, Sibai A, Kharroubi S. Identification of dietary patterns associated with elevated blood pressure among Lebanese men: a comparison of principal component analysis with reduced rank regression and partial least square methods. PLoS One. (2019) 14:e0220942. doi: 10.1371/journal.pone.0220942
13. Melaku Y, Gill T, Taylor A, Adams R, Shi Z. A comparison of principal component analysis, partial least-squares and reduced-rank regressions in the identification of dietary patterns associated with bone mass in ageing Australians. Eur J Nutr. (2018) 57:1969–83. doi: 10.1007/s00394-017-1478-z
14. Biesbroek S, van der A, Brosens M, Beulens J, Verschuren W, van der Schouw Y, et al. Identifying cardiovascular risk factor-related dietary patterns with reduced rank regression and random forest in the EPIC-NL cohort. Am J Clin Nutr. (2015) 102:146–54. doi: 10.3945/ajcn.114.092288
15. Mirmiran P, Esfahani F, Mehrabi Y, Hedayati M, Azizi F. Reliability and relative validity of an FFQ for nutrients in the Tehran lipid and glucose study. Public Health Nutr. (2010) 13:654–62. doi: 10.1017/S1368980009991698
16. Poustchi H, Eghtesad S, Kamangar F, Etemadi A, Keshtkar A, Hekmatdoost A, et al. Prospective epidemiological research studies in Iran (the PERSIAN cohort study): rationale, objectives, and design. Am J Epidemiol. (2018) 187:647–55. doi: 10.1093/aje/kwx314
17. Mirzaei M, Salehi-Abargouei A, Mirzaei M, Mohsenpour M. Cohort profile: the yazd health study (YaHS): a population-based study of adults aged 20–70 years (study design and baseline population data). Int J Epidemiol. (2018) 47:697–8h. doi: 10.1093/ije/dyx231
18. Zimorovat A, Moghtaderi F, Amiri M, Raeisi-Dehkordi H, Mohyadini M, Mohammadi M, et al. Validity and reproducibility of a semiquantitative multiple-choice food frequency questionnaire in Iranian adults. Food Nutr Bull. (2022) 43:171–88. doi: 10.1177/03795721221078353
19. Ghaffarpour M, Houshiar-Rad A, Kianfar H. The manual for household measures, cooking yields factors and edible portion of foods. Tehran Nashre Olume Keshavarzy. (1999) 7:42–58.
20. Haytowitz DB, Ahuja JKC, Wu X, Somanchi M, Nickle M, Nguyen QA, et al. USDA National Nutrient Database For Standard Reference, Legacy Release. Beltsville, MD: Nutrient Data Laboratory, Beltsville Human Nutrition Research Center, ARS, USDA (2019).
21. Moghaddam M, Aghdam F, Jafarabadi M, Allahverdipour H, Nikookheslat S, Safarpour S. The Iranian version of international physical activity questionnaire (IPAQ) in Iran: content and construct validity, factor structure, internal consistency and stability. World Appl Sci J. (2012) 18:1073–80.
22. Pate R, Pratt M, Blair S, Haskell W, Macera C, Bouchard C, et al. Physical activity and public health: a recommendation from the centers for disease control and prevention and the American college of sports medicine. JAMA. (1995) 273:402–7. doi: 10.1001/jama.1995.03520290054029
23. Prasad K, Mishra M. Do advanced glycation end products and its receptor play a role in pathophysiology of hypertension? Int J Angiol. (2017) 26:1–11. doi: 10.1055/s-0037-1598183
24. Uribarri J, Cai W, Pyzik R, Goodman S, Chen X, Zhu L, et al. Suppression of native defense mechanisms, SIRT1 and PPARγ, by dietary glycoxidants precedes disease in adult humans; relevance to lifestyle-engendered chronic diseases. Amino Acids. (2014) 46:301–9. doi: 10.1007/s00726-013-1502-4
25. Du R, Cao H. Dietary pattern associated with C-reactive protein and trajectories of blood pressure in Chinese adults: evidence from the china health and nutrition survey. J Hum Nutr Diet. (2022) 35:605–12. doi: 10.1111/jhn.12978
26. Currenti W, Godos J, Alanazi A, Grosso G, Cincione R, La Vignera S, et al. Dietary fats and cardio-metabolic outcomes in a cohort of italian adults. Nutrients. (2022) 14:4294. doi: 10.3390/nu14204294
27. Key T, Appleby P, Bradbury K, Sweeting M, Wood A, Johansson I, et al. Consumption of meat, fish, dairy products, and eggs and risk of ischemic heart disease. Circulation. (2019) 139:2835–45. doi: 10.1161/CIRCULATIONAHA.118.038813
28. Sakurai M, Stamler J, Miura K, Brown I, Nakagawa H, Elliott P, et al. Relationship of dietary cholesterol to blood pressure: the INTERMAP study. J Hypertens. (2011) 29:222–8.
29. Ohara Y, Peterson T, Harrison D. Hypercholesterolemia increases endothelial superoxide anion production. J Clin Investig. (1993) 91:2546–51. doi: 10.1172/JCI116491
30. Baumann M. Role of advanced glycation end products in hypertension and cardiovascular risk: human studies. J Am Soc Hypertens. (2012) 6:427–35. doi: 10.1016/j.jash.2012.10.005
31. Di Pino A, Currenti W, Urbano F, Scicali R, Piro S, Purrello F, et al. High intake of dietary advanced glycation end-products is associated with increased arterial stiffness and inflammation in subjects with type 2 diabetes. Nutr Metab Cardiovasc Dis. (2017) 27:978–84. doi: 10.1016/j.numecd.2017.06.014
32. Clifton P, Keogh J. Role of dietary advanced glycation end products. Curr Opin Lipidol. (2017) 28:514–5.
33. Koeth R, Wang Z, Levison B, Buffa J, Org E, Sheehy B, et al. Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. (2013) 19:576–85. doi: 10.1038/nm.3145
34. Micha R, Wallace S, Mozaffarian D. Red and processed meat consumption and risk of incident coronary heart disease, stroke, and diabetes mellitus: a systematic review and meta-analysis. Circulation. (2010) 121:2271–83.
35. He F, Li J, MacGregor G. Effect of longer term modest salt reduction on blood pressure: cochrane systematic review and meta-analysis of randomised trials. BMJ. (2013) 346:f1325.
36. Tsubota-Utsugi M, Ohkubo T, Kikuya M, Metoki H, Kurimoto A, Suzuki K, et al. High fruit intake is associated with a lower risk of future hypertension determined by home blood pressure measurement: the OHASAMA study. J Hum Hypertens. (2011) 25:164–71. doi: 10.1038/jhh.2010.48
37. Godos J, Vitale M, Micek A, Ray S, Martini D, Del Rio D, et al. Dietary polyphenol intake, blood pressure, and hypertension: a systematic review and meta-analysis of observational studies. Antioxidants. (2019) 8:152. doi: 10.3390/antiox8060152
38. Harasym J, Oledzki R. Effect of fruit and vegetable antioxidants on total antioxidant capacity of blood plasma. Nutrition. (2014) 30:511–7.
39. Barona J, Aristizabal J, Blesso C, Volek J, Fernandez M. Grape polyphenols reduce blood pressure and increase flow-mediated vasodilation in men with metabolic syndrome. J Nutr. (2012) 142:1626–32. doi: 10.3945/jn.112.162743
40. Schwingshackl L, Schwedhelm C, Hoffmann G, Knüppel S, Iqbal K, Andriolo V, et al. Food groups and risk of hypertension: a systematic review and dose-response meta-analysis of prospective studies. Adv Nutr. (2017) 8:793–803.
41. Sawicki C, Jacques P, Lichtenstein A, Rogers G, Ma J, Saltzman E, et al. Whole-and refined-grain consumption and longitudinal changes in cardiometabolic risk factors in the framingham offspring cohort. J Nutr. (2021) 151:2790–9. doi: 10.1093/jn/nxab177
42. Tang G, Meng X, Li Y, Zhao C, Liu Q, Li H. Effects of vegetables on cardiovascular diseases and related mechanisms. Nutrients. (2017) 9:857.
43. Marques F, Nelson E, Chu P, Horlock D, Fiedler A, Ziemann M, et al. High-fiber diet and acetate supplementation change the gut microbiota and prevent the development of hypertension and heart failure in hypertensive mice. Circulation. (2017) 135:964–77. doi: 10.1161/CIRCULATIONAHA.116.024545
44. Houston M, Harper K. Potassium, magnesium, and calcium: their role in both the cause and treatment of hypertension. J Clin Hypertens. (2008) 10:3–11.
45. Fekete Á, Givens D, Lovegrove J. Casein-derived lactotripeptides reduce systolic and diastolic blood pressure in a meta-analysis of randomised clinical trials. Nutrients. (2015) 7:659–81. doi: 10.3390/nu7010659
46. Xiao-Ge N, Can C, Fang-Chao L, Jian-Xin L, Huang K, Xue-Li Y, et al. Associations of tea consumption with blood pressure progression and hypertension incidence. J Geriatr Cardiol. (2021) 18:645.
47. Maleki S, Crespo J, Cabanillas B. Anti-inflammatory effects of flavonoids. Food Chem. (2019) 299:125124.
48. Shridhar K, Satija A, Dhillon P, Agrawal S, Gupta R, Bowen L, et al. Association between empirically derived dietary patterns with blood lipids, fasting blood glucose and blood pressure in adults-the India migration study. Nutr J. (2018) 17:15. doi: 10.1186/s12937-018-0327-0
Keywords: hypertension, dietary pattern, dietary analysis methods, principal component analysis, partial least-squares, reduced-rank regressions
Citation: Beigrezaei S, Jambarsang S, Khayyatzadeh SS, Mirzaei M, Mehrparvar AH and Salehi-Abargouei A (2023) A comparison of principal component analysis, partial least-squares, and reduced-rank regressions in the identification of dietary patterns associated with hypertension: YaHS-TAMYZ and Shahedieh cohort studies. Front. Nutr. 9:1076723. doi: 10.3389/fnut.2022.1076723
Received: 21 October 2022; Accepted: 21 December 2022;
Published: 12 January 2023.
Edited by:
Yong Zhu, General Mills, United StatesReviewed by:
Maryam Mahdavi, Shahid Beheshti University of Medical Sciences, IranErika Martinez-Lopez, University of Guadalajara, Mexico
Luisa Lampignano, National Institute of Gastroenterology S. De Bellis Research Hospital (IRCCS), Italy
Copyright © 2023 Beigrezaei, Jambarsang, Khayyatzadeh, Mirzaei, Mehrparvar and Salehi-Abargouei. 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: Amin Salehi-Abargouei, 
abargouei@ssu.ac.ir, abargouei@gmail.com
†Present Address: Amin Salehi-Abargouei, Department of Nutrition, School of Public Health, Shahid Sadoughi University of Medical Sciences, Yazd, Iran