Stress Overload and DNA Methylation in African American Women in the Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure Study.

Introduction: Experiencing psychosocial stress is associated with poor health outcomes such as hypertension and obesity, which are risk factors for developing cardiovascular disease. African American women experience disproportionate risk for cardiovascular disease including exposure to high levels of psychosocial stress. We hypothesized that psychosocial stress, such as perceived stress overload, may influence epigenetic marks, specifically DNA methylation (DNAm), that contribute to increased risk for cardiovascular disease in African American women.
Methods: We conducted an epigenome-wide study evaluating the relationship of psychosocial stress and DNAm among African American mothers from the Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure (InterGEN) cohort. Linear mixed effects models were used to explore the epigenome-wide associations with the Stress Overload Scale (SOS), which examines self-reported past-week stress, event load and personal vulnerability.
Results: In total, n = 228 participants were included in our analysis. After adjusting for known epigenetic confounders, we did not identify any DNAm sites associated with maternal report of stress measured by SOS after controlling for multiple comparisons. Several of the top differentially methylated CpG sites related to SOS score (P < 1 × 10-5), mapped to genes of unknown significance for hypertension or heart disease, namely, PXDNL and C22orf42. Conclusions: This study provides foundational knowledge for future studies examining epigenetic associations with stress and other psychosocial measures in African Americans, a key area for growth in epigenetics. Future studies including larger sample sizes and replication data are warranted.

Introduction

Evidence indicates that stress is an important risk factor in regards to cardiovascular disease.[1,2] For example, stress has been significantly associated with hypertension, a key risk factor for cardiovascular disease.[3-5] Black women have the unique intersectional reality of experiencing stressors such as discrimination on the basis of gender and race also known as “gendered racism,”[6] that has been shown to negatively affect health outcomes among African American women.[7,8] Other stressors such as parenting stress[9] and financial stress disproportionately affect Black women[10] and have been associated with adverse health outcomes.[11]

Women report higher overall levels and different types of stress compared with men [12]; and sex differences in pathophysiological and psychosocial mechanisms related to stress have been reported.[13-15] It is likely that genomic and social factors, including pyschosocial stress, are driving hypertension disparities.[16] Changes in DNA methylation can be influenced by environmental stressors, and has been identified as a potential mechanism by which stress may predispose individuals to chronic health conditions such as hypertension, obesity, and diabetes.[17-19] Previous research reported that parenting stress was associated with loss of methylation at over 90 cytosine-phosphate-guanine (CpG) sites, notably in genes associated with the stress signaling pathway including poly (ADP-ribose) polymerase-1 (PARP-1).[9] These findings indicated that higher levels of parenting stress were significantly associated with DNA methylation (DNAm) changes among young African American mothers.

To explore whether psychosocial stress is associated with DNA methylation, we conducted an epigenome-wide association study (EWAS) among young African American mothers who participated in the Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure (InterGEN) Study. We hypothesized that there would be associations between perceived stress and DNA methylation, given prior literature that have reported the linkages between stress and DNA methylation.[20-22]

Materials and Methods

For the present secondary analysis, we examined data from the InterGEN study, which was a prospective cohort study that examined the effects of genetic, epigenetic, and psychosocial factors on blood pressure in 250 mother/child dyads. Eligibility criteria included the following characteristics: mothers (⩾21 years old) who self-identified as African American or Black, spoke English and had no diagnosis of mental illnesses that could prevent collection of psychosocial measures. The mother-child dyads were recruited from childcare centers in Connecticut over a 5-year period between 2014 and 2019, and participated in 4 interviews over a period of 2 years. Sociodemographic measures as well as other self-reported psychosocial measures (eg, parenting stress, perceived racism and discrimination, and depression) were collected through using Audio Computer-Assisted Self-Interview (ACASI) software. Additional recruiting measures and practices have been previously described.[23,24] Institutional Review Board (IRB) approval was obtained from Yale University, New York University and Columbia University.

Survey measures

