Data are available via figshare: <https://figshare.com/s/d1e437c57838f5e63eee>.

Introduction {#sec005}
============

An ever-growing body of research suggests that early-life adversity (ELA) can program biological systems, which predispose individuals to later-life physical and mental-health problems \[[@pone.0221310.ref001], [@pone.0221310.ref002]\]. Empirical evidence exist for associations between ELA and elevated risk of depression, cardiovascular disease, diabetes, autoimmune diseases and cancer, to name a few (reviewed in \[[@pone.0221310.ref003]\]). Despite the salient role of ELA on disease risk, the biological mechanisms that play a downstream role in increased disease susceptibility are not well understood.

Mechanistic research on the biological embedding of ELA has emphasized maladaptive programming of the hypothalamic-pituitary-adrenal (HPA) axis with the associate release of cortisol through processes of allostasis \[[@pone.0221310.ref003], [@pone.0221310.ref004]\]. Specifically, studies have documented a shift in HPA axis function with hyper- or hypo-secretion of cortisol in depression and post-traumatic stress disorder (PTSD), respectively \[[@pone.0221310.ref005], [@pone.0221310.ref006]\]. Similar findings have been reported in individuals exposed to ELA without such diagnoses \[[@pone.0221310.ref007], [@pone.0221310.ref008]\]. This programming, in turn, can result in mitochondrial dysfunction, failure to down-regulate the inflammatory response and overall metabolic stress, thereby increasing circulatory levels of lipids, glucose, oxidants, and pro-inflammatory cytokines \[[@pone.0221310.ref009], [@pone.0221310.ref010]\]. Further mechanistic research on the biological embedding of ELA suggests physiological dysregulation may be mediated at the genetic level via epigenetic modifications that can persist over long periods of time \[[@pone.0221310.ref011]\], including evidence linking ELA and cortisol responses via methylation levels in the glucocorticoid receptor (*NR3C1*) gene \[[@pone.0221310.ref012], [@pone.0221310.ref013]\]. Other research suggest the involvement of telomere biology in mediating the longer-term link between ELA and disease risk \[[@pone.0221310.ref014]\]. What is less clear, however, is how target immune cells respond to stress *in vivo* as a consequence of ELA, via rapid gene expression regulation \[[@pone.0221310.ref015], [@pone.0221310.ref016]\]. This new knowledge can provide insights into an integrated and dynamic cellular regulatory system whose signal profiles could forecast disease risk associated with early adversity \[[@pone.0221310.ref017]--[@pone.0221310.ref019]\].

Cells show remarkable flexibility in response to stimuli by regulating gene expression in a transient manner \[[@pone.0221310.ref016], [@pone.0221310.ref020]\]. In one of the first studies relating peripheral blood mononuclear cells (PBMC) gene expression to trauma, basal gene expression signatures, both immediately following trauma and four months later, distinguished survivors who met diagnostic criteria for PTSD from those who did not \[[@pone.0221310.ref021]\]. Follow-up studies provided further support for associations between chronic stress and basal expression levels of genes related to glucocorticoid signaling and pro- and anti-inflammatory pathways \[[@pone.0221310.ref022]--[@pone.0221310.ref024]\]. Importantly, several studies have also provided evidence of rapid (e.g., from 30 minutes to 8 hours) gene expression activation in response to *in vitro* stimulation \[[@pone.0221310.ref025]\], psychological stress \[[@pone.0221310.ref026], [@pone.0221310.ref027]\], physical stress \[[@pone.0221310.ref028]\] and stress-reduction methods \[[@pone.0221310.ref029]\]. Notably, these response patterns were recently dubbed the "conserved transcriptional response to adversity" \[[@pone.0221310.ref030]\]. Taken together, theory and evidence suggests that programmed immune cells of individuals exposed to early adversity may show compromised adaptation in response to acute stress, which, if repeated, may play a downstream role in disease risk.

For example, the glucocorticoid-immune signaling pathway has been implicated as a key mechanism in relation to chronic stress (i.e., caregiving, poverty \[[@pone.0221310.ref023], [@pone.0221310.ref031]\]), through reduced glucocorticoid receptor (*NR3C1*) availability, ligand binding affinity, and functional capacity to regulate gene expression. Specifically, chronic stress, via extended exposure to cortisol, is associated with reduced *NR3C1* expression, leading to glucocorticoid resistance and impaired negative feedback inhibition of the HPA axis \[[@pone.0221310.ref032]\]. Reduced levels of glucocorticoid receptors, in turn, bind less cortisol. This effectively decreases the number of ligand-bound receptor complexes available to translocate to the nucleus and regulate the expression of genes, including anti-inflammatory genes. Thus, reduced levels of *NR3C1* expression can lead to impaired immune function \[[@pone.0221310.ref033]\]. Further, evidence exists for elevated levels of pro-inflammatory cytokines in children and adults exposed to ELA \[[@pone.0221310.ref034], [@pone.0221310.ref035]\]. Dysregulation of the immune system in the context of ELA, as well as the resulting increase of pro-inflammatory cytokines, can increase risk for a host of diseases, from autoimmune to atherosclerosis and cancer \[[@pone.0221310.ref036]\].

Here, we delineate a program of research to study ELA-related programming of biological systems with emphasis on processes related to glucocorticoid signaling and inflammation. Experimental studies usually employ either a between- or within-person designs \[[@pone.0221310.ref037]\], and measure physiological or gene expression changes over relatively short time frames. Our novel program of research uses a within-person, between-groups experimental design, and aims to test whether ELA leads to dysregulation of physiological, gene expression and pro-inflammatory cytokines in response to a canonical laboratory stressor, compared with a resting control condition, and compared with individuals without exposure to ELA. The within-person control condition enables us to disentangle the effects of acute stress from noisy gene expression measurements in the same individuals. The between-person acute stress condition enables tests of potential programming of biological systems between ELA and control individuals, which theory suggests in more likely to manifest during times of duress \[[@pone.0221310.ref026], [@pone.0221310.ref027]\]. Measuring stress-induced physiological, immune, and gene expression changes enable tests across a wide array of biological and cellular pathways that may play a downstream role in disease susceptibility. This conceptual framework integrates environmental, psychological, and biological elements of what Cohen et al., termed 'a stage model' of stress and disease \[[@pone.0221310.ref038]\]. Finally, to substantiate group differences between ELA and control individuals, we compared individuals who reported at least three major traumatic events up to age 18 years \[[@pone.0221310.ref039]\], while the control group were selected only if they reported none.

Specifically, in this pilot feasibility study, we applied a validated screening instrument \[[@pone.0221310.ref040]\] to recruit 12 men, 6 of whom who endorsed at least three significant adverse events ('ELA group') \[[@pone.0221310.ref039]\], and 6 who confirmed zero ('controls'). In a randomized within-subjects, between-groups experimental design, we induced acute stress in the lab (Trier Social Stress Test, TSST) and included a no-stress control condition separated by one week. During these sessions, we obtained repeated measurements of physiological reactivity, plasma levels of pro-inflammatory cytokines, and PBMC gene expression over a 4-hour window post-test. Primary analyses aimed to replicate previously known associations between physiological measures (i.e., blood-pressure and salivary cortisol) in response to acute laboratory stress. Our secondary analyses examined how young adults' exposure to ELA influence neuroendocrine and inflammatory responses to acute stress compared with non-exposed individuals by testing physiological, *NR3C1* gene expression, and pro-inflammatory cytokines responses. We focused on four pro-inflammatory cytokines, including interleukin-1β (IL-1β**)**, interleukin-6 (IL-6), interleukin-8 (IL-8), and tumor necrosis factor alpha (TNF-α). Based on prior literature, we hypothesized that individuals exposed to ELA will evince dysregulated physiological, *NR3C1* gene expression, and pro-inflammatory changes to an acute laboratory stressor, a response pattern that may reveal dynamic biological signatures of early-life programming with implications for life-long health.

Materials and methods {#sec006}
=====================

Participants {#sec007}
------------

Participants were healthy male college students at the Pennsylvania State University, recruited by word of mouth and advertisements on campus bulletin boards. We focused on men in this exploratory study due to known sex differences in the stress response \[[@pone.0221310.ref041]\] and the small sample size for stratified analyses. To obtain the sample who were exposed to ELA, a trained clinical interviewer conducted a phone interview to screen over 100 eligible men using the Stressful Life Events Screening Questionnaire (SLESQ) \[[@pone.0221310.ref040]\], a 13-item self-report measure that assesses lifetime exposure to traumatic events. We asked respondents 11 specific and two general categories of events, such as death of a parent or sibling, life-threatening accident, and sexual and physical abuse. Based on evidence that three or more traumatic events confers higher risk for disease \[[@pone.0221310.ref039]\], and considering the severity of the traumatic events, participants who responded to at least 3 incidents up to age 18 years (independently reviewed and reached consensus by MZ and IS) were invited to participate in the ELA group. Respondents' examples for adverse exposures in this study included (unsubstantiated) child abuse and neglect, severe violence exposure, parental loss, suicide of a close friend or a family member, severe illness of an immediate family member or car accidents. In addition, the SLESQ was used to screen participants without a history of traumatic exposures to serve as the control group. Selection criteria stipulated that subjects were between 18--25 years, without current medical illness or endocrine illness (for example, asthma, diabetes, thyroid disease or pituitary gland disorders confirmed by self-report and physical examination), were currently non-smokers and were not using medication on a regular basis, including psychiatric medication. The final sample included 12 men, 6 of whom experienced early adversity (i.e., 'ELA group') and 6 who did not (i.e., 'controls') (mean age = 21.25, SD = 2.3). Demographics of the sample are presented in **[Table 2](#pone.0221310.t002){ref-type="table"}**. The study was approved by the Ethics Committee at the Pennsylvania State University, registered at ClinicalTrials.gov (Identifier: NCT03637751), and all participants provided written informed consent. Participants received a modest monetary incentive for participation.

General procedure {#sec008}
-----------------

Testing was carried out at the Pennsylvania State's Clinical Research Center (CRC). Participants made two visits to the CRC during weekdays, one week apart, on the same day. Testing was scheduled to begin at 11:00am and end by 4:15pm. We used a randomized counter-balanced order for the two sessions (i.e., TSST and no-stress control conditions) blind to participants and lab personnel. Lab personnel were also blind to group status. Participants were given specific instructions to refrain from excessive physical activity on the day of the testing, consuming alcohol for 12 hours before their arrival, and eating and drinking (besides water) for 2 hours prior to the testing session. After arrival and consent, trained nurses completed a physical examination and inserted an IV catheter into the antecubital vein 30 minutes after arrival (30 minutes prior to testing). The TSST session was scheduled to begin at 12:00pm to minimize the effects of circadian changes in cortisol, and was carried out as described previously \[[@pone.0221310.ref042]\]. Briefly, the TSST consists of a free speech and a mental arithmetic task of 10 minutes duration performed in front of a panel of two committee members (mixed gender) with a camera and microphone situated between the interviewers. Participants were told that they would play the role of an interviewee for a job and have 5 minutes to make an argument for their candidacy. After 5 minutes, the second task emphasizing cognitive load commenced. In this task, participants were asked to count backwards from 1,687 in multiples of 13. If a mistake was made, they were instructed to start again from the beginning. In the no-stress control condition, participants were instructed to sit in a room, read magazines, and to refrain from any stressful activities (e.g., cell-phone use was restricted). After the second blood draw, approximately 60 minutes after the TSST session and 90 minutes after the first baseline measure in the no-stress control condition, participants were administered a set of questionnaires. These questionnaires were administered in both sessions and the average score was calculated before analyses (see below for details). Considering the long time-frame of the study and the repeated collection of multiple blood samples, a standardized low-calorie meal was provided after the third blood draw (approximately at 1:45pm). [Fig 1](#pone.0221310.g001){ref-type="fig"} outlines the study design.

![Study design for both sessions (TSST and no-stress control condition), separated by one week.](pone.0221310.g001){#pone.0221310.g001}

Primary endpoints: Physiological reactivity {#sec009}
-------------------------------------------

Salivary cortisol was repeatedly assessed from 7 saliva samples at the following time-points: 30 minutes after arrival (30 minutes prior to testing), 1 minute prior to testing, immediately after testing (15 minutes after last sample in the control condition), and 15, 30, 60 and 90 minutes post-test. Saliva samples were kept at room temperature throughout the session, were immediately centrifuged at the end of the session at 3000 rpm at 24°C for 15 minutes, and then stored at -80°C until assayed. Systolic and diastolic blood-pressure were measured at the same time points as salivary cortisol.

Salivette swabs (Sarstedt, Germany) were used to collect saliva. Salivary cortisol was assessed, in duplicate, through an enzyme immunoassay protocol (Salimetrics) with known controls. The lower detection limit of the assay is \<0.007 ug/dL. Intra-assay CV was 9.88% across all samples and inter-assay CV 5.79% across four plates. Participants' blood pressure was measured, while seated, using an automatic monitor (Omron HEM-712C).

Secondary endpoint: RNA extraction and gene expression assays {#sec010}
-------------------------------------------------------------

Gene expression changes were measured repeatedly from the four blood samples at each session at the following time-points: 30 minutes after arrival (30 minutes prior to testing), and at 30 (75 minutes after the first sample in the no-stress condition), 90 and 240 minutes post-test ([Fig 1](#pone.0221310.g001){ref-type="fig"}). Given known changes in immune cell redistribution and composition in response to acute stress \[[@pone.0221310.ref020]\], complete blood count with differential was measured within 24 hours by Quest Diagnostics using additional 4 ml EDTA collection tubes.

Whole blood samples were collected in 10 mL EDTA blood tubes via an IV catheter into the antecubital vein, and immediately centrifuged for 10 minutes at 1500g prior to collection of plasma. PBMCs were immediately isolated through density-gradient centrifugation using Ficoll. Immediately following isolation, cells were suspended in RNAlater solution (Ambion) before being stored at 4°C overnight. The duration from blood sampling to stabilization of RNA never exceeded 55 minutes. RNA extraction and cDNA synthesis were performed the following day using QIAamp RNA Blood Mini Kit and cDNA Synthesis Kit respectively (Qiagen), and then stored at -80°C until assayed. RNA purity was verified using Nanodrop 2000 spectrophotometer (Thermo Scientific).

All assays were performed on a real-time PCR (Rotor Gene Q, Qiagen). PCR reactions were set-up using the complementary QIAgility robotic pipettor (Qiagen) to ensure maximum pipetting accuracy. Samples were assayed in duplicate. All repeated, within-subject samples were run on the same plate. The reaction mix for gene expression assays consists of 5 uL TaqMan Gene Expression Master Mix (Thermo Fisher Scientific), 1x TaqMan gene expression primer, UltraPure Water (Rockland), and 100ng DNA in a 10 uL reaction. The cycling profile consists of an initial denaturing at 95°C for 15 seconds and annealing/extending at 60°C for 1 minute followed by fluorescence reading, 55 cycles. One hypothesis-driven gene (*NR3C1*: Hs00353740_m1) was normalized to a housekeeping gene (*GADD45A*: Hs00169255_m1). Expression of the *NR3C1* and the housekeeping genes were assessed on the same plate in two independent PCR reactions using cDNA from the same sample aliquot.

Sample normalization was done using the ΔΔCt method \[[@pone.0221310.ref043]\]. Briefly, a cycle threshold (Ct) is defined as the cycle number at which a sample's fluorescence reaches a defined threshold. The same threshold was used for reactions assessing housekeeping and *NR3C1* genes. Thus, each sample on a given plate has two Ct values (e.g. Ct~NR3C1~ and Ct~GADD45A~). The ΔCt is calculated as the difference between the Ct of the gene of interest and the Ct of the housekeeping gene (e.g. ΔCt = Ct ~NR3C1~ -Ct~GADD45A~). The ΔΔCt represents the within-subject normalization of the three post-test samples to expression levels at baseline. That is, ΔΔCt = ΔCt~POST-TEST~---ΔCt~BASELINE~. Thus, the ΔΔCt for the baseline sample for each session is always equal to zero. Lastly, fold change is calculated by exponentiating 2 by -ΔΔCt (i.e. Fold Change = 2^-ΔΔCt^). It follows that the fold change for each baseline sample is always equal to one (i.e. 2^0^).

Secondary endpoint: Pro-inflammatory cytokines {#sec011}
----------------------------------------------

Inflammatory assays were performed on plasma isolated from whole blood. Plasma samples were stored at -80°C prior to use. Plasma levels of IL-1β, IL-6, IL-8, and TNF-α were quantified using Meso Scale Discovery's Multi-Array technology (MSD, V-PLEX Human Proinflammatory Panel II) and analyzed on a Meso QuickPlex SQ 120 instrument (Meso Scale Discovery, Rockville, MD, USA). Sample concentrations were determined relative to standard curves generated by fitting electrochemiluminenscent signal from stock calibrators with known concentrations using MSD Discovery Workbench® software. Samples were run in duplicate. Intra-assay variability was 8.02% across all samples and inter-assay variability was 3.87% across the three plates. The lower limits of detection for inflammatory markers were 0.646 pg/mL (IL-1β), 0.633 pg/mL (IL-6), 0.591 pg/mL (IL-8), and 0.690 pg/mL (TNF-α). Samples with concentrations below the curve fit range were assigned a value of 0 for analyses considering those analytes. This occurred for 13 samples (13.5%) for IL-1β. Samples for all other analytes were within detection ranges. Summary information for inflammatory cytokines (pg/mL) by time, session, and group status is presented in [S1 Table](#pone.0221310.s001){ref-type="supplementary-material"}.

Self-reported measures and other covariates {#sec012}
-------------------------------------------

We administered several questionnaires to assess levels of adverse exposures in the past 12 months and mental health symptoms. Specifically, participants completed the following questionnaires at both sessions; the Life-Event Stress Scale (LESS) \[[@pone.0221310.ref044]\], which consists of 42 common events associated with some degree of disruption of an individual\'s life and provide a standardized measure of the impact of a wide range of common stressors in the past year; and the Life Events Questionnaire (LEQ) \[[@pone.0221310.ref045]\], an 82-item inventory-type questionnaire for the measurement of life changes during the past year. The LEQ consists of items that are designed primarily for use with students. We further assessed levels of anxiety and depressive symptoms using the Beck anxiety inventory \[[@pone.0221310.ref046]\], Beck depression inventory \[[@pone.0221310.ref047]\], and State-Trait Anxiety Inventory \[[@pone.0221310.ref048]\], as well as perceived stress levels using the 10-item Perceived Stress Scale \[[@pone.0221310.ref049]\].

As noted above, given gene expression changes may depend on specific cell populations \[[@pone.0221310.ref020]\], we measured complete blood cell counts during both experimental sessions, as well as PBMC counts, in duplicate, using a Countess automated cell counter (Invitrogen). Other potential covariates included; age, body mass index, and socioeconomic status (i.e., parental education and income).

Data reduction and final measures {#sec013}
---------------------------------

Statistical analyses of cortisol data used log transformed cortisol values at 7 time-points and area under the curve with respect to increase (AUCi) \[[@pone.0221310.ref050]\]. The variables were examined for outliers (\>3 SD) and none were detected. Blood pressure values were reduced to 4 measures, from 30 minutes prior to testing to 15 minutes after (samples 1--4) to evaluate the fast sympathetic response. Moreover, systolic and diastolic blood pressures were combined to derive a measure of the mean arterial pressure (MAP) to describe the average response in blood pressure (i.e., MAP = \[(2 x diastolic) + systolic\] / 3). Raw gene expression data was analyzed based on the 2\^-(ΔΔCt) method, with normalization to a housekeeping gene, and compared to the first baseline measure in each session \[[@pone.0221310.ref043]\]. AUCi was computed for *NR3C1* to assess overall responses from baseline. \[[@pone.0221310.ref051]\]

For the four pro-inflammatory cytokines, considering high correlations \[[@pone.0221310.ref052]\] (Pearson correlations ranged from .30 to .72), principal component analysis (PCA) indexing *systemic inflammation* of IL-1β, IL-6, IL-8 and TNF-α measures was conducted for the four repeated measures using data from both sessions. In each instance, the first component was extracted for use in subsequent analyses. The four repeated items mapped to components with eigenvalues of 2.55 for the first time point, which explained 63.81% of the variance across all four cytokines, 2.29 for the second time point (57.25% of variance), 2.46 for the third time point (61.54% of variance), and 2.80 for the fourth time point (69.88% of variance). PCA of AUCi for all four cytokines yielded two components with eigenvalues 1.93 and 1.01, which explained 48.27% and 25.37% of the variance respectively. Closer inspection of the factor loading scores for the PCA of AUCi revealed that the first component was largely representative of three cytokines (IL1-β, IL-8, TNF-α) with the second representing IL-6 ([Table 1](#pone.0221310.t001){ref-type="table"}). PCA was also conducted on the four repeated measures independently within each session (TSST and no-stress). The four repeated items in the TSST session mapped onto components with eigenvalues 2.02--2.56, which explained 50.41%-64.01% of variance at each time point. The four repeated items in the no-stress session mapped onto components with eigenvalues 1.65--2.29, which explained 41.28%-57.16% of variance at each time point. The components mapped using data from both sessions were used to investigate within-person differences across sessions, while the components mapped within each session independently were used to investigate between-person differences (i.e. ELA status) within each session. Scores in all repeated questionnaires for both sessions were averaged to increase reliability (Pearson correlations ranged from .72 to .93). None of the demographics measures differed significantly between the ELA and control groups ([Table 2](#pone.0221310.t002){ref-type="table"}), and thus were not included as covariates in the analysis.

10.1371/journal.pone.0221310.t001

###### PCA of pro-inflammatory cytokine AUCi: Factor loading scores.

![](pone.0221310.t001){#pone.0221310.t001g}

  Cytokine   Factor 1   Factor 2
  ---------- ---------- ----------
  IL1-β      0.597      0.224
  IL-6       0.174      0.942
  IL-8       0.908      -0.080
  TNF-α      0.849      -0.265

10.1371/journal.pone.0221310.t002

###### Sample characteristics.

![](pone.0221310.t002){#pone.0221310.t002g}

  Variable, mean (SD)   Total (N = 12)   Control (N = 6)   ELA (N = 6)     P value diff
  --------------------- ---------------- ----------------- --------------- --------------------------------------------
  Age                   21.25 (2.3)      20.83 (1.6)       21.67 (2.9)     0.56
  SES (average)         2.83             2.83              2.83            0.23[^1^](#t002fn001){ref-type="table-fn"}
  1\. Working class     1                0                 1               
  2\. Lower middle      2                1                 1               
  3\. Middle            7                5                 2               
  4\. Upper middle      2                0                 2               
  BMI                   25.40 (3.7)      26.26 (3.7)       24.55 (4.0)     0.46
  BAI                   8.83 (8.2)       6.00 (5.8)        11.67 (9.8)     0.25
  BDI                   5.79 (6.3)       3.00 (1.3)        8.58 (8.2)      0.16
  STAI                  70.38 (20.4)     61.58 (9.0)       79.17 (25.5)    0.14
  LESS                  156.50 (84.9)    114.92 (60.9)     198.08 (89.4)   0.09
  LEQ                   20.21 (12.5)     13.75 (5.0)       26.67 (14.8)    0.07
  PSS- TSST             16.08 (8.6)      11.00 (5.3)       21.67 (8.6)     0.03
  PSS- no-stress        15.75 (10.0)     9.50 (4.9)        22.00 (10.1)    0.02

^1^p-value from Chi-square

SES- socioeconomic status; BMI- body mass index; BAI- Beck anxiety inventory; BDI- Beck depression inventory; STAI- state-trait anxiety inventory; LESS- life-event stress scale; LEQ- life events questionnaire; PSS- perceived stress scale.

Statistical analysis {#sec014}
--------------------

All statistical tests were carried out using SPSS version 25 (Windows). Repeated measures general linear models (GLMs), ordinary least squares multiple regression analyses, Pearson product-moment correlations, and t-tests were carried out as appropriate. Statistical analyses of changes in gene expression, physiological responses and cytokines levels were subjected to multivariate GLMs, with salivary cortisol, cytokines and gene expression as the repeated measure, condition (stress/no stress) as a within-subjects factor, and status (risk/control) as between-subjects factors. In addition to these analyses, univariate tests were applied to summary cortisol measures (AUCi, \[[@pone.0221310.ref050]\]) to ascertain reliability of findings, as well as blood pressure (MAP), gene expression, and pro-inflammatory cytokine measures. Huynh-Feldt corrections were applied if sphericity (significant differences in variance between groups) was significant, and only adjusted results are reported.

Adjustment for multiple comparisons was performed using Bonferroni correction in accordance with FDA guidelines for multiplicity in clinical trials \[[@pone.0221310.ref053]\]. As the primary endpoints, analyses for MAP and cortisol were controlled for collectively, and analyses within each secondary endpoint (i.e. gene expression, cytokine PCA, and raw cytokines) were adjusted independently. For the purposes of multiplicity, analyses within subgroups defined by ELA status were each considered additional tests to be accounted for. Results following adjustment and domains within which analyses were adjusted are summarized in [S2 Table](#pone.0221310.s002){ref-type="supplementary-material"}.

Results {#sec015}
=======

Sample characteristics and self-reported measures {#sec016}
-------------------------------------------------

The ELA and control group did not differ in demographics measures (i.e., age, socioeconomic status and body max index) ([Table 2](#pone.0221310.t002){ref-type="table"}). While group differences existed in exposure to traumatic events up to age 18 (i.e. SLESQ≥3 for ELA; SLESQ = 0 for controls), the ELA and control groups did not differ significantly in their experiences of stressful life events during the past year (univariate ANOVA between-subjects effect: F = 3.55, p = 0.089 for LESS, estimated effect size η^2^ = 0.26; F = 4.10, p = 0.070 for LEQ, estimated effect size η^2^ = 0.29). There were no significant differences in levels of anxiety \[[@pone.0221310.ref048]\] (F = 2.54, p = 0.142, estimated effect size η^2^ = 0.20) or depressive symptoms \[[@pone.0221310.ref047]\] (F = 2.73, p = 0.129, estimated effect size η^2^ = 0.22) between the two groups, although levels tended to be higher in the ELA group relative to controls for all measures of stress, anxiety, and depression ([Table 2](#pone.0221310.t002){ref-type="table"}). Further, the ELA group self-reported significantly greater perceived stress in the TSST session compared with controls (F = 6.07, p = 0.033, effect size η^2^ = 0.38), as well as in the no-stress condition (F = 7.49, p = 0.021, effect size η^2^ = 0.43), and tended to report more stress in response to the TSST (Likert scale from 1--10) (F = 4.02, p = 0.080).

Physiological, gene expression, and pro-inflammatory cytokines response to acute psychosocial stress compared with a no-stress control condition {#sec017}
------------------------------------------------------------------------------------------------------------------------------------------------

### Primary endpoints: Physiological {#sec018}

In the whole sample, repeated measures GLMs indicated significant within-subjects effect for salivary cortisol in response to the TSST, compared with a no-stress condition (Time x Session, F = 4.47, p = 0.003, estimated effect size η^2^ = 0.17), as well as for mean arterial pressure (Time x Session, F = 5.31, p = 0.003, η^2^ = 0.20). Compared with controls, the ELA group exhibited significantly higher mean arterial pressure response to the TSST (Time x Status, F = 8.59, p\<0.001, η^2^ = 0.46), and a non-significant trend towards a higher cortisol response in the TSST relative to no-stress (ΔAUCi: F = 3.58, p = 0.088) ([Fig 2](#pone.0221310.g002){ref-type="fig"}). Notably, no significant differences were observed between the ELA and control groups in the no-stress condition (Time x Status, F = 1.38, p = 0.257 for salivary cortisol; F = 1.01, p = 0.402 for MAP). Overall, these findings confirm some \[[@pone.0221310.ref054]\], but not all studies \[[@pone.0221310.ref007]\], indicating increased physiological reactivity to acute stress in young adults exposed to early adversity, compared with non-exposed individuals.

![Normalized change score for physiological, gene expression, and pro-inflammatory cytokine response to the TSST relative to the no-stress session for ELA group, control group and full sample.\
Change scores were calculated by standardizing summary AUCi using data from both sessions and subtracting participant values from the no-stress session from those in the TSST session. Error bars represent standard error of the mean. Change scores are expressed for the full sample (grey), ELA group (red), and control group (green).](pone.0221310.g002){#pone.0221310.g002}

### Secondary endpoint: Gene expression {#sec019}

In the whole sample, there was a significant within-subjects effect of TSST vs. no-stress condition on *NR3C1* gene expression with increased levels in the TSST (Time x Session, F = 4.85, p = 0.006, estimated effect size η^2^ = 0.19). Group analysis revealed increased levels in the control group (Time x Session, F = 5.09, p = 0.013, estimated effect size η^2^ = 0.36), but not in the ELA group, which had a blunted response to the TSST (Time x Session, F = 1.00, p = 0.406, estimated effect size η^2^ = 0.09) ([Fig 3A](#pone.0221310.g003){ref-type="fig"}), suggestive of *NR3C1* expression resistance and lower levels to inhibit the HPA axis. Notably, no differences were observed in *NR3C1* expression between the ELA and control groups in the no-stress condition (Time x Session, F = 0.78, p = 0.491) ([Fig 3B](#pone.0221310.g003){ref-type="fig"}).

![**Fold change in NR3C1 for ELA (dashed lines) and control groups (solid lines) in response to the TSST (left) and during the no-stress sessions (right).** Error bars represent standard error of the mean.](pone.0221310.g003){#pone.0221310.g003}

### Secondary endpoint: Pro-inflammatory cytokines {#sec020}

In the whole sample, PCA for the four repeated measures of IL-1β, IL-6, IL-8 and TNF-α did not reveal a significant within-subjects effect of TSST vs. no-stress condition using repeated measures GLM analysis (Time x Session, F = 0.29, p = 0.831). Similarly, an analysis of the first AUCi component did not reveal significant differences in systemic inflammation between the TSST and no-stress conditions (F = 2.09, p = 0.165). However, an analysis of the second AUCi component (largely representing IL-6) showed significantly greater pro-inflammatory responses to the TSST relative to the no-stress condition (F = 5.85, p = 0.026) ([Fig 2](#pone.0221310.g002){ref-type="fig"}). No differences were observed between the ELA and control groups in response to the TSST relative to no-stress using either AUCi PCA components.

Exploratory analyses of each pro-inflammatory cytokine revealed a significant within-subjects effect for IL-6 in response to the TSST (Time x Session, F = 2.97, p = 0.044, AUCi: F = 7.70, p = 0.018), but not for the other three cytokines (IL-1β, Time x Session, F = 0.80, p = 0.500, AUCi, F = 1.37, p = 0.257; IL-8, Time x Session, F = 0.85, p = 0.470, AUCi, F = 0.93, p = 0.434; TNF-α, Time x Session, F = 0.25, p = 0.863, AUCi, F = 0.40, p = 0.537).

Adjustment for multiple comparisons & sensitivity analyses {#sec021}
----------------------------------------------------------

Bonferroni correction for multiple comparisons was performed across analyses of primary endpoints and within analyses of each secondary endpoint in accordance with FDA guidelines for multiplicity in clinical trials \[[@pone.0221310.ref053]\]. All significant findings within primary endpoints passed correction for multiple testing at a false discovery rate of 0.05, as did observed differences in *NR3C1* gene expression between sessions and between ELA status groups within the stress session. By contrast, results of analyses with inflammatory principal components, as well as those with raw inflammatory cytokines, did not pass correction for multiple testing. A summary of results and Bonferroni-corrected p-values is provided in [S2 Table](#pone.0221310.s002){ref-type="supplementary-material"}.

Sensitivity analyses were conducted using the 'leave-one-out' method. Overall, results were robust to the removal of any individual participant. Differences in sample characteristics and self-report measures remained consistent upon removal of any given participant, as did physiological, inflammatory, and gene expression responses to the TSST relative to the no-stress session. Likewise, the ELA group continued to display increased MAP responses to the TSST relative to the control group. Differences between ELA and control groups in cortisol response to the TSST relative to no-stress (ΔAUCi) were modestly attenuated by removal of any given participant, but not appear to be driven by a single individual. [S3](#pone.0221310.s003){ref-type="supplementary-material"}--[S5](#pone.0221310.s005){ref-type="supplementary-material"} Tables provide the full results from 'leave one out' sensitivity analyses.

Discussion {#sec022}
==========

This pilot study delineates a program of research to study ELA-related programming of biological systems using a within-person, between-groups experimental design. To our knowledge, this is the first investigation of stress-induced gene expression and pro-inflammatory cytokines changes within-individuals, comparing stratified groups of ELA-exposed and control individuals. By comparing a validated laboratory-based stressor to a no-stress condition within the same individuals, we are able to disentangle the effects of acute stress from noisy measurements in the same individuals. Further, this design allow us to distinctly identify if/when differences between ELA-exposed and control individuals are context dependent (i.e. manifesting only during stress). Preliminary results provide evidence in humans of a dysregulated pattern of *NR3C1* gene expression activation as a consequence of ELA. Importantly, these changes manifest more acutely in the presence of stress-induced cortisol release as compared to a no-stress resting condition.

As predicted by previous research, the ELA group evince higher cortisol response and lower *NR3C1* gene expression in response to the TSST compared with controls, with no difference between groups in the no-stress condition. For pro-inflammatory cytokines, only IL-6 increased significantly in response to the laboratory-induced stressor. Overall, this pilot study provides a framework to investigate the biological embedding of ELA via a dynamic and dysregulated pattern spanning multiple levels of analysis (genomic and physiological), and affords the opportunity to parse out whether such dysregulation is context dependent (i.e. manifesting only during duress). This conceptual framework integrates multiple definitions of stress and disease at the environmental (e.g., ELA), psychological (e.g., psychosocial stress response), and biological levels (e.g., physiological, gene expression and pro-inflammatory cytokine responses) \[[@pone.0221310.ref038]\].

These preliminary findings concur with the receptor-mediated model of glucocorticoid signaling resistance \[[@pone.0221310.ref055]\]. First, ELA was associated with increased physiological response to the TSST compared with controls, confirming some \[[@pone.0221310.ref054]\], but not all studies \[[@pone.0221310.ref007]\], indicating increased physiological reactivity in young adults exposed to early adversity. Second, chronic exposure to cortisol, as a consequence of ELA, can lead to a compensatory response whereby glucocorticoid sensitivity decreases (e.g. via decreased receptor availability). Here we replicated prior evidence of reduced *NR3C1* expression levels in ELA-exposed individuals, but only in response to acute laboratory stress. Third, in vitro studies have established a connection between glucocorticoid exposure and diminished capacity of immune cells to inhibit pro-inflammatory cytokines in individuals exposed to psychological stress. Here, only pro-inflammatory cytokine IL-6 increased significantly in response to acute stress, replicating previous studies \[[@pone.0221310.ref056]\].

Methodological strengths of this program of research include a laboratory-based within-subjects experimental design, which allows stronger causal inferences. We collected repeated measurements over a relatively long time scale to document changes in gene expression and pro-inflammatory cytokines. Our within-subjects, between-groups design, combined with four repeated measurements in each session, reduced biological variability and increased power to detect true associations. Although this preliminary study lacked sufficient power to investigate three-way interactions (i.e. Time x Session x Status) or perform multivariate analyses across outcomes, the strength of the design affords the opportunity to test for such associations at comparably smaller sample sizes than studies without as many repeated measures. Furthermore, inclusion of a no-stress control session allows researchers the chance to identify differences that are context (i.e. stress) dependent.

We acknowledge limitations. First, this was a pilot feasibility study with a small sample size. Although comparable to similar prior investigations \[[@pone.0221310.ref026], [@pone.0221310.ref031], [@pone.0221310.ref057]\], the results from this study still need to be interpreted with caution. Notwithstanding, the strength of the within-subjects experimental design combined with the leave-one-out sensitivity analysis and adjustment for multiple comparisons alleviate concerns about spurious findings of the results reported herein. Second, we focused on men in this exploratory study due to known sex differences in the stress response \[[@pone.0221310.ref041]\] and the small sample size for stratified analyses. Future studies with larger sample size including both males and females are needed to replicate these findings. Third, we focused on a single hypothesis-driven gene. There are multiple biological pathways that are activated in response to stress that may play a downstream role in disease susceptibility, such as the conserved transcriptional response to adversity pathway \[[@pone.0221310.ref030]\]. Prior research has investigated multiple genes using microarrays \[[@pone.0221310.ref023], [@pone.0221310.ref026]--[@pone.0221310.ref028], [@pone.0221310.ref031]\]. Future studies with adequate sample size will benefit by testing larger groups of genes/pathways. Fourth, this study did not consider specific types of ELA, or timing of exposure. Here, we focused specifically on severity of multiple (i.e., minimum of three) ELA exposures up to age 18 years. Future research can explore specific types of ELA in different populations and settings. Fifth, while this controlled laboratory experimental design allows stronger causal inferences, there are other alternative explanatory variables that can explain individual differences in response to stress (e.g., temperament), which can be considered in future studies. Further, our study included non-Hispanic white males and thus future research need to test whether the association generalizes to other populations. Finally, although we included a no-stress condition to control for the higher degree of noise associated with gene expression measurements, the control session did include the stress of venipuncture. However, this is unavoidable technical limitation for collecting sufficient immune cells of high quality for gene expression research.

In conclusion, ELA may program physiological systems in a maladaptive manner more likely to manifest during times of duress, predisposing individuals to the negative health consequences of everyday stressors. Although increased activation of the glucocorticoid-immune signaling in response to acute stress is considered adaptive in the short-term, persistent activation can increase risk for mental and physical health problems. This program of research and preliminary analyses provides a framework to investigate new targets for therapeutic interventions mitigating the negative effects of early adversity, such as pharmacological agents acting on the glucocorticoid receptor \[[@pone.0221310.ref058]\]. Further, while previous risk factors and biomarkers of stress contributed to our understanding of biological embedding processes, these are nevertheless static characteristics that have not explained health outcomes very well. For example, considering high failure rates for depression treatments, and in order to tailor individual interventions, identifying objective changes in stress-induced gene expression may help to predict short-term intervention efficacy in clinical and non-clinical settings. An example for such an effort could be to utilize models of dynamic cellular markers as individual-level factors to account for variation in intervention response and clinical outcomes \[[@pone.0221310.ref017]--[@pone.0221310.ref019]\]. Thus, future research in this area can have a range of impacts for basic science, intervention studies and clinical practice that will influence treatments to match the specific cellular processes operating within an individual.

Supporting information {#sec023}
======================

###### Summary information for inflammatory cytokines (pg/mL) by time, session, and group status.

Data presented as sMean (Standard Deviation).

(DOCX)

###### 

Click here for additional data file.

###### False discovery rate adjustment across study primary and secondary endpoints.

Adjustment was performed across both primary outcomes (cortisol and MAP) and within each secondary outcome. Significant results are bolded. \*Univariate analyses of raw cytokine AUCi were considered confirmatory following results of repeated measures models, and thus were not included in the count of total tests within raw cytokines secondary endpoint.

(DOCX)

###### 

Click here for additional data file.

###### Leave one out sensitivity analyses of Sample Characteristics (p-value).

Original results shown at top. Results which vary in significance from main findings are bolded.

(DOCX)

###### 

Click here for additional data file.

###### Leave one out sensitivity analyses of stress-induced changes in physiological, gene expression, and pro-inflammatory cytokines (p-value).

Original results shown at top. Results which vary in significance from main findings are bolded.

(DOCX)

###### 

Click here for additional data file.

###### Leave one out sensitivity analyses of group differences in stress-induced changes in physiological, gene expression, and pro-inflammatory cytokines (p-value).

Original results shown at top. Results which vary in significance from main findings are bolded.

(DOCX)

###### 

Click here for additional data file.

We thank all nurses at the CRC and the participants in this study.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.
