|Home | About | Journals | Submit | Contact Us | Français|
The period homolog genes Per1, Per2 and Per3 are important components of the circadian clock system. In addition to their role in maintaining circadian rhythm, these genes have been linked to mood disorders, stress response and vulnerability to addiction and alcoholism. In this study, we combined high-resolution sequence analysis and quantitative trait locus (QTL) mapping of gene expression and behavioral traits to identify Per3 as a compelling candidate for the interaction between circadian rhythm, alcohol and stress response. In the BXD family of mouse strains, sequence variants in Per3 have marked effects on steady-state mRNA and protein levels. As a result, the transcript maps as a cis-acting expression QTL (eQTL). We found that an insertion/deletion (indel) variant in a putative stress response element in the promoter region of Per3 causes local control of transcript abundance. This indel results in differences in protein binding affinities between the two alleles through the Nrf2 transcriptional activator. Variation in Per3 is also associated with downstream differences in the expression of genes involved in circadian rhythm, alcohol, stress response and schizophrenia. We found that the Per3 locus is linked to stress/anxiety traits, and that the basal expression of Per3 is also correlated with several anxiety and addiction-related phenotypes. Treatment with alcohol results in increased expression of Per3 in the hippocampus, and this effect interacts with acute restraint stress. Our data provide strong evidence that variation in the Per3 transcript is causally associated with and also responsive to stress and alcohol.
Numerous physiological processes and behaviors in mammals exhibit an inherent circadian oscillation that is vital for normal function. Disruption in the diurnal cycle is a common feature of different types of psychiatric illnesses including stress and anxiety disorders, bipolar disorder and depression.1 The circadian rhythm is under the coordinated regulation of clock genes such as Arntl, Dbp and Csnk1d, and the period homologs Per1, Per2 and Per3.2, 3 There is now growing evidence that implicates allelic variants in these genes in the development of neuropsychiatric diseases.4, 5 Additionally, clock genes have also been linked to responses to neuropharmacological agents in model organisms, as well as alcoholism and addiction in humans.2, 6 In mice, inactivation of Per1 and Per2 leads to changes in sensitivity to cocaine and ethanol,7 and in humans, gene variants in PER2 have been linked to alcoholism.8
There are close ties between stressful life experiences, mood disorders, substance abuse and sleep disturbances. Stressors have often been found to precede depressive episodes in humans.9 Studies in animals also show that stress can lead to depression-like behavior and sleep disruption with correlated changes in neural activity.10, 11 There is also evidence that both chronic stress and exposure to ethanol can have a far reaching impact on expression of circadian genes.12, 13 We have recently shown that prolonged restraint stress significantly upregulates expression of Dbp, Per1, Per2 and Per3 in multiple brain regions of C57BL/6J (B6) and DBA/2J (D2) inbred mouse strains.13 An important question then is whether or not natural variation in the expression of these genes also has a causal impact on stress and ethanol response in these mice?
Studies using inbred mouse strains have already demonstrated the importance of heritable genetic factors in stress-response and stress-induced alterations in sleep pattern.14 Here we investigate a specific period gene, Per3, and evaluate how sequence variants in this gene alter transcription factor binding, gene expression and correlated behavior. In humans, variants in PER3 have been associated with differences in sleep homeostasis,15, 16 susceptibility and onset age of bipolar disorder, and different aspects of mood disorder including differences in novelty seeking.17, 18, 19, 20 In mice, expression of Per3 in the brain is highly variable in strains derived by crossing the B6 and D2 parental strains (the BXD family). We have exploited this normal endogenous variation in transcript abundance to explore the association between a circadian gene and complex stress and alcohol related behaviors.
The BXDs are a genetically highly divergent set of inbred strains that show significant variation in alcohol preference and stress sensitivity.21, 22 They have been extensively phenotyped for behavioral traits and are particularly suitable for systems genetics analysis of sequence variants, expression differences and linkage to behavioral traits. In this study, we identify an insertion/deletion (indel) mutation in a promoter motif located in the 5′ UTR of Per3 that may underlie this variation in transcript abundance. We also examine the functional consequence of this sequence variant on downstream gene expression in the brain and behavioral traits related to stress and alcohol. As the hippocampus is a site vulnerable to both stress and alcohol,23, 24, 25 we also examined the effects of stress and alcohol on Per3 expression in the hippocampus of the BXDs. We provide evidence that expression of Per3 is both causally associated with behavioral differences in stress response and is influenced by stress and ethanol treatment.
The BXD strains were derived by inbreeding the F2 progeny of B6 and D2 strains.26 Newer sets of BXDs were created by inbreeding advanced intercross progeny.27, 28 Animals were weaned at 25 days-of-age and housed in same-sex cages, 2–5 mice per cage at the University of Tennessee Health Science Center. Animals were maintained on a 12:12 light/dark cycle at room temperature that ranged from 20–24°C and had free access to standard laboratory chow and water. All animal work was conducted according to an approved animal use protocol (UTHSC680) and in accordance with procedures approved by the Institutional Animal Care and Use Committee.
We have previously generated large array data sets that survey gene expression in different brain regions of the BXD strains. These data sets profile up to 71 members of the BXD family and have been used extensively for eQTL mapping. In this study, we used the Hippocampus Consortium M430v2 (June06) RMA as the primary mapping data set.29 Additional data sets from the whole brain, prefrontal cortex and cerebellum were used to verify the strong cis-effect for Per3 (see Table 2 of Mozhui et al.30 for an overview of these data sets). A detailed description of strain, sex, tissue preparation and microarray method, normalization and correction procedures for each individual data set is available at www.genenetwork.org.
We used the average strain expression signal detected by each probe set as the quantitative trait. QTL mapping was done using QTL Reaper.31, 32 A set of 3795 informative markers were used for mapping,28, 33 and up to a million permutations were performed to calculate the genome-wide empirical P-values and the associated likelihood ratio statistic score.
The following analytical steps were taken to identify transcripts that may be modulated by the variant in Per3 locus and are potential downstream targets of Per3: (1) We rank ordered all transcripts with expression levels >8 units that had significant point-wise linkage with the marker (rs4140148) closest to Per3. In all, 3880 transcripts were then selected that had trans QTLs near Per3 at a nominal linkage P-value <0.01. (2) Pearson's product correlation was then computed between Per3 and these potential downstream transcripts. Only transcripts that had significant expression correlations with Per3 (P<0.01) were selected for further analysis. (3) We performed partial correlation analysis controlling for the Per3 transcript and removed the transcripts that mapped back to the marker near Per3. (4) Finally, we used prior literature based information on these genes to nominate those transcripts most plausibly linked to Per3 function.
We correlated the first principal component (PC1) of Per3 with 2837 published phenotypes in GeneNetwork, which include numerous behavioral traits related to stress and alcohol response. To evaluate the false discovery rate, the adjusted P-value was computed by Westfall and Young's multiple testing procedure.34 Briefly, we randomly permuted the PC1 data 1000 times. For each permutation, we computed the adjusted P-value of Pearson's correlation between randomized PC1 and 2837 phenotypes. The adjusted P-value was determined by ranking the correlation coefficient.
To evaluate expression differences between the B6 and D2 alleles (B and D, respectively), we measured allelic expression in F1 hybrids using a combination of RT-PCR and a single-base primer extension technique (SNaPshot, Applied Biosystems, www.appliedbiosystems.com). Primer 335 was used to design a pair of PCR primers that target sequences on the same exon and flanking an informative single-nucleotide polymorphism (SNP).
Four pools of RNA from the hippocampus and four pools of genomic DNA from the spleen of F1 hybrids (male and female B6 × D2 F1 and D2 × B6 F1 hybrids) were prepared. To avoid contamination by genomic DNA, RNA pools were treated with Turbo DNase (Ambion, Austin, TX, USA), and then first strand cDNA was synthesized (GE Healthcare, Waukesha, WI, USA). The genomic DNA samples were used as controls, and both cDNA and genomic DNA samples were tested concurrently using the same assay to compare expression from the B and D alleles of Per3. We amplified the cDNA and genomic DNA using GoTaq Flexi DNA polymerase (Promega Corporation, www.promega.com). PCR products were purified using ExoSap-IT (USB Corporation, www.usbweb.com) followed by SNaPshot to extend primers by a single fluorescently labeled ddNTPs. Fluorescently labeled products were purified using calf intestinal phosphatase (CIP, New England BioLabs, Beverly, MA, USA) and separated by capillary electrophoresis on ABI3130 (Applied Biosystems, Foster City, CA, USA). Quantification was done using GeneMapper v4.0 software (Applied Biosystems), and transcript abundance was measured by peak intensities associated with each allele. Ratios of B and D allele expression in both cDNA and genomic DNA pools were computed, and a t-test was done to validate expression differences and polarity of parental alleles.
The transcriptional start site of Per3 was determined using the Database of Transcriptional Start Sites (DBTSS; http://dbtss.hgc.jp).36 Potential transcription factor-binding sites were then identified using the TRANSFAC database and P-Match software37 by screening the 2-kb upstream region of the Per3 sequence. All binding sites were found using default parameters in P-Match.
Two oligonucleotides with 5′-biotin modification representing both the B and D alleles of the Per3 promoter indel were synthesized by MWG Operon (Supplementary Table S1). These oligonucleotides were annealed to form double-stranded probes. Nuclear proteins were isolated from the hippocampus of B6 and D2 mice. EMSA was performed using the Lightshift Chemiluminescent EMSA Kit (Pierce, Rockford, IL, USA). Each 20μl binding reaction contained 50fmol of biotin end-labeled target DNA, 1 × binding buffer, 1μg of Poly(dI.dC) and 25μg of nuclear protein extract. These reactions were incubated for 20min at room temperature. In the supershift experiments, 3μl ATF3 or NRF2 antibody were added, and samples were electrophoresed for 90min at 100 volts in a 6% polyacrylamide gel (0.5% TBE running buffer). They were then transferred to a Biodyne B membrane (Pierce) and crosslinked in a Stratagene CrossLinker (Stratagene, La Jolla, CA, USA). The Protein–DNA complexes were detected by the ThermoShift Luminescent Nucleic Acid Detection Kit (Pierce). The membranes were exposed with X-ray film.
Fresh hippocampus was dissected from the following strains at 2 months-of-age: BXD2, BXD22, BXD28 and BXD42. Extraction of nuclear and soluble protein was performed following a standard protocol (see Supplementary Methods and Li et al.38). Immunoblotting was performed using the WesternDot 625 Goat Anti-Mouse Western Blot Kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. Briefly, 25μg of total protein (nuclear and soluble fractions combined) was loaded into each well of a 4–12% NuPAGE gel (Invitrogen). Anti-Per3 primary antibody was used at a concentration of 0.2μgml−1 and immunopositive bands were visualized using SuperSignal West Femto Maximum Sensitivity Substrate (Pierce) and the HRP-conjugated secondary antibody provided with the WesternDot kit. Blots were scanned using AlphaImager HP, and the intensity of each immunopositive band was quantified using AlphaView Imaging Software (version 18.104.22.168, AlphaInnotech Corp, Santa Clara, CA, USA). Unpaired t-tests were performed to determine if significant differences in protein levels were present between low (BXD2 and BXD22) and high (BXD28 and BXD42) strains.
The genome of DBA/2J strain has recently been sequenced by our group with 100 × average coverage using Illumina GAII and ABI SOLiD technologies.39 Reads with 25 and 50bp from SOLiD and reads with 100bp from Illumina GAII and HiSeq 2000 were mapped to the C57BL/6J reference genome using Corona-Lite v 4.0.2 (http://solidsoftwaretools.com/gf/) and the MAQ program v0.7.1,40 respectively. A threshold of three or more supporting reads and a consensus quality score >30 (Illumina) or a confidence value >0.5 (SOLiD) was used to declare SNPs. A complete compendium of new polymorphisms between the D2 and B6 (the source of the reference mouse genome) parental strains are freely available at GeneNetwork (www.genenetwork.org), and original sequence data are available at www.ncbi.nih.gov/sra.
To predict the potential functional impact of amino-acid substitutions we used two different methods—Sorting Intolerant From Tolerant (SIFT)41 and Polymorphism Phenotyping (PolyPhen).42 SIFT uses sequence homology to predict whether a substitution is likely to affect protein function, and PolyPhen predicts protein function on the basis of 3D structure and multiple alignments of homologous sequences.
The non-coding SNPs (ncSNPs) impact score is the product of the functional element and mutation effect scores. The functional element score reflects how likely it is that the SNP is located in a functional region of the genome. The mutation score reflects the likelihood that the SNP affects the function of an overlapping element based on the evolutionary conservation of the position in alignments with other placental mammals. Thus, all ncSNPs and indels with non-zero impact scores have some evidence of functional relevance.
Forty-three BXD strains were subjected to acute restraint stress (15min) followed by alcohol (1.8gkg−1 i.p.) or saline injection according to the standard operating procedure for Integrative Neuroscience Initiative on Alcoholism (INIA) (https://www.iniastress.org/sop.html). Mice were allocated into five groups: RSE (restraint stress followed by ethanol), RSS (restraint stress followed by saline), NOE (no stress followed by ethanol), NOS (no stress followed by saline) and NON (no stress and no injections, baseline control). For the RSE and RSS groups, mice were placed in an immobilization tube for 15min and then injected with either ethanol or saline (12.5% v/v, 1.8gkg−1), and returned to home cages. For the NOE and NOS groups, mice were injected with ethanol or saline control without restraint treatment. Four hours after treatment, mice were killed and tissues were processed as described below.
Mice from the normative naive set (NON) and restraint/ethanol treatment sets (NOS, NOE, RSE and RSS) were killed at ~2–3 months-of-age by cervical dislocation. Details on the strain, sex and age of each animal can be accessed from http://www.genenetwork.org/dbdoc/UT_ILM_BXD_hipp_NON_0909.html. Brains were removed quickly and stored in RNAlater (Ambion) overnight at 4°C, then kept at −80°C until use. The hippocampus was removed as described in Lu et al.,43 and left and right hippocampi were pooled and stored at −80°C. RNA was isolated from these tissues using RNA STAT-60 (Tel-Test, www.tel-test.com) as per the manufacturer's instructions. RNA purity and concentration were evaluated with a spectrophotometer (NanoDrop Technologies, Rockland, DE, USA), and RNA integrity was assessed on the Agilent 2100 Bioanalyzer (Palo Alto, CA, USA). We required an RNA integrity number of >8. This RNA integrity number value is based on the intensity ratio and amplitude of 18s and 28s rRNA signals.
Hippocampal gene expression was assayed using Illumina Sentrix Mouse-6 v1.1 arrays according to manufacturer's protocols (http://www.illumina.com/). Microarray runs for the five groups (RSE, RSS, NOE, NOS and NON) were carried out in a counterbalanced manner to avoid batch effects. Raw microarray data were normalized using the rank invariant method and background subtraction protocols provided by Illumina as part of the BeadStation software suite. We log-transformed the expression values and stabilized the variance of each array to a mean of 8 and a s.d. of 2 per array. A single unit on this scale (log2) corresponds to a twofold difference in expression.
Restraint and ethanol-treated data sets were used to evaluate the effect of these treatments on Per3. We used a linear nested-ANOVA model to evaluate the effect of Per3 genotype, sex and treatment (RSE, RSS, NOE, NOS and NON) on expression. The ANOVA model was analyzed using the R software (http://www.R-project.org). Differences were considered significant when P<0.05.
Two Per3 KO female mice (C57BL/6J background) were obtained from Vanderbilt University. Creation of the KOs is described in Pendergast et al.44 The original Per3 KOs were created on a 129/Sv background.45 The KOs were then backcrossed to the C57BL/6J strain.
We measured gene expression in the hippocampus of these KOs and wildtype using Illumina Sentrix Mouse-6 v1.1 arrays according to the manufacturer's protocols. Array data were normalized using the rank invariant method. Probes were then filtered by expression and only probe sets with above background expression were retained. Differential expression was then assessed by unpaired Student's t-test. Differences are considered significant when a P-value is <0.05.
The expression of Per3 varies widely among the BXD family with differences that range from approximately twofold to fivefold depending on the specific probe set (Supplementary Table S2). Higher expression is invariably associated with the D allele; an effect seen most prominently in the 3′ UTR (1460662_at and 1421087_at), and the 3′ UTR and last three exons of Per3 (1421086_at) (Figure 1a). None of the probe sets overlap SNPs between D and B genotypes (SNPs will tend to increase the apparent expression of the B genotype; see ref. 46) Additional probe sets targeting introns are either not abundantly expressed (1442243_at) or overlap SNPs (1458176_at and 1441445_at) and were excluded from analysis (Figure 1a).
To generate a single consensus estimate of expression across the panel of BXD strains, we computed the PC1 of three probe sets (1421087_at, 1460662_at and 1421086_at). We used this consensus Per3 eigengene to identify BXD strains with the lowest and highest expression. In the hippocampus, BXD2 and BXD22 have the lowest, and BXD28 and BXD42 have the highest expression of Per3 (Figure 2a). This genetic difference in transcript levels between the extreme strains was confirmed at the protein level by western blot analysis that showed lower PER3 level in BXD2 and BXD22 compared with BXD28 and BXD42 (Figure 2b,c).
All three Per3 probe sets (1421087_at, 1460662_at and 1421086_at) are associated with a significant QTL on distal chromosome 4, close to the location of the Per3 gene itself (150.38Mb). The maximum likelihood ratio statistic is 162 at SNP marker rs6358921 (149.929Mb). The D allele is associated with a positive-additive effect and a 1.66-fold (20.735) increase in expression per allele (Figure 1b). This cis-eQTL for Per3 was confirmed in five other brain regions (Supplementary Figure S1). This highly significant and robust cis-QTL is strong evidence that sequence variants in or near Per3 cause differences in transcript expression levels.
To validate the cis-regulation of Per3, we measured allele specific expression in B6D2F1 and D2B6F1 hybrids.46, 47 Three SNPs in the 3′ UTR and a SNP in exon 15 of Per3 were used to capture expression differences between the D and B alleles. The allele-specific expression analysis confirms that the D allele has approximately twofold higher expression than the B allele (Table 1) and that the effect is almost entirely due to one or more local polymorphisms.
To identify cis-acting polymorphisms in Per3 that cause expression variation between B and D alleles, we examined Per3 sequence differences using the recently sequenced genome of DBA/2J.39 There are a total of 163 SNPs and 30 indels in Per3 between the B6 and D2 genomes. Of these, four are non-synonymous mutations in exons 1 (Pro10Leu), 9 (Arg372Gln), 13 (Phe500Leu) and 15 (Glu590Gly) (Table 2). We used the SIFT and PolyPhen algorithms to predict the impact of these four non-synonymous mutations on protein function. The mutations Pro10Leu in exon 1 and Phe500Leu in exon 9 are predicted to be damaging and deleterious by SIFT and PolyPhen 2, respectively. The mutations Phe500Leu in exon 13 and Glu590Gly in exon 15 are predicted to be deleterious by PolyPhen 2 only. These amino-acid changes can potentially affect the structure and/or function of Per3 and have an influence on the expression of downstream genes.
To determine the potential impact of non-coding SNPs on Per3 transcription, we predicted impact scores for these SNPs based on overlap with conserved functional elements. There are 159 SNPs and 18 indels (Supplementary Table S3) within 5kb upstream or downstream of the transcription start site of Per3. Although the majority of these sequence variants have low-impact scores, two SNPs—rs32260283 and rs32042445—and an indel at 150.418813Mb are all within 1kb of the Per3 transcription start site and are located in conserved regulatory elements (Open Regulatory Annotation Database ID: OREG0048389 and OREG0033366). We searched the TRANSFAC database for DNA motifs in the upstream region of Per3 that may be disrupted by these variants.37 We found that an indel (-/T) disrupts a putative stress response core element (STRE: CCCCT/AGGGG) in BXD strains that inherit the B allele. This site is a potential candidate for the cis-QTL underlying the expression difference and the massive cis-eQTL associated with Per3.
To determine whether the indel in the STRE motif alters transcription factor binding, we performed EMSA using oligonucleotides from B6 (49bp) and D2 (50bp). The EMSA detected two DNA-protein complexes that were unique to the D2 promoter sequence and absent in B6 (Figure 3a). Super-shift assay with the PER3 antibody revealed no difference and no evidence of autoregulation (data not shown). The EMSA result supports the hypothesis that the -/T polymorphism in STRE affects binding of transcription factors.
The STRE motif is known to regulate gene expression in response to cellular stress, possibly by affecting the binding of transcription factors such as NRF2 and ATF3.48, 49 We examined whether NRF2 or ATF3 binds to the STRE site by performing a super-shift assay using antibodies against these two transcription factors. A new slow migrating complex containing NRF2 was present in both B6 and D2 protein extracts but was absent when ATF3 antibody was used (Figure 3b,c). This result suggests that NRF2 mediates the expression of Per3 through the STRE element and is likely to be involved in the differential expression of Per3 in the BXD family.
We next examined whether the marked variation in Per3 expression has an impact on putative downstream phenotypes. We carried out covariance analyses of PC1 of the three probe sets (1421087_at, 1460662_at and 1421086_at) with an extensive database of BXD phenotypes to identify traits that are strongly linked with variation in Per3 transcript abundance. The first component of Per3 expression is significantly correlated with five anxiety-related behaviors and five ethanol-responsive traits (P<0.05; Figure 4; Supplementary Table S4). Each of these 10 behavioral traits was assayed by different research groups using independent BXD cohorts and highlights a robust correlation between Per3 and anxiety- and alcohol-related traits.
To extract and evaluate intermediate molecular and genetic networks interposed between the Per3 polymorphism and higher order phenotypes, we extracted transcripts with trans QTLs that map to the Per3 locus in the neocortex (Figure 5). We found a set of 48 transcripts with likelihood ratio statistic scores above 10 that map near Per3 (nominal P<0.01). This set includes genes associated with RNA processing (Pabpc1, Rbm14, Rps11, Sfrs1 and Sfrs12), ion channel activity (Cacna1c, Trpm2 and Fxyd6) and synaptic function (Anks1b and Doc2a). We next asked if any of these 48 genes that are putatively trans-regulated by Per3 have functional associations with (1) circadian rhythm, (2) stress or (3) alcoholism. For this, we mined the literature for such associations using the Chilibot web resource (http://www.chilibot.net/).50 Four genes (Rbm14, Pabpc1, Irak2 and Il2rb) are associated with circadian rhythm, six genes (Trpm2, Erp29, Sepp1, Cacna1c, Smpd2 and Trim29) are associated with responses to physiological and cellular stress, two genes (Pdlim5 and Fxyd6) are associated with schizophrenia and one gene (Cast) is associated with alcoholism.
To confirm potential Per3 downstream targets, we compared expression levels between wild type and homozygous Per3 KO mice. We found 2217 differentially expressed genes (P<0.05), seven of which (Cast, Trpm2, Anks1b, Gpr81, Hsbp1, Pdlim5 and 4930528J18Rik) are also among the 48 transcripts with trans QTLs near Per3.
Given the associations between Per3 and stress- and alcohol-related traits, we evaluated the expression of Per3 under five different conditions: (1) RSE, (2) RSS, (3) NOS, (4) NOE and (5) NON. Four-way nested ANOVA revealed significant effects of strain and an interaction between stress and ethanol on Per3 expression (P<0.05; Figure 6a). Post hoc analysis showed significantly higher Per3 expression (P<0.05) in the ethanol and the stress with ethanol groups (NOE and RSE) compared with the naive (NON) and saline controls (NOS). However, there was no significant effect of restraint alone—that is, no difference between RSS and NON/NOS groups on Per3 expression (P<0.08) (Figure 6b).
Our work suggests that Per3, a circadian regulator, has a pleiotropic influence on response to alcohol and stress, and provides a potential causal link between stress, sleep disruption and addiction. Expression of Per3 varies up to fivefold among strains of mice and is under strong cis-acting control in the BXD family. In addition to these prominent genotypic effects on expression, there are also effects caused by alcohol and an interaction with restraint-induced stress. Our study is the first to attempt to systematically investigate the association between Per3 and behavioral phenotypes. The implication of our study is that Per3 serves as a new candidate gene for alcohol- and stress-responsive behavior, and supports a role for Per3 in circadian clock function.45
We have identified an indel in the stress response promoter motif that alters DNA–protein complexes formed by Nrf2, a transcriptional activator. This indel is an excellent candidate contributing to strain differences in Per3 expression and is also likely a molecular locus regulating interaction between stress and alcohol responses. The indel in the STRE promoter element is an especially strong candidate for the cis-QTL. EMSA showed two DNA–protein complexes that were unique to the D allele promoter variant associated with higher gene expression. In contrast, low expression of Per3 is linked to the B allele promoter variant. Interestingly, STREs modulate transcription of stress-responsive genes under various cellular challenges, including oxidative stress.49, 51 Because alcohol is known to induce oxidative stress in the brain,52, 53 the STRE variant may be an important regulator of basal and alcohol/stress-responsive expression of Per3.
The transcriptome and allele-specific expression analyses consistently detected higher expression for the D allele of Per3 relative to the B allele. This difference is accompanied by a large number of sequence differences between strains that inherit these alternative B and D alleles and haplotypes. By resequencing the Per3 gene we detected 163 SNPs that differentiated the B and D genotypes. The majority of these SNPs—142 of 163—are in introns. However, we detected four missense mutations, two of which are predicted by both SIFT and PolyPhen 2 to produce significant functional changes in the PER3 protein.
The set of mutations we have detected in both the promoter and untranslated regions are of significant relevance in explaining the high level of variation in Per3 mRNA levels. Complex transcription regulatory loops and posttranscriptional decay mechanisms are known to modulate the normal diurnal dynamics of Per3 expression.54, 55, 56 The sequence of the 3′ UTR is important in decay kinetics of Per3 mRNA, and mutations in a U-rich stretch can radically alter message stability.54 However, none of the 3′ UTR SNPs between B and D alleles overlap this crucial U-rich region. Per3 also has long and short 3′ UTR isoforms. We therefore cannot dismiss differences in 3′ UTR structure and mRNA decay rates as a possible source of variation in Per3 message levels.
The hippocampus is a key region that helps mediate response to stressors. A hallmark is feedback inhibition of the hypothalamic–pituitary–adrenal axis. Chronic stress leads to changes in both neuronal activity and neuronal morphology.25 The hippocampus is also particularly vulnerable to the effects of alcohol. Both human and rodent studies have shown that prolonged exposure to both stress and alcohol leads to hippocampal atrophy.23, 24 The hippocampus is an especially suitable site to study the interaction of stressors and alcohol. In this study, we found a significant interaction of stress and ethanol on Per3 expression, which suggests that ethanol and stress affect the expression of Per3 in an opposing manner. The interaction of ethanol and stress could be mediated via the hypothalamic–pituitary–adrenal system.
We were initially motivated to study Per3 because chronic restraint upregulates the expression of this gene in both B6 and D2 parental strains,13 and because the expression of Per3 is associated with a highly significant cis-QTL in the BXD family. Per1 and Per2 were also upregulated by the same restraint stress protocol, but unlike Per3, expression of these two genes is not self regulated. The very strong cis-regulation of Per3 with a likelihood ratio statistic of ~160 (LOD ~35) in some brain regions made this gene a prime candidate for an analysis of downstream consequences of expression differences using reverse complex trait analysis.45, 56 We found that variation in expression of Per3 has downstream effects on stress response, anxiety-associated traits and addiction-related behavior (Figure 4). The Per3 downstream network also modulates 48 transcripts associated with circadian rhythm, stress response, alcoholism and schizophrenia. Of seven downstream genes that are also differentially expressed between KO and wildtype, three genes (Cast, Trpm2 and Pdlim5) are associated with stress response, alcoholism and schizophrenia. Cast was identified as the susceptibility gene for alcohol dependence through genome-wide association study in humans.57 Trpm2 is highly expressed in the brain and is involved in response to stress and in the genetic etiology of bipolar disorder.58 Pdlim5 has been identified as a potential susceptibility gene for schizophrenia in humans.59 These results not only demonstrate a link between sequence differences, variation in transcript abundance and phenotypes, but also imply that Per3 is an excellent candidate gene for alcohol and stress responses.
In this study, we used an integrative approach to relate gene expression variation of Per3 with systems level phenotypes. We found that a local sequence variant contributes to significant differences in the mRNA and protein level of Per3 among the BXD family as well as the parental B6 and D2 strains. Such quantitative differences in mRNA were correlated with physiological and behavioral traits related to stress and addiction, suggesting that natural variations in Per3 may be causally associated with these traits. Treatments with alcohol and acute restraint in turn had significant effects on Per3 expression. The combined analysis of downstream genes and differential expression between KO and wildtype showed that Per3 has functional implications in alcohol and stress response.
This work was supported by the National Institute on Alcohol Abuse and Alcoholism (NIAAA) and Integrative Neuroscience Initiative on Alcoholism (INIA) (U01AA13499, U01AA017590 and U01AA014425), and supported by National Natural Science Foundation of China (30771200, 30971591 and 30770666). This work was also supported by the UTHSC Center for Integrative and Translational Genomics. We thank the CHDI Foundation and Dr Glenn D. Rosen for providing access to neocortical gene expression data for the BXD strains of mice.60
The authors declare no conflict of interest.
Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp)