SIDEBAR
»
S
I
D
E
B
A
R
«
Yoga e meditazione possono disattivare i geni dello stress
Maggio 29th, 2011 by admin

surya-namaskara-phusattna-asanaI numerosi studi, molti dei quali a lungo termine, hanno evidenziato come il benessere che deriva dallo yoga non sia esclusivamente un fatto mentale, ma dipenda dalla disattivazione di alcuni geni, proprio come ci spiega il dottor Herbert Benson, presidente emerito del Benson-Henry Institute for Mind-Body Medicine al Massachusetts General Hospital e professore associato di medicina all’Harvard Medical School. Per proteggerci dall’accumulo di energia negativa e disattivare i geni coinvolti nello stress, il dottor Benson consiglia pertanto venti minuti di rilassamento e meditazione ogni giorno.

Il benessere non è soltanto nella testa. Abbiamo scoperto che quando si evoca una risposta di rilassatezza gli stessi geni che sono attivati e disattivati dallo stress, vengono coinvolti nel processo e “spenti”.
La mente riesce a disattivare e ad attivare i geni e non è affatto separata dal corpo.
La mente riesce a disattivare e ad attivare i geni e non è affatto separata dal corpo.

Genomic Counter-Stress Changes Induced by the Relaxation ResponseJeffery A. Dusek1,2,3,6#, Hasan H. Otu3,4#, Ann L. Wohlhueter1, Manoj Bhasin3,4, Luiz F. Zerbini3,4, Marie G. Joseph4, Herbert Benson1,3,5*, Towia A. Libermann3,4*

1 Benson-Henry Institute for Mind Body Medicine at Massachusetts General Hospital, Chestnut Hill, Massachusetts, United States of America, 2 Department of Psychiatry, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, United States of America, 3 Department of Medicine, Division of Interdisciplinary Medicine and Biotechnology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, United States of America, 4 BIDMC Genomics Center, Beth Israel Deaconess Medical Center, Boston, Massachusetts, United States of America, 5 Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, United States of America, 6 Institute for Health and Healing, Abbott Northwestern Hospital, Minneapolis, Minnesota, United States of America

Citation: Dusek JA, Otu HH, Wohlhueter AL, Bhasin M, Zerbini LF, et al. (2008) Genomic Counter-Stress Changes Induced by the Relaxation Response. PLoS ONE 3(7): e2576. doi:10.1371/journal.pone.0002576

Editor: Philip Awadalla, University of Montreal, Canada

Received: February 12, 2008; Accepted: May 15, 2008; Published: July 2, 2008

Abstract

Background

Mind-body practices that elicit the relaxation response (RR) have been used worldwide for millennia to prevent and treat disease. The RR is characterized by decreased oxygen consumption, increased exhaled nitric oxide, and reduced psychological distress. It is believed to be the counterpart of the stress response that exhibits a distinct pattern of physiology and transcriptional profile. We hypothesized that RR elicitation results in characteristic gene expression changes that can be used to measure physiological responses elicited by the RR in an unbiased fashion.

Methods/Principal Findings

We assessed whole blood transcriptional profiles in 19 healthy, long-term practitioners of daily RR practice (group M), 19 healthy controls (group N1), and 20 N1 individuals who completed 8 weeks of RR training (group N2). 2209 genes were differentially expressed in group M relative to group N1 (p<0.05) and 1561 genes in group N2 compared to group N1 (p<0.05). Importantly, 433 (p<10−10) of 2209 and 1561 differentially expressed genes were shared among long-term (M) and short-term practitioners (N2). Gene ontology and gene set enrichment analyses revealed significant alterations in cellular metabolism, oxidative phosphorylation, generation of reactive oxygen species and response to oxidative stress in long-term and short-term practitioners of daily RR practice that may counteract cellular damage related to chronic psychological stress. A significant number of genes and pathways were confirmed in an independent validation set containing 5 N1 controls, 5 N2 short-term and 6 M long-term practitioners.

Conclusions/Significance

This study provides the first compelling evidence that the RR elicits specific gene expression changes in short-term and long-term practitioners. Our results suggest consistent and constitutive changes in gene expression resulting from RR may relate to long term physiological effects. Our study may stimulate new investigations into applying transcriptional profiling for accurately measuring RR and stress related responses in multiple disease settings.

Introduction

