Abstract

Posttraumatic stress disorder (PTSD) is characterized by the occurrence of a traumatic event that is beyond the normal range of human experience. The future of PTSD treatment may specifically target the molecular mechanisms of PTSD. In the US, approximately 20% of adults report taking herbal products to treat medical illnesses. L-theanine is the amino acid in green tea primarily responsible for relaxation effects. No studies have evaluated the potential therapeutic properties of herbal medications on gene expression in PTSD. We evaluated gene expression in PTSD-induced changes in the amygdala and hippocampus of Sprague-Dawley rats. The rats were assigned to PTSD-stressed and nonstressed groups that received either saline, midazolam, L-theanine, or L-theanine + midazolam. Amygdala and hippocampus tissue samples were analyzed for changes in gene expression. One-way ANOVA was used to detect significant difference between groups in the amygdala and hippocampus. Of 88 genes examined, 17 had a large effect size greater than 0.138. Of these, 3 genes in the hippocampus and 5 genes in the amygdala were considered significant () between the groups. RT-PCR analysis revealed significant changes between groups in several genes implicated in a variety of disorders ranging from PTSD, anxiety, mood disorders, and substance dependence.

1. Background

Posttraumatic stress disorder (PTSD) is characterized by the occurrence of a traumatic event that is beyond the normal range of human experience, for example, rape, torture, genocide, and war zone stress [1]. The range of symptoms of PTSD necessitates many different therapeutic approaches to treating the disorder. The most successful treatments are cognitive-behavioral therapy (CBT) and medications [1]. The selective serotonin reuptake inhibitors (SSRIs), sertraline (Zoloft) and paroxetine (Paxil), are approved by the FDA for treatment of PTSD [1]. The future of PTSD treatment may include medications specifically targeted for the molecular mechanisms responsible for PTSD symptoms. Zhang and colleagues used microarray technology to compare the gene expression profile of stressed maladapted rats (PTSD) and stressed well adapted rats in the areas of the brain associated with PTSD, the hippocampus, and the amygdala [2]. In the hippocampus they found about 110 genes were upregulated by a two-fold change or greater, and about 240 genes were downregulated by less than a half-fold change [2]. They found that multiple two signaling pathways were upregulated and multiple one signaling pathways were downregulated in the amygdala [2].

In the US, nearly 1 in 5 adults reports taking herbal products to treat a medical illness; in fact, as recently as 1980, the US Pharmacopeia consisted of 59% herbal medicines [3]. Green tea is an herbal product that has been consumed since ancient times to promote relaxation. L-theanine (marketed by the trade name Suntheanine) is the amino acid in green tea primarily responsible for its flavor and pharmacologic effect [4]. It was discovered in 1949 and is found almost solely in green tea leaves [4]. In a small-scale human study by Junega and colleagues, Suntheanine was found to promote the generation of alpha brain waves, indicative of a relaxed state. L-theanine is a derivative of L-glutamic acid, which is one of the major excitatory neurotransmitters in the brain [5].

Although little is known about the pharmacokinetics of L-theanine in humans, animal studies show that it is absorbed in the small intestine via sodium-coupled active transport, crosses the blood brain barrier via a leucine-preferring transport system, and is metabolized by the kidneys [5]. In rats, an intraperitoneal injection of L-theanine is taken up by the brain within thirty minutes without metabolic changes [4].

In addition to the study by Junega and colleagues, other researchers have investigated the ability of L-theanine to attenuate anxiety. Lu and colleagues compared the acute effects of L-theanine, alprazolam (a prototype anxiolytic), and a placebo on the behavioral measures of anxiety in healthy human subjects using an anticipatory anxiety model. Lu and colleagues found that L-theanine may have a relaxing effect under resting conditions, but neither L-theanine nor alprazolam had anxiolytic effects in the anticipatory anxiety model [6]. However, in a rat study by Heese et al., L-theanine concomitantly administered with midazolam was shown to have an interactive anxiolytic effect, and a decrease in fine and basic motor movements in a rodent model [7].

