Religious Involvement, Inflammatory Markers and Stress Hormones in Major Depression and Chronic Medical Illness ()
1. Introduction
Depressive disorder and religious involvement are common in persons with chronic medical illness. Depressive disorder has been reported in 19% to 45% of those with chronic illness depending on the setting and diagnostic method [1] -[3] , and physical disability is one of the strongest predictors of depressive symptoms in the general population [4] -[6] . Religious involvement is important to many individuals, especially those with chronic illness, since religion is often used to cope with the life changes caused by health problems. While over half of healthy persons in the U.S. indicate that religion is an important part of daily life [7] [8] , the importance of religion in- creases even further among those with medical illness. In some areas of the country, up to 90% of medically ill persons report that they depend on religion to cope with, with over 40% saying that it is the most important fac- tor that keeps them going [9] .
Depression is also often accompanied by physiological changes that can adversely affect the course of medi- cal illness over time, including increased levels of pro-inflammatory cytokines, decreased levels of anti-in- flammatory cytokines, and increased stress hormones. Alterations in immune and endocrine function associated with depression may adversely affect health by increasing risk of infection [10] , inflammatory disorders [11] , and even malignancy [12] -[14] . The etiological relationship between depression and these physiological changes, however, is a complex one that is likely bi-directional in nature [15] . In fact, depression is known to stimulate some components of the immune system and suppress others. Furthermore, certain immune elements (such as pro-inflammatory cytokines) can lead to sickness behaviors that resemble depression, which has led to the pos- sibility that altered immune and endocrine functions may be etiologically related to depression, especially when it develops in a setting of chronic stress [16] .
Regardless of direction of effect, major depressive disorder has been associated with a host of immune [17] , endocrine [18] , and inflammatory measures [19] , including an altered balance in the T-helper (Th)1/Th2 cyto- kine ratio, i.e., higher levels of pro-inflammatory Th1 cytokines (IL-6, IL-12, IFN-γ) and monocytic cytokines (IL-1, IL-6, TNF-α) [20] -[22] and lower levels of anti-inflammatory Th2 cytokines (IL-4, IL-10) [23] [24] . Im- portantly, these altered physiological functions associated with depression have been shown to normalize in re- sponse to treatment with electroconvulsive therapy (TNF-α), [25] antidepressant drug therapy (TNF-α, C-reac- tive protein, IL-1β, IL-6) [26] -[28] , and psychological interventions (IL1-RA, IL-6, IFN-γ) [29] . Coping beha- viors and psychological interventions that increase positive emotions and neutralize negative ones may also be effective in reducing pro-inflammatory markers such as IL-6 [30] [31] , TNF-α [31] , C-reactive protein (CRP) [31] , and IFN-γ [31] , and in normalizing catecholamine [31] and cortisol levels [32] -[33] .
Religious beliefs and practices may prevent the development of depression, promote the resolution of depres- sion, and/or help persons with depressive disorder cope with the illness [34] -[36] . Given the physiological alte- rations that occur in depression, religious beliefs and behaviors may help to normalize those changes. In fact, a number of studies have reported lower levels of pro-inflammatory markers (IL-6) [37] -[39] and stress hormones (specifically cortisol) among those who are more religious [40] -[44] . Furthermore, spiritual interventions have been shown to reduce the pro-/anti-inflammatory cytokine ratios [45] [46] , reduce cortisol [47] -[51] , and de- crease catecholamine [52] [53] levels.
2. Hypotheses
The present report examines cross-sectional relationships between religious involvement, depressive symptoms, indicators of inflammation, and stress hormones in persons with major depressive disorder and chronic medical illness. We hypothesize that:
1) Pro-inflammatory cytokines (CRP, IL-6) and stress hormone (cortisol, epinephrine, norepinephrine) in our sample with major depressive disorder will be higher than in a community-dwelling sample drawn from the Midlife in the United States (MIDUS) Survey;
2) After adjustment for demographic factors and physical functioning, depressive symptom severity will be greater among those with high levels of pro-inflammatory cytokines, high pro-inflammatory/anti-inflammatory cytokine ratios, and high stress hormone levels; as a specific hypothesis, depressive symptom severity will be greater among those with high levels of CRP, IL-6, and urinary cortisol;
3) After adjustment for demographics, depressive symptoms, and physical functioning, religious practices and experiences will be lower among those with high levels of pro-inflammatory cytokines, high pro-inflammatory/ anti-inflammatory cytokine ratios, and high stress hormone levels, as a specific hypothesis, overall religiosity will be lower among those with high serum CRP, serum IL-6, and urinary cortisol.
3. Methods
3.1. Procedures
Participants aged 18 to 85 were recruited into a randomized clinical trial conducted at two sites, one in Durham, North Carolina, and the other in Los Angeles County, California. Inclusion criteria were: 1) the presence of at least one chronic medical condition for 6 months or longer; 2) a DSM-IV diagnosis of major depressive disorder; 3) mild to moderately severe depressive symptoms (10 to 40 on the Beck Depression Inventory); and 4) religion or spirituality at least somewhat important (since the clinical trial involved a religious intervention). Exclusion criteria were: 1) significant cognitive impairment; 2) receipt of psychotherapy within the past two months; 3) a diagnosis of psychotic disorder, substance abuse, or posttraumatic stress disorder (PTSD) within the past year; 4) bipolar disorder; 5) active suicidal thoughts; and 6) human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS), autoimmune diseases, dementia, endocrine disorders, a prognosis of less than 6 months, or taking immunosuppressant drugs. Study coordinators screened potential participants by telephone and then arranged a face-to-face interview, when written informed consent was obtained. Those who met the inclusion criteria were then enrolled into the study and completed a baseline evaluation. The Duke University Medical Center institutional review board (protocol #26533) and Glendale Adventist Medical Center (3/17/11) approved the study.
3.2. Measures
Physical and mental. The 12-item Duke Activity Status Index (DASI) [54] was used to measure physical func- tioning across domains for physical and instrumental activities of daily living, with higher scores indicating bet- ter physical function (range 12 - 36). Cognitive functioning was assessed by the brief Mini-Mental State Exam [55] , which measures cognitive function on a range from 0 to 18, with significant cognitive impairment defined as scores below 15. Major depressive disorder was diagnosed using corresponding modules of the Mini-Inter- national Neuropsychiatric Interview (MINI) [56] , a structured psychiatric interview that follows standard DSM- IV criteria. Depression severity was assessed using the Beck Depression Inventory (BDI-II) [57] . The BDI is a 21-item self-report depression inventory with scores ranging from 0 to 63, and is a most widely used instrument for assessing depressive symptoms in primary care settings. Diagnoses of psychotic disorder, alcohol abuse, drug abuse or PTSD within the past year, lifetime history of bipolar disorder, and presence of active suicidal thoughts were determined using corresponding modules of the MINI.
Religious. Multiple domains of RPE were assessed including denomination, self-rated religiosity, public and private religious activity, intrinsic religiosity, daily spiritual experiences, and religious coping. To determine study eligibility, a single item was used to measure self-rated religiosity/spirituality and asked, “How important is religion/spirituality in your daily life?” (those indicating “not important” were excluded). Two items from the Duke University Religiosity Index [58] were used to measure public and private religious activity. For public re- ligious activity (adapted to chronically ill disabled patients), participants were asked, “When you were physic- ally able, how often did you attend religious services or other religious meetings?” (responses ranged from “never” [1] to “more than once a week” [6] ). For private religious activity, participants were asked: “How often do you spend time in private religious activities, such as prayer, meditation or scripture study?” (responses ranged from “rarely or never” [1] to “more than once a day” [6] ). Intrinsic religious motivation was assessed using a validated 10-item intrinsic religiosity scale that has been used in both community and medical popula- tions (range 10 - 50) [59] . Religious and spiritual experiences were examined using the 16-item Daily Spiritual Experiences Scale [60] (range 16 - 94). Finally, religious coping was measured using the 14-item Brief RCOPE [61] , which is made up of two 7-item subscales assessing positive (PRC) and negative religious coping (NRC) (0 - 21 range for each subscale). We reverse scored items on the NRC subscale and combined them with PRC scores to form an overall religious coping measure (range 0 - 42). To increase the power of the analysis for the primary hypotheses, religious variables were combined into an overall religiosity measure by summing public and private religious activity, daily spiritual experiences, intrinsic religiosity, and religious coping. In order to provide equal weighting for each religiosity variable, given varied metrics, we rescaled each measure to a 0 - 1 scale and summed them to create an overall religiosity variable (range 0 to 5) (alpha = 0.88).
Demographic. Demographic variables included age, gender, race, and education. Age and education (years of education) were left as continuous, whereas gender was categorized as female (1) vs. male (0) and race as white (1) vs. non-white (Black, Hispanic, Asian) (0).
Biomarkers. Serum was collected from venous blood collected in serum separator Vaccutainer tubes. Levels of serum inflammatory markers (TNF-α, IFN-γ, IL-1β, IL-4, IL-6, IL-10, IL-12(p70)) were measured using Millipore’s multiplexed high sensitivity cytokine magnetic bead-based immunoassay kits (Milliplex cat #HSTCMAG-28SK, EMD Millipore, Billerica, MA) according to the manufacturer’s instructions. IL-1RA was run using the Milliplex Human Cytokine kit (Hcytomag-60K according to manufacturer’s instruction with the following manufacturer’s recommendations to increase sensitivity: sample volume was doubled, from 25 μl to 50 μl and an additional standard dilution was added. The mean fluorescence intensity values were then divided by two in the data analyses to adjust for greater sample volume. Minimal detectable levels in pg/ml were 0.15 for TNF-α, 0.32 for IFN-γ, 0.12 for IL-1β, 0.2 for IL-1RA, 1.24 for IL-4, 0.13 for IL-6, 0.58 for IL-10, and 0.15 for IL-12(p70). Intra- and inter -assay coefficient of variance (CV) were <6% and <20% for all cytokines, re- spectively. Plates were read on the MAGPIX (Luminex Corp., Austin, TX), and the data analyzed on MasterP- lex 2010 (Hitachi Solutions America, San Bruno, CA). All samples were run in duplicate. Samples were re- peated if the CV between the duplicates was greater than 15%. Pro-inflammatory-to-anti-inflammatory ratios were calculated from mean values.
Serum CRP was measured using enzyme-linked immunosorbent assay (ELISA) kits from Assaypro, St. Charles, MO), which had a minimal detection of ~0.25 ng/ml and an intra- and inter-assay CV of 5.0% and 7.1%, respectively. Cortisol concentrations were determined in 12-h urine samples using ELISA kits (Enzo Life Sciences International Inc., Plymouth Meeting, PA) that had a lower limit of detection of 333 pg/ml, and intra- and inter-assay CV of 10.5% and 13.4%, respectively. Cortisol levels were normalized for urine volume using creatinine levels determined by parametric kits that employ the Jaffe reaction (R&D Systems, Minneapolis, MN; minimal detection of 0.01 mg/dL and intra- and inter-assay CV of 3.5% and 4.0%, respectively). Plates were read using a plate reader (µQuant, Biotek Instruments, Inc., Winooski, VT) set at the appropriate wavelength. All samples were run in duplicate along with duplicate standards that were used to generate a standard curve. The amount of analyte in the unknown samples was calculated from the standard curve, and expressed as a mean ± SEM of the two samples in pg or mg/ml as appropriate. If the coefficient of variance between the duplicates was greater than 15%, then another aliquot of the sample was thawed, and the assay repeated for that sample. Samples were assayed in batches to further limit any potential of inter-assay variability.
Twelve-hour urinary catecholamines (epinephrine and norepinephrine) were determined by high performance liquid chromatography with coulochem detection (HPLC-CD). All samples were run in batches along with standards and quality control samples (Biorad Lyphocheck I and II; City, State). Urine samples (1.0 ml) containing 1.0 ml of phosphate buffer (pH 7.0), 1.0 ml of 1.5 M (pH 8.6) Tris buffer and 200 µl 3,4- dihydroxybenzylamine (DHBA; internal standard) were vortexed. Next, 50 mg of acid washed alumina was added, and the tubes were agitated for 15 min. The samples were centrifuged for 1 min at 200 rpm, and liquid in the samples was aspirated. The samples were washed with 1 ml nanopure H2O three times and centrifuged after each aspiration to settle the alumina. To elute the catecholamines, 200 µl of 0.1 M HClO4 were added to each sample, and the samples were vortexed. The tubes were centrifuged for 5 min at 5000 rpm, and the filtrates were collected. Samples were placed in an ESA Model 542 autosampler with a Model 582 isocratic pump with a 3 um Atlantis T3 4.5 mm × 150 mm column (Waters, Canton, Massachusetts) and detected using a CouloChem III set at +200 mV and −200 mV. The mobile phase was delivered at a flow rate of 1.0 ml/min. Urinary catecholamine were analyzed using EZChrom Elite Software (Scientific Software Inc., Pleasanton, CA). Concentrations were determined based on standards of known concentrations (200 ng/ml) of norepinephrine, epinephrine, and dopamine and expressed as μg/g creatinine to normalize for urine volume.
3.3. MIDUS Sample
In order to determine how the ranges of biomarkers in the present study compared to those obtained in a com- munity-dwelling sample of largely healthy individuals, we identified a subsample of 572 participants from Project 4 of the MIDUS II biomarker study [62] . These participants were drawn from the overall dataset (n = 2500) by including those with at least one chronic medical illness who had available biomarker data, but excluding twins and other family members to avoid participants who were related. Participants completed the 20-item Center for Epidemiological Studies Depression (CES-D), whose scores were available for comparison with participants in our study (assessed by BDI). The presence of chronic medical illness was determined by the presence of self-reported or physician-diagnosed circulation problems, blood clots, heart murmur, TIA/stroke, anemia or other blood disease, hypercholesterolemia, diabetes, asthma, chronic obstructive pulmonary disease or emphysema, tuberculosis, thyroid disease, peptic ulcer disease, cancer, arthritis, glaucoma, cirrhosis/liver dis- ease, or other medical condition. The study was approved by the Institutional Review Board at each of the par- ticipating MIDUS centers, and informed written consent was obtained for all participants [63] . Biological sam- ples were collected from participants of MIDUS II after they had completed the telephone and mail surveys during a 2-day clinic visit. A 12-h urine sample was collected the day before the first clinic visit, and fasting blood samples were drawn in the morning of day 2 (see Dienberg et al. for more details [64] ).
Serum IL-6 levels from the MIDUS II study were determined using Quantikine® high-sensitivity ELISA kits (R&D Systems, Minneapolis, MN) with a reference range of 0.45 - 9.96 pg/ml, and intra- and inter-assay CV of 4.1% and 13%, respectively [65] . Serum CRP levels were measured using the Behring Nephelometer II (BNII) Analyzer (Dade Behring) using an immunonepholometric assay specific for CRP (Seimens Healthcare Diagnos- tics, Deer-242 field, IL), which provides equivalent results to the high-sensitivity ELISA kits for CRP [66] . Uri- nary cortisol levels were measured using an enzymatic colorimetric assay and liquid chromatography-tandem spectrometry (LC-MS/MS). Samples were spiked with 3H-cortisol as an internal standard. After extraction by on-line turbulent flow HPLC, samples were analyzed using LC-MS/MS with multiple-reaction monitoring in positive mode. Inter-assay CV was 5.7% - 8.8%, intra-assay CV was 4.7% - 5.0%, and the reference range was 29.1 - 63.4 ng/ml [67] .
3.4. Statistical Analysis
Descriptive statistics were used to describe the characteristics of the sample and visually compare the median and range values of inflammatory markers and stress hormones in our sample with those obtained in the MIDUS study (Table 1). Because inflammatory markers and stress hormone levels in our sample were non-normally distributed (and attempts to transform the data were not successful in achieving normality), we dichotomized levels into the bottom two-thirds (0) and top one-third (1) for analyses; values were excluded from analyses if they exceeded five standard deviations above the median. In bivariate analyses, the average level of depressive symptoms, each religious measure and overall religiosity were compared between participants with low and high levels of cytokines and stress hormones using the Student T-test (Table 2). Logistic regression was used to ex- amine the relationship between depressive symptoms and biomarkers, controlling for demographic factors and physical functioning (DASI) (Table 3). Logistic regression was also used to examine relationships between reli- gious measures and biomarkers, controlling for demographic factors, physical functioning, and depressive symptoms (Table 4). Given the highly exploratory nature of these analyses and relatively small sample size
Table 1. Characteristics of current sample compared to MIDUS dataset.
MIDUS = Midlife in the United States Study, CES-D = Center for Epidemiologic Studies Depression, CRP = C-reactive protein, TNF-α = tumor necrosis factor-α, IL = interleukin, IFN = interferon, CORT = cortisol, EPI = epinephrine, NE = norepinephrine, DA = dopamine, “---“ = no comparison available, 1n = 122 - 132, 2n = 561 - 570, 3C-reactive protein is a “positive” acute phase protein.
limiting power, p values were not reduced as would ordinarily be appropriate for the multiple comparisons made here. Statistical analyses were done using SAS (version 9.3, SAS Institute Inc., Cary, North Carolina).
4. Results
Participants. Between June 2011 and June 2013, a total of 450 participants were screened by telephone for eligi- bility. Of those, 187 underwent in-person screening and 132 were enrolled in the study. Three subjects were in- cluded who did not fulfill inclusion/exclusion criteria but were randomized in the trial, and these were included
Table 2. Comparison of depressive symptoms and religious indicators across high vs. low levels of inflammatory markers and stress hormones.
#0.05 < p < 0.10, *p ≤ 0.05, **p ≤ 0.01 (student t-test). High = top one-third of biomarker variable (n = 43 - 45); Low = bottom two-thirds of biomarker variable (n = 82 - 87). CRP=C-reactive protein, TNF-α = tumor necrosis factor-α, IL = interleukin, IFN = interferon, CORT = cortisol, EPI = epinephrine, NE = norepinephrine.
Table 3. Logistic regression models examining depressive symptoms (Beck Depression Inventory) as a predictor of high levels of inflammation and stress hormones (dependent variables).
OR = odds ratio, CI = confidence interval, Biomarkers dichotomized: high = top 1/3, low = bottom 2/3, *p ≤ 0.05; n = 122 - 129, Controlled for age, education, gender, race, physical functioning (Duke Activity Status Index).
Table 4. Logistic regression models examining religious characteristics as predictors of high levels of inflammatory markers and stress hormones (dependent variables).
Controlling for demographics (age, education, gender, race), physical function (Duke Activity Status Index), depressive symptoms (BDI), OR = odds ratio, CI = confidence interval (biomarkers dichotomized: high = top 1/3 and low = bottom 2/3), #0.05 < p < 0.10, *p ≤ 0.05, **p ≤ 0.01, n = 122 - 129.
in the analyses. The average age of participants was 51.6 years (range 24 - 84), average education level was 15.1 years (range 4 - 31), and the majority were female (68.9%) and white (53.0%) (Table 1). The mean BDI score was 25.3 ± 8.5, confirming that our participants were moderately depressed. Cytokine profiling in serum sam- ples revealed detectable levels of the acute phase protein, CRP, pro-inflammatory cytokines (TNF-α, IL-1β, IFN-g, IL-6, and IL-12), and anti-inflammatory cytokines (IL-4, IL-1RA, IL-10). Pro-inflammatory-to-anti-in- flammatory ratios were highest for IFN-g/IL-10 (1.3), TNF-α/IL-10 (1.1) and IFN-g/IL-4 (0.8), and lowest for IL-1β/IL-1RA (0.1), IL-1β/IL-4 (0.1), IL-1β/IL-10 (0.2), and IL-12/IL-4 (0.2).
Comparison to MIDUS. When MIDUS participants were compared to the present sample, the former were older (59.2 vs. 51.6 years, p < 0.05), less likely to be female (53.7% vs. 68.9%, p < 0.05), more likely to be white (92.0% vs. 53.0%, p < 0.05), but had similar education (14.8 vs. 15.1, p = 0.34) (Table 1). The average CES-D score (9.0, SD 10.3) was considerably below the usual cut-off of 16 or higher for significant depression, indicating that participants in the MIDUS study were largely not depressed. We did not statistically compare biomarkers due to their non-normal distributions. Instead, we provide the median value and range for the two pro-inflammatory markers (CRP, IL-6) and three stress hormones (cortisol, epinephrine, norepinephrine) com- mon to both samples. Although levels of IL-6 were similar between samples, the median CRP in our participants with major depressive disorder was nearly three times that of MIDUS participants (4.4 vs. 1.3 ng/ml). Likewise, median urinary cortisol level was three times higher (30.2 vs. 11.0 mg/L creatinine) in our depressed sample compared to MIDUS participants, urinary epinephrine was almost three times higher (5.4 vs. 1.8 mg/L), and median urine norepinephrine level was about 50% higher (39.5 vs. 26.0 mg/L). Given known relationships be- tween depressive disorder, inflammatory markers, and stress hormones, these findings are expected and support the accurate determination of biomarker levels reported in the present study.
Biomarkers and depressive symptoms. In the present sample, we hypothesized that greater depression severity would be related to higher levels of pro-inflammatory biomarkers, lower anti-inflammatory cytokines, higher pro-/anti-inflammatory cytokine ratios, and higher urinary stress hormones. Table 2 describes the average level of depressive symptoms across low and high levels of biomarkers, where high represents the top one-third of values. Table 3 illustrates these relationships in multivariate analyses controlling for demographics and physical functioning. Neither bivariate nor multivariate analyses supported our hypothesis. If anything, the opposite was true. In bivariate analyses, severity of depressive symptoms tended to be lower among those with high levels of the pro-inflammatory cytokines TNF-α (22.9 vs. 26.5 pg/ml, t = −2.29, p = 0.02) and IL-1β (23.4 vs. 26.3 pg/ml, t = −1.78, p = 0.08) and in those with high levels of urinary cortisol (23.3 vs. 26.2 mg/L creatinine, t = −1.84, p = 0.07). Multivariate models controlling for demographics and physical functioning (Table 3) largely confirmed these findings.
Biomarkers and religious characteristics. We hypothesized that religious activities and attitudes would be lower among those with 1) high levels of pro-inflammatory markers, 2) low anti-inflammatory cytokines, 3) high pro-/anti-inflammatory cytokine ratios, and 4) high stress hormones. Table 2 presents the average level of religious involvement across low and high levels of biomarkers, and Table 4 shows these relationships con- trolled for demographics, depressive symptoms, and physical functioning. Again, neither bivariate nor multiva- riate analyses supported our hypotheses (including our specific hypothesis), and most associations were non- significant. While there were some significant and near-significant values, no discernible pattern in one direction or the other was identified (see Figures 1-3).
Bivariate analyses. For religious attendance, while frequency of attendance as expected tended to be lower in those with high levels of TNF-α (3.5 vs. 4.0 pg/ml, t = −1.75, p = 0.08) and high levels of the pro-/anti-inflam- matory ratio TNF-α to IL-10 (3.4 vs. 4.1 pg/ml, t = −2.14, p = 0.03), attendance tended to be higher in those with high levels of the pro-/anti-inflammatory ratio IL-1β to IL-1RA (4.2 vs. 3.7, t = 1.79, p = 0.08). Likewise, private religious activities such as prayer and scripture reading were more frequent among those with high levels of the pro-/anti-inflammatory ratio IL-6 to IL-4 (4.1 vs. 3.4, t = 2.31, p = 0.02), high levels of the pro-/anti-in- flammatory ratio IFN-γ to IL-4 (4.0 vs. 3.4, t = 2.0, p = 0.04), and high levels of the pro-/anti-inflammatory ratio
Figure 1. C-reactive protein (median, interquartile range) by overall religiosity (ρ = 0.0004, p = ns, n = 129).
Figure 2. Interleukin-6 (median, interquartile range) by overall religiosity (ρ = −0.015, p = ns, n = 126).
Figure 3. Urine cortisol (median, interquartile range) by overall religiosity (ρ = 0.054, p = ns, n = 127).
IL-12 to IL4 (4.1 vs. 3.4, t = 2.1, p = 0.04). Similarly, intrinsic religiosity tended to be higher in those with high levels of the pro-inflammatory marker IFN-γ (36.7 vs. 34.1, t = 1.68, p = 0.09) and high levels of the pro-/anti- inflammatory ratio IFN-γ to IL-10 (37.1 vs. 34.0, t = 2.07, p = 0.04). One of the few findings consistent with our hypothesis, intrinsic religiosity was higher in those with high levels of the anti-inflammatory cytokine IL-1RA (37.5 vs. 33.5, t = 2.6, p = 0.01). Inconsistent with our hypothesis, however, daily religious and spiritual expe- riences tended to be more common among those with high levels of the pro-inflammatory IFN-γ (61.6 vs. 55.7, t = 1.97, p = 0.05) and high levels of the pro-/anti-inflammatory ratio IFN-γ to IL-10 (61.5 vs. 55.9, t = 1.85, p = 0.07). Religious coping was also higher among those with high levels of the pro-inflammatory IFN-γ (31.3 vs. 28.5, t = 2.50, p = 0.01), high levels of the pro-/anti-inflammatory ratio IFN-γ to IL-10 (31.9 vs. 28.3, t = 3.15, p = 0.002), and high levels of the pro-/anti-inflammatory ratio IL-12 to IL-10 (31.0 vs. 28.7, t = 1.91, p = 0.06). No significant associations were found with overall religiosity, although religiosity tended to be higher in those with high levels of IFN-γ (pro-inflammatory), IL-1RA (anti-inflammatory), and IFN-γ to IL-10 ratio (pro-in- flammatory).
Multivariate analyses. Multivariate analyses (Table 4) essentially confirmed bivariate analyses, with the vast majority of associations non-significant and a mixed pattern of expected and unexpected findings. Religious at- tendance tended to be lower in those with high levels of the pro-inflammatory TNF-α (OR = 0.80, 95% CI = 0.63 - 1.02) and high TNF-α to IL-10 ratio (OR = 0.75, 95% CI = 0.58 - 0.97), but was higher in those with high pro-/anti-inflammatory ratio IL-1β to IL-1RA (OR = 1.26, 95% CI = 0.98-1.64). Private prayer and scripture reading were also higher in those with high pro-/anti-inflammatory ratio IL-6 to IL-4 (OR = 1.30, 95% CI = 1.02 - 1.67) and high pro-/anti-inflammatory ratio IL-12 to IL-4 (OR = 1.26, 95% CI = 0.99 - 1.61). Likewise, intrinsic religiosity was higher among those with high levels of pro-inflammatory cytokines IL-1β (OR = 1.05, 95% CI = 1.00 - 1.11), IFN-γ (OR = 1.06, 95% CI = 1.01 - 1.12), and high pro-/anti-inflammatory ratio IFN-γ to IL-10 (OR = 1.05, 95% CI = 1.00 - 1.11), all contrary to our hypothesis. As expected, though, intrinsic religios- ity was higher among those with high levels of the anti-inflammatory cytokine IL-1RA (OR = 1.07, 95% CI = 1.02 - 1.13) and was lower in those with high levels of urinary norepinephrine (OR = 0.94, 95% CI = 0.89 - 0.99). But, contrary to expectations, daily spiritual experiences were significantly higher among those with high levels of IFN-γ (OR = 1.04, 95% CI = 1.01 - 1.07), and religious coping was more common among those with high levels of IFN-γ (OR = 1.11, 95% CI = 1.03 - 1.19), high levels of IFN-γ to IL-10 ratio (OR = 1.11, 95% CI = 1.03 - 1.19), and high levels of IL-12 to IL-10 ratio (OR = 1.06, 95% CI = 0.99 - 1.14). Similarly, overall reli- giosity tended to be higher among those with high levels of IFN-γ (OR = 1.54, 95% CI = 1.01 - 2.34) and high levels of the IFN-γ to IL-10 ratio (OR = 1.47, 95% CI = 0.97 - 2.24). Thus, if any consistent pattern was found in these relationships, religious activities and attitudes were related to higher levels of pro-inflammatory markers, the exact opposite of our hypothesis.
5. Discussion
This is one of the first studies to provide a detailed examination of associations between religious practices/ex- periences and a wide range of stress biomarkers in persons with major depressive disorder and chronic medical illness, while controlling for demographics, depressive symptoms and physical functioning. As expected, inflammatory markers and stress hormones were higher in the present sample than in a community sample of relatively healthy middle-aged adults. Contrary to expectations, depressive symptoms were largely unrelated to stress biomarkers in the current sample, and when related, were actually lower among those with higher levels of pro-inflammatory cytokines (TNF-α, IL-1β). This contrasts with what is reported in most of the literature, where those with more severe depression are usually found to have higher levels of pro-inflammatory cytokines, espe- cially IL-6, TNF-α and IFN-γ [67] [68] . All participants in the current study, however, had major depressive disorder to start with and the range of depressive symptoms was limited due to its mild to moderate severity (BDI 10 to 42). The presence of co-morbid chronic medical illness may also have added to the variability seen in biomarkers found here.
Likewise, religious practices and experiences (RPE) were largely unrelated to these stress biomarkers, and there was no evidence in favor of our primary hypothesis that overall religiosity would be related to lower levels of CRP, IL-6, or urinary cortisol (Figures 1-3). In fact, no clearly discernible pattern of association was found between RPE and inflammatory cytokine or stress hormone markers among the 138 analyses done here. Consis- tent with our hypothesis, RPE were higher among those with high levels of the anti-inflammatory cytokine IL-1RA and lower in those with high levels of urinary norepinephrine. However, RPE were also higher in those with high pro-inflammatory cytokines, especially IFN-γ, IFN-γ/IL-4, and IFN-γ/IL-10, which is contrary to our hypothesis. To our knowledge, this is the first study to report a link between religious activity and an anti-in- flammatory cytokine (IL-1RA). Others, however, have reported an increase in IFN-γ with spiritual practices [69] -[71] , which may result from the complex relationship between IFN-γ, depressive symptoms, and stress. On the one hand, IFN-γ increases in major depression [72] [73] and decreases in response to treatment [74] [75] , at least among studies in men. On the other hand, IFN-γ is suppressed by cortisol, decreases during acute psychological stress [24] , and may increase in response to treatment, especially in depressed women who may have lower IFN-γ levels than depressed men [74] [76] .
This is also to our knowledge the first report of a relationship between greater intrinsic religiosity and lower urinary norepinephrine, a stress hormone known to increase in depression [77] and be associated with a worse prognosis in cardiovascular diseases [78] . In a 32-year prospective study that compared 144 Catholic sisters with 138 healthy laywomen from a community in Italy, investigators reported that urinary norepinephrine was sig- nificantly lower in sisters at every six-year evaluation throughout the study period [79] . Furthermore, at least two studies have found that spiritual interventions (meditation) lowered serum or urinary norepinephrine in healthy adults [80] and in those with congestive heart failure [81] .
6. Limitations
Numerous limitations exist that should be taken into account when generalizing or interpreting the findings re- ported here. First, this was a relatively small sample and a sample of convenience (volunteers), making it difficult to generalize results to other populations. Second, the sample was made up of people with chronic medical illnesses taking an assortment of medications that could influence levels of inflammatory markers and stress hormones and make it difficult to detect more subtle associations with either depressive symptoms or RPE. Third, all participants in this study had major depressive disorder, and most median biomarker levels were three times that greater than in the comparison community population, thus ceiling effects may have been an issue. Fourth, as noted above, given the multiple statistical comparisons made here, at least 1 in 20 comparisons would be expected to be significant based on chance alone, increasing the likelihood that significant findings reported here were due to Type I error. Nevertheless, the study also has several strengths. Participants were drawn from two different sites, the East coast and the West coast, increasing the likelihood of generalizability; analyses were carefully controlled for demographics and physical health; and this is the first study to examine such associa- tions in persons with comorbid major depressive disorder and chronic medical illness. Furthermore, the breadth of religious characteristics assessed in this study is unparalleled in the literature for studies examining relation- ships between religiosity and biomarkers [82] .
7. Conclusion
Although the median and range inflammatory markers and stress hormones measured here were higher than those from a healthy non-depressed community sample as expected, we found little evidence for a relationship between depressive symptoms and these biomarkers in persons with major depressive disorder and chronic ill- ness. Furthermore, we found no consistent pattern between religious practices or experiences and either inflam- matory markers or stress hormones in these cross-sectional analyses (with the exception of higher levels of the anti-inflammatory cytokine IL-1RA, higher pro-inflammatory cytokine IFN-γ, and lower urinary norepineph- rine). Future research, particularly the ongoing clinical trial that sample is the baseline for, is needed to help de- termine whether religious interventions for major depression among those with chronic medical illness can help to reverse the high levels of pro-inflammatory markers and stress hormone levels associated with depression.
Funding
Funding support provided by the John Templeton Foundation.
NOTES
*Corresponding author.