The relaxation response (RR) has been defined as a mind-body intervention that offsets the physiological effects caused by stress [1], [2]. The RR has been reported to be useful therapeutically (often as an adjunct to medical treatment) in numerous conditions that are caused or exacerbated by stress [3][6].

Mind-body approaches that elicit the RR include: various forms of meditation, repetitive prayer, yoga, tai chi, breathing exercises, progressive muscle relaxation, biofeedback, guided imagery and Qi Gong [7]. One way that the RR can be elicited is when individuals repeat a word, sound, phrase, prayer or focus on their breathing with a disregard of intrusive everyday thoughts [2]. The non-pharmacological benefit of the RR on stress reduction and other physiological as well as pathological parameters has attracted significant interest in recent years to decipher the physiological effects of the RR. In addition to decreased oxygen consumption [8][10], other consistent physiologic changes observed in long-term practitioners of RR techniques include decreased carbon dioxide elimination, reduced blood pressure, heart and respiration rate [1], [2], [11], prominent low frequency heart rate oscillations [12] and alterations in cortical and subcortical brain regions [13], [14].

Despite these observations and the well-established clinical effects of RR-eliciting practices [15], [16], the mechanisms underlying the RR have not been identified. Similarly, the impact of the RR on gene expression and signaling pathways has not yet been explored in detail, although a transcriptional profiling study of Qi Gong [17] practitioners, another RR method, revealed apparent distinct gene expression differences between Qi Gong practitioners and age matched controls. It is likely that differences in gene expression may be an underlying factor in the physiologic and psychologic changes noted above. Toward that end, we conducted a study to explore the gene expression profile of healthy long-term practitioners versus healthy age and gender matched controls. As a further evaluation, we provided 8-weeks of RR training to the control subjects and reassessed their gene expression.

Materials and Methods

Participants

Nineteen healthy practitioners of various RR eliciting techniques (including several types of meditation, Yoga, and repetitive prayer) participated (M group; n = 19). Years of practice averaged 9.4 years (5.0 sd) and ranged from 4 to 20 years. Twenty individuals without any prior RR eliciting experience served as controls (N group; n = 20).. As shown in Table 2, the M and N groups are matched with respect to age, gender, race, height, weight, and marital status, which do not exhibit significant difference between the groups (p>0.05, t- and chi-square test).

thumbnailTable 2. Demographics

doi:10.1371/journal.pone.0002576.t002

Protocol

The study protocol was approved by the Committee on Clinical Investigations at the Beth Israel Deaconess Medical Center (BIDMC), Boston MA. All subjects provided written informed consent and the study was conducted in the General Clinical Research Center (GCRC) of the BIDMC. After providing written informed consent, participants were screened by a physician and had blood drawn to ensure good health. All participants completed a testing session in the GCRC. N1 (novice) subjects had 8-weeks of RR training, listened to a 20-minute RR-eliciting CD daily and returned to the GCRC for a repeat testing session (hereafter classified as the N2 group).

Relaxation-Response Training

N subjects received 8 weeks of RR training. Training included information about reducing daily stress, and a 20-minute elicitation of the RR [35]. Subjects randomized to the RR group received 8 weekly individual RR-training sessions from an experienced clinician as per our manualized research protocol [35]. The first session provided an educational overview of the stress response and the RR, instructions on how to elicit the RR, and a 20-minute guided RR experience. Sessions 2 through 8 consisted of a review of the subject’s home practice card for compliance and a 20-minute guided RR experience.

During the RR elicitation in the weekly session, the subject was guided through a RR sequence including: diaphragmatic breathing, body scan, as well as mantra and mindfulness meditation, while subjects passively ignored intrusive thoughts. The specific CD guided the subject through the same sequence and has our clinical research studies and clinical practice for more than 15 years [35]. Subjects were asked to listen to the RR-eliciting CD once a day for 20 minutes at home.

To measure compliance, participants’ daily home practice logs were reviewed each week and at the end of the 8 week training. These logs indicate that N subjects listened for an average of 17.5 minutes per day (3.7 sd) over 8-weeks.

Microarray Analysis