No studies have evaluated the potential therapeutic properties of herbal medications in PTSD or changes in gene expression. Because PTSD is categorized as an anxiety or nervous disorder, it is important to study the potential benefits of L-theanine in the treatment of PTSD and the molecular pathways responsible. The purpose of this study was to investigate if L-theanine modulates gene expression in the hippocampus and amygdala, regions of the brain that have been identified as critically important areas involved in the pathogenesis of PTSD [810]. This study was part of a larger study that also investigated neurobehavioral changes in this rodent model. This current investigation provides a beginning framework to pursue targeted treatments of PTSD by identifying genes involved in the neurobiology of this devastating condition.

2. Methods

Eighty male Sprague-Dawley rats (Harlan Sprague-Dawley Laboratories) weighing 300–350 grams were used. They were housed in groups of three in polycarbonate “shoebox-sized” cages lined with bedding. The animals went through a 7-day adaptation period in a temperature-controlled environment (°C, 60% humidity) with a reverse light-dark cycle where they received 12 hours of light (6:00 PM to 06:00 AM) and 12 hours of darkness (6:00 AM to 06:00 PM). They were allowed free access to food and water. The animals were only handled for weighing, drug administration, and cleaning of cages and were naïve to all instruments. After acclimation, the investigators measured baseline body weight. Statistical analysis demonstrated equivalency between the control and PTSD-stressed groups. The use of laboratory rats in this protocol was in accordance with the NIH Guide for the Care and Use of Laboratory Animals and received approval from the Institutional Animal Care and Use Committee (IACUC) at the US Army Institute of Surgical Research, San Antonio, Texas.

The animals were randomly assigned to one of eight groups by use of a computer number generator. Rats were assigned to the nonstressed groups or the 3-day restraint shock PTSD rodent model groups. There were four groups in the nonstressed rats: control vehicle (C-V), L-theanine (C-L), midazolam (C-M), and L-theanine + midazolam (C-M + L). There were also four groups in the PTSD-induced rats: control vehicle (P-V), L-theanine (P-L), midazolam (P-M), and L-theanine + midazolam (P-M + L). Each animal received a subcutaneous injection of the following: (1) control (saline); (2) L-theanine 10 mg/kg [11]; (3) midazolam (Versed) 1.5 mg/kg; or (4) L-theanine 10 mg/kg + midazolam (Versed) 1.5 mg/kg (to evaluate interaction effects). The vehicle (saline) group served as a control. All animals received equivalent volumes of subcutaneous injection to total 2 mL/kg body weight. Additionally, all experimentation occurred on a timed schedule between 8 AM and 4 PM over 4 consecutive days to control for the circadian rhythm of the animals.

2.1. PTSD Stress Model

Stress exposure consisted of 2-hour sessions per day of immobilization and tail-shocks for 3 consecutive days. The animals were stressed in the morning (between 08:00 AM and 12:00 PM). They were restrained in a plexiglass tube, and 40 electric shocks (2 mA, 3-second duration) were applied at varying intervals (140–180 seconds). Commercial hardware and software controlled the timing and amplitude of the shock (Precision Animal Shocker, Coulbourn Instruments, Columbus, Ohio, USA). Rats were stressed for 3 consecutive days. The rationale for the approach was because it had been previously demonstrated that repeated stress sessions for 3 days were more effective than a single stress session in producing sustained physiological and behavioral abnormalities, such as exaggerated acoustic startle response and reduced body weight [12].

Several laboratory studies have established and tested the inescapable tail-shock model of stress in rats. These tests have verified the neurobiological and behavioral alterations induced in the rat model and are similar to those found in human PTSD patients [12, 13]. Rats were stressed for three consecutive days based on previous literature which indicated repeated stress sessions were more effective than a single stress session in producing physiological and behavioral abnormalities, including exaggerated ASR and reduced body weight (see Table 1) [12].

