CoChat
Sign Up
Cited by 318Mindfulness and Compassion InterventionsBiofield Effects and BiophysicsReligion, Spirituality, and Psychology

The DOI matches its Crossref record; 71 of the 71 references it cites were found in a registry or an index.

Relaxation Response Induces Temporal Transcriptome Changes in Energy Metabolism, Insulin Secretion and Inflammatory Pathways

Open access
Bhasin, Manoj K., Dusek, Jeffery A., Chang, Bei-Hung, Joseph, Marie G., Denninger, John W., Fricchione, Gregory L., Benson, Herbert, Libermann, Towia A.·2013·PLoS ONE
DOI 10.1371/journal.pone.0062817PMC PMC3641112
Open this paper in CoChat — free accountOpen PDF

Abstract

The relaxation response (RR) is the counterpart of the stress response. Millennia-old practices evoking the RR include meditation, yoga and repetitive prayer. Although RR elicitation is an effective therapeutic intervention that counteracts the adverse clinical effects of stress in disorders including hypertension, anxiety, insomnia and aging, the underlying molecular mechanisms that explain these clinical benefits remain undetermined. To assess rapid time-dependent (temporal) genomic changes during one session of RR practice among healthy practitioners with years of RR practice and also in novices before and after 8 weeks of RR training, we measured the transcriptome in peripheral blood prior to, immediately after, and 15 minutes after listening to an RR-eliciting or a health education CD. Both short-term and long-term practitioners evoked significant temporal gene expression changes with greater significance in the latter as compared to novices. RR practice enhanced expression of genes associated with energy metabolism, mitochondrial function, insulin secretion and telomere maintenance, and reduced expression of genes linked to inflammatory response and stress-related pathways. Interactive network analyses of RR-affected pathways identified mitochondrial ATP synthase and insulin (INS) as top upregulated critical molecules (focus hubs) and NF-κB pathway genes as top downregulated focus hubs. Our results for the first time indicate that RR elicitation, particularly after long-term practice, may evoke its downstream health benefits by improving mitochondrial energy production and utilization and thus promoting mitochondrial resiliency through upregulation of ATPase and insulin function. Mitochondrial resiliency might also be promoted by RR-induced downregulation of NF-κB-associated upstream and downstream targets that mitigates stress.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1371/journal.pone.0062817

Cites