Following previously described protocols, the transcriptional profile of samples were probed using Affymetrix HG-U 133 Plus 2.0 chips representing over 47,000 transcripts and variants using more than 54,000 probesets. Scanned image output files were visually examined for major chip defects and hybridization artifacts and then analyzed with Affymetrix GeneChip Microarray Analysis Suite 5.0 (MAS5) software (Affymetrix). The image from each chip was scaled such that the 2% trimmed mean intensity value for all arrays was adjusted to target intensity and reported as a non-negative quantity. Chips used for subsequent analysis consisted of 19 M, 19 N1 and 20 N2 samples (one chip from the N1 group had insufficient signal values). A hierarchical clustering technique was used to construct an Unweighted Pair Group Method with Arithmetic-mean (UPGMA) tree using Pearson’s correlation as the metric of similarity [36]. The expression data matrix was row-normalized for each gene prior to the application of average linkage clustering. When comparing 2 groups of samples to identify genes enriched in a given group, we used combination of three criteria. We considered genes with significantly different expression across the two groups using t-test (p<0.05) that further remained significant at a 5% false discovery rate (FDR) using permutation testing with 1,000 permutations [37], [38]. In order to finalize a set of genes significantly up-regulated in a given group compared to another group, among the genes that passed the aforementioned steps, we filtered the ones that are “present” in at least half of the samples in the enriched group using Affymetrix’ MAS5 Presence/Absence (P/A) calls. We used a paired t-test when comparing samples in groups N1 and N2.

Data deposition: All data sets have been deposited in the Gene Expression Omnibus, www.ncbi.nlm.nih.gov/geo (accession nos. GSE10041 and GSM253663-253734).

Gene Ontology and Gene Set Enrichment Analyses

Differentially expressed genes between the 3 groups (N2 vs. N1, M vs. N1 and M vs. N2) were separately analyzed using EASE to identify biologically relevant categories that are over-represented in the input set [21]. EASE analyses tested each list against all genes on the chip and overrepresentation describes a group of genes belonging to a certain GO category that appear more often in the given input list than expected to occur if the distribution were random. GO categories that had EASE scores of 0.05 or lower were selected as significantly over-represented. Gene Set Enrichment Analysis (GSEA 2.0 package http://www.broad.mit.edu/gsea/) was used to determine whether an a priori defined set of genes showed statistically significant, concordant differences between 2 groups (N2 vs. N1, and M vs. N1) in the context of known biological pathways. We tested expression values of all the genes in the relevant sample groups against 1687 gene sets obtained from the MSigDB2.0 for enrichment belonging to various metabolic pathways, chromosomal locations and functional sets (gene sets related to cancer/cancer cells are not included). The enriched gene sets have nominal p-value (NPV) less than 1% and False Discovery Rate (FDR) <50% after 100 random permutations. These criteria ensure that there is minimal chance of identifying false positives.

Supplementary Methods are located at http://bidmcgenomics.org/MIND_BODY_RR/(L​ogin: benson) (Password: test1).

Results

Patient characteristics

This study includes both cross sectional and an 8-week prospective design. Healthy adults were enrolled, comprising 2 groups: individuals with a long-term RR practice (group M; n = 19) or those with no prior RR experience (novice; group N1; n = 19). Group N1 novices, furthermore, underwent 8-weeks of RR training (Group N2; n = 20) for the prospective analysis. In the cross sectional study, we compare gene expression profiles (GEP) in whole blood between groups M and N1, whereas in the prospective study GEP is compared for each individual novice subject before and after RR experience, matched individuals of groups N1 versus N2 respectively.

Gene expression changes associated with the RR

Transcriptional differences between the different groups and within individuals before and after the RR are assessed by microarray analysis using Affymetrix HG-U133 Plus 2.0 genechips (www.affymetrix.com). This technology is a well established and reliable method to assess global gene expression differences [18]. Comparing group M (subjects with long term RR practice) to group N1 (subjects prior to RR training), we find statistically significant differential expression of 2209 genes; 1275 significantly up-regulated and 934 significantly down-regulated in M vs. N1. Additionally, 1561 genes are differentially expressed in novices after RR experience, N2 vs. N1; 874 significantly up-regulated and 687 significantly down-regulated. Comparison of gene lists from M vs. N1, N2 vs. N1 and M vs. N2 with Venn diagrams reveals significant overlap (Fig. 1a). Significance of overlaps is calculated using hypergeometric distributional assumption [19] and p-values are adjusted using Bonferroni correction for multiple comparisons [20]. Heatmaps were generated from genes in the intersecting areas of the Venn diagrams (Fig. 1b). We find 316 up-regulated and 279 down-regulated genes are differentially expressed in group M compared to both group N1 and N2; these changes in GEP are only observed in long-term RR practitioners. Similarly, 260 genes are up-regulated and 168 genes are down-regulated in both groups M and N2 compared to N1; they represent GEP changes characteristic of RR practice over at least 8 weeks.

thumbnailFigure 1. Gene Ontology Analysis.

Analysis of differentially expressed genes: a) Venn diagrams: * indicates significant overlaps (p<106); b) Heatmaps of the 595 differentially regulated genes in both M vs. N1 and M vs. N2 (left) and the 418 differentially regulated genes in both M vs. N1 and N2 vs. N1; c) Heatmap of 15 genes in the intersection of all three groups (gene symbols listed on the right). In heatmaps, rows represent genes and columns represent samples from N1, N2, and M groups. Genes are clustered using row-normalized signals and mapped to the [−1,1] interval (shown in scales beneath each heatmap). Red and green represent high and low expression values, respectively.