2.2. Gene Analysis

After the completion of neurobehavioral tests, the investigators anesthetized the male Sprague-Dawley rats using isoflurane in a bell jar, and the calvarium was removed exposing the whole brain. The whole brain was removed intactly and placed immediately on a brain block on ice [14]. Coronal slices (200 μm thickness) in the vicinity of interaural and bregma coordinates 6.96 mm and −2.04, respectively, were dissected using a brain block (see red shaded region in Figure 1(a)). Circular punches (80–100 μm diameter) were retrieved from the left and right amygdala (see lower left red circle in area 6 of Figure 1(b)). Following the dissection of the amygdala, the investigator grossly dissected the bilateral hippocampi (see upper medial red outlined area in area 6 of Figure 1(b)). Tissue was snap-frozen in liquid nitrogen, placed in prelabeled eppendorf tubes, and packed in dry ice for shipping to the QIAGEN Service Core for Genomics and Gene Expression in Frederick, MD.

RNA was isolated using the QIAGEN RNEasy Mini Kit (Cat. no. 74104) following the manufacturer’s protocol for tissue. RNA quality was determined using the Agilent Bioanalyzer (Agilent) with RNA 6000 Nano Kits (Agilent, Cat. no. 5067-1511). Total RNA yield, 260/280, and 260/230 ratios were measured using a NanoDrop spectrophotometer (Thermo). A reverse transcription reaction using 500 ng of total RNA was completed using the QIAGEN RT2 First Strand Kit (QIAGEN, Cat. no. 330401). cDNA samples were assayed using a modified QIAGEN RT2 PCR array (Cat. no. PARN-60) according to manufacturer’s instructions [15, 16]. This array, containing 84 assays related to rat neurotransmitter receptors and regulators, was modified to include four additional assays of genes from QIAGEN. These genes are related to neurotransmission pathways that have been identified in previous PTSD studies: p11 (S100a10) (Cat. no. PPR06766) [17]; 5H receptor (Htr2a) (Cat. no. PPR06850) [18]; alpha1 adrenergic receptor (Adra1a) (Cat. no. PPR43329) [1921]; and Egr1 (Zif/268) (Cat. no. PPR44272) an effector immediate early gene [15, 16].

2.3. Rat Neurotransmitter Receptors and Regulators PCR Array

The rat neurotransmitter receptors and regulators RT2 profiler PCR array profiles the expression of 84 genes involved in modulating the biological processes of neurotransmitter biosynthesis, uptake, transport, and signaling through neurotransmitter receptors. This array contains receptors for specific neurotransmitters, such as acetylcholine, benzodiazepine, dopamine, gamma-aminobutyric acid (GABA), glutamate, serotonin, somatostatin, and neuropeptides. Genes involved in the regulation of neurotransmitter levels were included as well. Analysis of the expression of a focused panel of genes related to neuronal systems with this array can be performed using real-time PCR.

2.4. Statistical Analyses

For this cross-sectional, randomized, prospective study, a one-way ANOVA was performed to compare the 8 groups, of which each unit was randomly allocated. All assumptions were examined (e.g., homogeneity of variance, univariate normality, etc.) and in the event that the assumptions are not tenable remedial (e.g., transformation) or alternative (e.g., nonparametric tests, such as Kruskal-Wallis) strategies were considered. Multiple comparison procedures (MCP) were performed in the event that significance was obtained, with the Tukey post hoc test. A variance explained statistic eta-squared (η2) was used as the reported effect size, and though ascertaining what constitutes a small/medium/large effect is primarily a function of the context/discipline, for this report .01/.059/.138 was small/medium/large [22]. All data were statistically analyzed based on gene cycle thresholds (CTs) normalized to five housekeeping genes per QIAGEN procedures [15, 16].

3. Results

