Aya Jamal Abusheikha1
Corbin S. C Johnson2
Noah Snyder-Mackler3,4,5
Kip D. Zimmerman6,7
Jacob D. Negrey1,8
Kenneth L. Chiou3,4
Brett M. Frye1,9
Timothy D. Howard6,10
Carol A. Shively1*
Thomas C. Register1,6*- 1Department of Pathology/Comparative Medicine, Wake Forest University School of Medicine, Winston-Salem, NC, United States
- 2Department of Neurology, University of Washington, Seattle, WA, United States
- 3Center for Evolution and Medicine, Arizona State University, Tempe, AZ, United States
- 4School of Life Sciences, Arizona State University, Tempe, AZ, United States
- 5School of Human Evolution and Social Change, Arizona State University, Tempe, AZ, United States
- 6Center for Precision Medicine, Wake Forest University Health Sciences, Winston-Salem, NC, United States
- 7Department of Molecular Medicine/Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC, United States
- 8School of Anthropology, University of Arizona, Tucson, AZ, United States
- 9Department of Biology, Emory and Henry College, Emory, VA, United States
- 10Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, NC, United States
Western diets and social subordination are associated with increased risk of cardiovascular disease. In this study, we investigated the impact of Western versus Mediterranean diets and social status on atherogenesis and arterial transcriptional profiles in a 30-month randomized study in middle-aged, cynomolgus monkeys (Macaca fascicularis). Atherosclerosis (intimal area) in the left anterior descending (LAD) coronary artery was higher in the Western diet group compared to the Mediterranean diet group (F = 5.25, p = 0.03). There was no effect of diet on intimal lesion size in the iliac and carotid arteries (p > 0.05). Diet altered the transcriptome in iliac arteries; at an FDR threshold of 0.05, seven transcripts were upregulated (WDR62, PKDCC, SLC29A2, MARS1, RAD21L1, MAMDC4, and ENSMFAG00000052859), and 13 transcripts were downregulated (PIK3R1, PABPC1, PAQR8, ZNF667, FGGY, EIF4B, ALDH3A2, ANP32A, KDM3B, XPO7, RPS20, TOMM20, and CHCHD7) in the Western compared to the Mediterranean diet cohort. These genes are associated with endothelial dysfunction, smooth muscle proliferation and migration, angiogenesis, and abnormal extracellular matrix (ECM) dynamics. In addition, two transcripts (ENSMFAG00000064154 [LncRNA] and ENSMFAG00000057515 [small nucleolar RNA U13]) were downregulated in subordinate monkeys relative to their dominant counterparts (FDR < 0.05). There was no effect of diet on the carotid artery transcriptome, but we did identify significant social status effects: Eleven transcripts were upregulated (KCNQ4, STIM1, TNKS1BP1, CSNK1D, INPPL1, PNPLA7, F10, RAD9A, KCNIP3, ENSMFAG00000059809 [LncRNA], and ENSMFAG00000053865 [secreted protein A0A7N9CS45]), and seven transcripts were downregulated (IRAK1BP1, KIAA0513, SMIM15, PSMD14, TOPORS, ARPC2, and ENSMFAG00000050714 [LncRNA]) in subordinate relative to dominant monkeys. These alterations were associated with dysregulated vascular tone and smooth muscle contractility, apoptosis, and abnormal ECM dynamics. These findings demonstrate differential effects of diet composition and social status depending on arterial sites. The effects of Western diet were observed primarily in the coronary and iliac arteries, whereas social status differences were observed primarily in the carotid arteries. Our results demonstrate that Western diets and social subordination have adverse, yet distinct and tissue-specific impacts on arterial atherogenesis and transcriptional profiles, highlighting the interplay between diet, social hierarchy, and vascular health.
1 Introduction
Dietary patterns play a pivotal role in human health, influencing various physiological processes and disease outcomes. Among dietary regimens, the Western and Mediterranean diets represent two contrasting paradigms with distinct compositions and health implications.
Western diets derive much of their protein and fat content from animal sources and are generally high in saturated fats, sodium, omega-6 (n-6) fatty acids, and processed simple sugars (1). Western diets have been shown to promote sympathetic nervous system arousal, insulin resistance, oxidative stress, and inflammation (2–4). Inflammation and oxidative stress contribute to the pathogenesis of atherosclerosis and poor cardiovascular health (5, 6). Western diets are also associated with increased risk of cardiometabolic diseases including obesity, hypertension, type 2 diabetes and coronary artery atherosclerosis (7–9). In contrast, Mediterranean diets derive most of their protein and fat from plant and fish sources, have higher proportions of monounsaturated, polyunsaturated omega-3 (n-3) fatty acids, and natural antioxidants; and have lower levels of refined simple sugars. Fruits and vegetables, which make a larger contribution to Mediterranean diets, are rich in antioxidants that help reduce oxidative stress and inflammation (6, 10, 11). Moreover, n-3 fatty acids (found in fish, nuts, and flaxseed) and polyphenols in olive oil possess anti-inflammatory and antioxidant properties (3, 11). The Mediterranean diet’s anti-inflammatory properties, attributed to its rich array of phytonutrients and n-3 fatty acids, may help protect against chronic diseases including cardiovascular disease (CVD) (12–15).
Diet may also adversely impact health through activation of stress mechanisms. In population studies, the Western diet pattern has been associated with greater perceived stress and higher sympathetic activity (8) and increased urinary cortisol levels (16); whereas the Mediterranean diet pattern has been associated with lower perceived stress (4, 17, 18) and cortisol levels (19). Psychological stress also promotes inflammation (20), and increases CVD risk (21–23). This raises the possibility that one of the pathways through which a healthy diet reduces CVD risk is by reducing psychological stress responses. This is especially pertinent in societies organized by social status hierarchies, as lower socioeconomic status (SES) has been associated with behavioral and physiological characteristics of stress (24, 25), as well as increased risk for atherosclerosis and CVD (26, 27).
Social determinants of health (SDHs), like SES and ethnicity, exhibit a significant impact on vascular diseases; in which SDHs contribute to disparities in disease prevalence, treatment accessibility, and health outcomes (28). Low SES correlates with higher disease prevalence and poorer postoperative outcomes, partly driven by limited healthcare access and lifestyle factors associated with lower-income settings (28). For example, individuals with lower SES face higher risks of peripheral artery disease and chronic venous disease, often attributed to occupational hazards prevalent in low-income employment sectors, that involves prolonged standing (29–31). Similarly, low SES populations in rural areas have limited access to specialized healthcare services, which further exacerbate health inequalities (32). Low SES is significantly associated with limited access to health care, as individuals with low incomes are more likely to be uninsured, experience poor-quality health care, and primarily seek medical attention during emergencies. Additionally, low-SES individuals often face limitations in care due to underfunded health facilities, contributing to poorer health outcomes (33).
SES significantly influences diet quality as well, in which individuals from low-SES groups typically consume poorer diet quality, characterized by lower amounts of fruits, vegetables, whole grains, fiber, and fish, and higher amounts of energy-dense, nutrient-poor processed foods (34, 35). These disparities are driven by economic limitations and reduced access to healthier food options, which are often more expensive and less available in low-income neighborhoods (36). Diet quality, assessed through the Dietary Quality Score (DQS), Nutrient Adequacy Ratio (NAR), and Dietary Diversity Score (DDS) were shown to consistently improve with increasing SES level (34). Another study showed that individuals of higher SES are significantly associated with greater food expenditure; thus, purchasing healthier food options (37).
On a physiological level, many studies explained the intricate relationships between biomarkers and SDHs in influencing the progression and outcomes of vascular diseases. Lower SES was shown to be significantly correlated with elevated levels of inflammatory biomarkers such as high-sensitivity C-reactive protein (hs-CRP) and interleukin-6 (IL-6) (38). Additionally, environmental factors like air pollution, prevalent in low-income areas, have been associated with increased matrix metalloproteinase levels, accelerating vascular damage and disease progression (39, 40). Moreover, social stressors linked to low SES, including limited education and healthcare access, exacerbate biological stress responses, further influencing disease severity through elevated cortisol and inflammatory markers (41). It is important to mention that it is difficult to separate the effects of low SES-related stress from the other factors associated with low SES in people, like air pollution, poor quality diet, and poor health care, that impact health. Therefore, a randomized preclinical trial eliminates most of these confounders.
Despite the substantial evidence supporting the health benefits of the Mediterranean diet, Western dietary patterns continue to prevail in many regions, contributing to the global burden of complex chronic diseases. Moreover, diet composition and its effect on physiological mechanisms have been mostly investigated in short term controlled trials which may not predict long term health effects (42). Long-term observational studies are confounded by other factors and behaviors (e.g., smoking behavior, alcohol consumption, physical activity levels), and rely on self-reported food intake (42–44). Long-term randomized feeding trials are thus required to better predict long term health effects. These relationships and uncertainties led us to design an experiment to address diet and behavior interactions in a rigidly controlled study in a well-established nonhuman primates (NHP) model of diet and psychosocial stress effects on health.
NHPs are useful models to study human health and complex diseases as they share many genetic, physiological, and behavioral phenotypes with humans, including key characteristics of energy metabolism, susceptibility to develop aging related diseases including metabolic and vascular disorders when consuming Western diets (45–49). Additionally, similar to humans, individual NHPs differ in their degree of social integration and isolation, in which multiple factors can affect sociality, including sex, age, social status, and kin networks (50, 51). NHPs can also participate reliably in long term trials with extensive phenotyping in controlled studies. In both wild and captive groups, female macaques exhibit social status hierarchies (1). Like humans, social hierarchies in macaques are associated with social inequalities in health. Similarities between NHPs and humans suggest that NHPs may share mechanisms linking social status and health with humans. The ability to directly record social behaviors in NHPs avoids the confounders found in human studies, thus supporting their use as a model system.
Subordinate monkeys (defined here as those in the lower half of the hierarchy) spend more time alone, receive more aggression, are groomed less, and are more fearful and anxious than their dominant counterparts (1, 47, 50, 52–54). Subordinates respond to a standardized stressor with a higher initial heart rate and recover more slowly than dominants (47, 50, 55). Subordinates also have higher visceral obesity, more inflammation, lower bone density (56, 57), and develop more diet-induced coronary and carotid artery atherosclerosis than their dominant counterparts (58, 59). With respect to social status and health, it is important to note that most of the available data from epidemiological studies and macaques has been derived from subjects consuming a Western diet (60, 61).
Here we report the results of our 30-month randomized trial, in which diet was manipulated, environment and confounders were controlled, and multiple cardiovascular tissues were collected at the end of the study. The purpose of the present study was to determine the effects of diet (Western versus Mediterranean) and social status on atherosclerosis in the coronary, iliac, and carotid arteries and to identify putative molecular mechanisms underlying the effects. We aimed to determine the transcriptome-related effects of diet and social status in the iliac arteries, which are a good surrogate for coronary arteries in macaques (49, 62, 63), and the carotid arteries, which are important in blood pressure control and associated with stroke risk (64). Our overall hypotheses were that the Mediterranean diet would promote overall vascular health and potentially protect against adverse effects of subordinate status.
2 Materials and methods
2.1 Subjects
Thirty-eight middle-aged (mean = 9.0, range = 8.2–10.4 years, estimated by dentition) female cynomolgus monkeys (Macaca fascicularis) were obtained from Shin Nippon Biomedical Laboratories US SRC (Alice, TX) and housed at the Wake Forest University School of Medicine (Winston-Salem, NC) (42). Based on previously published comparative aging models, a 9-year-old Macaca fascicularis is approximately equivalent to a 35–40-year-old human female in terms of physiological and reproductive aging (e.g., age at menarche, menopause, and onset of age-related pathologies) (65). Monkeys were quarantined in single cages for 1 month, then socially housed in groups of 3–4 animals in indoor enclosures (3 m x 3 m x 3 m), maintained on a 12/12 light/dark cycle with exposure to daylight, and water ad libitum. All animal manipulations were performed according to the guidelines of state and federal laws, the US Department of Health and Human Services, and the Animal Care and Use Committee of Wake Forest University School of Medicine.
2.2 Design
The study design was a preclinical randomized trial 39 months in length. During an 8-month baseline phase, the animals consumed standard monkey chow (Supplementary Table 1; monkey diet 5037/5038; LabDiet, St. Louis, Missouri) with ad libitum access to water. Monkeys were allowed to adapt to their social groups for 3 months followed by assessments of behavior, physiology, and social status during the last 4 months of the baseline phase. Individual social groups were then assigned to receive a Western diet (Western n = 21) or a Mediterranean diet (Mediterranean n = 17) by stratified randomization to balance body mass index, basal plasma cortisol and triglyceride concentrations (42). Comprehensive physical, behavioral, and physiological characteristics were assessed during the diet intervention phase.
2.3 Experimental diets
Western and Mediterranean diets were formulated to be matched on proportions of calories from macronutrients (protein, carbs, and fats), cholesterol content (~ 320 mg / 2000 calories/day), and caloric density (2.8 cal/g wet diet; Supplementary Table 1). Diet composition was designed and prepared at the Wake Forest School of Medicine Comparative Medicine Diet Lab. Western diet was designed to be similar to the nutritional composition of the typical North American diet, with protein and fat derived mainly from animal sources, the diet was high in saturated fats and low in mono- and poly-unsaturated (66). The Mediterranean diet derived most of the protein and fat from plant sources, and was relatively high in mono- and poly-unsaturated fatty acids (67). Further, the n-3: n-6 fatty acid ratio was approximately 5-fold higher in the Mediterranean diet compared to the Western diet (68).
2.4 Food consumption
Methods for monitoring food consumption have been previously described (42). Briefly, custom individual feeding cages were fabricated and placed inside group pens. Monkeys were trained to run into the feeding cages on voice command, consume their food, and then be released back into their group pen. Each monkey was offered 120 kcal diet/kg BW/day. Throughout the experiment, the animals finished feeding in the first 30 min. During the baseline phase and the first year of the treatment phase, the monkeys were fed for 2 h once a day at ~0730 h. During the last half of the treatment phase, beginning in month 16, one-half of the ration was fed at ~0730 h for 1 h, and the other half was fed at ~1,330 h for 1 h. Food consumption was estimated by weight before and after the feeding periods during month 5 of the baseline phase and months 3–6, 7–10, 17–20, and 27–30 of the experimental phase; average consumption was calculated for each monkey for each time period.
2.5 Quantification of social status
Social status rankings were determined throughout the experiment by recording the outcomes of agonistic interactions within monkeys of the same social group (47, 50, 69, 70). The highest-ranking monkey of each group was the female that defeated all others, as demonstrated by her ability to consistently elicit submissive responses. The female that defeated all but the first-ranking monkey was designated as the second-ranking monkey, etc. Monkeys formed relatively stable linear hierarchies, and social status relationships within groups were typically transitive, meaning that if monkey 1 was dominant over monkey 2, and monkey 2 was dominant over monkey 3; then monkey 1 was dominant over monkey 3 as well. We determined social status by defining the number of monkeys submissive to a given monkey, corrected for group size so the resulting relative rank could be compared across social groups (71). Accordingly, relative rankings varied from 0 to 1. For some analyses, we dichotomized social status, in which females with a relative rank greater than 0.5 were considered dominant and labeled as 1, and all others were considered subordinate (50, 72). Social subordination appears stressful, as subordinates in this study received more aggression, were groomed less, spent more time fearfully scanning, had higher systolic blood pressure, and had higher sympathetic nervous system activity than did dominants (73). Thus, we used social status as a proxy for psychosocial stress.
2.6 Behavioral characterization
Agonistic and affiliative behavior were recorded in 10 min focal animal observations (74) beginning in the third month of the experimental phase as previously described (73). Behavioral data were collected weekly during two 10 min focal observations, balanced for time of day, for 6 weeks during the baseline phase (2 h/monkey total) and for 14 months during the experimental phase (mean = 17.6 h/monkey total). Behaviors recorded included affiliative behaviors: percent of time spent alone (out of monkey’s arm’s reach), percent of time spent in body contact with other monkeys, and percent of time spent in close proximity (within monkey’s arm’s reach) (67, 75–77). Anxiety behavior observed is the frequency of scratching, itching, yawning, and self-grooming (67, 75–77). Percent of time spent depressed is defined by recording the time a monkey spends in a depressive behavior (slumped or collapsed body posture accompanied by a lack of responsivity to environmental events with eyes open) (78–81).
2.7 Arterial tissue collection and processing
At the end of the experimental phase, monkeys were anesthetized with pentobarbital (30–50 mg/kg) to obtain a surgical level of anesthesia, and we collected coronary, iliac and carotid arteries. The heart and coronary arteries were perfusion-fixed with 4% paraformaldehyde at 100 mm Hg for 1 h, followed by storage in 70% ethanol at 4°C. Twelve blocks (each 3 mm in length) were cut perpendicular to the long axis of the left circumflex (LCX), left anterior descending (LAD), and right coronary artery (RCA). Iliac and carotid arteries were subdivided perpendicular to long axis into 5 sections, sections 1,3, and 5 were fixed for assessment of lesion size, sections 2 and 4 were rinsed in ice-cold lactated Ringer’s solution and placed in RNAlater® for removal of adventitia and adherent connective tissue. Intima-media sections for transcriptomic analyses were then slam frozen and stored in a − 80°C freezer until processing.
2.8 Atherosclerosis measures
Blocks of coronary, carotid, and iliac arteries were embedded in paraffin, and 5 μm sections were stained with Verhoeff-van Gieson’s stain. Images were captured for analysis using ImageProPlus software, and cross-sectional area of the arterial intima was measured as described previously (63).
2.9 RNA extraction and processing
RNA was isolated from iliac and carotid arteries samples using the QIAGEN Quick-Start RNeasy Mini Kit (QIAGEN, Hilden, Germany). RNA integrity number (RIN) scores were determined to assess quality with an average score of 8 across samples. One iliac artery RNA sample with a RIN score of only 4.4 was excluded from all analyses, leaving a final sample size of 35 iliac arteries (n = 20 Western diet, n = 15 Mediterranean diet). We processed RNA for bulk-RNA sequencing by mapping 37,781 reads to the Macaca fascicularis reference genome (Macaca_fascicularis_6.0) using STAR-2.5.2a (82), then converted to a sample-by-gene read count matrix using featureCounts (83).
2.10 Quality control
Read counts data were filtered using the rpkm function from the edgeR R package (84), excluding genes with median reads per kilobase per million reads mapped (rpkm) of less than one, which resulted in 13,020 and 13,206 genes for downstream analyses in iliac and carotid arteries, respectively. Read counts were then normalized using the weighted trimmed mean of M-value (TMM) method (85) followed by log2 transformation using the voom function of the R package limma (86), which also applies normalization scaling factors that take into account the differences in sequencing depth. To control for technical variables related to RNA quality and reads mapping, residual gene expression was determined from a model of normalized gene expression as a function of the technical variables that significantly contribute to variation in gene expression. Principal component analysis (PCA) (87) was performed using the prcomp function in R to determine effects of technical variables. For the iliac artery, we controlled for library size, percent of uniquely mapped reads, and RNA concentration (ng/ul); for the carotid artery, we controlled for percent of uniquely mapped reads, RNA concentration (ng/ul), and RNA integrity. The resulting residual gene expression values (called the residual matrix) were used in the subsequent analyses. Hierarchal clustering analysis was carried out using the hclust function from the WGCNA R package (88) to identify similar groups and potential sample outliers.
2.11 Modeling differential gene expression by diet composition and social status
Differential expression analysis was employed to identify differences in transcript levels in response to diet, social status and potential interactions. The R package EMMREML was used to model the residual expression of each gene as a function of diet and/or social status with a linear mixed effects model controlling for relatedness among monkeys (89). We estimated genetic relatedness using the ngsRelate program (90) with SNP genotype inferred from the RNA-seq reads using bcftools mpileup (91). Social status contributions to outcomes were assessed using both a continuous (relative rank score) or dichotomous (dominant rank vs. subordinate rank) approach. We used the qvalue R package (92) to calculate the false discovery rates (FDR) for each gene, and reported genes that passed an FDR threshold of 0.05 as significant. To examine global patterns of variations in gene expression, principal component analysis was applied to the correlation matrix of normalized residual gene expression using the prcomp function in R.
2.12 Enrichment analysis and identifying biological pathways
Enrichment of biological processes and regulatory networks exhibiting diet-specific and social status-specific activation patterns were assessed using Pathway Enrichment Analyses. Pathway analysis was performed with QIAGEN’s Ingenuity Pathway Analysis (IPA, QIAGEN Redwood City).1 Given the limited number of differentially expressed transcripts (DETs) identified at FDR threshold of 0.05, we expanded our criteria to include genes with an FDR threshold of 0.2 (n = 174 for iliac arteries, n = 216 for carotid arteries).
2.13 Correlations between gene expression and phenotypic variables
During the last half of the experimental phase (13–31 months on experimental diet), study subjects underwent an extensive characterization of different phenotypes that represent categories of anatomy, physiology, and behavior. Detailed definitions and references describing these phenotypes are presented in Supplementary Table 6. Briefly, phenotype categories analyzed in this study were: stress physiology (hypothalamic–pituitary–adrenal function) (73, 79), 24 h heart rate via telemetry and heart rate variability (HRV) (73), blood pressure (73), morphometrics (42), insulin resistance (42), ovarian functioning (93), and coronary artery atherosclerosis (81). We measured cortisol levels for acute and ACTH challenge tests using RIA kits from DiaSource (IBL America) (73, 79). Electrocardiograms (ECG) were recorded over a 24 h period at 12 and 29 months of experimental diet consumption using a portable ECG telemetry unit (Life Sensing Instrument Co, Tullahoma, TN, United States) (73, 79, 81, 94). Dependent variables derived from the ECG included HRV in 2-h daytime (1600–1800) and nighttime (0100–0300) blocks derived using Nevrokard-HRV software (Nevrokard Kiauta, d. o. o. Izola Slovenia) (73, 95, 96). Blood pressure was measured using high definition oscillometry via tail cuff (73). Ovarian function was assessed during months 16–27 of experimental phase by measuring peak luteal-phase progesterone concentration (ng/ml), cycle length (days), and cycle regularity for each menstrual cycle. Blood sampling was performed 3X a week or progesterone assay while vaginal swabbing was conducted 6X a week for menses (93, 97). Progesterone was determined using a commercially available radioimmunoassay kit (DiaSource Progesterone RIA-CT KIP1458, Louvain-la-Neuve, Belgium). Plasma estradiol was assessed at the time of necropsy using a sensitive estradiol ELISA kit (Cat # IB78239) from IBL America (Minneapolis, MN, United States). Detailed definition and references describing the collected phenotypes are presented in Supplementary Table 6. Relationships between the phenotypes and the normalized transcript levels of the top differentially expressed genes (FDR < 0.05) were analyzed in iliac arteries as a function of diet, and in carotid arteries as a function of social status (relative rank score). To measure correlation coefficient and significance, we conducted a Spearman’s correlation test and adjusted for diet in iliac arteries, and for relative rank in carotid arteries, using the corr.test function in R.
3 Results
3.1 Western diet promoted intimal lesion formation in the left anterior descending coronary arteries
Atherosclerosis was assessed in the left anterior descending artery (LAD), left circumflex artery (LCX), right coronary artery (RCA), in the common iliac artery (LCI), and left carotid artery (LCC) as previously described (49). Intimal lesion areas were assessed for the main effects of diet (Western vs. Mediterranean) and social status (dominant vs. subordinate). We also tested the interaction between diet and social status. Western diet fed monkeys had significantly larger intimal lesions in the LAD coronary artery (0.04 ± 0.03 mm2; n = 21) compared to monkeys fed the Mediterranean diet (0.02 ± 0.008 mm2; n = 17; F (1, 36) = 5.25, p = 0.03; Table 1; Figure 1). There were no significant main effects of diet, social status, or diet by status interactions on lesions in any other arteries.
Table 1. Effects of diet and social status on atherosclerosis.
Figure 1. Western diet promoted intimal lesion formation in the proximal left anterior descending coronary compared to the Mediterranean diet. Tissue sections were stained with Verhoeff-van Gieson’s stain, and cross-sectional area of the arterial intima was measured. Two-way ANOVA showed a significant main effect of diet (Western 0.04 ± 0.03 mm2 vs. Mediterranean 0.02 ± 0.008 mm2; F(1, 36) = 5.25, p = 0.03).
3.2 Overview of diet and social status effects on the transcriptional profiles in iliac and carotid arteries
Diet was associated with significant differences in transcript levels in iliac but not carotid arteries; while social status was associated with significant differences in transcript levels mainly in the carotid artery, and to a smaller extent in the iliac artery. The social status effect on transcript levels was assessed using both a continuous (relative rank score) and dichotomous (dominant rank vs. subordinate rank) approach; effect size on transcript levels was larger when using the continuous value, thus we are only reporting differential expression in response to the relative rank score in this paper. There was no significant interaction between diet and social status in either artery. A summary of the findings is shown in Table 2.
Table 2. Summary of diet and social status differential effects on iliac and carotid arterial transcriptome.
3.2.1 Effects of diet and social status on the iliac artery transcriptome
For the iliac artery, 20 transcripts were differentially expressed between the two diets (FDR < 0.05; Figure 2). In the Western diet cohort, seven genes were upregulated (WDR62, PKDCC, SLC29A2, MARS1, RAD21L1, MAMDC4, and ENSMFAG00000052859 [uncharacterized coding protein]), and 13 genes were downregulated (PIK3R1, PABPC1, PAQR8, ZNF667, FGGY, EIF4B, ALDH3A2, ANP32A, KDM3B, XPO7, RPS20, TOMM20, and CHCHD7) relative to the Mediterranean diet cohort (Figure 2; Supplementary Table 2).
Figure 2. Transcript levels (TMM normalized) differentially expressed by diet in the iliac artery (FDR < 0.05).
Using a relaxed FDR threshold of 0.2 to identify a larger set of genes for pathway analysis, we identified 174 diet-associated DETs (Figure 3A; Supplementary Table 2). Regardless of the effect size direction (increased or decreased in Western diet compared to Mediterranean diet), these DETs were enriched for a number of pathways including (1) Translation Initiation, (2) mTOR Signaling, (3) Signaling by PDGF, (4) Autophagy, (5) Thrombin Signaling, (6) Insulin Secretion Signaling Pathway, (7) Integrin Signaling, (8) Neutrophil Extracellular Trap Signaling Pathway (Fisher’s Exact Test, p < 0.05; Figure 3B). A complete list of the enriched biological pathways is available in Supplementary Table 3. We also explored pathway enrichment taking into account the effect size direction (increased or decreased in Western diet compared to Mediterranean diet) of the analyzed DETs. We identified and predicted the activity of cardiovascular disease- associated causal networks using IPA, including angiogenesis and autophagy. Autophagy was predicted to be inhibited in the Western diet cohort relative to the Mediterranean cohort (Z-score = −0.572, FDR = 0.02) according to the effect size of the following genes: FKBP5, MAP1LC3A, MSH2, NDRG1, PIK3R1, PLK1, PRODH, RAC3, SREBF2, and CAMK1 (Figure 3C). Angiogenesis was also predicted to be inhibited in the Western cohort relative to the Mediterranean diet cohort (Z-score = −1.732, FDR = 0.05) according to the effect size of the following genes: PALD1, PLK1, RIPK3, SYNE2, AMIGO2, BMX, HMGB2, NDRG1, PDGFB, PIK3R1, PTAFR, PTGES, and RGS5 (Figure 3C).
Figure 3. Diet was associated with differential transcriptional profiles in iliac arteries. (A) Hierarchal clustering of normalized transcript levels of the top DETs (n = 174, FDR < 0.2) in response to diet. Calculated residual transcript level values were converted into z-scores before clustering. (B) Canonical pathways identified by Ingenuity Pathway Analysis from 174 DETs with FDRs < 0.2 in iliac arteries. Ratio indicates the number of target genes in the dataset divided by the total number of genes in the pathway. (C) Predicted activity of cardiovascular function- related pathways in response to the expression state of DETs (FDR < 0.2) using IPA. Arrows indicate directional relationships. Orange indicates predicted activation and blue indicates predicted inhibition of the biological pathway based on the expression levels of the shown genes in the Western cohort relative to the Mediterranean cohort. Red indicates upregulation and green indicates downregulation of the labeled gene in the Western cohort relative to the Mediterranean cohort.
Social status was significantly associated with the expression of two genes (FDR < 0.05), both of which are non-coding: (ENSMFAG00000064154 [LncRNA], and ENSMFAG00000057515 [small nucleolar RNA U13]; Table 2). At a relaxed FDR threshold of 0.2, we only identified 39 DETs, which is not enough to conduct a pathway analysis.
3.2.2 Effects of social status on the carotid artery transcriptome
For the carotid artery, 18 transcripts were differentially expressed between dominants and subordinates (FDR < 0.05; Figure 4). In the subordinate monkeys, 11 transcripts were upregulated (KCNQ4, STIM1, TNKS1BP1, CSNK1D, INPPL1, PNPLA7, F10, RAD9A, KCNIP3, ENSMFAG00000059809 [LncRNA], and ENSMFAG00000053865 [secreted protein A0A7N9CS45]), and seven transcripts were downregulated (IRAK1BP1, KIAA0513, SMIM15, PSMD14, TOPORS, ARPC2, and ENSMFAG00000050714 [LncRNA]) relative to the dominant monkeys (Figure 4; Supplementary Table 4).
Figure 4. Transcript levels (TMM normalized) differentially expressed by social status in the carotid artery (FDR < 0.05).
At a relaxed FDR threshold of 0.2, we identified 216 social status-associated DETs (Figure 5A; Supplementary Table 4). Regardless of the effect size direction (increased or decreased in subordinates compared to dominants), these DETs were enriched for a number of pathways including (1) Cardiac conduction, (2) ID1 Signaling Pathway, (3) Salvage Pathways of Pyrimidine Ribonucleotides, (4) Glioblastoma Multiforme Signaling, (5) Clathrin-mediated endocytosis, (6) TR/RXR Activation (Fisher’s Exact Test, p < 0.05; Figure 5B). A complete list of the enriched biological pathways is available in Supplementary Table 5. We also explored pathway enrichment taking into account the effect size direction (increased or decreased in subordinates compared to dominants) of the analyzed DETs. We identified and predicted the activity of cardiovascular disease- associated causal networks using IPA, including apoptosis and cytoskeleton formation. Cytoskeleton formation was predicted to be activated in subordinate compared to dominant monkeys (Z-score = 2.14, FDR = 0.03) according to the effect size of the following genes: ARHGAP18, ARHGEF10L, ARL2, CDC25B, CDC42EP2, COL5A3, NOTCH1, BIN1, and TACR1 (Figure 5C). Apoptosis was predicted to be activated in subordinate compared to dominant monkeys (Z-score = 1.4, FDR = 0.06) according to the effect size of the following genes: ACIN1, ARHGAP18, ARMC10, ATN1, BIN1, BRF1, CD44, CDIP1, GCLM, INPPL1, MDM2, MGP, MRE11, MST1R, PDGFC, PELI2, PIK3C3, PRDX4, SRCAP, SRI, TACR1, TAOK2, TGS1, TNFRSF10D, TSC2, TYMS, VPS13A, and YTHDF2 (Figure 5C).
Figure 5. Social status is associated with a differential transcriptional profile in carotid arteries. (A) Hierarchal clustering of normalized transcript levels of the top DETs (n = 216, FDR < 0.2) in response to social status. Calculated residual transcript level values were converted into z-scores before clustering. (B) Canonical pathways identified by Ingenuity Pathway Analysis from 216 DETs with FDRs < 0.2 in carotid arteries. Ratio indicates the number of target genes in the dataset divided by the total number of genes in the pathway. (C) Predicted activity of cardiovascular function- related pathways in response to the expression state of DETs (FDR < 0.02) using IPA. Arrows indicate directional relationships. Orange indicates predicted activation and blue indicates predicted inhibition of the biological pathway based on the expression levels of the shown genes in the subordinate cohort relative to the dominant cohort. Red indicates upregulation and green indicates downregulation of the labeled gene in the subordinate cohort relative to the dominant cohort.
There was no main effect of diet on the carotid artery transcriptome (all FDR > 0.05, Table 2). With a relaxed FDR threshold of 0.2, we identified 35 diet dependent transcripts in carotid arteries (Table 2).
3.3 Associations of iliac and carotid arteries transcriptomes with physical, behavioral, and physiological phenotypes
Differential transcriptional profiles in iliac and carotid arteries were compared with physical, behavioral, and physiologic phenotypes collected throughout the study using Spearman’s correlation test. A complete list of the analyzed phenotypes and the correlation coefficients with DETs is shown in Supplementary Tables 7, 8.
Iliac arteries: Relationships between the top diet associated DETs (FDR < 0.05) and multisystem phenotypes were assessed using Spearman’s test (Figure 6A). Coronary artery atherosclerosis (LAD-lesion area) was positively correlated with WDR62 and MARS1 transcript levels, and negatively correlated with FGGY, PAQR8, and RPS20 transcript levels (p < 0.05; Figure 6A). Telemetry heart rate variability measures were significantly correlated with the transcript levels of most of the analyzed DETs (p < 0.05; Figure 6A). Relative rank was positively correlated with WDR62, PIK3R1, PKDCC, ALDH3A2, SLC29A2, ZNF667, and CHCHD7 transcript levels; and negatively correlated with FGGY, PAQR8, MAMDC4, RAD21L1, and RPS20 transcript levels (p < 0.05; Figure 6A). The percentage of time spent depressed was positively correlated with PIK3R1, PABPC1, KDM3B, ANP32A, PAQR8, XPO7, ALDH3A2, TOMM20, CHCHD7, and ZNF667 transcript levels (p < 0.05; Figure 6A). Menstrual cycle length was positively correlated with PIK3R1, PABPC1, ANP32A, EIF4B, PAQR8, MARS1, RAD21L1, TOMM20, CHCHD7, and RPS20 transcript levels (p < 0.05; Figure 6A). Estradiol concentration at time of necropsy positively correlated with ANP32A and TOMM20 transcript levels (p < 0.05; Figure 6A).
Figure 6. Association of normalized transcript levels of differentially expressed genes in iliac and carotid arteries with physical, behavioral, and physiological phenotypes. (A) Diet adjusted Spearman’s correlation between DETs (FDR < 0.05) in iliac arteries with collected phenotypes. (B) Relative rank adjusted Spearman’s correlation between DETs (FDR < 0.05) in carotid arteries with collected phenotypes. Spearman’s test was conducted to identify significant correlations. Squares with colored circles inside indicate significant correlations (FDR < 0.05) between the indicated gene and the respective phenotype. Blue indicates positive correlation and red indicates negative correlation. RelRank: continuous value of social status, EstrogenLevels: plasma estrogen levels in pg./mL measured at necropsy, PulseWaveVelocity: pulse wave velocity measured by tail high definition oscillometry, DayHFHROscilation: percentage of high frequency (HF) oscillations during the day at 1600–1800 h, DayLF. HFHROscilation: ratio of LF/HF oscillations during the day at 1600–1800 h, NightLF. HFHROscilation: ratio of LF/HF oscillations during the night at 0100–0300 h, NightLFHROscilation: percentage of low frequency (LF) oscillations during the night at 0100–0300 h, NightHFHROscilation: percentage of high frequency (HF) oscillations during the night at 0100–0300 h, MeanCoronaryLesionArea: mean lesion area of the three coronary arteries collected, LADLesionArea: lesion area of the left anterior descending coronary artery, MeanCarotidLesionArea: mean lesion area of the three left common carotid arteries collected, MenstrualCycleLength: average length of cycles during experimental phase, measured in days, PerTimeDepressed: percent of time spent in depressive behavior.
Carotid arteries: Relationships between the top status-associated DETs (FDR < 0.05) with multisystem phenotypes were assessed using Spearman’s test (Figure 6B). LAD-lesion area (mm2) was positively correlated with ARPC2, IRAK1BP1, KIAA0513, PSMD14, SMIM15, and TOPORS transcript levels; and negatively correlated with CSNK1D, KCNQ4, PNPLA7, RAD9A, STM1, and TNKS1BP1 transcript levels (p < 0.05; Figure 6B). Telemetry heart rate variability measures were correlated with the transcript levels of many DETs (p < 0.05; Figure 6B). Diet was positively correlated with INPPL1, KCNQ4, STIM1, KCNIP3, and KIAA0513 transcript levels; and negatively correlated with ARPC2, PNPLA7, and TNKS1BP1 transcript levels (p < 0.05; Figure 6B). Percentage of time spent depressed was positively correlated with ARPC2, CSNK1D, F10, INPPL1, IRAK1BP1, PSMD14, RAD9A, SMIM15, STIM1, TNKS1BP1, KCNIP3, and TOPORS transcript levels (p < 0.05; Figure 6B). Menstrual cycle length was positively correlated with INPPL1, PSMD14, and RAD9A transcript levels (p < 0.05; Figure 6B).
4 Discussion
Our results demonstrate the associations of environmental influences (e.g., diet and social status) with atherosclerosis and vascular biology. Western diet promoted coronary artery atherogenesis and altered the arterial transcriptome in iliac but not carotid arteries relative to the Mediterranean diet. Social status was associated with 18 DETs in the carotid artery transcriptome and two DETs in the iliac artery. Transcriptional profiles in iliac and carotid arteries were significantly correlated with multiple physical, behavioral, and physiological phenotypes collected throughout the study, including coronary atherosclerosis, heart rate variability measures, menstrual cycle length, and depression behavior.
Atherosclerotic lesion development was increased in response to Western diet in the LAD coronary artery but not in any other coronary or peripheral arteries. Endothelial cell biology plays a key role in the initiation, development, and composition of atherosclerotic plaques (98, 99), and both local environmental factors and vascular bed origin influence atherosclerotic plaque development (100). Atherosclerosis susceptibility and lesion characteristics vary across different arterial sites (98, 101) with plaque formation more prevalent in coronary compared to femoral and carotid arteries (99). The LAD is the largest coronary artery, and supplies blood to the left side of the heart (102, 103), and is more atherosclerosis prone than the RCA and LCX (104), as observed in the present study. Intimal lesion area was measured in the sections close to branch points and bifurcations of the arteries; which are atherosclerosis-prone regions with disturbed blood flow promoting endothelial cell dysfunction (105, 106).
Atherosclerosis development was modest in the present study. Female macaques such as those in the present study are relatively protected from atherosclerosis, similar to the protection observed in premenopausal women (107). The Western and Mediterranean diets both contained a cholesterol content of 320 mg / 2000 calories/day with 30% of calories from fat, although the sources and types of fat in the two diets were significantly different, demonstrating that the sources of macronutrients and micronutrients contribute significantly to cardiovascular disease risk.
In the iliac arteries, diet-altered DETs had links to endothelial dysfunction, smooth muscle proliferation and migration, angiogenesis, and abnormal ECM dynamics (108–118). WDR62 (WD repeat-containing protein 62) is involved in the regulation of the cell cycle pathway and interacts directly with other cell cycle proteins to induce spindle formation and mitotic progression (109), and is a regulator of the cardiac myocyte cell cycle and a susceptibility gene for congenital heart defects (CHD) (109, 110). WDR62 upregulation in the Western cohort may be associated with hyperactive cell cycle and increased proliferation. WDR62 iliac transcripts were strongly positively correlated with plaque formation in the LAD coronary artery.
PIK3R1 (Phosphoinositide-3-Kinase Regulatory Subunit 1) mRNA expression was lower in the Western diet cohort relative to the Mediterranean diet cohort. PIK3R1 is involved in the regulation of insulin and PIK3R1/AKT pathways and exerts anti-apoptotic and protective effects on endothelial cells. PIK3R1 downregulation and/or inhibition is shown to promote atherosclerosis (112–115). Western diet also downregulated ANP32A (Acidic Nuclear Phosphoprotein 32 Family Member A), a translational repressor of the Wnt pathway known to impact cardiac disease (108). ANP32a-deficient mice have activated Wnt signaling and develop cardiac hypertrophy (108). The transcription factor ZNF667 (Zinc Finger Protein 667), can repress the expression of anti-angiogenesis genes and regulate expression of Wnt signaling genes in vascular endothelial cells, which in turn promote angiogenesis pathway (111). ZNF667 transcripts were downregulated in the Western diet cohort relative to the Mediterranean diet cohort, suggesting reduced angiogenesis.
Western diet downregulated transcripts for FGGY (FGGY carbohydrate kinase domain containing), and PAQR8 (Progestin and AdipoQ Receptor Family Member 8) a transmembrane progesterone receptor involved in steroid hormone response pathways. Both genes have been associated with vascular disease, cardiac structure and function, and ischemic stroke (116–118). Our analysis showed that transcript levels of both genes were strongly negatively correlated with intimal area in the LAD coronary artery, suggesting an association with plaque formation.
Ingenuity pathway analysis of the significant DETs in iliac arteries predicted an inhibition of autophagy in the Western diet cohort relative to the Mediterranean diet cohort. Autophagy is responsible for the degradation of proteins and organelles, to maintain metabolic homeostasis within the cell (119). Autophagy dysfunction is closely associated with many aging-related diseases including cardiovascular disorders (120). Multiple studies have demonstrated a protective effect of maintaining basal autophagy in atherosclerosis (121–123). Impaired autophagy leads to accumulation of cytotoxic and dysfunctional organelles and aggregates (124), contributing to endothelial dysfunction, monocyte/macrophage migration and adhesion that can lead to the development of atherosclerosis (125).
Collectively, this present study demonstrates that diet had a significant impact on coronary and iliac arteries that are associated with promoting atherosclerosis. In addition to this data, we have also shown before that diet induced numerous other effects spanning several tissues and phenotypes. Our previous findings demonstrated that compared to a Western diet, Mediterranean diet was protective against neuroinflammation (126) and pro-inflammatory monocyte polarization (127), reduced triglyceride levels (42), reduced social isolation and anxiety (67), protected against ovarian dysfunction (93) and insulin resistance 19, and enhanced stress resilience as indicated by lower sympathetic activity, overt heart rate responses to acute stress, and lower cortisol responses to acute psychological stress (73).
These findings highlight the multifaceted consequences of dietary composition on a wide range of physiological systems, suggesting that nutritional habits impact multiple health outcomes including cardiovascular, metabolic, endocrine, neurological, and behavioral endpoints, which shows the complexity of diet effects. This demonstrates the strength of this randomized preclinical approach in which it allows for comprehensive assessment of these widespread effects, providing strong evidence for the role of diet as a critical modulator of health and disease.
Diet did not induce major effects on the transcriptional profile of the carotid artery; on the other hand, social status- a proxy for psychosocial stress – impacted patterns of gene expression linked to dysregulated vascular tone and smooth muscle contractility, apoptosis, endothelial dysfunction, and abnormal ECM dynamics (128–132). Upregulation of KCNQ4 (voltage-gated potassium channel subunit Kv7.4) and KCNIP3 (potassium voltage-gated channel interacting protein 3) observed in the subordinates relative to dominants, has been shown previously in the literature to be associated with dysregulated vascular tone and contractility of carotid arteries. KCNQ4 forms a potassium channel that enables the selective passage of potassium ions across the cell membrane and play a critical role in the regulation of neuronal excitability, and the cardiac action potential (128). KCNQ4 is expressed in the central nervous system, as well as carotid, thoracic, and femoral arteries, specifically in the cell membranes of vascular smooth muscle cells (VSMC) (133). Activation of Kv7.2–5 channels results in vasodilation and relaxation of carotid arteries, while inhibition causes vasoconstriction and contraction, indicating that these channels play a significant role in controlling arterial tone of carotid arteries (128, 133). KCNIP3 plays a role in regulating potassium current across cell membrane and in turn regulates arterial tone. The transcript levels of both KCNQ4 and KCNIP3 are significantly correlated with heart rate variability measures, further suggesting an association between KCNQ4 and KCNIP3 transcript levels and dysregulated vascular tone in subordinate monkeys.
STIM1 (stromal interaction molecule 1), upregulated in subordinate carotids, functions as a calcium sensor that mediates Ca2+ influx after depletion of intracellular Ca2+ stores to regulate VSMC contractility. Increased activity of STIM1 can cause vasoconstriction and may affect the development or progression of hypertension (129, 130). STIM1-regulated Ca+2 homeostasis is crucial for induction of VSMC proliferation, and development. STIM1 mRNA and protein levels were reported to be upregulated in vivo in VSMC from balloon-injured rat carotid arteries compared to non-injured control vessels, and the knockdown of STIM1 suppressed neo-intimal hyperplasia in a rat carotid artery balloon injury model (134–136). STIM1 was also found to be upregulated in human coronary artery VSMCs and major regulator of their proliferation (137). Proliferation and migration of VSMC and ECM rearrangement are important early mechanisms in the pathogenesis of vascular disease.
Subordinate monkeys carotid arteries also showed increased INPPL1 (inositol polyphosphate phosphatase like1), also called Shc homology 2–containing inositol 5 phosphatase-2 (SHIP2), a lipid phosphatase that inhibits Akt phosphorylation and activation, inhibits insulin signaling, and induces apoptosis (131). Apoptosis induction plays a critical role in the pathogenesis of atherosclerotic cardiovascular disease (138). IPA pathway analysis of DETs predicted apoptosis to be activated in the subordinate monkeys relative to dominant monkeys. SHIP2 was also found to be associated with hypertension, obesity and type 2 diabetes (139). Carotid arteries of subordinate monkeys also showed upregulation of F10 (coagulation factor X), an important factor in forming blood clots at sites of injuries. Beyond the coagulation cascade, F10 can also affect inflammatory signaling by acting on endothelial, VSMCs and inflammatory cells, again associated with the pathogenesis of atherosclerosis. Coagulation protease inhibitors have been shown to attenuate atherogenesis in animal models 118. PNPLA7 and CSNK1D were upregulated in subordinate relative to dominant monkeys; PNPLA7 (patatin-like phospholipase), a regulator of adipocyte differentiation induced by metabolic stimuli, was shown to be upregulated in subjects with hypertension (140). CSNK1D (Casein kinase I delta) regulates cellular growth and survival including Wnt signaling, DNA repair and circadian rhythms. Large scale omics studies have shown that CSNK1D is associated with immune and inflammatory responses, LDL cholesterol concentrations, cardiometabolic risk, and myocardial infarction (141, 142). Subordinate monkey carotid arteries had reduced transcripts for TOPORS, an E3 ubiquitin-protein ligase that was shown in the literature to interact with syndecan-1 and inhibit cell growth in murine arterial smooth muscle cells (143).
It is interesting that social subordination had a higher impact on the transcriptional profile of carotid arteries than iliac arteries. Previous work demonstrated elevated carotid artery atherosclerosis in subordinate female monkeys (64). Additionally, Shively et al., showed previously that social subordination in macaques is stressful, evident by the heightened cortisol response to challenges (47, 50, 59), and elevated heart rate and blood pressure (50, 80, 144) which can strain the carotids. Psychosocial stress can also increase inflammation and cholesterol levels, contributing to plaque formation (52, 71). It was also shown that subordinate female long-tail macaques exhibit more depressive behavior than their dominant counterparts (80, 94). These characteristics of subordinates may explain in part the significant social status differences in gene expression of carotid arteries.
Given these data, and what is known previously in literature about SDH, it is critical to address SES-related disparities in vascular health through targeted interventions and policy reforms. This can be done by implementing comprehensive intervention strategies to enhance awareness of social disparities and their influence on disease progression and treatment outcomes. Additionally, the quality of diet of socially disadvantaged people must be improved, by implementing policies aimed at increasing the affordability and accessibility of healthy foods. Educational programs that promote nutrition awareness are critical. These efforts are vital for improving public health outcomes and ensuring equitable access to nutritious foods across all socioeconomic levels. Moreover, integrating social and biological data into clinical decision-making processes and risk assessment models is crucial for developing predictive strategies that incorporate both traditional biomarkers and socio-markers, thereby enhancing precision medicine.
This study provides evidence in support of differential effects of diet and psychosocial stress on vascular health. Strengths of the study include the use of stratified randomization to equalize the treatment groups, standardized and consistent environmental conditions across all participants, and the use of a well-established NHP model of cardiovascular and other health outcomes. Limitations include the use of bulk RNAseq to assess transcriptomic changes in arteries, which are complex tissues with cell types that are diverse in both provenance and function. Thus, we were unable to disentangle the contributions of various cell types, such as smooth muscle cells, fibroblasts, endothelial cells, and immune cell populations to the overall effects of diet and social stress. Future studies utilizing single-cell RNA sequencing of arterial tissues would facilitate a comprehensive analysis of cell specific responses to these environmental stressors. Our trial was also limited by relatively small sample size with limited statistical power to capture diet and social status interactions. This study did not include a group maintained on a standard laboratory chow diet, because “monkey chow” is not similar to any diet that humans or NHP eat. Furthermore, monkey chow has some very specific characteristics which may influence experimental results such as high levels of soy isoflavones which are selective estrogen receptor modulators. Our aim was to model human dietary patterns more realistically by directly comparing Western and Mediterranean diets. Finally, this study sample was focused on females, therefore additional study will be needed to fully understand the effects of diet and psychosocial stress on atherosclerosis and gene expression in male subjects.
In conclusion, our data demonstrates differential diet and social status effects dependent on arterial sites. Diet was associated with differences in coronary atherosclerosis and in the iliac artery transcriptome, while social status was associated with differences in the carotid artery transcriptome with a lesser effect in the iliac. In both iliac and carotid arteries, the resulting dysregulated pathways in response to diet and social status are all characteristics of atherosclerosis development and progression. These findings further demonstrate the detrimental effects of Western dietary habits and the stress from social subordination on vascular health. They also point toward the mechanisms through which Mediterranean diets may provide protection from the development of cardiovascular disease. Ongoing analyses are underway to assess relationships between arterial, adipose, and brain responses to diet and social status in this study.
Data availability statement
The original contributions in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Animal Care and Use Committee of Wake Forest University School of Medicine. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
AA: Data curation, Formal analysis, Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing. CJ: Methodology, Writing – review & editing, Formal analysis, Validation. NS-M: Methodology, Supervision, Writing – review & editing. KZ: Writing – review & editing, Formal analysis. JN: Writing – review & editing. KC: Writing – review & editing. BF: Writing – review & editing. TH: Methodology, Writing – review & editing. CS: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing. TR: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing, Formal analysis, Validation.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by grants from NIH: R01 HL087103 (CAS), R01 HL122393 (TCR), U24DK097748 (TCR), and the Wake Forest University School of Medicine Graduate School.
Acknowledgments
The authors thank Maryanne Post, JD Bottoms, Edison Floyd, Jason Lucas, Joshua Long, Rebecca Marcus, Diane Wood, and Sherri Samples for technical support.
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.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
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.2025.1564741/full#supplementary-material
Footnotes
References
1. Shively, CA, Frye, BM, Negrey, JD, Johnson, CSC, Sutphen, CL, Molina, AJA, et al. The interactive effects of psychosocial stress and diet composition on health in primates. Neurosci Biobehav Rev. (2023) 152:105320. doi: 10.1016/j.neubiorev.2023.105320
2. Kiecolt-Glaser, JK. Stress, food, and inflammation: psychoneuroimmunology and nutrition at the cutting edge. Psychosom Med. (2010) 72:365–9. doi: 10.1097/PSY.0b013e3181dbf489
3. O’Keefe, JH, Gheewala, NM, and O’Keefe, JO. Dietary strategies for improving post-prandial glucose, lipids, inflammation, and cardiovascular health. J Am Coll Cardiol. (2008) 51:249–55. doi: 10.1016/j.jacc.2007.10.016
4. Straznicky, NE, Louis, WJ, McGrade, P, and Howes, LG. The effects of dietary lipid modification on blood pressure, cardiovascular reactivity and sympathetic activity in man. J Hypertens. (1993) 11:427–37.
5. Maes, M, Song, C, Lin, A, De Jongh, R, Van Gastel, A, Kenis, G, et al. The effects of psychological stress on humans: increased production of pro-inflammatory cytokines and TH1-like response in stress-induced anxiety. Cytokine. (1998) 10:313–8.
6. Lopez-Garcia, E, Schulze, MB, Fung, TT, Meigs, JB, Rifai, N, Manson, JE, et al. Major dietary patterns are related to plasma concentrations of markers of inflammation and endothelial dysfunction. Am J Clin Nutr. (2004) 80:1029–35. doi: 10.1093/ajcn/80.4.1029
7. Clemente-Suárez, VJ, Beltrán-Velasco, AI, Redondo-Flórez, L, Martín-Rodríguez, A, and Tornero-Aguilera, JF. Global impacts of Western diet and its effects on metabolism and health: a narrative review. Nutrients. (2023) 15:2749. doi: 10.3390/nu15122749
8. Kopp, W. How western diet and lifestyle drive the pandemic of obesity and civilization diseases. Diabetes Metab Syndr Obes. (2019) 12:2221–36. doi: 10.2147/DMSO.S216791
9. Rakhra, V, Galappaththy, SL, Bulchandani, S, and Cabandugama, PK. Obesity and the western diet: how we got here. Mo Med. (2020) 117:536–8.
11. Giugliano, D, Ceriello, A, and Esposito, K. The effects of diet on inflammation. J Am Coll Cardiol. (2006) 48:677–85. doi: 10.1016/j.jacc.2006.03.052
12. Vahid, F, Wilk, P, and Bohn, T. Longitudinal effects of diet quality on healthy aging - focus on cardiometabolic health: findings from the Canadian longitudinal study on aging (CLSA). Aging Clin Exp Res. (2025) 37:157. doi: 10.1007/s40520-025-03058-9
13. Frölich, S, Lehmann, N, Weyers, S, Wahl, S, Dragano, N, Budde, T, et al. Association of dietary patterns with five-year degree and progression of coronary artery calcification in the Heinz Nixdorf recall study. Nutr Metab Cardiovasc Dis. (2017) 27:999–1007. doi: 10.1016/j.numecd.2017.09.002
14. Martínez-Beamonte, R, Barranquero, C, Gascón, S, Mariño, J, Arnal, C, Estopañán, G, et al. Effect of virgin olive oil as spreadable preparation on atherosclerosis compared to dairy butter in Apoe-deficient mice. J Physiol Biochem. (2024) 80:671–83. doi: 10.1007/s13105-024-01029-8
15. Lian, Z, Perrard, XYD, Peng, X, Raya, JL, Hernandez, AA, Johnson, CG, et al. Replacing saturated fat with unsaturated fat in Western diet reduces foamy monocytes and atherosclerosis in male Ldlr−/− mice. Arterioscler Thromb Vasc Biol. (2020) 40:72–85. doi: 10.1161/ATVBAHA.119.313078
16. Laugero, KD, Falcon, LM, and Tucker, KL. Relationship between perceived stress and dietary and activity patterns in older adults participating in the Boston Puerto Rican health study. Appetite. (2011) 56:194–204. doi: 10.1016/j.appet.2010.11.001
17. Luciano, M, Corley, J, Cox, SR, Valdés Hernández, MC, Craig, LCA, Dickie, DA, et al. Mediterranean-type diet and brain structural change from 73 to 76 years in a Scottish cohort. Neurology. (2017) 88:449–55. doi: 10.1212/WNL.0000000000003559
18. Romagnolo, DF, and Selmin, OI. Mediterranean diet and prevention of chronic diseases. Nutr Today. (2017) 52:208–22. doi: 10.1097/NT.0000000000000228
19. Alufer, L, Tsaban, G, Rinott, E, Kaplan, A, Meir, AY, Zelicha, H, et al. Long-term green-Mediterranean diet may favor fasting morning cortisol stress hormone; the DIRECT-PLUS clinical trial. Front Endocrinol. (2023) 14:910. doi: 10.3389/fendo.2023.1243910
20. Meyer-Lindemann, U, and Sager, HB. Neuroimmune crosstalk: how mental stress fuels vascular inflammation. Herz. (2024) 49:249–53. doi: 10.1007/s00059-024-05254-1
21. Ebong, IA, Quesada, O, Fonkoue, IT, Mattina, D, Sullivan, S, Oliveira, GMMde, et al. The role of psychosocial stress on cardiovascular disease in women. J Am Coll Cardiol (2024);84:298–314. doi: 10.1016/j.jacc.2024.05.016
22. Vaccarino, V, and Bremner, JD. Stress and cardiovascular disease: an update. Nat Rev Cardiol. (2024) 21:603–616. doi: 10.1038/s41569-024-01024-y
23. Dimsdale, JE. Psychological stress and cardiovascular disease. J Am Coll Cardiol. (2008) 51:1237–46. doi: 10.1016/j.jacc.2007.12.024
24. Sherman, GD, and Mehta, PH. Stress, cortisol, and social hierarchy. Curr Opin Psychol. (2020) 33:227–32. doi: 10.1016/j.copsyc.2019.09.013
25. Nurius, PS, Uehara, E, and Zatzick, DF. Intersection of stress, social disadvantage, and life course processes: reframing trauma and mental health. Am J Psychiatr Rehabil. (2013) 16:91–114. doi: 10.1080/15487768.2013.789688
26. Franks, P, Winters, PC, Tancredi, DJ, and Fiscella, KA. Do changes in traditional coronary heart disease risk factors over time explain the association between socio-economic status and coronary heart disease? BMC Cardiovasc Disord. (2011) 11:1–6. doi: 10.1186/1471-2261-11-28
27. Schultz, WM, Kelli, HM, Lisko, JC, Varghese, T, Shen, J, Sandesara, P, et al. Socioeconomic status and cardiovascular outcomes: challenges and interventions. Circulation. (2018) 137:2166–78. doi: 10.1161/CIRCULATIONAHA.117.029652
28. Costa, D, Ielapi, N, Bevacqua, E, Ciranni, S, Cristodoro, L, Torcia, G, et al. Social determinants of health and vascular diseases: a systematic review and call for action. Sociol Sci. (2023) 12:214. doi: 10.3390/socsci12040214
29. Tüchsen, F, Krause, N, Hannerz, H, Burr, H, and Kristensen, TS. Standing at work and varicose veins. Scand J Work Environ Health. (2000) 26:414–20. doi: 10.5271/sjweh.562
30. Demsas, F, Joiner, MM, Telma, K, Flores, AM, Teklu, S, and Ross, EG. Disparities in peripheral artery disease care: a review and call for action. Semin Vasc Surg. (2022) 35:141–54. doi: 10.1053/j.semvascsurg.2022.05.003
31. Davies, HO, Popplewell, M, Bate, G, Kelly, L, Koutsoumpelis, A, and Bradbury, AW. Impact of UK NICE clinical guidelines 168 and social deprivation on access to interventional treatment for symptomatic varicose vein and specialist referral for leg ulceration. Phlebology. (2017) 32:548–52. doi: 10.1177/0268355516677874
32. Barshes, NR, and Minc, SD. Healthcare disparities in vascular surgery: a critical review. J Vasc Surg. (2021) 74:6S–14S.e1. doi: 10.1016/j.jvs.2021.03.055
33. Becker, G, and Newsom, E. Socioeconomic status and dissatisfaction with health care among chronically ill African Americans. Am J Public Health. (2003) 93:742–8. doi: 10.2105/AJPH.93.5.742
34. Gómez, G, Kovalskys, I, Leme, A, Quesada, D, Rigotti, A, Cortés Sanabria, L, et al. Socioeconomic status impact on diet quality and body mass index in eight Latin American countries: ELANS study results. Nutrients. (2021) 13:2404. doi: 10.3390/nu13072404
35. Merhout, F, and Doyle, J. Socioeconomic status and diet quality in college students. J Nutr Educ Behav. (2019) 51:1107–12. doi: 10.1016/j.jneb.2019.06.021
36. Raffensperger, S, Fanelli Kuczmarski, M, Hotchkiss, L, Cotugna, N, Evans, MK, and Zonderman, AB. Effect of race and predictors of socioeconomic status on diet quality in the HANDLS study sample. J Natl Med Assoc. (2010) 102:923–30. doi: 10.1016/S0027-9684(15)30711-2
37. Pechey, R, and Monsivais, P. Socioeconomic inequalities in the healthiness of food choices: exploring the contributions of food expenditures. Prev Med. (2016) 88:203–9. doi: 10.1016/j.ypmed.2016.04.012
38. Scalise, E, Costa, D, Gallelli, G, Ielapi, N, Turchino, D, Accarino, G, et al. Biomarkers and social determinants in atherosclerotic arterial diseases: a scoping review. Ann Vasc Surg. (2025) 113:41–63. doi: 10.1016/j.avsg.2024.12.076
39. Alimohamadi, Y, Jalilvand, H, and Sepandi, M. The impact of air pollutants on the hospitalization of cardiovascular patients in the emergency room: a ten-year study. Med Sci. (2024) 28:1–12. doi: 10.54905/disssi.v28i147.e45ms3344
40. Serra, R, Abramo, A, Ielapi, N, Procopio, S, and Marino, P. Environmental pollution and peripheral artery disease. Risk Manag Healthc Policy. (2021) 14:2181–90. doi: 10.2147/RMHP.S307150
41. Miller, KG, Gianaros, PJ, Kamarck, TW, Anderson, BA, Muldoon, MF, and Manuck, SB. Cortisol activity partially accounts for a relationship between community socioeconomic position and atherosclerosis. Psychoneuroendocrinology. (2021) 131:105292. doi: 10.1016/j.psyneuen.2021.105292
42. Shively, CA, Appt, SE, Vitolins, MZ, Uberseder, B, Michalson, KT, Silverstein-Metzler, MG, et al. Mediterranean versus Western diet effects on caloric intake, obesity, metabolism, and hepatosteatosis in nonhuman primates: Mediterranean and Western diet in nonhuman primates. Obesity. (2019) 27:777–84. doi: 10.1002/oby.22436
43. Suchanek, P, Poledne, R, and Hubacek, JA. Dietary intake reports fidelity--fact or fiction? Neuro Endocrinol Lett. (2011) 32:29–31.
44. Stice, E, and Durant, S. Elevated objectively measured but not self-reported energy intake predicts future weight gain in adolescents. Appetite. (2014) 81:84–8. doi: 10.1016/j.appet.2014.06.012
45. Jarczok, MN, Koenig, J, Shively, CA, and Thayer, JF. Behavioral depression is associated with increased vagally mediated heart rate variability in adult female cynomolgus monkeys (Macaca fascicularis). Int J Psychophysiol. (2018) 131:139–43. doi: 10.1016/j.ijpsycho.2017.11.004
46. Kromrey, SA, Czoty, PW, Nader, SH, Register, TC, and Nader, MA. Preclinical laboratory assessments of predictors of social rank in female cynomolgus monkeys: predictors of social status in female cynomolgus monkeys. Am J Primatol. (2016) 78:402–17. doi: 10.1002/ajp.22514
47. Shively, CA. Social subordination stress, behavior, and central monoaminergic function in female cynomolgus monkeys. Biol Psychiatry. (1998) 44:882–91.
48. Clarkson, TB, Koritnik, DR, Weingand, KW, and Miller, LC. Nonhuman primate models of atherosclerosis: potential for the study of diabetes mellitus and hyperinsulinemia. Metabolism. (1985) 34:51–9.
49. Register, TC. Primate models in women’s health: inflammation and atherogenesis in female cynomolgus macaques (Macaca fascicularis). Am J Primatol. (2009) 71:766–75. doi: 10.1002/ajp.20722
50. Shively, CA, and Day, SM. Social inequalities in health in nonhuman primates. Neurobiol Stress. (2015) 1:156–63. doi: 10.1016/j.ynstr.2014.11.005
51. Willard, SL, and Shively, CA. Modeling depression in adult female cynomolgus monkeys (Macaca fascicularis): adult female monkey model of depression. Am J Primatol. (2012) 74:528–42. doi: 10.1002/ajp.21013
52. Shively, CA, Laber-Laird, K, and Anton, RF. Behavior and physiology of social stress and depression in female cynomolgus monkeys. Biol Psychiatry. (1997) 41:871–82.
53. Abbott, DH, Keverne, EB, Bercovitch, FB, Shively, CA, Mendoza, SP, Saltzman, W, et al. Are subordinates always stressed? A comparative analysis of rank differences in cortisol levels among primates. Horm Behav. (2003) 43:67–82. doi: 10.1016/S0018-506X(02)00037-5
54. Karasek, RA, Theorell, T, Schwartz, JE, Schnall, PL, Pieper, CF, and Michela, JL. Job characteristics in relation to the prevalence of myocardial infarction in the US health examination survey (HES) and the health and nutrition examination survey (HANES). Am J Public Health. (1988) 78:910–8.
55. Kaplan, JR, Chen, H, Appt, SE, Lees, CJ, Franke, AA, Berga, SL, et al. Impairment of ovarian function and associated health-related abnormalities are attributable to low social status in premenopausal monkeys and not mitigated by a high-isoflavone soy diet. Hum Reprod. (2010) 25:3083–94. doi: 10.1093/humrep/deq288
56. Walker, SE, Adams, MR, Franke, AA, and Register, TC. Effects of dietary soy protein on iliac and carotid artery atherosclerosis and gene expression in male monkeys. Atherosclerosis. (2008) 196:106–13. doi: 10.1016/j.atherosclerosis.2007.02.007
57. Walker, SE, Register, TC, Appt, SE, Adams, MR, Clarkson, TB, Chen, H, et al. Plasma lipid-dependent and -independent effects of dietary soy protein and social status on atherogenesis in premenopausal monkeys: implications for postmenopausal atherosclerosis burden. Menopause. (2008) 15:950–7. doi: 10.1097/gme.0b013e3181612cef
58. Shively, CA, Register, TC, Grant, KA, Johnson, JL, and Cline, JM. Effects of social status and moderate alcohol consumption on mammary gland and endometrium of surgically postmenopausal monkeys. Menopause. (2004) 11:389–99. doi: 10.1097/01.GME.0000109312.11228.62
59. Shively, CA, Register, TC, and Clarkson, TB. Social stress, visceral obesity, and coronary artery atherosclerosis: product of a primate adaptation. Am J Primatol. (2009) 71:742–51. doi: 10.1002/ajp.20706
60. Snyder-Mackler, N, Sanz, J, Kohn, JN, Brinkworth, JF, Morrow, S, Shaver, AO, et al. Social status alters immune regulation and response to infection in macaques. Science. (2016) 354:1041–5. doi: 10.1126/science.aah3580
61. Kohn, JN, Snyder-Mackler, N, Barreiro, LB, Johnson, ZP, Tung, J, and Wilson, ME. Dominance rank causally affects personality and glucocorticoid regulation in female rhesus macaques. Psychoneuroendocrinology. (2016) 74:179–88. doi: 10.1016/j.psyneuen.2016.09.005
62. Register, TC, Appt, SE, and Clarkson, TB. Atherosclerosis and vascular biologic responses to estrogens: histologic, immunohistochemical, biochemical, and molecular methods In: KM Eyster, editor. Estrogen receptors. New York, NY: Springer New York (2016). 517–32.
63. Clarkson, TB, Anthony, MS, and Morgan, TM. Inhibition of postmenopausal atherosclerosis progression: a comparison of the effects of conjugated equine estrogens and soy phytoestrogens 1. J Clin Endocrinol Metab. (2001) 86:41–7. doi: 10.1210/jcem.86.1.7151
64. Shively, CA, Kaplan, JR, and Clarkson, TB. Carotid artery atherosclerosis in cholesterol-fed female cynomolgus monkeys. Effects of oral contraceptive treatment, social factors, and regional adiposity. Arterioscler Thromb Vasc Biol. (1990) 10:358–66.
65. Mattison, JA, and Vaughan, KL. An overview of nonhuman primates in aging research. Exp Gerontol. (2017) 94:41–5. doi: 10.1016/j.exger.2016.12.005
66. Drescher, HK, Weiskirchen, R, Fülöp, A, Hopf, C, de San Román, EG, Huesgen, PF, et al. The influence of different fat sources on steatohepatitis and fibrosis development in the western diet mouse model of non-alcoholic steatohepatitis (NASH). Front Physiol. (2019) 10:1–17. doi: 10.3389/fphys.2019.00770
67. Johnson, CSC, Frye, BM, Register, TC, Snyder-Mackler, N, and Shively, CA. Mediterranean diet reduces social isolation and anxiety in adult female nonhuman Primates. Nutrients. (2022) 14:2852. doi: 10.3390/nu14142852
68. Pontzer, H, Wood, BM, and Raichlen, DA. Hunter-gatherers as models in public health: hunter-gatherer health and lifestyle. Obes Rev. (2018) 19:24–35. doi: 10.1111/obr.12785
69. Shively, CA, and Kaplan, JR. Stability of social status rankings of female cynomolgus monkeys, of varying reproductive condition, in different social groups. Am J Primatol. (1991) 23:239–45.
70. Shively, CA, Silverstein-Metzler, M, Justice, J, and Willard, SL. The impact of treatment with selective serotonin reuptake inhibitors on primate cardiovascular disease, behavior, and neuroanatomy. Neurosci Biobehav Rev. (2017) 74:433–43. doi: 10.1016/j.neubiorev.2016.08.037
71. Shively, CA, Kaplan, JR, and Adams, MR. Effects of ovariectomy, social instability and social status on female Macaca fascicularis social behavior. Physiol Behav. (1986) 36:1147–53.
72. Shively, C. The evolution of dominance hierarchies in nonhuman primate society In: SL Ellyson and JF Dovidio, editors. Power, dominance, and nonverbal behavior. New York, NY: Springer New York (1985). 67–87.
73. Shively, CA, Appt, SE, Chen, H, Day, SM, Frye, BM, Shaltout, HA, et al. Mediterranean diet, stress resilience, and aging in nonhuman primates. Neurobiol Stress. (2020) 13:100254. doi: 10.1016/j.ynstr.2020.100254
75. Troisi, A. Displacement activities as a behavioral measure of stress in nonhuman primates and human subjects. Stress. (2002) 5:47–54. doi: 10.1080/102538902900012378
76. Troisi, A, Belsanti, S, Bucci, AR, Mosco, C, Sinti, F, and Verucci, M. Affect regulation in alexithymia: an ethological study of displacement behavior during psychiatric interviews. J Nerv Ment Dis. (2000) 188:13–8. doi: 10.1097/00005053-200001000-00003
77. Maestripieri, D, Schino, G, Aureli, F, and Troisi, A. A modest proposal: displacement activities as an indicator of emotions in primates. Anim Behav. (1992) 44:967–79.
78. Suomi, SJ, Eisele, CD, Grady, SA, and Harlow, HF. Depressive behavior in adult monkeys following separation from family environment. J Abnorm Psychol. (1975) 84:576–8.
79. Shively, CA, Register, TC, Friedman, DP, Morgan, TM, Thompson, J, and Lanier, T. Social stress-associated depression in adult female cynomolgus monkeys (Macaca fascicularis). Biol Psychol. (2005) 69:67–84. doi: 10.1016/j.biopsycho.2004.11.006
80. Shively, CA, Musselman, DL, and Willard, SL. Stress, depression, and coronary artery disease: modeling comorbidity in female primates. Neurosci Biobehav Rev. (2009) 33:133–44. doi: 10.1016/j.neubiorev.2008.06.006
81. Shively, CA, Register, TC, Appt, SE, and Clarkson, TB. Effects of long-term sertraline treatment and depression on coronary artery atherosclerosis in premenopausal female Primates. Psychosom Med. (2015) 77:267–78. doi: 10.1097/PSY.0000000000000163
82. Dobin, A, Davis, CA, Schlesinger, F, Drenkow, J, Zaleski, C, Jha, S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635
83. Liao, Y, Smyth, GK, and Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. (2014) 30:923–30. doi: 10.1093/bioinformatics/btt656
84. Robinson, MD, McCarthy, DJ, and Smyth, GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. (2010) 26:139–40. doi: 10.1093/bioinformatics/btp616
85. Robinson, MD, and Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. (2010) 11:R25. doi: 10.1186/gb-2010-11-3-r25
86. Ritchie, ME, Phipson, B, Wu, D, Hu, Y, Law, CW, Shi, W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47–7. doi: 10.1093/nar/gkv007
87. Maćkiewicz, A, and Ratajczak, W. Principal components analysis (PCA). Comput Geosci. (1993) 19:303–42.
88. Langfelder, P, and Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. (2008) 9:1–13. doi: 10.1186/1471-2105-9-559
89. Akdemir, D. EMMREML: fitting mixed models with known covariance structures. (2015); Available online at: https://CRAN.R-project.org/package=EMMREML (Accessed January 10, 2025).
90. Hanghøj, K, Moltke, I, Andersen, PA, Manica, A, and Korneliussen, TS. Fast and accurate relatedness estimation from high-throughput sequencing data in the presence of inbreeding. Gigascience. (2019) 8:1–9. doi: 10.1093/gigascience/giz034/5481763
91. Li, H, Handsaker, B, Wysoker, A, Fennell, T, Ruan, J, Homer, N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352
92. Storey, J, Bass, A, Dabney, A, and Robinson, D. Qvalue: Q-value estimation for false discovery rate control. R package version 2.36.0. (2024); Available online at: http://github.com/jdstorey/qvalue (Accessed January 10, 2025).
93. Frye, BM, Register, TC, Appt, SE, Vitolins, MZ, Uberseder, B, Chen, H, et al. Differential effects of western versus mediterranean diets and psychosocial stress on ovarian function in female monkeys (Macaca fascicularis). Psychoneuroendocrinology. (2023) 153:106107. doi: 10.1016/j.psyneuen.2023.106107
94. Shively, CA. Depression and coronary artery atherosclerosis and reactivity in female Cynomolgus monkeys. Psychosom Med. (2002) 64:699–706. doi: 10.1097/01.PSY.0000021951.59258.C7
95. Shaffer, F, and Ginsberg, JP. An overview of heart rate variability metrics and norms. Front Public Health. (2017) 5:258. doi: 10.3389/fpubh.2017.00258
96. Shively, CA, Mietus, JE, Grant, KA, Goldberger, AL, Bennett, AJ, and Willard, SL. Effects of chronic moderate alcohol consumption and novel environment on heart rate variability in primates (Macaca fascicularis). Psychopharmacology. (2007) 192:183–91. doi: 10.1007/s00213-007-0709-z
97. Shively, CA, and Clarkson, TB. Social status and coronary artery atherosclerosis in female monkeys. Arterioscler Thromb Vasc Biol. (1994) 14:721–6.
98. Ross, R, and Epstein, FH. Atherosclerosis — an inflammatory disease. N Engl J Med. (1999) 340:115–26.
99. Dalager, S, Paaske, WP, Bayer Kristensen, I, Marsvin Laurberg, J, and Falk, E. Artery-related differences in atherosclerosis expression: implications for Atherogenesis and dynamics in intima-media thickness. Stroke. (2007) 38:2698–705. doi: 10.1161/STROKEAHA.107.486480
100. Debakey, ME, Lawrie, GM, and Glaeser, DH. Patterns of atherosclerosis and their surgical significance. Ann Surg. (1985) 201:132.
101. Aird, WC. Phenotypic heterogeneity of the endothelium: I. Structure, function, and mechanisms. Circ Res. (2007) 100:158–73. doi: 10.1161/01.RES.0000255691.76142.4a
102. Chatzizisis, YS, Giannoglou, GD, Parcharidis, GE, and Louridas, GE. Is left coronary system more susceptible to atherosclerosis than right? Int J Cardiol. (2007) 116:7–13. doi: 10.1016/j.ijcard.2006.03.029
103. Ahmed, IAM, Khalid, NHM, Abd-Elmagid, AEM, Abdullah, MAM, Musa, AMI, and AL-Qarni, NO. Common coronary artery occlusions in patients with myocardial infarction. Pan Afr Med J. (2022) 42:1–13. doi: 10.11604/pamj.2022.42.254.34762
104. Geerlings-Batt, J, and Sun, Z. Evaluation of the relationship between left coronary artery bifurcation angle and coronary artery disease: a systematic review. J Clin Med. (2022) 11:5143. doi: 10.3390/jcm11175143
106. Burridge, KA, and Friedman, MH. Environment and vascular bed origin influence differences in endothelial transcriptional profiles of coronary and iliac arteries. Am J Physiol-Heart Circ Physiol. (2010) 299:H837–46. doi: 10.1152/ajpheart.00002.2010
107. Nathan, L, and Chaudhuri, G. Estrogens and atherosclerosis. Annu Rev Pharmacol Toxicol. (1997) 37:477–515.
108. Monteagudo, S, Cornelis, FMF, Wang, X, de Roover, A, Peeters, T, Quintiens, J, et al. ANP32A represses Wnt signaling across tissues thereby protecting against osteoarthritis and heart disease. Osteoarthr Cartil. (2022) 30:724–34. doi: 10.1016/j.joca.2022.02.615
109. Hao, L, Ma, J, Wu, F, Ma, X, Qian, M, Sheng, W, et al. WDR62 variants contribute to congenital heart disease by inhibiting cardiomyocyte proliferation. Clin Transl Med. (2022) 12:1–18. doi: 10.1002/ctm2.941
110. Zaki, MS, Salam, GMHA, Saleem, SN, Dobyns, WB, Issa, MY, Sattar, S, et al. New recessive syndrome of microcephaly, cerebellar hypoplasia, and congenital heart conduction defect. Am J Med Genet A. (2011) 155:3035–41. doi: 10.1002/ajmg.a.34078
111. Wang, W, Shang, W, Zou, J, Liu, K, Liu, M, Qiu, X, et al. ZNF667 facilitates angiogenesis after myocardial ischemia through transcriptional regulation of VASH1 and Wnt signaling pathway. Int J Mol Med. (2022) 50:1–14. doi: 10.3892/ijmm.2022.5185
112. Li, H, Song, D, Liu, Q, Li, L, Sun, X, Guo, J, et al. Mir-351 promotes atherosclerosis in diabetes by inhibiting the ITGB3/PIK3R1/Akt pathway and induces endothelial cell injury and lipid accumulation. Mol Med. (2022) 28:1–21. doi: 10.1186/s10020-022-00547-9
113. Eyileten, C, Wicik, Z, De Rosa, S, Mirowska-Guzel, D, Soplinska, A, Indolfi, C, et al. MicroRNAs as diagnostic and prognostic biomarkers in ischemic stroke—a comprehensive review and Bioinformatic analysis. Cells. (2018) 7:249. doi: 10.3390/cells7120249
114. Jamshidi, Y, Gooljar, SB, Snieder, H, Wang, X, Ge, D, Swaminathan, R, et al. SHP-2 and PI3-kinase genes PTPN11 and PIK3R1 may influence serum apoB and LDL cholesterol levels in normal women. Atherosclerosis. (2007) 194:e26–33. doi: 10.1016/j.atherosclerosis.2006.12.013
115. Zhou, QB, Chen, Y, Zhang, Y, Li, DD, Wang, HQ, Jia, ZJ, et al. Hypermethylation effects of Yiqihuoxue decoction in diabetic atherosclerosis using genome-wide DNA methylation analyses. J Inflamm Res. (2022) 15:163–76. doi: 10.2147/JIR.S335374
116. Ren, P, Wang, JY, Chen, HL, Lin, XW, Zhao, YQ, and Guo, WZ. Diagnostic model constructed by nine inflammation-related genes for diagnosing ischemic stroke and reflecting the condition of immune-related cells. Front Immunol. (2022) 13:966. doi: 10.3389/fimmu.2022.1046966
117. Chai, T, Tian, M, Yang, X, Qiu, Z, Lin, X, and Chen, L. Association of circulating cathepsin B levels with blood pressure and aortic dilation. Front Cardiovasc Med. (2022) 9:1–11. doi: 10.3389/fcvm.2022.762468
118. Wild, PS, Felix, JF, Schillert, A, Teumer, A, Chen, MH, Leening, MJG, et al. Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function. J Clin Invest. (2017) 127:1798–812. doi: 10.1172/JCI84840
119. Mizushima, N, Levine, B, Cuervo, AM, and Klionsky, DJ. Autophagy fights disease through cellular self-digestion. Nature. (2008) 451:1069–75. doi: 10.1038/nature06639
120. Xu, J, Kitada, M, Ogura, Y, and Koya, D. Relationship between autophagy and metabolic syndrome characteristics in the pathogenesis of atherosclerosis. Front Cell Dev Biol. (2021) 9:1–14. doi: 10.3389/fcell.2021.641852
121. Kim, KH, and Lee, MS. Autophagy—a key player in cellular and body metabolism. Nat Rev Endocrinol. (2014) 10:322–37. doi: 10.1038/nrendo.2014.35
122. Kim, J, Lim, YM, and Lee, MS. The role of autophagy in systemic metabolism and human-type diabetes. Mol Cell. (2018) 41:11–7. doi: 10.14348/molcells.2018.2228
123. Nussenzweig, SC, Verma, S, and Finkel, T. The role of autophagy in vascular biology. Circ Res. (2015) 116:480–8. doi: 10.1161/CIRCRESAHA.116.303805
124. Zhang, Y, Whaley-Connell, AT, Sowers, JR, and Ren, J. Autophagy as an emerging target in cardiorenal metabolic disease: from pathophysiology to management. Pharmacol Ther. (2018) 191:1–22. doi: 10.1016/j.pharmthera.2018.06.004
125. Ma, M, Song, L, Yan, H, Liu, M, Zhang, L, Ma, Y, et al. Low dose tunicamycin enhances atherosclerotic plaque stability by inducing autophagy. Biochem Pharmacol. (2016) 100:51–60. doi: 10.1016/j.bcp.2015.11.020
126. Frye, BM, Negrey, JD, Johnson, CSC, Kim, J, Barcus, RA, Lockhart, SN, et al. Mediterranean diet protects against a neuroinflammatory cortical transcriptome: associations with brain volumetrics, peripheral inflammation, social isolation, and anxiety in nonhuman primates (Macaca fascicularis). Brain Behav Immun. (2024) 119:681–92. doi: 10.1016/j.bbi.2024.04.016
127. Johnson, CS, Shively, CA, Michalson, KT, Lea, AJ, DeBo, RJ, Howard, TD, et al. Contrasting effects of Western vs Mediterranean diets on monocyte inflammatory gene expression and social behavior in a primate model. eLife. (2021) 10:293. doi: 10.7554/eLife.68293
128. Yeung, SYM, Pucovský, V, Moffatt, JD, Saldanha, L, Schwake, M, Ohya, S, et al. Molecular expression and pharmacological identification of a role for K v 7 channels in murine vascular reactivity. Br J Pharmacol. (2007) 151:758–70. doi: 10.1038/sj.bjp.0707284
129. Guo, R w, Yang, L x, Li, M q, Pan, X h, Liu, B, and Deng, Y l. Stim1- and Orai1-mediated store-operated calcium entry is critical for angiotensin II-induced vascular smooth muscle cell proliferation. Cardiovasc Res. (2012) 93:360–70. doi: 10.1093/cvr/cvr307
130. Xu, F, Zhu, J, Chen, Y, He, K, Guo, J, Bai, S, et al. Physical interaction of tropomyosin 3 and STIM1 regulates vascular smooth muscle contractility and contributes to hypertension. Biomed Pharmacother. (2021) 134:111126. doi: 10.1016/j.biopha.2020.111126
131. Watt, NT, Gage, MC, Patel, PA, Viswambharan, H, Sukumar, P, Galloway, S, et al. Endothelial SHIP2 suppresses Nox2 NADPH oxidase–dependent vascular oxidative stress, endothelial dysfunction, and systemic insulin resistance. Diabetes. (2017) 66:2808–21. doi: 10.2337/db17-0062
132. Maeda, M, Tsuboi, T, and Hayashi, T. An inhibitor of activated blood coagulation factor X shows anti-endothelial senescence and anti-atherosclerotic effects. J Vasc Res. (2019) 56:181–90. doi: 10.1159/000499975
133. Jackson, WF. Ion channels and vascular tone. Hypertension. (2000) 35:173–8. doi: 10.1161/01.HYP.35.1.173
134. Bisaillon, JM, Motiani, RK, Gonzalez-Cobos, JC, Potier, M, Halligan, KE, Alzawahra, WF, et al. Essential role for STIM1/Orai1-mediated calcium influx in PDGF-induced smooth muscle migration. Am J Physiol-Cell Physiol. (2010) 298:C993–C1005. doi: 10.1152/ajpcell.00325.2009
135. Mancarella, S, Potireddy, S, Wang, Y, Gao, H, Gandhirajan, RK, Autieri, M, et al. Targeted STIM deletion impairs calcium homeostasis, NFAT activation, and growth of smooth muscle. FASEB J. (2013) 27:893–906. doi: 10.1096/fj.12-215293
136. Dai, F, Zhang, Y, Wang, Q, Li, D, Yang, Y, Ma, S, et al. Overexpression of SARAF ameliorates pressure overload–induced cardiac hypertrophy through suppressing STIM1-Orai1 in mice. Cell Physiol Biochem. (2018) 47:817–26. doi: 10.1159/000490036
137. Aubart, FC, Sassi, Y, Coulombe, A, Mougenot, N, Vrignaud, C, Leprince, P, et al. RNA interference targeting STIM1 suppresses vascular smooth muscle cell proliferation and Neointima formation in the rat. Mol Ther. (2009) 17:455–62. doi: 10.1038/mt.2008.291
138. DeKroon, R, Robinette, JB, Hjelmeland, AB, Wiggins, E, Blackwell, M, Mihovilovic, M, et al. APOE4 -VLDL inhibits the HDL-activated phosphatidylinositol 3-kinase/Akt pathway via the Phosphoinositol phosphatase SHIP2. Circ Res. (2006) 99:829–36. doi: 10.1161/01.RES.0000245479.03190.9f
139. Marcano, ACB, Burke, B, Gungadoo, J, Wallace, C, Kaisaki, PJ, Woon, PY, et al. Genetic association analysis of inositol polyphosphate phosphatase-like 1 (INPPL1, SHIP2) variants with essential hypertension. J Med Genet. (2007) 44:603–5. doi: 10.1136/jmg.2007.049718
140. Zhao, XC, Yang, SH, Yan, YQ, Zhang, X, Zhang, L, Jiao, B, et al. Identification of differential gene expression profile from peripheral blood cells of military pilots with hypertension by RNA sequencing analysis. BMC Med Genet. (2018) 11:1–10. doi: 10.1186/s12920-018-0378-2
141. Scott, WA, and Mahoney, E. Defects of glucose-6-phosphate and 6-phosphogluconate dehydrogenases in Neurospora and their pleiotropic effects. Curr Top Cell Regul. (1976) 10:205–36. doi: 10.1016/B978-0-12-152810-2.50012-7
142. Ahmadi, SA, Tranmer, JE, Ritonja, JA, Flaten, L, Topouza, DG, Duan, QL, et al. DNA methylation of circadian genes and markers of cardiometabolic risk in female hospital workers: an exploratory study. Chronobiol Int. (2022) 39:735–46. doi: 10.1080/07420528.2022.2032729
143. Braun, KR, DeWispelaere, AM, Bressler, SL, Fukai, N, Kenagy, RD, Chen, L, et al. Inhibition of PDGF-B induction and cell growth by Syndecan-1 involves the ubiquitin and SUMO-1 ligase, Topors. PLoS One. (2012) 7:e43701. doi: 10.1371/journal.pone.0043701
144. Shively, CA, Register, TC, and Clarkson, TB. Social stress, visceral obesity, and coronary artery atherosclerosis in female Primates. Obesity. (2009) 17:1513–20. doi: 10.1038/oby.2009.74
Glossary
CVD - Cardiovascular Disease
DDS - Dietary Diversity Score
DET - Differentially Expressed Transcripts
DQS - Dietary Quality Score
ECG - Electrocardiogram
ECM - Extracellular Matrix
IPA - Ingenuity Pathway Analysis
LAD - Left Anterior Descending
LCX - Left Circumflex
NHP - Nonhuman Primate
PCA - Principal Component Analysis
RCA - Right Coronary Artery
RIN - RNA Integrity Number
SDH - Social Determinants Of Health
SES - Socioeconomic Status
TMM - Trimmed Mean Of M-Value
WGCNA - Weighted Gene Co-Expression Network Analysis
ACIN1 - Apoptotic Chromatin Condensation Inducer 1
ALDH3A2 - Aldehyde Dehydrogenase 3 Family Member A2
ANP32A - Acidic Nuclear Phosphoprotein 32 Family Member A
ARHGAP18 - Rho Gtpase Activating Protein 18
ARHGEF10L - Rho Guanine Nucleotide Exchange Factor 10 Like
ARL2 - ARF Like Gtpase
ARMC10 - Armadillo Repeat Containing 10
ARPC2 - Actin-Related Protein 2/3 Complex Subunit 2
ATN1 - Atrophin 1
BIN1 - Bridging Integrator 1
BRF1 - BRF1 General Transcription Factor IIIB Subunit
CD44 - CD44- Receptor For Hyaluronic Acid
CDC25B - Cell Division Cycle 25
CDC42EP2 - CDC42 Effector Protein 2
CDIP1 - Cell Death Inducing P53 Target 1
CHCHD7 - Coiled-Coil-Helix-Coiled-Coil-Helix Domain Containing 7
COL5A3 - Collagen Type V Alpha 3 Chain
CSNK1D - Casein Kinase 1 Delt
EIF4B - Eukaryotic Translation Initiation Factor 4B
ENSMFAG00000050714 - [Lncrna]
ENSMFAG00000052859 - Unnamed Gene
ENSMFAG00000053865 - [Secretedproteina0a7n9cs45]
ENSMFAG00000057515 - Secreted Protein A0A7N9CS45
ENSMFAG00000059809 - [Lncrna]
ENSMFAG00000064154 - [Lncrna]
F10 - Coagulation Factor X
FGGY - FGGY Carbohydrate Kinase Domain-Containing Protein
GCLM - Glutamate-Cysteine Ligase Modifier Subunit
INPPL1 - Inositol Polyphosphate-5 Phosphatase-Like 1
IRAK1BP1 - Interleukin 1 Receptor Associated Kinase 1 Binding Protein 1
KCNIP3 - Potassium Voltage-Gated Channel Interacting Protein 3
KCNQ4 - Potassium Voltage-Gated Channel Subfamily Q Member 4
KDM3B - Lysine Demethylase 3B
KIAA0513 - Undefined Gene
MAMDC4 - MAM Domain Containing 4
MARS1 - Methionyl-tRNA Synthetase 1
MDM2 - Mouse Double Minute 2 Homolog
MGP - Matrix Gla Protein
MRE11 - MRE11 Homolog, Double Strand Break Repair Nuclease
MST1R - Macrophage Stimulating 1 Receptor
NOTCH1 - Notch Receptor 1
PABPC1 - Polyadenylate-Binding Protein 1
PAQR8 - Progestin and AdipoQ Receptor Family Member 8
PDGFC - Platelet Derived Growth Factor C
PELI2 - Pellino E3 Ubiquitin Protein Ligase Family Member 2
PIK3C3 - Phosphatidylinositol 3-Kinase Catalytic Subunit Type 3
PIK3R1 - Phosphatidylinositol 3-Kinase
PKDCC - Protein Kinase Domain Containing, Cytoplasmic
PNPLA7 - Patatin Like Domain 7, Lysophospholipase
PRDX4 - Peroxiredoxin 4
PSMD14 - Proteasome 26S Subunit, Non-Atpase 14
RAD21L1 - RAD21 Cohesin Complex Component Like 1
RAD9A - RAD9 Checkpoint Clamp Component A
RPS20 - Ribosomal Protein S20
SLC29A2 - Solute Carrier Family 29 Member 2
SMIM15 - Small Integral Membrane Protein 1
SRCAP - Snf2 Related CREBBP Activator Protein
SRI - Sorcin
STIM1 - Stromal Interaction Molecule 1
TACR1 - Tachykinin Receptor 1
TAOK2 - Thousand And One Amino Acid Protein Kinase
TGS1 - Trimethylguanosine Synthase 1
TNFRSF10D - TNF Receptor Superfamily Member 10d
TNKS1BP1 - Tankyrase 1 Binding Protein 1
TOMM20 - Translocase Of Outer Mitochondrial Membrane 20
TOPORS - TOP1 Binding Arginine/Serine Rich Protein, E3 Ubiquitin Ligase
TSC2 - TSC Complex Subunit 2 (Tuberin)
TYMS - Thymidylate Synthetase
VPS13A - Vacuolar Protein Sorting 13 Homolog A
WDR62 - WD Repeat-Containing Protein 62
XPO7 - Exportin 7
YTHDF2 - YTH N6-Methyladenosine RNA Binding Protein F2 2
ZNF667 - Zinc Finger Protein 667
Keywords: iliac artery, carotid artery, atherosclerosis, Mediterranean diet, Western diet, nonhuman primates, social status, arterial transcriptome
Citation: Abusheikha AJ, Johnson CSC, Snyder-Mackler N, Zimmerman KD, Negrey JD, Chiou KL, Frye BM, Howard TD, Shively CA and Register TC (2025) Differential effects of Mediterranean vs. Western diets on coronary atherosclerosis and peripheral artery transcriptomics. Front. Nutr. 12:1564741. doi: 10.3389/fnut.2025.1564741
Edited by:
Gaytán-Martínez Marcela, Autonomous University of Queretaro, MexicoReviewed by:
Davide Costa, Magna Græcia University, ItalyMaría De La Luz Reyes, Autonomous University of Coahuila, Mexico
Copyright © 2025 Abusheikha, Johnson, Snyder-Mackler, Zimmerman, Negrey, Chiou, Frye, Howard, Shively and Register. 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: Thomas C. Register, cmVnaXN0ZXJAd2FrZWhlYWx0aC5lZHU=; Carol A. Shively, Y3NoaXZlbHlAd2FrZWhlYWx0aC5lZHU=