doi:10.1371/journal.pone.0002576.g001

Heatmaps generated using these genes exhibit consistent GEP changes across the three groups with a few samples in each group resembling the GEP of another group. To determine if any demographic characteristics (e.g. age, ethnicity , etc.) influences this observation, we clustered each group separately using the same set of genes. For each cluster analysis, we calculated the significance of observing a characteristic among the samples in the subgroups formed. We found that number of times M subjects reported eliciting the RR per week was significantly associated with the subgroups formed when M samples were clustered using genes differentially expressed in long term RR practitioners only. Specifically, there were 316 up-regulated and 279 down-regulated genes differentially expressed in group M compared to both group N1 and N2; (Fig. 1b). All remaining cluster analyses revealed no such significant influence of demographic characteristics (see online supplementary data).

Finally, the intersection of all 3 areas (M vs. N1, N2 vs. N1 and M vs. N2) identifies genes with expression behavior that is monotonically changed between N1 to N2 to M (Fig. 1c). These results clearly demonstrate that short term as well as long term RR practice lead to distinct and consistent gene expression changes in hematopoietic cells.

Signaling pathways modulated by the RR

We performed Expression Analysis Systematic Explorer (EASE) analysis [21] using M vs. N1, and N2 vs. N1 data-sets, to identify Gene Ontology (GO) categories where specific genes in these data occur more often than would be expected by random distribution of genes. These findings (and those of the validation data-set below) are summarized in Table 1, where select over-represented GO categories are listed along with specific genes differentially expressed in our data-sets. These categories include oxidative phosphorylation, ubiquitin-dependent protein catabolism, nuclear messenger RNA (mRNA) splicing, ribosomes, metabolic processes, regulation of apoptosis, NF-κB pathways, cysteine-type endo-peptidase activity and antigen processing. Most are significant in both long-term (M vs. N1) and short-term (N2 vs. N1) practitioners of daily RR practice (see Table).

thumbnailTable 1. Gene Ontology Categories

doi:10.1371/journal.pone.0002576.t001

Even though our analyses of differentially expressed genes and GO categories associated with RR practice meet widely accepted criteria for statistical significance, we were concerned about the relatively small fold changes that were observed (see Supplementary Methods). To address this issue we employed Gene Set Enrichment Analysis (GSEA). GSEA has proven to be useful for capturing subtle expression changes in complex gene signatures based on predefined gene sets or pathways [22]. As described above, we examined expression data for 2 comparisons, M vs. N1 and N2 vs. N1. The selected pathways or gene sets that are significantly enriched (False Discovery Rate (FDR)<50%, nominal p-value (NPV)< = 0.02) are shown in Figure 2, with gene sets for N2 vs. N1 and M vs. N1 in Fig 2A and Fig 2B respectively.

thumbnailFigure 2. GSEA Analysis.

The analysis has been performed for >1200 predefined datasets using GSEA 2.0 software. Signal values for each gene are obtained by collapsing the probe values using max_probe algorithm. Representative datasets, significantly enriched (FDR<50%, or NPV< = 0.01) between any two groups and corresponding heatmaps (depicting relative gene expression changes of core enrichment) are shown in a) N2 vs. N1 and b) M vs. N1. Datasets that are enriched in both the original and validation analyses are marked with *. c) Heatmaps of ribosomal proteins and ubiquitin mediated proteolysis illustrate transitional trends in gene expression across the N1, N2 and M groups.

doi:10.1371/journal.pone.0002576.g002