71 of 71 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Individual gene-based differential expression analysis. A) Differentially expressed genes identified by 3 across-group comparisons (N1 vs. N2, N1 vs. M, and N2 vs. M) at T0, T1 and T2. Venn diagrams depict the overlap of genes identified by these 3 comparisons at each time point. B) Heat map of genes that were significantly differentially expressed comparing N1 vs. N2 and N1 vs. M at T1 and T2 (marked with arrow in Venn diagrams). Gene expression is shown with a pseudocolor scale (−1 to 1) with red color denoting increased and green color denoting decreased fold change in gene expression. The rows represent the genes and columns represent subjects in N1, N2 and M groups at T0, T1 and T2. C) Differentially expressed genes identified by 3 within-group comparisons at different time points (T0 vs. T1, T0 vs. T2 and T1 vs. T2). Venn diagrams depict the overlap of genes identified by the 3 comparisons within each group.
Figure 1 Individual gene-based differential expression analysis. A) Differentially expressed genes identified by 3 across-group comparisons (N1 vs. N2, N1 vs. M, and N2 vs. M) at T0, T1 and T2. Venn diagrams depict the overlap of genes identified by these 3 comparisons at each time point. B) Heat map of genes that were significantly differentially expressed comparing N1 vs. N2 and N1 vs. M at T1 and T2 (marked with arrow in Venn diagrams). Gene expression is shown with a pseudocolor scale (−1 to 1) with red color denoting increased and green color denoting decreased fold change in gene expression. The rows represent the genes and columns represent subjects in N1, N2 and M groups at T0, T1 and T2. C) Differentially expressed genes identified by 3 within-group comparisons at different time points (T0 vs. T1, T0 vs. T2 and T1 vs. T2). Venn diagrams depict the overlap of genes identified by the 3 comparisons within each group.CC BY · PubMed Central
Figure 2 Temporal genomic expression patterns during one session of RR elicitation. Genes that were differentially expressed either across or within groups comparisons at different time points were used as the seed set of genes for Self-Organizing Map (SOM) analysis. These differentially expressed genes were partitioned to 18 separate maps according to Pearson correlation coefficient based distance metrics (Figure S2). Selected biologically interesting SOM maps were manually clustered into 4 biologically relevant categories based on the gene expression of N1, N2 and M groups at the 3 time points in one session of RR elicitation: Long-term Downregulation; Long-term Upregulation; Progressive Upregulation; and Progressive Downregulation. One representative pattern for each of these 4 biologically relevant categories is shown in the figure. The figure displays the box plot of the gene expression with X-axis representing time points and groups, and Y-axis representing scaled gene expression data from −1 to +1.
Figure 2 Temporal genomic expression patterns during one session of RR elicitation. Genes that were differentially expressed either across or within groups comparisons at different time points were used as the seed set of genes for Self-Organizing Map (SOM) analysis. These differentially expressed genes were partitioned to 18 separate maps according to Pearson correlation coefficient based distance metrics (Figure S2). Selected biologically interesting SOM maps were manually clustered into 4 biologically relevant categories based on the gene expression of N1, N2 and M groups at the 3 time points in one session of RR elicitation: Long-term Downregulation; Long-term Upregulation; Progressive Upregulation; and Progressive Downregulation. One representative pattern for each of these 4 biologically relevant categories is shown in the figure. The figure displays the box plot of the gene expression with X-axis representing time points and groups, and Y-axis representing scaled gene expression data from −1 to +1.CC BY · PubMed Central
Figure 3 Significantly enriched pathways with progressive patterns identified using gene set enrichment analysis. A) Upregulated Pathways B) Downregulated Pathways. The solid dots indicate significantly affected pathways (False Discovery Rate <25%) identified from across group comparisons (N1 vs. N2, N1 vs. M and N2 vs. M) at a particular time point (T0, T1 and T2). The asterisks represent significance and directionality of enrichment (P value<0.09 *, P value<0.05 **, P value<0.01 ***) identified from within group comparisons at different time points (T0 vs. T1, T0 vs. T2, T1 vs. T2). The red and green color asterisks indicate up- and down-regulated enrichment of pathways respectively. The heatmaps depicting relative expression of selected genes from representative pathways are shown in panels on the right side. Gene expression is shown with a pseudo color scale (−3 to 3) with red and green colors denoting increased and decreased relative expression respectively. Pathways with progressive patterns were enriched (up- or down- regulated) in N2 and M groups with greater significance of enrichments in M group. Furthermore, increasing enrichment over time within one session of RR elicitation was observed in M group.
Figure 3 Significantly enriched pathways with progressive patterns identified using gene set enrichment analysis. A) Upregulated Pathways B) Downregulated Pathways. The solid dots indicate significantly affected pathways (False Discovery Rate <25%) identified from across group comparisons (N1 vs. N2, N1 vs. M and N2 vs. M) at a particular time point (T0, T1 and T2). The asterisks represent significance and directionality of enrichment (P value<0.09 *, P value<0.05 **, P value<0.01 ***) identified from within group comparisons at different time points (T0 vs. T1, T0 vs. T2, T1 vs. T2). The red and green color asterisks indicate up- and down-regulated enrichment of pathways respectively. The heatmaps depicting relative expression of selected genes from representative pathways are shown in panels on the right side. Gene expression is shown with a pseudo color scale (−3 to 3) with red and green colors denoting increased and decreased relative expression respectively. Pathways with progressive patterns were enriched (up- or down- regulated) in N2 and M groups with greater significance of enrichments in M group. Furthermore, increasing enrichment over time within one session of RR elicitation was observed in M group.CC BY · PubMed Central
Figure 4 Significantly enriched pathways with long-term patterns identified using gene set enrichment analysis. A) Upregulated Pathways B) Downregulated Pathways. The solid dots indicate significantly affected pathways (False Discovery Rate <25%) identified from across group comparisons (N1 vs. N2, N1 vs. M and N2 vs. M) at a particular time point (T0, T1 and T2). The asterisks represent significance and directionality of enrichment (P value<0.09 *, P value<0.05 **, P value<0.01 ***) identified from within group comparisons at different time points (T0 vs. T1, T0 vs. T2, T1 vs. T2). The red and green color asterisks indicate up- and down-regulated enrichment of pathways respectively. The heatmaps depicting relative expression of selected genes from representative pathways are shown in panels on the right side. Gene expression is shown with a pseudo color scale (−3 to 3) with red and green colors denoting increased and decreased relative expression respectively. Pathways with long-term patterns were enriched (up- or down- regulated) only in M group. Furthermore, increasing enrichment over time within one session of RR elicitation was observed in M group.
Figure 4 Significantly enriched pathways with long-term patterns identified using gene set enrichment analysis. A) Upregulated Pathways B) Downregulated Pathways. The solid dots indicate significantly affected pathways (False Discovery Rate <25%) identified from across group comparisons (N1 vs. N2, N1 vs. M and N2 vs. M) at a particular time point (T0, T1 and T2). The asterisks represent significance and directionality of enrichment (P value<0.09 *, P value<0.05 **, P value<0.01 ***) identified from within group comparisons at different time points (T0 vs. T1, T0 vs. T2, T1 vs. T2). The red and green color asterisks indicate up- and down-regulated enrichment of pathways respectively. The heatmaps depicting relative expression of selected genes from representative pathways are shown in panels on the right side. Gene expression is shown with a pseudo color scale (−3 to 3) with red and green colors denoting increased and decreased relative expression respectively. Pathways with long-term patterns were enriched (up- or down- regulated) only in M group. Furthermore, increasing enrichment over time within one session of RR elicitation was observed in M group.CC BY · PubMed Central
Figure 5 Interactive network and top focus gene hubs identified from significantly affected pathways. The figure represents the top focus genes. A) Progressive upregulated Pathways, B) Progressive downregulated Pathways, and C) Integrated network of Long-term and Progressive affected pathways. The top focus hubs were identified from complex interactive networks generated from pathways with progressive and long-term patterns. The focus gene hubs were identified using the bottleneck algorithm for identification of the most interactive molecules with tree like topological structure. The bottleneck algorithm ranks genes on the basis of significance level with smaller rank indicating increasing confidence. The pseudocolor scale from red to green represents the bottleneck ranks from 1 to 20.
Figure 5 Interactive network and top focus gene hubs identified from significantly affected pathways. The figure represents the top focus genes. A) Progressive upregulated Pathways, B) Progressive downregulated Pathways, and C) Integrated network of Long-term and Progressive affected pathways. The top focus hubs were identified from complex interactive networks generated from pathways with progressive and long-term patterns. The focus gene hubs were identified using the bottleneck algorithm for identification of the most interactive molecules with tree like topological structure. The bottleneck algorithm ranks genes on the basis of significance level with smaller rank indicating increasing confidence. The pseudocolor scale from red to green represents the bottleneck ranks from 1 to 20.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 71 references