The Stress Overload Scale (SOS) is a validated instrument that was developed to measure psychosocial stress.[25] The SOS is a practical, multiculturally sensitive, and brief questionnaire that captures perception of overload, which is a state of feeling overwhelmed or overextended. It was derived using community-based samples and has been used in diverse populations. The SOS is comprised of 2 subscales that total 24 items (min: 24; max: 119). The 2 subscales include a 12-item personal vulnerability scale (min: 12; max 59) and a 12-item event load scale (min: 12; max: 60). Both measures capture levels of self-reported vulnerability (ie, feelings of powerlessness and frailty) in the past week as well as self-reported past-week event-load (ie, external demands and pressures).[25] Items are scored on a 5-point Likert scale from 1 (not at all) to 5 (a lot). Subscales were calculated by adding totals in each subscale, and the total SOS score was calculated by summing the total subscale scores.

Potential confounding variables

Mothers self-reported their age, if they were currently smokers (yes/no) and other demographic data at the initial interview. We adjusted for age and smoking, which are known confounders in epigenetic studies.[26] We also adjusted for batch effects and potential heterogeneity in cell proportions from saliva using the reference-free EWAS method by Houseman et al.[27] This is also a method frequently utilized in previous research and literature.[28]

DNA collection, processing and DNA methylation profiling

Saliva was collected from participants for DNA analysis using Oragene (OG)−500 format tubes.[29] This procedure requires that participants fill a tube with 2 mL of saliva. DNA extraction and processing was performed based on protocol and standard DNA processing for ReliaPrep, which has been described previously.[23] To maintain confidentiality and ensure accuracy of analyses, tubes were barcoded and stored in a laboratory freezer. A barcode scanning system was used for DNA pipetting with robotic workstations to maintain the accuracy and integrity during the transfer process from tubes to plates. This process ensured that accurately identified participants’ DNA for correct merging to genotype calls.

To quantify DNA methylation, the Illumina Infinium Methylation EPIC (850K) BeadChip was used, and quantile-normalization of beta values was completed for autosomal CpG sites. Laboratory-based quality control procedures (missing rate <10% and no sex mismatch) were implemented for each saliva sample. CpG sites were excluded if they had a missing rate greater than 10% (n = 3343), overlapped with single nucleotide polymorphisms (n = 87 074), or were listed in the recent Illumina quality notice (n = 977). Ultimately, a total of 756 448 autosomal sites were included in the association analyses.

Statistical analyses

An epigenome-wide association study (EWAS) analysis of total stress overload scores was conducted among InterGEN mothers (N = 228). We modeled DNAm beta-values as the dependent variable and linear mixed effects models with 2 random effects (batch and chip) were applied to ascertain epigenetic associations with stress overload scores adjusted for age, smoking status and cell-type heterogeneity. We used the reference-free EWAS method to account for cell proportions for saliva samples.[27] Principal components (PCs) were also derived using DNAm data[30] and we included the top 10 PCs as fixed-effect covariates to account for potential population stratification. False discovery rate (FDR) was used to correct for multiple comparisons and FDR-corrected P < .05 was defined as EWAS significance threshold. We utilized this FDR given its use in other studies in InterGEN that has been previously discussed.[31] We conducted analyses using the R statistical package environment[32] and a selection of packages from Bioconductor.

Results

In our study sample of 250 mothers, 16 were missing data for the SOS total score. Of the 234 participants with complete SOS data, 6 participants’ DNAm samples did not pass quality assessment and were excluded from analysis. A total of N = 228 mothers were included in the EWAS analysis of stress overload, and characteristics of the study sample are summarized in Table 1. The mean maternal age was 31.3 years. Fifty-eight (25.4%) mothers self-reported that they were smokers. Over half of the sample had completed at least some college (59.2%). With regards to income, nearly half of the sample (46.1%) earned less than or equal to $15 000. More than two-thirds (66.2%) of the sample designated single as their marital status. Mothers reported a mean score of 59.9 (s.d. = 24.2) for stress overload. The mean scores for personal vulnerability and event load were 26.2 (s.d. = 12.4) and 33.7 (s.d. = 13.3) respectively.

Table 1.

Descriptive statistics of study sample (N = 228).