GSEA analysis of N2 vs. N1 showed highly significant enrichment in gene sets related to various cellular stressors/stress responses and metabolism. To a pronounced degree these observations complement the results of GO analysis presented in the Table, also depicting significant alterations in cellular response to stress, oxidative and primary metabolism. The transition effect of the RR from novice to short term (8 weeks) to long term RR practitioners has been denoted through a colorgram of ribosomal proteins and ubiquitin mediated proteolysis gene sets (Fig. 2C). Whereas expression of ribosomal genes is significantly upregulated in RR practitioners at 8 weeks and more pronounced in long term practitioners, ubiquitin mediated proteolysis gene expression in general shows an opposite trend. Closer inspection of the colorgrams for ribosomal proteins and Ub proteolysis gene sets shows some variation in the GEP in each subgroup (N1, N2 M). The GEP of a few N1 or N2 subjects resembles the GEP of M subjects and vice versa. To elucidate the association between GEP and subject characteristics (Race, Age, etc), we performed clustering of each subgroup separately (N1, M) using the enriched gene sets (Ribosomal and Ub Proteolysis). This analysis identified a subcluster in the N2 subgroup that has significant over-representation of Asian subjects (P value <0.05) when the clustering was performed using the ribosomal protein gene set. This observation needs further validation on a larger dataset as the current study contains only five Asian subjects. No other characteristic exhibited significant association with the ribosomal protein gene set. No significant association between the GEP profiles and subject characteristic was found when clustering was performed using the Ub proteolysis gene set. This analysis provides further insight into the stress response related genes that are influenced by RR practice.

Independent validation set analysis

As a validation of our results, we repeated the experimental and analysis procedures defined in the “Methods” section on a new set of samples consisting of 5 N1, 5 N2 and 6 M subjects. We found 1846 and 2390 probe sets differentially expressed between M vs. N1, and N2 vs. N1 groups. The validation data-set showed a significant (p<10−5) number of genes in common with the original analysis of 58 samples. We also found that 70–75% of all GO categories from the original analysis were retained in the validation set (p ~0), and 30–65% of significantly over-represented GO categories were shared. Of note, biologically relevant GO categories such as oxidative phosphorylation, regulation of apoptosis, and antigen presentation, come up as significantly over-represented in both the original and validation analyses. Results of the validation set and comparison analyses can be found in the online supplementary data. In addition, validation GSEA analysis on N2 vs. N1 subjects shows enrichment of ribosomal proteins and platelet expressed gene sets and enrichment of ribosomal proteins, oxidative phosphorylation and electron transport chain gene sets in M vs. N1 subjects (Fig. 2A and 2B). The similarities between the original and validation results from GSEA analysis argues against random chance accounting for the observed enrichment of these gene sets.

Discussion

Results from our study indicate that there are distinct differences in the GEPs between individuals with many years of RR practice (group M) and those without such experience (group N1). Furthermore we find significant GEP changes within the same individuals before (N1) and after 8 weeks of RR training (N2). Finally, the changes in GEP found in M vs. N1, and those of N2 vs. N1, are to a great degree similar when assessed by analysis of differentially expressed genes, GO analysis and GSEA.

It is becoming increasingly clear that psychosocial stress can manifest as system-wide perturbations of cellular processes, generally increasing oxidative stress and promoting a pro-inflammatory milieu [23][25]. Stress associated changes in peripheral blood leukocyte expression of single genes have been identified [26][28]. More recently, chronic psychosocial stress has been associated with accelerated aging at the cellular level. Specifically, shortened telomeres, low telomerase activity, decreased anti-oxidant capacity and increased oxidative stress are correlated with increased psychosocial stress [29] and with increased vulnerability to a variety of disease states [30]. Stress-related changes in GEP have been demonstrated by microarray analysis in healthy subjects, including up-regulation of several cytokines/chemokines and their receptors [31], and in individuals suffering from post-traumatic stress disorder, including inflammation, apoptosis and stress response [32] as well as metabolism and RNA processing pathways [33]. The pro-inflammatory transcription factor NF-kappa B (NF-κB) which is activated by psychosocial stress has been identified as a potential link between stress and oxidative cellular activation [34].