This study was designed to investigate the effects of L-theanine on PTSD-induced changes in rat brain gene expression. Eighty-eight genes were investigated in this study. Data analysis showed a number of significant differences in gene expression in both the hippocampus and amygdala. Volcano plots were developed to compare the gene expression of the 40 nonstressed to the 40 PTSD-stressed rats for both the hippocampus and amygdala as demonstrated in Figure 2. The volcano plots graphically display the relationship of –log 10 value and the Log 2 fold change for the combined treatment and control groups for all of the genes. For the 8-group one-way ANOVA, a large effect size was chosen to serve as a springboard for discussion and interpretation of relative between-group differences on gene expression.

3.1. Hippocampus

In the hippocampus, the genes with an effect size > 0.138 are displayed in Figure 3. In accordance with Figure 3, Table 2 summarizes between-group changes in gene expression in the hippocampus. The Egr1 gene was found to be significant and between-group differences were obtained: , . There was one significant post hoc test with the P-M group having a significantly higher mean than the P-V group . A significant difference was found between groups for the Maoa gene: , . The post hoc test found one significant pairwise comparison: the C-M + L group had the highest mean and the P-V group had the lowest . The Anxa9 gene was found to have significance between groups: , . Multiple comparison tests resulted in one significant pairwise result: the P-M group had a significantly higher mean than the P-V group .

The following genes in the hippocampus did not show significant differences between groups; however, due to the large effect size (variance explained) these genes are presented in Table 2. The S100a10 gene did not have significant between-group differences: , . The P-M group had the highest mean and the P-V group had the lowest . Significant between-group differences were not obtained for the Gabrb2 gene: , . The P-M group had the highest mean and the P-V group had the lowest . Finally, the Chrm2 gene did not show significance between groups. A statistically significant difference was not obtained: , . Though not significant, the C-M + L group had the highest mean and the P-V group had the lowest .

3.2. Amygdala

Genes in the amygdala with an effect size greater than 0.138 are displayed in Figure 4. In accordance with Figure 4, Table 3 summarizes between-group differences in gene expression within the amygdala. The Gabra4 gene was , . Significant post hoc pairwise comparisons were such that the C-L group had a higher mean than each of the following: P-L , (2) P-M , and C-M group (M = 3.09). A significant difference for the Slc5a7 gene was found: , . The C-V group had the highest mean and the P-L group had the lowest ; however as per the post hoc test, this pairwise comparison was not significant. Significant between-group differences were also found for the Drd2 gene: , . There was a significant difference between P-V and the P-M group . The Drd1a gene was also significant: , . A significant difference was found between the C-M + L and the P-M group (M = 4.83). Lastly, the Glra2 showed significance: , . There was one significant post hoc test between the C-M + L group with the higher mean and the P-M group with the lower mean .

Again the following genes were not found to be significant; however, based on the large effect size (variability) used in the study these genes are included in Table 3. No difference was found for the Chrna6 gene: , . The C-M group had the highest mean and the C-M + L group had the lowest . The Htr3a gene showed no significant difference between the groups: , . The C-L group had the highest mean and the P-M + L group had the lowest . The Chrnb4 also showed no significant difference: , . The C-V group had the highest mean and the C-M + L group had the lowest . No significant difference was found for the Gabrd gene between the groups: , . The C-M + L group had the highest mean and the P-M group had the lowest . The Prima1 gene did not have a significant difference between the groups: , . The P-V group had the highest mean and the P-M group had the lowest . No difference between the groups was found for the Chrna3 gene: , . The P-L group had the highest mean and the C-M + L group had the lowest .

4. Discussion

L-theanine has previously been demonstrated to produce anxiolytic effects alone or in combination with concomitant administered benzodiazepines [7]. No studies have examined the administration of L-theanine’s effects on PTSD-induced changes in brain gene expression. We explored the modulation of gene expression in PTSD-induced changes in the hippocampus and amygdala and the effects of L-theanine and midazolam on these changes.

4.1. Hippocampus