Characteristic	Mean (SD) or N (%)
Age (years)	31.3 (5.7)
Stress Overload Scale total score	59.9 (24.2)
Personal Vulnerability subscale	26.2 (12.4)
Event Load subscale	33.7 (13.3)
Current smokers	58 (25.4)
Education level
Less than high school	12 (5.3)
High school	81 (35.5)
Some college or higher	135 (59.2)
Annual household income
Less than $5000	51 (22.4)
$5000 to $9999	29 (12.7)
$10 000 to $14 999	25 (11.0)
$15 000 to $19 999	18 (7.9)
$20 000 or above	101 (44.3)
Marital Status
Married	52 (22.8)
Single	151 (66.2)
Divorced	11 (4.8)
Separated	3 (1.3)
Living with significant other	11 (4.8)

We did not find significant associations between SOS score and DNA methylation among InterGEN mothers. (Figures 1 and 2) However, several of the top (P < 1 × 10−5) differentially methylated CpGs were positionally mapped to PXDNL, C22orf42, and SND1 genes in mothers (Table 2). Two genes, PXDNL and C22orf42, were identified as top signals for the total stress overload scores in mothers and were also top hits for 2 of the stress overload subscales: personal vulnerability and event load (Supplemental Tables 1 and 2). Hypomethylation of cg20043066 (mapped to intron 4 of PXDNL) was associated with increased total stress overload score (P = 7 × 10−7) as well as increased event load (P = 9.5 × 10−6). Additionally, hypomethylation of cg24060571 (mapped to 927 bp upstream of C22orf42) was associated with increased total stress score (P = 4.6 × 10−6) and increased personal vulnerability (P = 1 × 10−5).

Figure 1.

Quantile-quantile plot of epigenome-wide associations with self-reported stress for mothers. Associations of DNAm level at each CpG site with total stress score was tested using linear mixed model adjusted for age, smoking, cell-type proportions and top 10 principal components. The total sample size is n = 228, and the overall inflation factor is 0.983 showing well-controlled type I error rate. Observed -log10(p) (y-axis) was plotted against expected -log10(p) derived from a uniform distribution (x-axis). Red straight line is y = x and the red curved lines indicate the 95% confidence interval for the expected −log10(p).

Figure 2.

Manhattan plot of epigenome-wide associations with self-reported stress for mothers. Associations of DNAm level at each CpG site with total stress score was tested using linear mixed model adjusted for age, smoking, cell-type proportions and top 10 principal components. The total sample size is n = 228. −log10(p) of each CpG site is plotted against its genomic position (hg19). Red line indicates the P-value cutoff for epigenome-wide significance.

Table 2.

Top epigenetic associations with self-reported stress for mothers adjusted for age, smoking, cell-type proportions and 10 principal components.

CpG	CHR	BP*	Gene	Beta	SE	P-value
cg20043066	8	52467264	PXDNL	−4.9 × 10−4	9.4 × 10−5	7.0 × 10−7
cg24060571	22	32556170	C22orf42	−1.5 × 10−4	3.1 × 10−5	4.6 × 10−6
cg17866650	2	37798344		−6.4 × 10−4	1.3 × 10−4	5.6 × 10−6
cg08575835	3	170448125		1.4 × 10−4	2.9 × 10−5	6.3 × 10−6
cg10057227	7	127607004	SND1	−1.4 × 10−4	3.0 × 10−5	1.0 × 10−5

Base-pair position using human genome build hg19.

Discussion

This study explored associations between perceived stress overload and epigenome-wide DNA methylation in African American mothers. In this cohort of African American women, perceived stress overload was not associated with differences in epigenome-wide DNA methylation. However, we did identify differentially methylated sites based on stress overload scores that warrant follow-up testing in larger cohorts with improved power to detect significant differences across the epigenome.

PXDNL, the gene that encodes the peroxidasin-like protein, has been found to be exclusively expressed in the heart[33] and previous research suggests that mutations in the PXDNL gene may lead to ventricular arrhythmias.[34] Although the exact function of peroxidasin-like protein remains unknown, it has been implicated in cardiovascular health.[35-37] Other peroxidases in the heart are responsible for breaking down hydrogen peroxide (H2O2), a common reactive oxygen species generated during aerobic respiration with known roles in the etiology of hypertension.[38-40] It is possible that if hypomethylation of this gene is associated with increased expression, it may be an adaptive marker associated the cellular responsiveness to experiencing increased stress. Further research is required to understand the exact consequences of differential methylation in the expression of PXDNL, and how that could contribute to the etiology of hypertension.

Another gene that exhibited a large difference in methylation associated with stress was C22orf42, a validated open reading frame that has been found to be expressed profusely in testis, prostate, some brain cells and adrenal gland.[41,42] There is growing evidence to suggest that there is a relationship between stress in early life and the HPA axis. Some studies have linked the stress and HPA activation with mental health outcomes[43] and alcohol dependence.[44] However, there is still a critical research gap regarding the role of the protein produced by C22orf42 and its potential role in cardiovascular health.

The construct of stress overload among African American women is complex and may have been underreported as given cultural nuances and norms that normalize portraying strength and suppressing stress.[45] Others have explained how maladaptive coping mechanisms such as John Henryism, which involves “high effort coping” may lead to adverse health outcomes; particularly as it pertains to cardiovascular disease risk.[46,47] Other constructs such as the “Strong Black Woman Schema” or “Superwoman Syndrome” posit that despite exposure to stress and discrimination, African Americans may feel compelled to relay strength whilst not acknowledging or suppressing their perceptions of these stressors.[48,49] In turn, this may have deleterious health effects and may elevate risk for hypertension and cardiovascular risk.[50] Steinhardt et al.[51] found that psychosocial intervention that explored stress, resilience, and diabetes self-management among diabetic African Americans found improvements in physical risk factors such as BMI and systolic and diastolic blood pressure that were statistically significant as well as in psychosocial factors such as empowerment and diabetes self-management. However, in regards to stress, there were no statistically significant changes. These findings suggest that that while individuals’ perception of their stress may not change, the underlying effects of stress could be influencing epigenetic changes to occur.

A prior study examining stress, coping and methylation of blood pressure-related genes in InterGEN, showed that stress overload and coping strategies were associated with epigenetic changes, though these relationships did not hold after accounting for multiple testing adjustment.[52] Additionally, stress and coping were not found to be significantly associated with blood pressure candidate genes.[52] This present analysis expanded on the previous study by employing an epigenome-wide approach to assess potential associations of stress overload with DNAm. Other InterGEN analyses have also demonstrated the importance of chronic stressors such as neighborhood stress and experiencing racial discrimination and how exposure to such stressors negatively impact mental and physical health.[53-55]

It is imperative that future research explore culturally-sensitive ways to capture stress and examine possible changes to stress levels over time. While there are several validated measures to capture stress such as the Perceived Stress Scale (PSS) (past 30 days stress),[56] Weekly Stress Inventory (WSI) (past week stress),[57] and others, there is no consensus for measuring stress/ideal timeframes for measurement. Similarly, this is also apparent with measures of stress in childhood,[58,59] Longitudinal data regarding the examination and dynamic nature of stress is needed.

This study addressed a gap in the literature on epigenomic associations with stress among a cohort of African American women, a group that is often underrepresented in epigenomic studies.[60] However, there are some limitations to consider. First, we do not have a replication analysis for our results. Future studies can expand on our sample size using multiple cohorts with meta analyses to increase power for gene discovery. Given our sample size, it is possible that we were underpowered to detect these associations. Also, our measurement of stress may not have captured the full spectrum of stressors faced by African American women. Studies have shown that African American women face a variety of stressors, including stress related to parenting,[9] as well as racism and gender discrimination.[61] Additionally, the SOS assesses past week levels of stress and may not be representative of cumulative stress, which may have a considerable influence on DNAm. Lastly, it is also important to consider that there are many other social determinants outside the scope of this analysis that could have a considerable effect on stress and/or DNAm such as neighborhood stress/safety[62] and social support,[63] for example.

Conclusion

There has been a significant interest in expanding knowledge on how social epigenomics may be influencing disparities; particularly among African American women who are disproportionately burdened by health disparities. Despite increased interest, there remains a dearth of epigenetic studies that include African Americans participants. This analysis provides a foundation for future studies examining epigenetic associations with stress and other psychosocial measures among African Americans.

Click here for additional data file.

Supplemental material, sj-docx-1-gae-10.1177_25168657221126314 for Stress Overload and DNA Methylation in African American Women in the Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure Study by Jolaade Kalinowski, Yunfeng Huang, Martin A Rivas, Veronica Barcelona, Michelle L Wright, Cindy Crusto, Tanya Spruill, Yan V Sun and Jacquelyn Y Taylor in Epigenetics Insights