The RR is clinically effective for ameliorating symptoms in a variety of stress-related disorders including cardiovascular, autoimmune and other inflammatory conditions and pain [15]. We hypothesize that RR elicitation is associated with systemic gene expression changes in molecular and biochemical pathways involved in cellular metabolism, oxidative phosphorylation/generation of reactive oxygen species and response to oxidative stress and that these changes to some degree serve to ameliorate the negative impact of stress. Genome-wide evaluation of PBL GEP is a reasonable approach to survey the transcriptional changes that are involved in elicitation of the RR. The GEP of RR practitioners presented here reveals altered gene expression in specific functional groups which suggest a greater capacity to respond to oxidative stress and the associated cellular damage. Genes including COX7B, UQCRB and CASP2 change in opposite direction from that in the stress response [31], [32].

Our findings are relatively consistent with those found in a study of Qi Gong [17], a practice that elicits the RR. In their study of 6 Qi Gong practitioners and 6 aged matched controls, practitioners had down-regulation of ubiquitin, proteasome, ribosomal protein and stress response genes and mixed up- and down-regulation of genes involved in apoptosis and immune function. We find a similar pattern of GO categories that are significantly over-represented in GO or enriched in GSEA in our cross sectional comparison, M vs. N1. However, in our data-set ribosomal proteins were up-regulated.

Overall, similar genomic pattern changes occurred in practitioners of a specific mind body technique (Qi Gong) as well as in our long-term practitioners who utilized different RR practices including Vipassna, mantra, mindfulness or transcendental meditation, breath focus, Kripalu or Kundalini Yoga, and repetitive prayer. This indicates there is a common RR state regardless of the techniques used to elicit it.

Our study is the first to prospectively evaluate GEP changes in individuals before and after a short-term (8 week) RR training which consequently enables an appreciation of the parallel GEP changes that occur with short- and/or long-term RR practice. Replications in larger cohorts are warranted. Future investigations could better define the therapeutic value and required duration of RR training to counter stress-related disorders.