We found a downregulation of the Egr1 gene in the PTSD-vehicle group compared to the PTSD-midazolam group. Egr1 is an early growth response nuclear transcription factor that plays a role in regulating the expression of genes required for learning and memory [23, 24]. In particular, previous studies have demonstrated a role for Egr1 in hippocampal long-term potentiation, a process critical to learning and memory [25]. Expression of the Maoa gene was decreased in the PTSD-vehicle group compared to the control group that received L-theanine and midazolam.

Decreased Maoa enzymatic activity leads to increased concentrations of its substrate biogenic amines, including norepinephrine and serotonin. As reviewed by Bortolato and colleagues, previous research has demonstrated strong links between decreased Maoa activity and a number of reactive social behaviors, including specific links to violence and aggression [26]. Similarly, decreased expression of Maoa in humans predicted, among other changes, amygdala and hippocampus hyperreactivity during aversive recall [27].

Although no significant difference was found between groups, the PTSD-vehicle group had the lowest mean expression of the S100a10 gene and the PTSD-midazolam group the highest. The S100a10 gene and its associated protein p11 are important in the regulation of serotonin receptors [28], and decreases in S100a10 gene expression have been linked to PTSD, major depressive disorder, and bipolar disorder and proposed as a likely biomarker for suicidality [29, 30].

A downregulation of the Gabrb2 gene in the PTSD-vehicle group compared to the PTSD-midazolam group was similar to Zhao and colleagues findings of a significantly decreased expression of the Gabrb2 gene in the anterior cingulate cortex in subjects suffering from major depressive disorder or bipolar disorder [31].

A downregulation of the Anxa9 gene was found in the PTSD-vehicle group compared to the PTSD-midazolam group. The Anxa9 gene, part of the calcium-dependent phospholipid-binding annexin family, may play a role in suppression of apoptosis and is currently under investigation as a therapeutic target for the treatment of breast cancer [32]. In a study using the Flinders rat model, Chrm2 gene expression levels were found to be significantly lower in the Flinders depression sensitive line [33].

4.2. Amygdala

According to Ponomarev and colleagues, downregulation of the Gabra4 gene can be correlated with stress enhanced fear learning (SEFL) in rats. These data suggest a significant role of GABA signaling in SEFL. The method of inducing SEFL was similar to the method of inducing PTSD in our study [34]. Interestingly, our findings suggest an upregulation of this gene in the nonstressed-L-theanine group compared to PTSD-L-theanine, PTSD-midazolam, and nonstressed-midazolam groups.

PCR analysis of tissue samples from the amygdala showed a 2.28-fold increase in Drd2 gene expression in the PTSD-vehicle group compared to the PTSD-midazolam group, which is inconsistent with the findings of a study that investigated the effects of maternal deprivation and chronic mild stress in rats [35]. Zhu and colleagues showed that mRNA levels of the Drd2 gene were significantly decreased in all stressed groups compared to controls, which enabled the investigators to correlate a downregulation of the Drd2 gene with depressive phenotypes [35].

Drd1a genes were significantly downregulated in the PTSD-midazolam group compared to the nonstressed-midazolam + L-theanine group concurrent with the results of another study exploring the effects of stress on the mouse prefrontal cortex transcriptome, demonstrating a deregulation of the Drd1a gene in chronic mild stress exposed subjects [36].

The SLc5a7 gene was significantly downregulated in the PTSD-L-theanine group compared to the nonstressed-vehicle group. In comparison the SLc5a7 gene was found to be upregulated in tasks that required higher levels of attention, and it was found that decreased levels of cholinergic activity are involved in ADHD [37].

Prima1 DNA methylation increased by 10–15% in subjects suffering from major depressive disorder [38] which interestingly shows some correlation with our finding of the highest level of Prima1 gene expression in the PTSD-vehicle group, which is a 0.69-fold increase in expression compared with the PTSD-midazolam group with the lowest level of expression.

