- Research note
- Open Access
Integrated microRNA and mRNA gene expression in peripheral blood mononuclear cells in response to acute psychosocial stress: a repeated-measures within-subject pilot study
BMC Research Notes volume 14, Article number: 222 (2021)
The impact of psychosocial stress on a variety of negative health outcomes is well documented, with current research efforts directed at possible mechanisms. Here, we focused on a potential mechanism involving differential expression of mRNA and microRNA in response to acute psychosocial stress. We utilized a validated behavioral paradigm, the Trier Social Stress Test (TSST), to induce acute psychosocial stress in a cohort of volunteers. Stress reactivity was assessed repeatedly during the TSST using saliva samples that were analyzed for levels of cortisol. Peripheral blood mononuclear cells were extracted from blood drawn at baseline and at two time points following the stress paradigm. Total RNA was extracted, and mRNA and microRNA microarrays were utilized to assess within-subject changes in gene expression between baseline and the two post-stressor time points.
For microarray gene expression analysis, we focused on 12 participants who showed a robust cortisol response to the task, as an indicator of robust HPA-axis activation. We discovered a set of mRNAs and miRNAs that exhibited dynamic expression change in response to the TSST in peripheral blood mononuclear cells, further characterizing the link between psychosocial stress and cellular response mechanisms.
Organismal stress responses are associated with widespread changes in the transcriptome of circulating peripheral blood mononuclear cells (PBMCs), with changes in expression of genes regulating cell signaling, adhesion, differentiation, and immune pathways . Such stress-associated transcriptional changes in PBMCs are mediated, in part, through microRNAs (miRNA) [2, 3], such as miRs 7b, 16, 20b, 21, 26b, 29a, 29c, 126 and 144 [4,5,6,7].
A popular laboratory-based paradigm for the study of psychosocial stress physiology is the Trier Social Stress Test (TSST) . A previous study reported TSST-associated changes in salivary levels of miR-20b, -21 and 26b , but their expression levels in PBMCs and their association with gene expression remain unknown. Here, we report results from a TSST pilot study examining within-subject repeated-measures of mRNA and miRNA in PBMCs.
Participants were Caucasian males recruited via flyers posted on the Stony Brook University campus and through advertisements in local newspapers in surrounding communities on Long Island, NY, USA. Study exclusion criteria included: age under 18 years, significant current or prior presence of psychiatric illness, smoking, presence of substance or alcohol abuse in the past six months, hormonal, mood-altering, heart disease and/or cardiovascular conditions related, or psychoactive drug medication, BMI ≥ 30. Participants were pre-screened for exclusion criteria by phone. Fifty-five participants were part of the study, which was approved by the Institutional Review Board (IRB) and Committee on Research Involving Human Subjects (CORIHS) of Stony Brook University, Stony Brook, USA. We selected 12 participants who showed a robust cortisol response to the TSST (“Cortisol Responders”; based on a criterion suggested by Miller et al. ) as a discovery cohort to generate repeated-measures, within-subject genome-wide mRNA and miRNA expression data.
Psychosocial stress procedure: TSST
All participants were scheduled for a 4-h weekday laboratory session conducted at the General Clinical Research Center (GCRC) at the Stony Brook University Hospital from 1400 to 1800 h. Upon arrival, participants were informed about the study, completed the consent procedure, and underwent insertion of an intravenous catheter into the arm for the first blood draw. Participants then rested for a period of at least 45 min before they were introduced to the psychosocial stress paradigm. The Trier Social Stress Test (TSST) is a standardized social stress protocol that was used to induce psychosocial stress and that has been found to elicit a strong and reliable physiological response to laboratory stress . Briefly, the standardized laboratory stressor consists of a 5-min preparation period, a 5-min free speech task during which participants give a mock job interview, and a 5-min mental arithmetic task wherein participants perform serial subtraction by seventeen all in front of a non-responsive audience.
Biological measures: cortisol
To assess cortisol levels, ten salivary samples were obtained from each participant throughout the duration of the study. Samples were collected using Salivettes® (Sarstedt) at arrival, 2 min before beginning the TSST, and again at 2, 10, 20, 30, 45, 60, 90- and 105-min post stress task. Participants were instructed to gently chew on the cotton swab for approximately one minute during each collection. Cotton swabs were then transferred to plastic containers and immediately stored at − 20 °C until shipment to the analyzing laboratory (N. Rohleder, Brandeis University, Waltham, MA, USA). Salivary cortisol concentrations were measured using a commercially available chemiluminescence immunoassay (CLIA; RE62019; IBL International, Toronto, Canada). Inter-assay variability was 3.3%, and intra-assay variability was 3.4%. Salivary cortisol concentrations have been shown to highly correlate with serum cortisol concentrations, making salivary cortisol a reliable indicator of HPA activity . Participants had no oral intake except for water for 3 h before beginning the study session. Cortisol response to the stress paradigm was assessed by subtracting the highest mean cortisol concentration following the stressor minus the baseline cortisol concentration (two minutes before the TSST). Using a criterion suggested by Miller et al. , those participants who showed a cortisol increase of at least 1.5 nmol/l from baseline were assigned in the “Cortisol Responder” group (n = 33, mean cortisol increase of 8.55 nmol/l (± 5.63), while those who did not show an increase of at least 1.5 nmol/l from baseline were placed in the “Cortisol Non-Responder” group (N = 22, mean cortisol increase of − 3.889 nmol/l (± 6.45). The two groups differed in age (responders: 34.58 years ± 17.09; non-responders: 26.18 years ± 13.44; t-test p < 0.05). There were no significant differences in psychological variables of interest, as the BAI, BDI, TiCS, STAI-S and CTQ sum scores were not significantly different between the two groups (all p > 0.05). There was no significant difference in the number of stressful life events (SLEs) between the two groups (p > 0.05). Subsequent gene expression analyses with microarray profiling were conducted on the first 12 Cortisol Responders.
Biological measures: blood collection
Directly upon arrival and following the informed consent procedure, an intravenous (IV) catheter was inserted into each participant’s arm for the first blood draw (Time Point A). Participants then rested for a period of at least 45 min before they were introduced to the psychosocial stress paradigm. The second blood draw occurred 45 min after the stress paradigm (Time Point B) and the third blood draw occurred 105 min post stress (Time Point C). Ten ml of EDTA blood (VacuetteR Greiner Bio-one) was collected at each of the three time points. Blood was immediately processed for PBMC extraction according to the Leucosep® Instruction Manual (Greiner Bio-one), and cell pellets were frozen and stored at − 80 °C.
RNA isolation, quantification, quality control
Total RNA was extracted from PBMCs at each time point using the Qiagen All-Prep kit with company-recommended modification to extract small RNA (Qiagen). RNA quantity was assessed using a Nanodrop Technologies ND-1000 instrument (NanoDrop Technologies). Total RNA quality was measured using the Agilent RNA 6000 Pico kit and the Agilent Bioanalyzer 2100A (Agilent Technologies), while miRNA content was resolved using the Small RNA Kit (Agilent). RNA samples were stored at − 80 °C.
MiRNA microarray profiling
miRNA microarray profiling was performed at the Microarray Core Facility of Cold Spring Harbor Laboratory using the Affymetrix GeneChip miRNA 2.0 Array according to the manufacturer’s protocol. Briefly, 400 ng of total RNA was labelled using the 3DNA Array Detection FlashTag™ Biotin HSR Kit following the manufacturer's recommendations (Genisphere). Arrays were scanned on the Affymetrix GeneChip Scanner 3000, and feature extraction was conducted using Affymetrix Command Console software. Microarray background correction, quantile normalization, probe set summarization, and log2 transformation were performed with the miRNA QC Tool Version 188.8.131.52 (Affymetrix).
MRNA microarray profiling
mRNA microarray profiling was performed at the ShanghaiBio Corporation using Affymetrix Human Gene ST 1.0 Arrays according to the manufacturer’s protocol. Briefly, 500 ng of total RNA was used in complementary DNA synthesis and labeling, followed by hybridization, washing, scanning, and image processing of GeneChips according to standard Affymetrix protocols. Quantitative image files were loaded into the Bioconductor Affymetrix linear modeling Graphical User Interface software  and background adjustment, quantile normalization, probe summarization, and log2 transformation was performed using the Robust Multi-Array Average method .
Statistical analysis of mRNA and miRNA microarray data
Time-course microarray data analysis was conducted using the Linear Models for Microarray Data (limma) Bioconductor package . Linear models were fitted to distinguish differentially expressed mRNAs between the three time points (A, B, C) on the Affymetrix platform. mRNA showing limma p < 0.01 for at least one of the time point comparisons (B to A, C to A, B to C) correcting for multiple comparisons using false discovery rate (FDR) correction, as well as with a fold change (FC) |FC|≥ 1.3 were highlighted as genes that show greater variation in expression in response to the TSST. The fold change threshold of |FC|≥ 1.3 was chosen as literature has conventionally utilized fold changes ≥ 1.3–2 to suggest differential expression in mRNA microarray data . Also, fold change of gene expression in blood has usually been less than two-fold [3, 16] and multiple studies have used a cut-off of FC > 1.2 for identifying changes of gene expression specifically in blood [16,17,18]. miRNAs generally show more subtle changes in expression than mRNAs, but even small expression changes can have significant biological effect as each miRNA has multiple gene targets and can influence multiple genes in a given pathway. Therefore, miRNAs showing a |FC|≥ 1.2 and limma p value < 0.01 for at least one of the time point comparisons (B to A, C to A, B to C) were highlighted as showing change in expression in response to the TSST.
MiRNA/mRNA integration and functional profiling
Integration of mRNA and miRNA expression was performed by the computational gene network prediction tool Ingenuity Pathway Analysis (IPA, Ingenuity® Systems). Our study showed a limited number of mRNAs that were significant according to limma p < 0.01 and |FC|≥ 1.3, possibly in part due to the small study sample size and conservative statistical inclusion criteria. It is worthwhile noting that a single miRNA can subtly regulate multiple genes in a given pathway or process, having an additive biological effect. Therefore, mRNAs that showed significant changes in expression at any time contrast according to limma p < 0.01 were input into IPA. In parallel, human orthologues of differentially expressed miRNAs (limma p < 0.01 and |FC|≥ 1.2) were input into the IPA program as well. To integrate the two datasets, identification of negatively correlated miRNA/mRNA pairs was performed using the Target Filter Analysis in IPA. The algorithm uses experimentally validated interactions as well as predicted microRNA–mRNA interactions from TargetScan. Fold change at each miRNA and mRNA was integrated at Time B vs Time A, Time C vs Time A, and Time C vs Time B. The expression pairing analysis was used to identify negatively correlated miRNA/mRNA pairs that had experimentally demonstrated targeting or a moderate to high degree of likelihood of a targeting relationship based on prediction algorithms.
Participant characteristics and gene expression analysis
Table 1 lists participant characteristics and Table 2 lists 9 mRNAs that were differentially expressed according to the criteria of limma p < 0.01 and |FC|≥ 1.3. mRNAs that showed significant changes in expression at any time contrast according to limma p < 0.01 were put into IPA for expression pairing and are presented in Additional file 1: Table S1. Limma analysis revealed 33 miRNA probe sets on the Affymetrix array that showed significant expression changes (p < 0.01 and |FC|> 1.2) in Cortisol Responders over the course of the TSST. Only the miRNAs with identified human orthologues that readily mapped into IPA were considered in further analysis. In the instance that stem-loops were identified by the microarray, both mature products were put into IPA when as per verbal instructions from IPA Technical Support. The list of resulting miRNAs is presented in Table 3.
Negatively correlated miRNA/mRNA pairs were identified at contrasts B vs A, C vs A, or C vs B using IPA’s Target Filter Analysis. Since miRNAs often result in subtle changes in gene expression, mRNAs that showed significant change in expression according to limma p < 0.01 were used in an expression-pairing analysis along with miRNAs significant according to criteria limma p < 0.01 and |FC|≥ 1.2. Resultant miRNA/mRNA expression pairs at each time contrast are listed in Additional file 2: Table S2A–C.
This pilot study examined genome-wide changes in mRNA and miRNA gene expression in PBMCs in response to TSST-induced psychosocial stress within-subjects in a repeated-measures design. By removing between-subject genomic variance that would otherwise obscure any relationship between stress and gene regulation, we could identify differentially expressed miRs and target mRNAs at FDR-corrected levels, although we replicated none of the previously reported stress-related miRs. Non-replication is common is small-scale studies, but the aggregate of our results and previous reports suggests that properly powered large-scale studies of stress-related gene regulation by miRNAs are warranted to better understand the molecular mechanisms of stress-related gene dysregulation.
Limitations of this study include a small sample size consisting of healthy Caucasian males only. Thus, cautious interpretation is advised. Replication in a larger and more diverse samples that are more broadly representative of the variation in human stress-related phenotypes should be considered in future work.
In addition, microarray results should be confirmed by an orthologous method such as RT-qPCR in order to confirm changes in gene expression.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Nater UM, Whistler T, Lonergan W, Mletzko T, Vernon SD, Heim C. Impact of acute psychosocial stress on peripheral blood gene expression pathways in healthy men. Biol Psychol. 2009;82:125–32.
Leung AK, Sharp PA. MicroRNA functions in stress responses. Mol Cell. 2010;40:205–15.
Mendell JT, Olson EN. MicroRNAs in stress signaling and human disease. Cell. 2012;148:1172–87.
Gidron Y, De Zwaan M, Quint K, Ocker M. Influence of stress and health-behaviour on miRNA expression. Mol Med Rep. 2010;3:455–7.
Katsuura S, Kuwano Y, Yamagishi N, Kurokawa K, Kajita K, Akaike Y, Nishida K, Masuda K, Tanahashi T, Rokutan K. MicroRNAs miR-144/144* and miR-16 in peripheral blood are potential biomarkers for naturalistic stress in healthy Japanese medical students. Neurosci Lett. 2012;516:79–84.
Honda M, Kuwano Y, Katsuura-Kamano S, Kamezaki Y, Fujita K, Akaike Y, Kano S, Nishida K, Masuda K, Rokutan K. Chronic academic stress increases a group of microRNAs in peripheral blood. PLoS ONE. 2013;8:e75960.
Vaisvaser S, Modai S, Farberov L, Lin T, Sharon H, Gilam A, Volk N, Admon R, Edry L, Fruchter E, et al. Neuro-epigenetic indications of acute stress response in humans: the case of MicroRNA-29c. PLoS ONE. 2016;11:e0146236.
Kirschbaum C, Pirke KM, Hellhammer DH. The ‘Trier Social Stress Test’—a tool for investigating psychobiological stress responses in a laboratory setting. Neuropsychobiology. 1993;28:76–81.
Wiegand C, Heusser P, Klinger C, Cysarz D, Bussing A, Ostermann T, Savelsbergh A. Stress-associated changes in salivary microRNAs can be detected in response to the Trier Social Stress Test: an exploratory study. Sci Rep. 2018;8:7112.
Miller R, Plessow F, Kirschbaum C, Stalder T. Classification criteria for distinguishing cortisol responders from nonresponders to psychosocial stress: evaluation of salivary cortisol pulse detection in panel designs. Psychosom Med. 2013;75:832–40.
Kirschbaum C, Hellhammer DH. Salivary cortisol in psychoneuroendocrine research: recent developments and applications. Psychoneuroendocrinology. 1994;19:313–33.
Wettenhall JM, Simpson KM, Satterley K, Smyth GK. affylmGUI: a graphical user interface for linear modeling of single channel microarray data. Bioinformatics. 2006;22:897–9.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:1544–6115.
Hartmann B, Castelo R, Blanchette M, Boue S, Rio DC, Valcarcel J. Global analysis of alternative splicing regulation by insulin and wingless signaling in Drosophila cells. Genome Biol. 2009;10:R11.
Sharp FR, Jickling GC, Stamova B, Tian Y, Zhan X, Liu D, Kuczynski B, Cox CD, Ander BP. Molecular markers and mechanisms of stroke: RNA studies of blood in animals and humans. J Cereb Blood Flow Metab. 2011;31:1513–31.
Jickling GC, Xu H, Stamova B, Ander BP, Zhan X, Tian Y, Liu D, Turner RJ, Mesias M, Verro P, et al. Signatures of cardioembolic and large-vessel ischemic stroke. Ann Neurol. 2010;68:681–92.
Stamova B, Xu H, Jickling G, Bushnell C, Tian Y, Ander BP, Zhan X, Liu D, Turner R, Adamczyk P, et al. Gene expression profiling of blood for the prediction of ischemic stroke. Stroke. 2010;41:2171–7.
We thank Drs. Leighann Litcher-Kelly and Arthur Stone for their suggestions in designing the TSST, and the other undergraduate and graduate students of the Canlilab for their assistance in conducting the TSST, and Renee Wu for her statistical advice.
Open Access funding enabled and organized by Projekt DEAL. This research was supported by funding from NSF BCS-0843346 to T.C. The work of A. M.A. was supported by a stipend from the German Research Foundation (MU-2847).
Ethics approval and consent to participate
Ethical approval and a written informed consent to participate of every subject was obtained and approved by the IRB (IRB Approved by CORIHS Stony Brook University on 4/15/2011). None of the participants were under the age of 16 at the time of our study. This study complied with the most recent version of ethical principles regarding research with human participants (Declaration of Helsinki) in terms of being rigorous, using all the available good-quality data, answering the most important question(s) about the topic, presenting a sound statistical analysis, presenting the results clearly, and aiming to lead to beneficial outcomes.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Table S1A–C. Differentially Expressed mRNAs according to limma p < 0.01. These tables list mRNAs found to be differentially expressed according to limma p < 0.01 for one or more of the contrasts (Time Point B vs Time Point A, Time Point C vs Time Point A, and Time Point C vs Time Point B). This is the complete list of genes used in IPA pathway analysis. mRNAs highlighted in yellow have fold change ≥ 1.3.
Additional file 2: Table S2A–C
. mRNA/miRNA expression pairings. Negatively correlated mRNA/miRNA pairs identified by IPA Target Filter analysis at Time Point B vs A, Time Point C vs A and Time Point C vs B are presented in Table S2A-C, respectively.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jurkiewicz, M.M., Mueller-Alcazar, A., Moser, .A. et al. Integrated microRNA and mRNA gene expression in peripheral blood mononuclear cells in response to acute psychosocial stress: a repeated-measures within-subject pilot study. BMC Res Notes 14, 222 (2021). https://doi.org/10.1186/s13104-021-05635-3
- Psychosocial stress
- Gene expression