Psychological stress is a negative affective condition in which environmental demands tax or exceed one’s adaptive capacity(Reference Cohen and Gordon1). A recent study of 329 adults in Saudi Arabia found that 47·7 % were experiencing stress, and rates were significantly higher among females(Reference El-Gamal, Nadhrah and Abdulghani2). In a larger study (n 1597) focusing on stress in Saudi Arabia during the coronavirus disease 2019 pandemic, 12 % of respondents reported moderate to severe stress levels, and it was significantly higher among females and younger participants(Reference Alamri, Algarni and Shehata3). Similarly, among young adults, particularly those attending university in Saudi Arabia, the prevalence of stress is evident, with results again indicating that female students experience higher stress levels than male students(Reference Abdulghani, AlKanhal and Mahmoud4,Reference Al-Sowygh5) . Moreover, this period between adolescence and adulthood is considered stress arousing in general, given that college students are in a transitional time wherein they are obtaining independence and self-sufficiency(Reference Zirkel and Cantor6). This is concerning given that high levels of stress are linked to other negative affective states, such as depression and anxiety, and can influence biological processes that increase susceptibility to chronic disease(Reference Cohen and Gordon1).
It is thought that nutrition can play an essential role in mental health and that dietary shifts towards less healthy dietary patterns may be detrimental(Reference Freeman7). Diets in Saudi Arabia are transitioning towards higher intakes of saturated fats, cholesterol and refined carbohydrates and lower intakes of polyunsaturated fats and fibre(Reference Musaiger8). This is considered a Western dietary pattern and has been associated with higher levels of negative affective conditions(Reference Jacka, Pasco and Mykletun9). However, investigations into the relationship between plant-based diets, including vegan and vegetarian diets, and the negative affective conditions have resulted in conflicting results. For instance, according to some studies, plant-based diets are associated with decreased levels of negative affective conditions, including stress(Reference Beezhold, Radnitz and Rinne10,Reference Zamani, Daneshzad and Siassi11) , while other studies have not found such an association(Reference Lavallee, Zhang and Michalak12,Reference Michalak, Zhang and Jacobi13) . Two recent systematic reviews and meta-analyses attempted to address these inconsistencies(Reference Iguacel, Huybrechts and Moreno14,Reference Askari, Daneshzad and Darooghegi Mofrad15) . Iguacel and colleagues(Reference Iguacel, Huybrechts and Moreno14) concluded that there were no associations between consuming a vegetarian diet, continuous depression scores or stress. Likewise, Askari and colleagues(Reference Askari, Daneshzad and Darooghegi Mofrad15) found no relationship between a plant-based diet and either depression or anxiety and concluded that the available pooled data were insufficient to confirm the relationship between plant-based diets and psychological stress, highlighting the need for more research(Reference Askari, Daneshzad and Darooghegi Mofrad15). The adherence to plant-based diets in Saudi Arabia has not been thoroughly investigated. Only a couple of small cross-sectional studies related to disordered eating have examined the rates of vegetarianism among young women (between 7 % and 11 % of the study population)(Reference Fatima, Fatima and Anwar16,Reference Fatima and Ahmad17) .
However, not all plant-based diets are of comparable quality. The discrepancies in previous findings may be due to variations in the healthfulness and completeness of these diets in providing necessary nutrients. Only a couple of studies have investigated the relationship between psychological stress and diet while accounting for differences in the healthfulness of plant-based diets(Reference Zamani, Daneshzad and Siassi11,Reference Daneshzad, Keshavarz and Qorbani18) . These studies found that adherence to healthier plant-based diets was associated with reductions in psychological stress(Reference Zamani, Daneshzad and Siassi11), while unhealthy plant-based diets were associated with increases in psychological stress(Reference Daneshzad, Keshavarz and Qorbani18). Both studies used the Depression, Anxiety, and Stress Scale, a tool used to generate scores for general negative affective disorders(Reference Lovibond and Lovibond19). Although depression, anxiety and stress are invariably related, they are distinct diagnoses that warrant individual and focused investigations. Furthermore, the relationship between plant-based diets and mental health has been shown to vary between countries(Reference Lavallee, Zhang and Michalak12), emphasising the importance of regionally specific studies, especially given the noticeable lack of evidence on these topics in Saudi Arabia. The aim of this study was to examine the relationship between the healthfulness of plant-based diets and psychological stress in female students in Saudi Arabia.
Methods
Study design
The current cross-sectional study involved a total of 401 college students aged 19–35 years from the female campus of King Saud University, Riyadh, Saudi Arabia. The details of the study design were published previously(Reference Alfreeh, Abulmeaty and Abudawood20). Briefly, a pilot study was initially carried out on thirty volunteers to standardise the methodology; the pilot study volunteers and their data were not part of the final study sample. Students were invited to take part in the study through the university’s email directory, social media networks and ads. After students provided their consent, data were collected over a period of 4 months, with each participant providing demographic, dietary and anthropometric information at the clinics of the College of Applied Medical Sciences, King Saud University. Figure 1 shows the flow chart of the study. For biochemical analysis, blood samples were collected from the phlebotomy room of the same college. Students were excluded if they had been diagnosed with a chronic or inflammatory disease (e.g., CVD, diabetes mellitus, rheumatoid arthritis) or infection; were taking high dosages (> 300 mg/d) of any nonsteroidal anti-inflammatory medications (e.g., aspirin); were on hormonal contraceptives (e.g., birth control pills) or were pregnant (n 10). Participants had the right to withdraw consent at any time without offering any reason. This study was approved by the local institutional ethics committee of the specified site (reference number 19/0105/IRB).
Dietary assessment
A validated Saudi FFQ with high overall reproducibility(Reference Gosadi, Alatar and Otayf21) was administered by trained interviewers using a standardised protocol to collect data regarding individuals’ food and beverage intake during the prior year. The Goldberg equation(22) was used to define under- and over-reporters of energy. Based on this equation, twenty-seven students were excluded from the analysis (total energy intake was < 650 or > 4200 kcal/d), bringing the total number of participants to 401. To validate the FFQ data, two 24-h dietary recalls were additionally used on 25 % of the participants (n 100). The FFQ comprised 133 food items and also collected data on fast food intake and the use of artificial sweeteners and added fats and salt. The USA Department of Agriculture food composition tables(23) were used to identify the nutritional content of foods. Dietary data were entered using NutriBase version 20 (CyberSoft Inc.).
Plant-based diet indices
Plant-based diet indices (PDIs) have been defined previously(Reference Satija, Bhupathiraju and Spiegelman24). In brief, a PDI assigns positive scores to plant-based foods and negative scores to animal-based foods. The healthy PDI (hPDI) and unhealthy (uPDI) are versions of the overall PDI that are distinguished based on healthfulness. Therefore, for the hPDI, plant-based foods that are established as healthier were assigned positive scores (e.g. whole grains, fruits, vegetables, nuts and legumes, vegetable oils, tea and coffee). Less healthy plant foods were assigned negative scores (e.g. fruit juices, sugary drinks, refined grains, potatoes, sweets and desserts), as were animal-based food sources (e.g. dairy products, eggs, fish or shellfish and total meat) and miscellaneous food sources (e.g. hamburgers). Local foods were also added to each plant index, according to type. For the uPDI, less healthy plant-based food groups were given positive scores, and healthy plant-based and animal-based food sources were assigned negative scores.
The Saudi Dietary Guidelines(25) were used as references for each food group, where intakes equal to or above the guideline recommendations were identified as being in the highest quartile and given a score of 4 (online Supplementary Table S1). Afterwards, intakes below those of the guideline recommendations were divided into tertiles, and scores of 3 to 1 were assigned.
Scores from the plant-based indices were determined by summing the scores of all food groups (the maximum and minimum scores were 54 and 12, respectively). For the three indices, a high score means low intake of animal-based foods.
Blood biochemical test
Blood samples (2 ml) were withdrawn from each participant during the clinic visit. Plasma and serum were isolated via low-speed centrifugation and were frozen (–80°C). High-sensitivity C-reactive protein (hs-CRP) was measured using an enzyme-linked immunosorbent assay kit (Aviva Systems Biology, OKBA00016) according to the manufacturer’s protocol. The hs-CRP assays were calibrated using international standards crafted by the Centers for Disease Control and Prevention. The hs-CRP was measured for 289 randomly selected participants, and duplicate runs for samples and standards were conducted. The cut-off values of hs-CRP were set as > 3·00 mg/l for a high risk of cardiovascular disease; 1·00–3·00 mg/l for a moderate risk and < 1·00 mg/l for a low risk.
Perceived stress scale
The perceived stress scale (PSS-10; Arabic version) was used to assess each participant’s stress level. This is a validated stress-assessment tool with adequate psychometric properties of reliability and validity that can be applied to different populations(Reference Cohen and Gordon1). The PSS-10 is composed of ten items and measures the prevalence of stressful incidences over the previous month on a five-point Likert scale that includes never (0), almost never (1), sometimes (2), fairly often (3) and very often (4). To sum the PSS-10, scores are reversed for select questions (4, 5, 7 and 8) such that 0 = 4, 1 = 3, 2 = 2, 3 = 1 and 4 = 0. Low, moderate and high levels of stress are defined as score tertiles of 4 − 17, 18 − 23 and 24–36, respectively.
Assessment of confounding variables
When the participants visited the clinic, trained interviewers collected data on demographic characteristics including age (continuous variable), marital status (single, married), education level (bachelor’s, master’s and PhD degree), duration of physical activity (min/d) and intensity of physical activity (low, moderate and severe) defined according to the Centers for Disease Control and Prevention and the American College of Sports Medicine guidelines(Reference Pate, Pratt and Blair26), smoking (yes, no), family income per month (categorised according to the Saudi General Authority of Statistics(27), sleeping (h/d), medical history (presence of any medical condition refers to chronic constipation, gastroesophageal reflux disease, prolactinoma, vitamin D deficiency, polycystic ovary syndrome and hypothyroidism), previous use of a weight loss diet (yes, no) and medication use (yes, no).
Weight and height were recorded twice while participants were not wearing shoes or any heavy clothing. Weight was measured using a mechanical beam scale (Detecto) to the nearest 0·1 kg, and height was measured using a stadiometer (Detecto) to the nearest 0·1 cm. The BMI (kg/m2) was calculated using the average of the two measures.
Statistical analysis
SAS version 9.4 (SAS Institute Inc.) was used for the analyses. The sample size was calculated assuming a stress level ranging between 45 and 65 %(Reference Abdulghani, AlKanhal and Mahmoud4,Reference Al-Sowygh5,Reference Al-Daghri, Al-Othman and Albanyan28) , with a margin of missing data set at 5 % and a confidence level of at least 95 %. The normality of each quantitative variable was tested before analysis using the Kolmogorov–Smirnov test. Socio-demographic characteristics were presented by quartiles of PDI, hPDI and uPDI using the generalised linear age-adjusted model procedure in SAS for continuous variables, with type III SS P-value of 0·05 or less considered significant. Frequencies were used for categorical variables. The relation between PSS score, PDI indices, hs-CRP and anthropometrics was identified using Pearson partial correlation.
Multivariate linear regression models were applied to estimated the mean differences in the PSS-10 score with a six point (1 sd) higher score of PDI, hPDI and uPDI in participants. The models were adjusted for potential lifestyle confounders: age, marital status, education level, college, family income, medical condition, previous weight loss, sleeping h/d, duration of physical activity and BMI. Further adjustment for hs-CRP was included. Additional models were used to identify food groups that contribute the most to the relation between the PDI scores and PSS-10 by summing up the foods other than the intended variable to avoid over-adjusting (e.g. a variable summing up food groups other than fruits was added to the model when investigating the influence of fruits in the relation between PDI scores and PSS-10). P values < 0·05 were considered statistically significant.
Results
Characteristics
About (57 %) of participants were classified as moderately stressed and (30 %) as highly stressed with a mean PSS-10 score of 20·12 ± 6·13 (online Supplementary Table S2). Compared with participants with low hPDI scores, participants with high hPDI scores had higher incomes; spent more hours doing physical activity, especially high-intensity physical activity; had more body muscle; had lower BMIs; had lower body fat percentages; had lower hs-CRP (mg/l) levels and lower PSS-10 scores (Table 1).
hPDI, healthy plant-based diet index; hs-CRP, high-sensitivity C-reactive protein; PA, physical activity; PDI, plant-based diet index; PSS-10, perceived stress scale; SR, Saudi riyals; uPDI, unhealthy plant-based diet index.
* Values are mean (95 % CI) or %.
† Medical condition refers to (chronic constipation, gastroesophageal reflux disease, prolactinoma, vitamin D deficiency, polycystic ovary syndrome and hypothyroidism).
‡ hs-CRP was collected from 289 participants.
In contrast, compared with participants with low uPDI scores, participants with high uPDI scores were less physically active; had higher BMIs; had higher hs-CRP (mg/l) levels and had dietary intakes rich in saturated fat and low in dietary fibre.
According to the analysis of the nutrient compositions of the scores, compared with lower hPDI scores, higher hPDI scores were lower in dietary cholesterol and higher in vitamin E, vitamin C, vitamin B12, thiamine, dietary Ca, K and Mg (online Supplementary Table S3).
The partial correlation between the hPDI and uPDI scores was inverse (r = −0·29), while the scores on the overall PDI were positively correlated with hPDI scores (r = 0·28) and inversely correlated with uPDI scores (r = −0·30; Table 2). The PSS-10 scores were inversely correlated with the overall PDI scores (r = −0·13) and hPDI scores (r = −0·10) and were positively correlated with uPDI scores (r = 0·10) (Table 2).
* Correlations were significant between (r = 0·10, −0·10).
Associations between perceived stress scale-10 and overall plant-based diet index and healthy plant-based diet index scores
A regression analysis revealed no association between the PSS-10 score and a six-point higher overall PDI score (model 2; Table 3). However, a six-point higher hPDI score was associated with a 0·16-point lower PSS-10 score (95 % confidence interval, –0·24, –0·08) after controlling for lifestyle factors. These results remained significant after adjusting for BMI. Further adjustment for hs-CRP attenuated the association between hPDI and PSS-10 and was no longer significant (model 3b, Table 3). Additional individual adjustments for healthy food groups attenuated the association between the hPDI and PSS-10 scores for vegetables and fruits (models 3c and 3d; Table 3). The data revealed no evidence of potential effect modification by BMI when using stratified analyses and interaction terms.
hPDI, healthy plant-based diet index; PDI, plant-based diet index; PSS-10, perceived stress scale; sd, standard deviation; uPDI, unhealthy plant-based diet index.
* Values are presented as means (95 % CIs).
† Model 1 is a crude (unadjusted) model. Model 2 is adjusted for age, marital status, education level, college, family income, medical condition, previous weight loss, sleeping h/d and duration of physical activity. Model 3a is model 2 adjusted for BMI. Model 3b is model 2 adjusted for hs-CRP. Model 3c is model 2 adjusted for intakes of fruits and ‘sum of food groups other than fruits’. Model 3d is model 2 adjusted for intake of vegetables and ‘sum of food groups other than vegetables’.
Association between perceived stress scale-10 score and ulhealthy plant-based diet index score
A six-point higher uPDI score was not associated with the PSS-10 score. The results were similar with adjustments for individual unhealthy food groups (Table 3).
Discussion
To date, studies focused on the associations between plant-based diets and negative affective disorders have found mixed results(Reference Iguacel, Huybrechts and Moreno14,Reference Askari, Daneshzad and Darooghegi Mofrad15) . In this context, this study fills in an important and noted gap between the relationship between stress and plant-based diets(Reference Askari, Daneshzad and Darooghegi Mofrad15). This study’s key findings were an inverse relationship between hPDI and stress and lower hs-CRP with higher adherence to an hPDI. The food groups that contributed most to the inverse association between stress and hPDI were fruits and vegetables.
This study found that higher adherence to an hPDI was associated with lower levels of stress. This finding is similar to other studies investigating the healthfulness of plant-based diets and negative affective disorders(Reference Zamani, Daneshzad and Siassi11,Reference Daneshzad, Keshavarz and Qorbani18,Reference Baden, Kino and Liu29) . Using the Depression, Anxiety, and Stress Scale questionnaire in a sample of Iranian women, Zamani and colleagues(Reference Zamani, Daneshzad and Siassi11) recently found inverse associations between higher overall PDI and hPDI scores and depression, anxiety and psychological distress and a positive association between uPDI scores and depression. Similarly, also using the Depression, Anxiety, and Stress Scale questionnaire in a sample of women with diabetes, Daneshzad and colleagues(Reference Daneshzad, Keshavarz and Qorbani18) found that those with the highest adherence to a uPDI had greater risks of depression, anxiety and stress. Likewise, using mental component scores, Baden and colleagues(Reference Baden, Kino and Liu29) found that higher overall PDI and hPDI scores were associated with mental health improvements. Results of the present study are also similar to those from investigations of plant-based diets more generally, which have demonstrated the benefits of vegan and vegetarian diets on symptoms of stress, depression and/or anxiety(Reference Beezhold, Radnitz and Rinne10,Reference Beezhold and Johnston30) . In addition, findings from other studies have indicated that the relationship between plant-based diets and stress may be moderated by gender(Reference Beezhold, Radnitz and Rinne10), but whether this stems from biology or society remains unclear(Reference Davis, Matthews and Twamley31). Unlike most previous literature investigating the relationship between stress and the healthfulness of plant-based diets, this study utilised the PSS-10 questionnaire, which is explicitly designed to measure stress. PSS-10 scores have been correlated with biomarkers of stress(Reference Van Eck and Nicolson32–Reference Beezhold, Johnston and Daigle34), and the use of the PSS-10 has often been evaluated among college students(Reference Lee35).
Other research has demonstrated rates of stress similar to those observed in the present study. Studies conducted in Saudi Arabia suggest that 47·7 % of participating adults experienced stress, with 12 % experiencing moderate to severe stress(Reference El-Gamal, Nadhrah and Abdulghani2,Reference Alamri, Algarni and Shehata3) . Other studies of university students in Saudi Arabia have found that rates of stress among females were as high as 48·6 % to 75·7 % among medical students and 38·7 % among non-medical students(Reference Abdulghani, AlKanhal and Mahmoud4,Reference Al-Sowygh5) . Research investigating the prevalence of stress among university students in other countries has found similar results. For instance, 59·6 % and 51·9 % of the female medical and non-medical university students participating in studies in India and Turkey, respectively, experienced stress(Reference Iqbal, Gupta and Venkatarao36,Reference Bayram and Bilgel37) .
Hs-CRP is considered an independent risk factor for CVD and is an important health measure to consider(Reference Lagrand, Visser and Hermens38), and associations between hPDI and stress attenuated when adjusted for hs-CRP. It is thought that psychosocial factors, such as psychological stress, can influence CVD processes, including inflammation(Reference Wirtz and von Känel39). The present study found that hs-CRP levels were lower in participants with higher adherence to an hPDI and lower adherence to a uPDI. These findings concur with that of Baden and colleagues, who found hPDI scores were inversely related to hs-CRP in cross-sectional and longitudinal analyses, and uPDI was positively associated with changes in hs-CRP in longitudinal analyses(Reference Baden, Satija and Hu40).
Plant-based foods considered healthy in the PDI score include whole grains, fruits, vegetables, nuts, legumes, vegetable oils, tea and coffee. The present study found that fruits and vegetables positively contributed to the inverse association between stress and hPDI scores. This finding is similar to results from other investigations into plant-based diets and stress(Reference Beezhold, Radnitz and Rinne10) and concurs with studies demonstrating that fruits and vegetables are associated with lower odds of or improved psychological distress(Reference Richard, Rohrmann and Vandeleur41–Reference McMartin, Jacka and Colman43). Several mechanisms are thought to mediate the relationship between fruit and vegetable intakes and psychological well-being(Reference Rooney, McKinley and Woodside44), including a slow release of glucose, which results in the stimulation of serotonin, and the provision of B6, folate and certain minerals, which have been shown to impact neurotransmitter synthesis and improve brain function and processing(Reference Rooney, McKinley and Woodside44,Reference Parletta, Milte and Meyer45) . Additionally, plant-based foods, including fruits and vegetables, are an important source of antioxidants, which are posited to reduce oxidative stress and inflammation, that can be detrimental to psychological well-being(Reference Rooney, McKinley and Woodside44,Reference Parletta, Milte and Meyer45) .
Nutrients such as vitamin B12 (Reference Coppen and Bolander-Gouaille46) and vitamin D(Reference Anglin, Samaan and Walter47) have also been studied regarding psychological well-being, with low levels and deficiencies being associated with depression. However, animal products are a major source of these nutrients, with vitamin B12 found in fish, meats, eggs and milk(48) and vitamin D found in fish, meats, eggs and fortified foods(49). Thus, there are some concerns that consuming an entirely plant-based diet might lead to nutrient deficiencies of dietary components associated with mental health(Reference McEvoy, Temple and Woodside50,Reference Craig and Mangels51) . However, when examined in terms of stress levels, the present study’s findings do not suggest deficiencies in those participants following a plant-based diet. The nutrient compositions of the scores revealed that a higher hPDI was higher in vitamin B12 and a higher uPDI was lower in vitamin D. Furthermore, the nature of the PDI scoring does not require complete abstention from animal products, although a higher PDI score does result from a higher proportion of plant-based sources in the diet.
Meanwhile, those with higher uPDI scores had higher intakes of sugary drinks and sweets. Similarly, sweet intake in other studies was positively related to stress and anxiety in females(Reference Mohamed, Mahfouz and Badr52–Reference Yannakoulia, Panagiotakos and Pitsavos54). In fact, according to a recent cross-sectional study in Saudi Arabia, the perceived stress among females led to an increased preference for sweets(Reference Mohamed, Mahfouz and Badr52). Furthermore, Western dietary patterns, characterised in part by the consumption of high-energy, low-nutrient dense foods, such as sweets and sugary drinks, have been associated with higher levels of negative affective conditions(Reference Musaiger8).
The current study has several strengths. It used assessment tools, such as the FFQ and PSS-10, that are validated, interview-based, and follow thorough and rigorous protocols. Measurement of hs-CRP provided an objective standard for evaluating inflammation. As per WHO guidelines, all measurements were recorded twice and the average was used for analyses; however, residual confounding due to measurement error is possible. In addition, the temporal nature of the association between plant-based diets and stress cannot be determined with a cross-sectional study. It could be that stress leads to the consumption of less healthy diets or that healthier diets lead to reduced stress. Moreover, misreporting of energy is unavoidable, although efforts were made to ensure accurate dietary data, including through the use of the Goldberg cut-off values and validation with 24-h recalls of 20 % of the total FFQ. On the other hand, the use of hs-CRP as a marker of inflammation has its limitations. For instance, hs-CRP measures were taken only once, which does not account for variations in measurements over time(Reference Koc, Karaarslan and Abali55,Reference Bower, Lazo and Juraschek56)
In conclusion, this study found that healthy plant-based diets are associated with lower levels of psychological stress in young Saudi Arabian students. This finding underscores the importance of following diets that are not only plant-based but are also healthy, especially among female students, given that the relationship with stress is thought to be mediated by gender.
Acknowledgements
The authors extend their appreciation to the ‘Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project no (IFKSURG-1441-545)’.
G. S. A.’s roles were conceptualisation, data curation, formal analysis and investigation, writing and reviewing the final version.
I declare no conflict of interest.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114521001689