References

  • Wallace RK, Benson H, Wilson AF (1971) A wakeful hypometabolic physiologic state. Am J Physiol 221: 795–799. Find this article online

  • Sternberg EM (1997) Emotions and disease: from balance of humors to balance of molecules. Nat Med 3: 264–267. Find this article online

  • Benson H, Kotch JB, Crassweller KD (1978) Stress and hypertension: interrelations and management. Cardiovasc Clin 9: 113–124. Find this article online

  • Nakao M, Myers P, Fricchione G, Zuttermeister PC, Barsky AJ, et al. (2001) Somatization and symptom reduction through a behavioral medicine intervention in a mind/body medicine clinic. Behav Med 26: 169–176. Find this article online

  • Benson H, Goodale IL (1981) The relaxation response: your inborn capacity to counteract the harmful effects of stress. J Fla Med Assoc 68: 265–267. Find this article online

  • Benson H (1983) The relaxation response: Its subjective and objective historical precedents and physiology. TINS 6: 281–284. Find this article online

  • Benson H, Steinert RF, Greenwood MM, Klemchuk HM, Peterson NH (1975) Continuous measurement of O2 consumption and CO2 elimination during a wakeful hypometabolic state. J Human Stress 1: 37–44. Find this article online

  • Kesterson J, Clinch NF (1989) Metabolic rate, respiratory exchange ratio, and apneas during meditation. Am J Physiol 256: R632–638. Find this article online

  • Warrenburg S, Pagano RR, Woods M, Hlastala M (1980) A comparison of somatic relaxation and EEG activity in classical progressive relaxation and transcendental meditation. J Behav Med 3: 73–93. Find this article online

  • Beary JF, Benson H (1974) A simple psychophysiologic technique which elicits the hypometabolic changes of the relaxation response. Psychosom Med 36: 115–120. Find this article online

  • Peng CK, Henry IC, Mietus JE, Hausdorff JM, Khalsa G, et al. (2004) Heart rate dynamics during three forms of meditation. Int J Cardiol 95: 19–27. Find this article online

  • Lazar SW, Bush G, Gollub RL, Fricchione GL, Khalsa G, et al. (2000) Functional brain mapping of the relaxation response and meditation. Neuroreport 11: 1581–1585. Find this article online

  • Jacobs GD, Benson H, Friedman R (1996) Topographic EEG mapping of the relaxation response. Biofeedback Self Regul 21: 121–129. Find this article online

  • Astin JA, Shapiro SL, Eisenberg DM, Forys KL (2003) Mind-body medicine: state of the science, implications for practice. J Am Board Fam Pract 16: 131–147. Find this article online

  • Esch T, Fricchione GL, Stefano GB (2003) The therapeutic use of the relaxation response in stress-related diseases. Med Sci Monit 9: RA23–34. Find this article online

  • Li QZ, Li P, Garcia GE, Johnson RJ, Feng L (2005) Genomic profiling of neutrophil transcripts in Asian Qigong practitioners: a pilot study in gene regulation by mind-body interaction. J Altern Complement Med 11: 29–39. Find this article online

  • Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA, et al. (2002) A stem cell molecular signature. Science 298: 601–604. Find this article online

  • Hosack DA, Dennis G Jr., Sherman BT, Lane HC, Lempicki RA (2003) Identifying biological themes within lists of genes with EASE. Genome Biol 4: R70. Find this article online

  • Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550. Find this article online

  • Irie M, Asami S, Nagata S, Miyata M, Kasai H (2002) Psychological mediation of a type of oxidative DNA damage, 8-hydroxydeoxyguanosine, in peripheral blood leukocytes of non-smoking and non-drinking workers. Psychother Psychosom 71: 90–96. Find this article online

  • Yamaguchi T, Shioji I, Sugimoto A, Yamaoka M (2002) Psychological stress increases bilirubin metabolites in human urine. Biochem Biophys Res Commun 293: 517–520. Find this article online

  • Zheng KC, Ariizumi M (2007) Modulations of immune functions and oxidative status induced by noise stress. J Occup Health 49: 32–38. Find this article online

  • Glaser R, Kennedy S, Lafuse WP, Bonneau RH, Speicher C, et al. (1990) Psychological stress-induced modulation of interleukin 2 receptor gene expression and interleukin 2 production in peripheral blood leukocytes. Arch Gen Psychiatry 47: 707–712. Find this article online

  • Glaser R, Lafuse WP, Bonneau RH, Atkinson C, Kiecolt-Glaser JK (1993) Stress-associated modulation of proto-oncogene expression in human peripheral blood leukocytes. Behav Neurosci 107: 525–529. Find this article online

  • Platt JE, He X, Tang D, Slater J, Goldstein M (1995) C-fos expression in vivo in human lymphocytes in response to stress. Prog Neuropsychopharmacol Biol Psychiatry 19: 65–74. Find this article online

  • Epel ES, Blackburn EH, Lin J, Dhabhar FS, Adler NE, et al. (2004) Accelerated telomere shortening in response to life stress. Proc Natl Acad Sci U S A 101: 17312–17315. Find this article online

  • Epel ES, Lin J, Wilhelm FH, Wolkowitz OM, Cawthon R, et al. (2006) Cell aging in relation to stress arousal and cardiovascular disease risk factors. Psychoneuroendocrinology 31: 277–287. Find this article online

  • Morita K, Saito T, Ohta M, Ohmori T, Kawai K, et al. (2005) Expression analysis of psychological stress-associated genes in peripheral blood leukocytes. Neurosci Lett 381: 57–62. Find this article online

  • Zieker J, Zieker D, Jatzko A, Dietzsch J, Nieselt K, et al. (2007) Differential gene expression in peripheral blood of patients suffering from post-traumatic stress disorder. Mol Psychiatry 12: 116–118. Find this article online

  • Segman RH, Shefi N, Goltser-Dubner T, Friedman N, Kaminski N, et al. (2005) Peripheral blood mononuclear cell gene expression profiles identify emergent post-traumatic stress disorder among trauma survivors. Mol Psychiatry 10: 500–513, 425. Find this article online

  • Bierhaus A, Wolf J, Andrassy M, Rohleder N, Humpert PM, et al. (2003) A mechanism converting psychosocial stress into mononuclear cell activation. Proc Natl Acad Sci U S A 100: 1920–1925. Find this article online

  • Dusek JA, Chang BH, Zaki J, Lazar S, Deykin A, et al. (2006) Association between oxygen consumption and nitric oxide production during the relaxation response. Med Sci Monit 12: CR1–10. Find this article online

  • Sneath PHA (1973) Numerical taxonomy; the principles and practice of numerical classification (W. H. Freeman, San Francisco, CA).

  • Li C, Wong WH (2001) Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA 98: 31–36. Find this article online

  • Tusher VG, Tibshirani R, Chu G (2001) Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA 98: 5116–5121. Find this article online

http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0002576#s1


Comments are closed

»  Substance:WordPress   »  Style:Ahren Ahimsa