Research works printed but not linked: 1

  1. HB Benson, 1974, The relaxation response, Psychiatry: Journal for the Study of Interpersonal Processes, 37, 37, 10.1080/00332747.1974.11023785 · DOI

Cited by

Sorted by citations

Related

See all

Metadata

Authors
Bhasin, M. K., Dusek, J. A., Chang, B.-H., Joseph, M. G., Denninger, J. W., Fricchione, G. L., Benson, H., & Libermann, T. A.
Affiliations
  • Beth Israel Deaconess Medical Center (US)
  • Abbott Northwestern Hospital (US)
  • Boston University (US)
  • Harvard University (US)
Corresponding author
Libermann, Towia A.
ORCID
Venue
PLoS ONE
Publisher
Public Library of Science (PLoS)
Published
2013-05-01
Volume
8
Issue
5
Pages
e62817
ISSN
1932-6203 (ISSN-L)
DOI
10.1371/journal.pone.0062817
PMC
PMC3641112
OpenAlex
W2163249197
Cited by
304
References
71
Citations by year
2026: 8 · 2025: 15 · 2024: 13 · 2023: 38 · 2022: 34 · 2021: 14 · 2020: 24 · 2019: 16 · 2018: 33 · 2017: 25 · 2016: 20 · 2015: 26 · 2014: 26 · 2013: 11
FWCI
30.78
Type
journal article
Language
en
Update notices
correction (2017-02-21)
Open access
Licence
Locations
  • PLoS ONE (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, other-oa)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Digital Access to Scholarship at Harvard (DASH) (Harvard University) (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
  • Figshare (submittedVersion, cc-by)
Funders
  • Deutsches Krebsforschungszentrum
  • National Institutes of Health (RR 01032)
  • Centers for Disease Control and Prevention
  • U.S. Public Health Service (H75/CCH 123424)
  • National Center for Complementary and Alternative Medicine
  • National Center for Research Resources (RR 01032, M01-RR-01032, M01RR001032)
  • National Center for Chronic Disease Prevention and Health Promotion (R01DP000339)
  • National Center for Complementary and Integrative Health (R01AT006464)
Topics
  • Social Sciences › Psychology › Clinical Psychology › Mindfulness and Compassion Interventions
  • Health Sciences › Medicine › Complementary and alternative medicine › Biofield Effects and Biophysics
  • Social Sciences › Health › Religion, Spirituality, and Psychology
Keywords
relaxation response, meditation, gene expression, transcriptome, energy metabolism, mitochondrial function, insulin secretion, telomere maintenance, NF-κB signaling, inflammation, stress response, mitochondrial ATP synthase, peripheral blood, stress resilience
MeSH terms
Insulin Secretion, Cross-Sectional Studies, Energy Metabolism, Gene Expression Regulation, Humans, Immune System, Insulin, Mitochondria, Nitric Oxide, Prospective Studies, Relaxation, Religion, Stress, Physiological, Yoga, Inflammation Mediators, Meditation, Oligonucleotide Array Sequence Analysis, Gene Regulatory Networks, Transcriptome, Telomere Homeostasis

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review