Although no studies have investigated gene expression changes of Chrnb4 in PTSD or anxiety, variations in this gene were shown to increase sensitivity to nicotine’s effects on attention, mood, and vulnerability to tobacco dependence [39]. The Htr3a gene and serotonin have been linked to anxiety and moods as well as psychoses [40], and interactions with 5HT3 receptors are also involved with aggressive behavior, depression, and the stress response [41]. The Gabrd gene is a potential susceptibility gene for childhood-onset mood disorders [42], and variations in this gene are associated with a range of genetic epilepsy syndromes [43] including febrile seizures [44]. Chrna6 and Chrna3 genes are involved in dopaminergic neurotransmission, specifically noted in mesodiencephalic dopamine releasing neurons [45]. At this time, there is no sufficient research linking the Glra2 gene with PTSD, anxiety, or mood disorders.

5. Conclusions

RT-PCR analysis revealed significant changes between groups in several genes implicated in a variety of disorders ranging from PTSD, anxiety, mood disorders, and substance dependence. These data further elucidate the transcriptomic footprint of PTSD in the rodent amygdala and hippocampus as well as transcriptome changes effected by L-theanine and midazolam interventions.

These different neurotransmitter genes are significant biologically and clinically as they are associated with various central nervous system pathologies. The hippocampal genes related to long-term potentiation (Egr1), violence and aggression (Maoa), and major depressive and bipolar disorder (S100a10, serotonin receptors) were shown to be downregulated in PTSD animals and expression levels significantly changed. Other neurotransmitter genes Gabrb2 and Chrm2 were also modulated across treatment groups and warrant further exploration as they are associated with depression related clinical conditions. Amygdalar genes that were changed across treatment groups also were related to psychological conditions such as depression (Drd2), stress enhanced fear learning (Gabra4), chronic mild stress (Drd1a), and anxiety, mood disorders, and psychoses (Htr3a). These findings taken together demonstrate a relationship of possible gene expression changes in important and known neurotransmitter system genes and associated clinical pathologies. This study lays the foundation for future studies to interrogate and investigate the translation of these genes into proteins and the potential for these proteins to act as sites for pharmacologic therapeutic intervention.

This understanding of the effect of PTSD at the level of mRNA transcription contributes to a more complete understanding of the overall metabonomics of PTSD and may allow for more targeted pharmacologic intervention in the future.

The next step to consider in clarifying the molecular mechanism of PTSD may be to characterize PTSD-induced changes to the rat proteome. Changes in mRNA expression are often linked to changes in protein expression, but not in a predictable manner. Small changes in mRNA expression can lead to variable changes in corresponding protein expression. We propose repeating the current experiment using protein assay methods to correlate changes in mRNA expression with changes in protein expression induced by PTSD and the effects of L-theanine and midazolam on these changes. Multiple methods and computational approaches to protein assays are available [46].

We anticipate further investigation into the network of gene relationships between control and PTSD-induced rats using bioinformatic methods such as gene enrichment analysis and network enrichment analysis. Characterization of a comprehensive proteomic network in both control and PTSD-induced rats will allow for comparative proteomic analysis that may lead to the identification of biomarkers for PTSD susceptibility as well as potential targets for pharmacologic interventions [4749]. Additionally, we anticipate investigation into the chemical components of other herbals or pharmaceuticals that may have similar effects on gene expression at the level of mRNA transcription, with a particular emphasis on those with demonstrated efficacy in treating PTSD.

Disclaimer

The views expressed in this paper are those of the authors and do not reflect the official policy or position of the Department of the Army, Department of Defense, or the US Government.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The researchers would like to acknowledge their families and the support they provide, as well as the US Army Graduate Program in Anesthesia Nursing Faculty for their expertise and guidance. In addition, the researchers are very grateful to Jerry Merrill, Ph.D., Ms. Kathleen McKay, VSB Facility Manager at the US Army Institute of Surgical Research, the US Army Institute of Surgical Research, Dr. Dale Glaser, and the Department of Clinical Investigations at Brooke Army Medical Center for the use of their facilities and support. This study was funded by TriService Nursing Research Program.