Molecular laterality encodes stress susceptibility in the medial prefrontal cortex
Molecular Brain volume 14, Article number: 92 (2021)
Functional lateralization of the prefrontal cortex has been implicated in stress and emotional disorders, yet underlying gene expression changes remains unknown. Here, we report molecular signatures lateralized by chronic social defeats between the two medial prefrontal cortices (mPFCs). Stressed mice show 526 asymmetrically expressed genes between the mPFCs. This cortical asymmetry selectively occurs in stressed mice with depressed social activity, but not in resilient mice with normal behavior. We have isolated highly asymmetric genes including connective tissue growth factor (CTGF), a molecule that modulates wound healing at the periphery. Knockdown of CTGF gene in the right mPFC by shRNA led to a stress-resistant behavioral phenotype. Overexpression of CTGF in the right mPFC using viral transduction induces social avoidance while the left mPFC thereof prevent stress-induced social avoidance. Our study provides a molecular window into the mechanism of stress-induced socioemotional disorders, which can pave the way for new interventions by targeting cortical asymmetry.
Stress induces numerous physiological changes in the brain, but only some individuals develop emotional disorders. Studies have sought to identify key mechanisms and molecules that explain the relationship between stress and emotional disorders. Brain imaging studies have revealed that patients with post-traumatic stress disorder (PTSD) [1, 2] or depression [3,4,5] show asymmetric activity between the two prefrontal cortices (PFC). The subgenual PFC (sgPFC) is a common site of the asymmetric changes observed in patients with bipolar and unipolar depression . Chronic stress has been shown to induce asymmetric changes in the structure [7, 8] and activity  of the medial PFC (mPFC), the rodent counterpart of the sgPFC, supporting the relevance of cortical asymmetry in stress-induced emotional disorders.
According to the functional divergence hypothesis, the two mPFC hemispheres play different roles in stress perception and resilience. In mice  and rats , a lesion in the right mPFC has been shown to cause a stress-resistant phenotype, whereas the effects of such lesions in the left mPFC are not significant . Chronic social defeat stress specifically depresses activity in the left mPFC but not the right mPFC; consistent with this, optogenetic stimulation of activity in the depressed left mPFC, but not the depressed right mPFC, restores stress-induced social avoidance, which is a pro-resilience effect . Despite the data derived from these lesion and gain-of-function studies, however, the physiological mechanisms and related molecular markers associated with stress-induced emotional changes remain unclear.
Several studies have sought to identify molecular effectors that mediate stress responses by exploring stress-induced changes of gene expression in the mPFC [11,12,13]. Comparative analyses have found that certain genes associated with the vascular system, brain injury, stress hormone responses, epigenetic modulation, and other phenomena are differentially expressed between mice that are susceptible and resilient to chronic social defeat stress [14,15,16], suggesting possible mechanisms underlying the resilience to chronic stress. However, whether the mPFC shows hemisphere-specific molecular changes under chronic stress was previously unknown. In the present study, we compared gene expression profiles between the two mPFC hemispheres in stressed mice that were classified into susceptible and resilient groups according to whether they showed social avoidance. We found that stress-induced gene expression profiles are lateralized between the left and right mPFC hemispheres, and that this can contribute to explaining some socio-emotional behaviors.
Methods and materials
C57BL/6 J male mice were used for all experiments. Before stress exposure, mice were housed in groups of 5 or 6 under a 12:12 h light:dark cycle (lights on, 07:00) with food and water available ad libitum. All animal care and experimental procedures were performed in accordance with protocols approved by the Animal Care and Use Committee of Korea Advanced Institute of Science and Technology (KAIST).
Five-week-old mice were anesthetized with Avertin (2,2,2-tribomoethanol; Sigma, USA) and placed in a stereotaxic frame (Neurostar, Germany) for injection of viral vectors. For small hairpin RNA (shRNA)-mediated knockdown of CTGF, the viral vector AAV2/9-GFP-U6-mCTGF-shRNA (Vector Biolabs, USA, 0.5 ul of 1013 GC/ml titer) was delivered to the mPFC (anteroposterior, + 1.95 mm; mediolateral, ± 0.35 mm; dorsoventral, − 1.80 mm) using a Nanofil syringe with a 33-gague injection needle (World Precision Instruments, Inc., USA). As a control, AAV2/5-Scramble shRNA-CMV-mCherry-hGH (Penn Vector Core, USA, 0.5 μl of 2.4 \(\times\) 1012 GC/ml titer) was injected into the same location. Overexpression of CTGF in the mPFC was achieved by delivering AAV2/9-CamKIIα-mCTGF-IRES-mCherry (Vector Biolabs, USA, 0.5 μl of 2.5 \(\times\) 1013 GC/ml titer) to the mPFC; as a control, AAV2/5-CamKIIα-mCherry (Penn Vector Core, USA, 0.5 μl of 4 \(\times\) 1012 GC/ml titer) was injected into the mPFC in the same manner. After surgery, mice were returned to their home cages and housed for 3 weeks to allow recovery and viral expression.
Chronic social defeat stress and social interaction test
Retired CD-1 breeder mice (> 15 weeks old) were screened for aggression as previously described [15, 16]. Each 8-week-old experimental mouse was introduced into the home cage of an unfamiliar, aggressive CD-1 mouse for a 10-min physical defeat exposure. After the attack, the CD-1 and C57BL/6 J mice were separated by a transparent acryl panel containing multiple holes that enabled sensory contact. Control, non-defeated mice were housed in identical cages and rotated in a similar manner. This social defeat stress paradigm was performed in the few hours before the onset of the dark phase (17:00–18:00). After 10 days of social defeat, all mice were singly housed.
The test for social avoidance was performed in a dark room in a white acrylic open-field box (42 × 42 cm) with a removable iron mesh cage (10 × 6.5 × 42 cm) that was used to secure the social target in the middle of one side of the box. Each experimental mouse was placed in the center of the open field and allowed to move freely for 150 s in the absence and presence of a social target (i.e., an unfamiliar CD-1 mouse). All mouse behavior was recorded under infrared lights using a digital video camera, and the time spent in the interaction zone was analyzed using EthoVision XT (Noldus, Netherland). The interaction zone was defined as previously described (14 × 24 cm rectangular area, 8 cm around iron mesh cage) . A social preference index (or interaction ratio) was calculated as the ratio of the time spent in the interaction zone in the presence or absence of the target CD-1 aggressor mouse. Mice with interaction ratios > 1 were defined as resilient, whereas mice with interaction ratios < 1 were deemed susceptible.
Forced swim test (acute stress)
Experimental mice were subjected to a forced swim stress as previously reported . Briefly, mice were placed in a transparent plastic cylinder (diameter 10 cm, height 30 cm) filled with 20 cm of water at a temperature of 24 °C ± 1 °C for 6 min. Mouse behavior was recorded using a video camera, and the onset and total duration of immobility during the last 4 min were analyzed using EthoVision XT (Noldus, Netherland).
Sucrose preference test
The sucrose preference test (SPT) was adapted from the previously reported 8-day sucrose preference protocol , which is a highly reliable method of testing chronic stress-induced anhedonia. Mice subject to social defeat stress were introduced to the sucrose preference paradigm on the day of the social interaction test. The utilized SPT apparatus was designed to have 10 chambers of the same size. Adaptation was performed over day 1 to day 4, with day 1 to day 3 being used for sucrose solution adaptation (1% (wt/vol)) and day 3 and day 4 being used for apparatus adaptation. Two rounds of baseline measurements were collected from day 4 to day 6, after which the mice were deprived of food and water for 1 day. On day 7 to day 8, the preference test was carefully done for 12 h. The sucrose preference calculated by the equation, sucrose intake/total intake (sucrose + regular water) \(\times\) 100%. Anhedonia was defined as a reduction of the sucrose preference compared to the control.
Microarray experiments and data analysis
The raw microarray data have been deposited in the NCBI Gene Expression Omnibus at accession number GSE114224. Mice exposed to chronic social defeat stress were subdivided into two groups according to their sociability index (SI). The average of SI of mice for microarray analysis is 1.76 (control; n = 8), 1.61 (resilient; n = 8), and 0.5 (susceptible; n = 7). Non-defeated control mice were also selected from those who underwent a social interaction test . One day after the social-interaction tests, the mice were sacrificed and their left and right mPFCs were dissected separately on an ice-chilled plate and immediately frozen in liquid nitrogen (LN2). Total RNA was prepared from each of the representative mice using the TRIzol reagent (Invitrogen, USA) and further purified using an RNeasy kit (Qiagen, USA). RNAs were isolated from the left and right mPFC tissues of susceptible, resilient, and control mice and used for triplicate microarray analyses on a MouseRef-8_V2 Expression BeadChip (Illumina, USA).
The prepared microarrays were scanned with an Illumina BeadArray reader and preprocessed with the Illumina GenomeStudio software (version 1.0.6). Raw signals of at least two biological replicates per condition were transformed to the log2 scale and normalized using the Bioconductor lumi package . In cases where the probes for a given gene yielded a p-value (detection p-value) greater than 0.05, the gene was excluded from further analysis. Differentially expressed genes (DEGs) were identified using the Bioconductor limma package , with a false discovery rate (FDR)-adjusted p-value cutoff of 0.05. The expression of significant DEGs was compared in seven pairwise groups: susceptible left vs. susceptible right; resilient left vs. resilient right; control left vs. control right; susceptible left vs. control left; susceptible right vs. control right; resilient right vs. control right; and resilient left vs. control left.
The functions of DEGs with an FDR-adjusted p-value cutoff of 0.05 were predicted by gene ontology (GO) analysis using GeneMANIA . Briefly, GeneMANIA receives a set of input genes as well as other genes related to the input and compares them to a large set of functional association data, including protein and genetic interactions, pathways, co-expression, co-localization, and protein domain similarity.
Real time quantitative reverse transcription-polymerase chain reaction (RT-qPCR)
Our microarray gene expression results were confirmed by RT-qPCR using three biological replicates of the four independent samples applied to the microarray analyses. Commercially available TaqMan Gene Expression Assay primers for the target genes, Cux2 (Mm00500377_m1), Wfs1 (Mm00495979_m1), Mbp (Mm01266402_m1), and Rprm (Mm00469773_s1), Ctgf (Mm01192933_g1), and the reference gene, Gapdh (Mm99999915_g1), were obtained from Applied Biosystems (USA).
Each mouse brain mPFC was harvested, snap-frozen in LN2, and stored at − 70 °C until use. The brain tissue was homogenized with T-Per buffer (tissue protein extraction buffer, #78510; Thermo Scientific, USA) containing 1 × protease inhibitor cocktail (cOmplete, EDTA-free; Roche Diagnostics, USA) and 1 × phosphatase inhibitor cocktail (PhosSTOP; Roche Diagnostics). Lysates were incubated on ice for 1 h with frequent vortexing, after which cleared lysates (supernatants) were collected by centrifugation at 13,000 rpm for 20 min at 4 °C. For detection of CTGF and β-actin, equal amounts of total protein (30 μg) were first incubated at 95 °C for 5 min in sodium dodecyl sulfate (SDS) sample buffer and then resolved by 12% SDS-PAGE (polyacrylamide gel electrophoresis). The resolved proteins were transferred to nitrocellulose (NC) membranes, which were incubated with primary anti-CTGF (sc-14939, sc-365970; Santa Cruz Biotechnology, USA) or anti-β-actin (sc47778; Santa Cruz Biotechnology) antibodies, and then with horseradish peroxidase (HRP)-conjugated secondary antibodies. Immunoreactive proteins were visualized using a SuperSignal West Pico System (Thermo Scientific). The amount of protein expressed was calculated using the ImageJ program (NIH, USA).
After all behavioral experiments were completed, mice were overdosed with Avertin (2,2,2-tribomoethanol; Sigma) and perfused first with heparin solution (66.5 μg/ml) in phosphate-buffered saline (PBS) and then with 4% formaldehyde in PBS. Brains were removed, post-fixed by an overnight incubation in PBS containing 4% formaldehyde, and cut into 40-μm-thick coronal sections using a vibratome (VT1200S; Leica, Germany). The resulting brain slices were permeabilized with PBS containing 0.5% Triton-X for 30 min and then transferred to blocking solution containing 10% normal donkey serum (D9663; Sigma) for 1 h. Brain sections were incubated overnight at 4 °C with goat anti-CTGF primary antibody (1:200, sc-14939; Santa Cruz), and then incubated for 2 h with Cy5-conjugated anti-goat (1:200, 305–175-003; Jackson ImmunoResearch, USA) or fluorescein isothiocyanate (FITC)-conjugated anti-goat (1:200, 305-095-003; Jackson ImmunoResearch, USA) secondary antibodies plus 4’,6-diamidino-2-phenylindole (DAPI; 1:1000). Finally, brain sections were mounted on glass slides and imaged with a confocal laser-scanning microscope (LSM780; Zeiss, Germany).
The numbers of mice used (n) are indicated in the Results and Figure Legends. The graphed values present the mean ± SEM. Comparison between two groups were analyzed using the two-tailed Student’s t-test. We used the two-way RM ANOVA when significant differences between paired two measurements in the same individual analyzed (e.g., no-target/target). When the data did not pass normality, we used the Mann–Whitney Rank Sum Test instead of t-test. Comparisons across more than two groups and between multiple variables were made using two-way ANOVA, followed by post-hoc tests when any of the main effect or interaction was significant at p < 0.05. All statistical analyses were performed using SigmaPlot 12.0 (Systat Software Inc.) with a significance threshold of p < 0.05. (Indicated as n.s. > 0.05, * < 0.05, ** < 0.01) (Additional file 9: Table S7).
Chronic social defeat stress induces asymmetric gene expression in the mPFC
To determine whether stress induces gene expression asymmetry in the mPFC, we performed chronic social defeat experiments using C57BL/6 J mice, as previously described  (Fig. 1a). We first divided stressed mice into two groups—resilient and susceptible—according to their sociability index (SI), which is a measure of the presence or absence of a social preference for a novel encounter (CD-1 mice) (Additional file 1: Figure S1a). The average SI of non-stressed control, resilient, and susceptible mice were 1.5, 1.5, and 0.8, respectively (Additional file 1: Figure S1a), which is consistent with previous descriptions [23,24,25]. We then excluded mice with intermediate SI values between 1.1 ~ 0.9, as susceptible and resilience traits may potentially overlap in this range. For the mice that were finally selected in this way for our microarray analysis, the average SI were 1.76 (control; n = 8), 1.61 (resilient; n = 8), and 0.5 (susceptible; n = 7) (Additional file 1: Figure S1b, c). One day after the sociability tests, we collected the left and right mPFC, including the infralimbic and prelimbic areas, from the three groups of mice. We extracted and pooled equal amounts of RNAs from individual tissues and performed a microarray analysis (Fig. 1a, Additional file 1: Figure S1c).
To assort DEGs according to functions, we calculated a laterality score, defined as the log ratio of gene expression intensity between the left (L) and right (R) mPFC (log2 L/R) of in each of control, resilient and susceptible mice. The laterality of genes between hemispheres of control and stressed mice were analyzed. A correlation analysis using 526 DEGs with meaningful log2 L/R values (FDR-adjusted p-value < 0.05) revealed that susceptible mice showed a higher laterality score than resilient and non-stressed control mice (Fig. 1b, Additional file 4: Table S2). Resilient mice showed a smaller number of lateralized genes (Fig. 1c, 1d, Additional file 5: Table S3) despite changes in a significant number of genes compared with non-stressed controls (Additional file 3: Table S1), suggesting that these genes are changed in a similar direction in both hemispheres.
Highly asymmetric genes in susceptible mice
We then performed a rank-order analysis of the 526 genes according to their Log2 L/R values to identify individual genes that showed extreme laterality in the mPFC (Fig. 1e, Additional file 6: Table S4). Genes showing greater expression in the left than right mPFC were Cux2 (cut-like homeobox 2), Wfs1 (wolframin ER transmembrane glycoprotein), Tnnc1 (troponin C, cardiac/slow skeletal) and Igf2 (insulin-like growth factor 2). Genes with higher expression in the right mPFC were Mbp (myelin basic protein), Trf (transferrin), Rprm (reprimo, TP53-dependent G2 arrest mediator candidate) and Ctgf (connective tissue growth factor). The laterality of these genes was confirmed by RT-qPCR analysis of individual RNA samples obtained from susceptible mice (Additional file 2: Figure S2). The right or left lateralized genes were checked for their transcript expressions compare to control and resilient mice. The normalized expression (by control left) and laterality ratio (L/R) are shown in several highly lateralized genes. The left oriented laterality genes (Cux2, Wfs1) (Additional file 2: Figure S2a, c) and the expression of right oriented genes (Ctgf, Mbp, Rprm) (Additional file 2: Figure S2b, d) are analyzed and its expression is confirmed.
To check the possibility that the DEGs results from difference in the amount of total RNA, we compared the expression levels of the housekeeping genes, Gapdh (glyceraldehyde-3-phosphate dehydrogenase), Actb (β-actin), and B2m (β2 microglobin). We found that these genes showed no difference between the two hemispheres in pairwise group comparisons (p > 0.1 for Actb; p > 0.4 for Gapdh; p > 0.1 for B2m), indicating that DEGs are screened from the equivalent amounts of RNA from each hemisphere (Additional file 7: Table S5). These results suggest that the gene expression laterality between the two mPFC hemispheres involves stress susceptibility.
Functional annotation of lateralized genes in susceptible mice
To better understand the functions of lateralized genes in susceptible mice, we categorized them based on their functional properties in the gene ontology (GO) database (Fig. 2, Additional file 8: Table S6). The GO terms for genes with negative left log2 FC (fold change) values (reflecting lower expression in the left mPFC) included ensheathment of neurons (Cldn11, Fa2h, Gal3st1, Mal, Mbp, Mtmr2, Plp1, Trf, Ugt8a) (Fig. 2a), myelination (Fa2h, Gal3st1, Mal, Mbp, Mtmr2, Plp1, Trf, Ugt8a) (Fig. 2b), action potentials (Cldn11, Drd1a, Gal3st1, Mal, Mbp, Plp1, Scn4b, Tac1, Ugt8a) (Fig. 2c), and gliogenesis (Drd1a, Fa2h, Fgf10, Nfib, Plp1, Sox5, Trf) (Fig. 2d). The GO terms for genes with positive log2 FC values (reflecting increased expression in the left or decreased expression in the right mPFC) included neuronal synaptic plasticity (Bdnf, Camk2a, Egr2, Rasgrf1, Vgf) (Fig. 2e), cognition (Bdnf, Casp1, Igf2, Mef2c, Pde4d, Reln, Serpinf1, Vip) (Fig. 2f), response to metal ions (Anxa11, Fos, Fosb, Junb, Mef2c, Mt3, Tnnc1) (Fig. 2g), and synaptic transmission (Adcyap1, Bdnf, Camk2a, Egr2, Itpka, Mef2c, Rasgrf1, Reln, Rims3, Slc1a3, Sphk1, Vgf) (Fig. 3h). These functional clusters were prominent in susceptible mice and have previously been associated with various neuronal responses to stress, including neural activity depression [26, 27], neuronal degeneration with gliosis , and cognitive alterations . Our results are also consistent with the previous finding that depression of neural activity in the left mPFC contributes to social depression in susceptible mice [5, 9].
Differential expression of CTGF in social defeat stress
CTGF is a key growth factor involved in wound healing at the periphery [30, 31] and in the brain . It is known to mediate interneuron degeneration , demyelination, and gliosis [28, 34, 35], and thereby contributes to effects similar to those induced by chronic stress in the mPFC. We thus, selected CTGF as a functional case study. A recent study of postmortem samples showed that CTGF is increased in the hippocampus and amygdala of patients with major depressive disorder (MDD) . The most highly lateralized gene on our microarray results was Ctgf (Fig. 1). We validated this finding by examining Ctgf expression in susceptible mice, as assessed by laterality ratio (L/R) in susceptible mice compared to those in control and resilient mice (Fig. 3a). We found that the three groups show different asymmetry patterns of Ctgf gene expression, and that the transcript of Ctgf recapitulated the microarray results. We observed a significant difference in RNA expression as laterality score (Left/Right) (Fig. 3a; H = 9.620, df = 2, **p = 0.008, Kruskal–Wallis One Way Analysis of Variance on Ranks; between control and susceptible, *p < 0.05, between control and resilient, p > 0.05. between resilient and susceptible, *p < 0.05, Turkey test post-hoc). We checked whether the functional protein level was consistent with the transcript profile. Importantly, we found that the CTGF protein level increased in the right hemisphere of susceptible mice and both hemispheres of resilient mice under stress (Fig. 3b). Our analysis revealed that CTGF proteins showed a significant right-oriented expression in susceptible mice (Fig. 3b, upper panel) and a significant difference in expression between the left and right hemispheres (Fig. 3b, lower panel: For group F(2, 38) = 14.445, **p < 0.001, laterality F(1, 38) = 3.634, p = 0.064, group \(\times\) laterality interaction F(2, 38) = 4.489, *p = 0.018, laterality within susceptible **p = 0.002, two-way ANOVA with Holm-sidak post-hoc analysis). In the resilient mice, we did not detect any significant between-hemisphere difference in CTGF expression (Fig. 3b, lower panel; laterality within resilient p = 0.357; two-way ANOVA). The CTGF expression in resilient and susceptible mice showed a significant difference compare to control (within left mPFC, control and resilient, **p < 0.001, control and susceptible, p = 0.120, resilient and susceptible, *p = 0.016; within right mPFC, control and resilient, *p = 0.025, control and susceptible, **p < 0.001, resilient and susceptible, p = 0.166, two-way ANOVA with Holm-sidak post-hoc analysis). These findings indicate that the RNA and protein expression levels of CTGF are lateralized in susceptible mice but balanced in resilient and control mice, suggesting that the laterality score could be a relevant biomarker for stress susceptibility.
CTGF is functionally lateralized in the mPFC hemispheres
To induce hemisphere-specific knockdown of CTGF, we injected a viral vector expressing an shRNA targeting CTGF (AAV2/9-GFP-U6-mCTGF-shRNA) into the left or right mPFC. As a control for the effects of viral transfection, we injected a non-targeting shRNA-expression vector (AAV2/5-Scramble shRNA-CMV-mCherry-hGH) into the contralateral mPFC (Fig. 3c). After 3 weeks, we subjected the same mice to 10 days of chronic social defeat stress followed by social interaction tests (Fig. 3c, lower panel). Using immunostaining and Western blotting, we found that the CTGF shRNA effectively reduced the CTGF protein level in the mPFC (Fig. 3c, d; *p = 0.0464, t4 = 2.850, t-test). After the 10-day social defeat protocol, control mice with bilateral mPFC injections of scrambled shRNA showed a significant increase in social avoidance of CD-1 mice (Fig. 3e; For stress, F(1, 28) = 1.731, p = 0.199, for target F(1, 28) = 7.721, *p = 0.010, stress \(\times\) target interaction, F(1, 28) = 5.126, *p = 0.032, two-way RM ANOVA; within non-stressed, effect of target, **p = 0.002; within stressed, effect of target, p = 0.709; within non-target, effect of stress, p = 0.739, within target, effect of stress, *p = 0.019, Post-hoc analysis with Holm-Sidak method). Mice with left mPFC CTGF depletion showed levels of social avoidance similar to those of control mice (Fig. 3f; For stress, F(1, 28) = 2.514, p = 0.124, for target F(1, 28) = 20.319, **p < 0.001, stress \(\times\) target interaction, F(1, 28) = 4.478, *p = 0.043, two-way RM ANOVA; within non-stressed, effect of target, **p < 0.001; within stressed, effect of target, p = 0.080; within non-target, effect of stress, p = 0.949, within target, effect of stress, *p = 0.014, Post-hoc analysis with Holm-Sidak method). Mice with right mPFC CTGF depletion showed no social avoidance (Fig. 3g; For stress, F(1, 24) = 1.968, p = 0.173, for main target F(1, 24) = 58.986, **p < 0.001, stress \(\times\) target interaction, F(1, 24) = 2.480, p = 0.128, two-way RM ANOVA; within non-stressed, effect of target, **p < 0.001; within stressed, effect of target, **p < 0.001; Post-hoc analysis with Holm-Sidak method). Analysis of representative heatmaps showed that mice with right mPFC CTGF knockdown exhibited reduced social avoidance when the CD-1 target was present (Fig. 3h). These results are highly consistent with a previous report that a lesion in the right mPFC leads to stress resistance [9, 10] and suggest that the upregulation of CTGF in the right mPFC plays an important role in stress perception of all stressed mice.
CTGF overexpression in the right mPFC facilitates social avoidance
Our analysis of CTGF gene expression revealed that CTGF protein increased in the right mPFC of all socially stressed mice, while the corresponding changes in the left mPFC differed by the stress susceptibility of the mice (Fig. 3b). To test the role of increased CTGF in the right mPFC, we expressed AAV2/9-CamKIIα-mCTGF in the right mPFC (Fig. 4a). To generate control mice, we injected the virus, AAV2/5-CamKIIα-mCherry, into the ipsilateral right mPFC. In this experiment, the left mPFC was left intact to check the right hemisphere-specific role of CTGF when social stress is experienced. The expression of CTGF was confirmed by immunohistochemistry (Fig. 4b).
After 3 weeks of expression and 10 consecutive days of social defeat stress, we measured the sociability value by measuring each test mouse’s interaction time with a CD-1 mouse within the interaction zone (Fig. 4c). Heatmaps were generated for the sociability of control mice and those with CTGF expressed in the right-side mPFC (Right-mPFCCTGF; R-CTGF) in the absence and presence of a target (for representative maps, see Fig. 4a). Reflecting their normal perception of social stress, control mice showed increased social avoidance when the target CD-1 mouse was present (Fig. 4c left panel; For target, F(1, 17) = 0.00147, p = 0.970, stress, F(1, 17) = 4.008, p = 0.061, target \(\times\) stress interaction, F(1, 17) = 17.763, **p < 0.001, two-way RM ANOVA; within non-stress, effect of target, *p = 0.012, within stress the effect of target, **p = 0.005, within target presence, effect of stress, **p < 0.001; Post-hoc analysis with Holm-Sidak method). However, right-mPFCCTGF mice showed no more social avoidance when socially defeated (Fig. 4c right panel; main target effect, F (1, 19) = 18.637, **p < 0.001, stress, F(1, 19) = 3.371, p = 0.082, target \(\times\) stress interaction, F(1, 19) = 3.965, p = 0.061, two-way RM ANOVA; within non-stress, effect of target p = 0.125, within stress, effect of target, **p < 0.001; Post-hoc analysis with Holm-Sidak method).
Furthermore, right-mPFCCTGF mice showed social avoidance under the non-stressed condition, to the point that there was a significant difference in sociability between control and right-mPFCCTGF mice that had not been subjected to stress (Fig. 4d, left panel; For group, F(1, 16) = 1.266, p = 0.277, for target, F(1, 16) = 0.143, p = 0.710, group \(\times\) target interaction, F(1, 16) = 8.921, **p = 0.009, two-way RM ANOVA; within control, effect of target, *p = 0.038, within R-CTGF, effect of target, p = 0.068, within target presence, the effect of R-CTGF, *p = 0.025; Post-hoc analysis with Holm-Sidak method). This increase of social avoidance under the non-stressed condition indicates that CTGF expression in the right mPFC led mice to exhibit a socioemotionally depressive state. Our analysis of sociability index (SI) values clearly showed that the right-mPFCCTGF mice were socially avoidant when non-stressed probably due to increased stress sensitivity (Fig. 4d, right panel; t16 = 3.020, **p = 0.00814, t-test).
To confirm that the depressive phenotype of right-mPFCCTGF mice is relevant to other depression measures, we used the sucrose preference and forced swim tests as measures of anhedonia and despair, respectively (Fig. 4e). As expected from the results of our social interaction test (Fig. 4d), we observed a significant difference in anhedonia between control and right-mPFCCTGF mice (Fig. 4e left panel; *p = 0.033, T = 54, Mann–Whitney Rank Sum Test). Similarly, the forced swim test showed that right-mPFCCTGF mice exhibited more despair during a forced swim than control mice, as represented by their immobility time (Fig. 4e right panel, *p = 0.0256, t12 = -2.547, t-test). The representative video shows that right-mPFCCTGF mice exhibited dramatically more immobile behavior compared to control mice under the stress-naïve condition (Additional file 10: Video 1). These results suggest that CTGF mediates the stress-induced emotional changes in the right mPFC (Fig. 6).
CTGF overexpression in the left mPFC facilitates social resilience
After social defeat, resilient mice showed no laterality in CTGF expression, but instead exhibited balanced and increased expression in both mPFC hemispheres (Fig. 3a, b). Considering that the knockout of CTGF in the left mPFC shows stress avoidance to social defeats (Fig. 3f), the increased CTGF in the left mPFC of resilient mice may contribute to enhancing resilience. To test the hypothesis that increased CTGF in the left mPFC made mice socially resilient under social stress, we injected a CTGF-overexpressing viral vector (AAV2/9-CamKIIα-mCTGF) into the left mPFC (Left-mPFCCTGF; L-CTGF, Fig. 5). To control for stress-induced social response- and lesion-related effects, we injected a control viral vector (AAV2/5-CamKIIα-mCherry) to the right mPFC (Fig. 5a). As control mice, we injected a control viral vector (AAV2/5-CamKIIα-mCherry) bilaterally into the left and right mPFC. Western blotting confirmed that the CTGF protein level was increased in the left mPFC of L-CTGF mice compared with control mice (Fig. 5b; **p < 0.01, t6 = − 6.142, t-test).
After a 10-day social defeat protocol, we measured the sociability value by monitoring the time spent within the interaction zone of a CD-1 target mouse (Fig. 5c). We assessed the sociability of control mice and left-mPFCCTGF mice in the absence (left) and presence (right) of a target mouse (for representative heatmaps, see Fig. 5a, upper panel; non-stressed, lower panel; stressed). Whereas control mice showed increased social avoidance of CD-1 mice (Fig. 5c, left panel; For stress, F(1, 32) = 2.517, p = 0.122, target, F(1, 32) = 29.886, **p < 0.001, stress \(\times\) target interaction, F(1, 32) = 9.537, **p = 0.004, two-way RM ANOVA; within non-stressed, the effect of target, **p < 0.001; within stressed, the effect of target, p = 0.063; within target presence, effect of stress, **p = 0.004; Post-hoc analysis with Holm-Sidak method), Left-mPFCCTGF mice showed no social avoidance (Fig. 5c, right panel; For stress, F(1, 27) = 0.0540, p = 0.818, main target effect, F(1, 27) = 20.560, **p < 0.001, stress \(\times\) target interaction, F(1, 27) = 0.266, two-way RM ANOVA, within L-CTGF non-stressed, the effect of target,**p < 0.001; within L-CTGF stressed, the effect of target, *p = 0.021; Post-hoc analysis with Holm-Sidak method).
To evaluate the stress coping role of CTGF in acute stress, we measured the latency to immobility after forced swim stress as an acute stress paradigm. Under these conditions, we failed to find any difference in the latency (Fig. 5d left panel, p = 0.110, t12 = 1.728, t-test) or duration (Fig. 5d right panel, p = 0.758, t12 = 0.315, t-test) of immobility in Left-mPFCCTGF mice. These results suggest that CTGF expression is specifically involved in the adaptation to chronic stress, which typically leads to emotional and behavioral changes.
These results show the increase of CTGF expression in the left mPFC can critically prevent the development of stress-induced social avoidance (Fig. 5). Without the increase of CTGF in the left mPFC, stressed mice may show depression-like phenotypes as shown in CTGF knockout experiments (Fig. 3).
Hemispheric lateralization has been implicated in emotional disorders. Transcranial magnetic stimulation (TMS) of the left PFC, which is known to reduce anxiety by retrieving positive memory , has been used empirically to treat depression [38,39,40]. However, the mechanism of PFC lateralization has been poorly studied because of the lack of robust animal models. Our study reveals that chronic social defeat stress induces hemisphere-specific gene expression in the mPFC, opening a new avenue for studying the molecular mechanisms of hemispheric lateralization and stress-induced mental illness (Fig. 6).
Screening for stress adaptation-related genes by hemisphere-specific analysis
Consistent with previous reports that patients with emotional disorders show lateralization of neuronal activity in the PFC [9, 10], we herein report that the two mPFC hemispheres undergo differential molecular changes under chronic stress conditions (Fig. 1). We identified key molecules that modulate various neuronal functions of stress adaptation, including neuronal degeneration, excitability, synaptic transmission, and cognition (Fig. 2), all of which have been associated with changes induced by chronic stress. In addition, we found that mice exhibiting stress-induced behavioral changes show greater molecular laterality than resilient or control mice, suggesting that laterality score may represent a relevant biomarker for screening stress-modulating genes (Fig. 1).
We identified three groups of stress-related genes. The genes of the first group were altered by stress but did not show any evidence of lateralization. These include Hba-a1 (hemoglobin alpha, adult chain 1) and Hbb-b1 (hemoglobin, beta adult major chain), which were increased in both hemispheres of mice exposed to chronic social defeat stress, but not in those exposed to acute stress . Also included in this group is Arc (activity-regulated cytoskeleton-associated protein), an immediate early gene that is induced by neural activity and was herein found to be reduced in both hemispheres of stressed mice (Additional file 3: Table S1). This finding is consistent with changes that have been found in postmortem brain tissue of depressed patients and mice exposed to social defeat stress . The second group of genes showed a strong left dominance, as measured by log2 L/R values, in response to chronic stress (Fig. 1e, Additional file 6: Table S4). This group included Adcyap1 (adenylate cyclase activating polypeptide 1, also known as PACAP), which is associated with post-traumatic stress disorder (PTSD) . Notably, suppression of PACAP expression has been shown to reduce corticosterone secretion and depression-like behaviors . The third group of genes showed right dominance, characterized by negative log2 L/R values. The genes in this group include Adra2a (alpha-2A adrenergic receptor), which was reduced in the left mPFC (Additional file 6: Table S4). Consistent with this, loss of ADRA2A is associated with socially withdrawn behaviors , depression and autism [45, 46]; moreover, an increase in ADRA2A has been shown to prevent the withdrawal phenotype and increase novelty-seeking behavior, as observed in ADHD . Increased right mPFC expression of Ctgf (Fig. 1e, Additional file 6: Table S4), another group 3, right-dominant gene, was recently observed in the amygdala of an MDD patient . In addition, CTGF levels are reportedly increased in the hippocampus of rats exposed to chronic social defeat stress  and in the hippocampus and amygdala of rats exposed to predator-scent stress .
Studies have shown that there appear to be gender differences in the functional laterality of the brain [50,51,52]. While the right vmPFC is involved in socioemotional decision making in males, the left vmPFC is involved in females . Furthermore, females are known to be much more vulnerable to depression and PTSD  and show a higher sensitivity in their CRF response  than males. Thus, it is plausible that the molecular laterality revealed in our study may provide new insights into gender differences in stress perception and adaptation. Female mice may show more laterality between the two mPFC hemispheres compared to males. As a new social defeat paradigm was recently developed for females [55,56,57], we could test this hypothesis in a future study.
There have been numerous reports regarding the lateralization of hippocampus under stress conditions [58,59,60,61]. Here in this report, we focused on the medial prefrontal cortex, however, there are possibilities of molecular lateralization in other subcortical area such as the hippocampus. With network analysis of regions connected to the mPFC, we may obtain more information about the mechanism of functional laterality induced by emotional stress. Further efforts to study brain regions that are systemically connected to the mPFC will be a future target to consolidate the mechanism of laterality induced by emotional stress.
CTGF modulates the hemisphere-specific roles of mPFC
In line with the previous finding that right mPFC activity induces stress recognition and left mPFC activity facilitates resilience to chronic stress , CTGF was increased in the right mPFC of all stressed mice and in the left mPFC of resilient mice (Fig. 3b). This all-or-none stress information-coding role in the right hemisphere was demonstrated by knockdown and overexpression of CTGF in the mPFC (Fig. 3, 4, 5). Our results showed that right KD mice exhibited increased sociability (Fig. 3g) and right-mPFCCTGF mice exhibited increased stress vulnerability (Fig. 4). The role of the left hemisphere in detecting stress was supported by the decreased sociability of left KD mice (Fig. 3f) and the increased resilience of left-mPFCCTGF mice (Fig. 5c).
The stress sensitivity increment under non-stressed conditions is supported by the increased social avoidance, anhedonia and despair responses in right-mPFCCTGF mice (R-CTGF). (Fig. 4d, e) We can further test the stress sensitivity of R-CTGF mice by measuring corticosterone because this stress hormone is a reliable indicator of stress awareness [62, 63]. Actually, in a pilot study with R-CTGF mice, we found a meaningful stress hormonal increment in the R-CTGF mice under non-stressed conditions indicating that R-CTGF mice were more vulnerable to stress (unpublished data). Further efforts to solve the mechanism of social avoidance of R-CTGF mice might provide answers to this question. We conclude that right-mPFCCTGF mice became depressive even during stress-naïve conditions and are more vulnerable to social stress.
CTGF is known to activate tumor growth factor (TGF)-β signaling through Smad2/3, which in turn, stimulates CTGF expression [64, 65], creating a positive-feedback loop that may contribute to long-term adaptation to chronic stress. Additional studies will be necessary to elucidate the details of this mechanism. Thus, the data from studies on the role of the PFC in emotional control should be interpreted with caution, and any plans for their application must consider the issues of cortical asymmetry and functional divergence. For example, blocking CTGF may alleviate stress perception through a right mPFC-dependent mechanism, but it could harmfully attenuate the pro-resilience effects mediated by the left mPFC. In sum, we herein provide evidence supporting the idea that molecular laterality is a critical mechanism of stress adaptation that encodes information key to the development of stress-related emotional disorders, providing a framework for the potential development of new strategies for treating stress-associated mental illnesses.
Availability of data and materials
The raw microarray data are available in NCBI Gene Expression Omnibus at accession number GSE114224.
Meyer T, Quaedflieg CW, Giesbrecht T, Meijer EH, Abiad S, Smeets T. Frontal EEG asymmetry as predictor of physiological responses to aversive memories. Psychophysiology. 2014;51:853–65.
Meyer T, Smeets T, Giesbrecht T, Quaedflieg CW, Smulders FT, Meijer EH, Merckelbach HL. The role of frontal EEG asymmetry in post-traumatic stress disorder. Biol Psychol. 2015;108:62–77.
Henriques JB, Davidson RJ. Left frontal hypoactivation in depression. J Abnorm Psychol. 1991;100:535–45.
Martinot JL, Hardy P, Feline A, Huret JD, Mazoyer B, Attar-Levy D, Pappata S, Syrota A. Left prefrontal glucose hypometabolism in the depressed state: a confirmation. Am J Psychiatry. 1990;147:1313–7.
Bruder GE, Stewart JW, McGrath PJ. Right brain, left brain in depressive disorders: clinical and theoretical implications of behavioral, electrophysiological and neuroimaging findings. Neurosci Biobehav Rev. 2017;78:178–91.
Drevets WC, Price JL, Simpson JR Jr, Todd RD, Reich T, Vannier M, Raichle ME. Subgenual prefrontal cortex abnormalities in mood disorders. Nature. 1997;386:824–7.
Perez-Cruz C, Muller-Keuker JI, Heilbronner U, Fuchs E, Flugge G. Morphology of pyramidal neurons in the rat prefrontal cortex: lateralized dendritic remodeling by chronic stress. Neural Plast. 2007;2007:46276.
Czeh B, Muller-Keuker JI, Rygula R, Abumaria N, Hiemke C, Domenici E, Fuchs E. Chronic social stress inhibits cell proliferation in the adult medial prefrontal cortex: hemispheric asymmetry and reversal by fluoxetine treatment. Neuropsychopharmacology. 2007;32:1490–503.
Lee E, Hong J, Park YG, Chae S, Kim Y, Kim D. Left brain cortical activity modulates stress effects on social behavior. Sci Rep. 2015;5:13342.
Sullivan RM, Gratton A. Lateralized effects of medial prefrontal cortex lesions on neuroendocrine and autonomic stress responses in rats. J Neurosci. 1999;19:2834–40.
Bondar N, Bryzgalov L, Ershov N, Gusev F, Reshetnikov V, Avgustinovich D, Tenditnik M, Rogaev E, Merkulova T. Molecular adaptations to social defeat stress and induced depression in mice. Mol Neurobiol. 2018;55:3394–407.
Stankiewicz AM, Goscik J, Swiergiel AH, Majewska A, Wieczorek M, Juszczak GR, Lisowski P. Social stress increases expression of hemoglobin genes in mouse prefrontal cortex. BMC Neurosci. 2014;15:130.
Farhang S, Barar J, Fakhari A, Mesgariabbasi M, Khani S, Omidi Y, Farnam A. Asymmetrical expression of BDNF and NTRK3 genes in frontoparietal cortex of stress-resilient rats in an animal model of depression. Synapse. 2014;68:387–93.
Bagot RC, Cates HM, Purushothaman I, Vialou V, Heller EA, Yieh L, LaBonte B, Pena CJ, Shen L, Wittenberg GM, Nestler EJ. Ketamine and imipramine reverse transcriptional signatures of susceptibility and induce resilience-specific gene expression profiles. Biol Psychiatry. 2017;81:285–95.
Berton O, McClung CA, Dileone RJ, Krishnan V, Renthal W, Russo SJ, Graham D, Tsankova NM, Bolanos CA, Rios M, et al. Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress. Science. 2006;311:864–8.
Krishnan V, Han MH, Graham DL, Berton O, Renthal W, Russo SJ, Laplant Q, Graham A, Lutter M, Lagace DC, et al. Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions. Cell. 2007;131:391–404.
Golden SA, Covington HE 3rd, Berton O, Russo SJ. A standardized protocol for repeated social defeat stress in mice. Nat Protoc. 2011;6:1183–91.
Can A, Dao DT, Arad M, Terrillion CE, Piantadosi SC, Gould TD. The mouse forced swim test. J Vis Exp JoVE. 2012. https://doi.org/10.3791/3638.
Liu MY, Yin CY, Zhu LJ, Zhu XH, Xu C, Luo CX, Chen H, Zhu DY, Zhou QG. Sucrose preference test for measurement of stress-induced anhedonia in mice. Nat Protoc. 2018;13:1686–98.
Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24:1547–8.
Wettenhall JM, Smyth GK. limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004;20:3705–6.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214-220.
Anacker C, Scholz J, O’Donnell KJ, Allemang-Grand R, Diorio J, Bagot RC, Nestler EJ, Hen R, Lerch JP, Meaney MJ. Neuroanatomic differences associated with stress susceptibility and resilience. Biol Psychiatry. 2016;79:840–9.
Jiang C, Lin WJ, Sadahiro M, Labonte B, Menard C, Pfau ML, Tamminga CA, Turecki G, Nestler EJ, Russo SJ, Salton SR. VGF function in depression and antidepressant efficacy. Mol Psychiatry. 2018;23:1632–42.
Wook Koo J, Labonte B, Engmann O, Calipari ES, Juarez B, Lorsch Z, Walsh JJ, Friedman AK, Yorgason JT, Han MH, Nestler EJ. Essential role of mesolimbic brain-derived neurotrophic factor in chronic social stress-induced depressive behaviors. Biol Psychiatry. 2016;80:469–78.
Wang M, Perova Z, Arenkiel BR, Li B. Synaptic modifications in the medial prefrontal cortex in susceptibility and resilience to stress. J Neurosci. 2014;34:7485–92.
Sinha R, Lacadie CM, Constable RT, Seo D. Dynamic neural activity during stress signals resilient coping. Proc Natl Acad Sci U S A. 2016;113:8837–42.
Lehmann ML, Weigel TK, Elkahloun AG, Herkenham M. Chronic social defeat reduces myelination in the mouse medial prefrontal cortex. Sci Rep. 2017;7:46548.
Uribe-Marino A, Gassen NC, Wiesbeck MF, Balsevich G, Santarelli S, Solfrank B, Dournes C, Fries GR, Masana M, Labermeier C, et al. Prefrontal cortex corticotropin-releasing factor receptor 1 conveys acute stress-induced executive dysfunction. Biol Psychiatry. 2016;80:743–53.
Igarashi A, Okochi H, Bradham DM, Grotendorst GR. Regulation of connective tissue growth factor gene expression in human skin fibroblasts and during wound repair. Mol Biol Cell. 1993;4:637–45.
Oemar BS, Luscher TF. Connective tissue growth factor. Friend or foe? Arterioscler Thromb Vasc Biol. 1997;17:1483–9.
Liu Y, Liu Z, Li X, Luo B, Xiong J, Gan W, Jiang M, Zhang Z, Schluesener HJ, Zhang Z. Accumulation of connective tissue growth factor+ cells during the early phase of rat traumatic brain injury. Diagn Pathol. 2014;9:141.
Khodosevich K, Lazarini F, von Engelhardt J, Kaneko H, Lledo PM, Monyer H. Connective tissue growth factor regulates interneuron survival and information processing in the olfactory bulb. Neuron. 2013;79:1136–51.
Ercan E, Han JM, Di Nardo A, Winden K, Han MJ, Hoyo L, Saffari A, Leask A, Geschwind DH, Sahin M. Neuronal CTGF/CCN2 negatively regulates myelination in a mouse model of tuberous sclerosis complex. J Exp Med. 2017;214:681–97.
Mendes FA, Coelho Aguiar JM, Kahn SA, Reis AH, Dubois LG, Romao LF, Ferreira LS, Chneiweiss H, Moura Neto V, Abreu JG. Connective-tissue growth factor (CTGF/CCN2) induces astrogenesis and fibronectin expression of embryonic neural cells in vitro. PLoS ONE. 2015;10:e0133689.
Turner CA, Sharma V, Hagenauer MH, Chaudhury S, O’Connor AM, Hebda-Bauer EK, Thompson RC, Myers RM, Bunney WE, Barchas JD, et al. Connective tissue growth factor is a novel prodepressant. Biol Psychiatry. 2018. https://doi.org/10.1016/j.biopsych.2018.04.013.
Balconi M, Ferrari C. Left DLPFC rTMS stimulation reduced the anxiety bias effect or how to restore the positive memory processing in high-anxiety subjects. Psychiatry Res. 2013;209:554–9.
Gershon AA, Dannon PN, Grunhaus L. Transcranial magnetic stimulation in the treatment of depression. Am J Psychiatry. 2003;160:835–45.
Pascual-Leone A, Rubio B, Pallardo F, Catala MD. Rapid-rate transcranial magnetic stimulation of left dorsolateral prefrontal cortex in drug-resistant depression. Lancet. 1996;348:233–7.
Klein E, Kreinin I, Chistyakov A, Koren D, Mecz L, Marmur S, Ben-Shachar D, Feinsod M. Therapeutic efficacy of right prefrontal slow repetitive transcranial magnetic stimulation in major depression: a double-blind controlled study. Arch Gen Psychiatry. 1999;56:315–20.
Covington HE 3rd, Lobo MK, Maze I, Vialou V, Hyman JM, Zaman S, LaPlant Q, Mouzon E, Ghose S, Tamminga CA, et al. Antidepressant effect of optogenetic stimulation of the medial prefrontal cortex. J Neurosci. 2010;30:16082–90.
Ressler KJ, Mercer KB, Bradley B, Jovanovic T, Mahan A, Kerley K, Norrholm SD, Kilaru V, Smith AK, Myers AJ, et al. Post-traumatic stress disorder is associated with PACAP and the PAC1 receptor. Nature. 2011;470:492–7.
Lehmann ML, Mustafa T, Eiden AM, Herkenham M, Eiden LE. PACAP-deficient mice show attenuated corticosterone secretion and fail to develop depressive behavior during chronic social defeat stress. Psychoneuroendocrinology. 2013;38:702–15.
Rubin DH, Althoff RR, Ehli EA, Davies GE, Rettew DC, Crehan ET, Walkup JT, Hudziak JJ. Candidate gene associations with withdrawn behavior. J Child Psychol Psychiatry. 2013;54:1337–45.
Goodwin RD, Fergusson DM, Horwood LJ. Early anxious/withdrawn behaviours predict later internalising disorders. J Child Psychol Psychiatry. 2004;45:874–83.
Ooi YP, Rescorla L, Ang RP, Woo B, Fung DS. Identification of autism spectrum disorders using the child behavior checklist in Singapore. J Autism Dev Disord. 2011;41:1147–56.
de Cerqueira CC, Polina ER, Contini V, Marques FZ, Grevet EH, Salgado CA, da Silva PO, Picon FA, Belmonte-de-Abreu P, Bau CH. ADRA2A polymorphisms and ADHD in adults: possible mediating effect of personality. Psychiatry Res. 2011;186:345–50.
Stankiewicz AM, Goscik J, Majewska A, Swiergiel AH, Juszczak GR. The effect of acute and chronic social stress on the hippocampal transcriptome in mice. PLoS ONE. 2015;10:e0142195.
Daskalakis NP, Cohen H, Cai G, Buxbaum JD, Yehuda R. Expression profiling associates blood and brain glucocorticoid receptor signaling with trauma-related individual differences in both sexes. Proc Natl Acad Sci. 2014;111:13529–34.
Martin-Soelch C, Szczepanik J, Nugent A, Barhaghi K, Rallis D, Herscovitch P, Carson RE, Drevets WC. Lateralization and gender differences in the dopaminergic response to unpredictable reward in the human ventral striatum. Eur J Neurosci. 2011;33:1706–15.
Sutterer MJ, Koscik TR, Tranel D. Sex-related functional asymmetry of the ventromedial prefrontal cortex in regard to decision-making under risk and ambiguity. Neuropsychologia. 2015;75:265–73.
Reber J, Tranel D. Sex differences in the functional lateralization of emotion and decision making in the human brain. J Neurosci Res. 2017;95:270–8.
Breslau N. Gender differences in trauma and posttraumatic stress disorder. J Gend Specif Med. 2002;5:34–40.
Bangasser DA, Wiersielis KR. Sex differences in stress responses: a critical role for corticotropin-releasing factor. Hormones (Athens). 2018;17:5–13.
Harris AZ, Atsak P, Bretton ZH, Holt ES, Alam R, Morton MP, Abbas AI, Leonardo ED, Bolkan SS, Hen R, Gordon JA. A novel method for chronic social defeat stress in female mice. Neuropsychopharmacology. 2018;43:1276–83.
Iniguez SD, Flores-Ramirez FJ, Riggs LM, Alipio JB, Garcia-Carachure I, Hernandez MA, Sanchez DO, Lobo MK, Serrano PA, Braren SH, Castillo SA. Vicarious social defeat stress induces depression-related outcomes in female mice. Biol Psychiatry. 2018;83:9–17.
Takahashi A, Chung JR, Zhang S, Zhang H, Grossman Y, Aleyasin H, Flanigan ME, Pfau ML, Menard C, Dumitriu D, et al. Establishment of a repeated social defeat stress model in female mice. Sci Rep. 2017;7:12838.
Woon FL, Sood S, Hedges DW. Hippocampal volume deficits associated with exposure to psychological trauma and posttraumatic stress disorder in adults: a meta-analysis. Prog Neuropsychopharmacol Biol Psychiatry. 2010;34:1181–8.
Tse YC, Montoya I, Wong AS, Mathieu A, Lissemore J, Lagace DC, Wong TP. A longitudinal study of stress-induced hippocampal volume changes in mice that are susceptible or resilient to chronic social defeat. Hippocampus. 2014;24:1120–8.
Liu X, Yuan J, Guang Y, Wang X, Feng Z. Longitudinal in vivo diffusion tensor imaging detects differential microstructural alterations in the hippocampus of chronic social defeat stress-susceptible and resilient mice. Front Neurosci. 2018;12:613.
Lee J, Song Y, Won E, Bang M, Lee SH. higher rightward laterality of the hippocampal tail and its association with early trauma in panic disorder. Clin Psychopharmacol Neurosci. 2020;18:311–21.
Parker KJ, Schatzberg AF, Lyons DM. Neuroendocrine aspects of hypercortisolism in major depression. Horm Behav. 2003;43:60–6.
Bale TL. Stress sensitivity and the development of affective disorders. Horm Behav. 2006;50:529–33.
Hong HH, Uzel MI, Duan C, Sheff MC, Trackman PC. Regulation of lysyl oxidase, collagen, and connective tissue growth factor by TGF-beta1 and detection in human gingiva. Lab Invest. 1999;79:1655–67.
Liu XH, Zhang QY, Pan LL, Liu SY, Xu P, Luo XL, Zou SL, Xin H, Qu LF, Zhu YZ. NADPH oxidase 4 contributes to connective tissue growth factor expression through Smad3-dependent signaling pathway. Free Radic Biol Med. 2016;94:174–84.
We thank Dr. P. Greengard at Rockefeller University for critical discussions about this work and Drs. K-W Lee, Y-S Oh and Y. Kim and Ms. Sehrish Javed for their assistance with experiments. Gene expression analytic and imaging devices and analysis were supported by the KAIST Institute for the Bio Core Center.
Design, analysis, materials and writing this manuscripts and all expenses used for this research was supported by grants from the Brain Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2017M3C7A1029612); the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2019M3E5D2A01066259); and the National R&D Program for Cancer Control, Ministry of Health & Welfare, Republic of Korea (1720100).
Ethics approval and consent to participate
All animal care and experimental procedures were performed in accordance with protocols approved by the Animal Care and Use Committee of Korea Advanced Institute of Science and Technology (KAIST).
Consent for publication
The authors declare no competing interest or potential financial disclosure.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sample preparation for microarrays. a. The distribution of control and socially defeated mice depending on their sociability index (SI). Resilient mice (defined as those exhibiting SI greater than 1) were highly sociable, whereas susceptible mice (defined as those exhibiting SI scores less than 1) were social avoidant. Non-stressed, control mice from social interaction tests showed similar distribution of SI with resilient as previously described . The average SI values for non-defeat control, resilient, and susceptible mice were 1.5, 1.5, and 0.8, respectively. b. The SI values of mice that were selected for the microarray analysis. The average SI of the selected mice were 1.76 (n = 8), 1.61 (n = 8), and 0.5 (n = 7) for the control, resilient, and susceptible groups, respectively. c. Full list of the SI values of the mice used in the microarray analysis.
RT-qPCR confirmation of lateralized genes, categorized according to hemispheric dominance. Changes in the expression levels of left- and right-dominant genes were confirmed by RT-qPCR analysis. a. The Gapdh-normalized expression of the left-dominant genes, Cux2 and Wfs1, in the mPFC, presented as a bar graph. The expression levels of Cux2 and Wfs1 were higher in the left mPFC of susceptible mice, while the levels of these genes were similar in the left and right mPFC of non-stressed and resilient mice. b. The Gapdh-normalized expression of the right-dominant genes, Ctgf, Mbp, and Rprm, in the mPFC. The expression levels of Ctgf, Mbp, and Rprm were higher in the right mPFC of susceptible mice, while the levels of these genes were similar between the left and right mPFC of non-stressed and resilient mice. c. The left/right (L/R) expression ratios of the left-dominant genes, Cux2 and Wfs1, in the mPFC, presented as a bar graph. The ratios of Cux2 and Wfs1 were higher in susceptible mice than in non-stressed and resilient mice, whose ratios were similar. d. The L/R expression ratios of the right-dominant genes, Ctgf, Mbp, and Rprm, in the mPFC. The ratios of Ctgf, Mbp, and Rprm were lower in susceptible mice than non-stressed and resilient mice, whose ratios were similar.
Analysis of DEGs in the two hemispheres of the mPFC in mice with social defeat stress versus non-stressed mice. Significant DEGs with a FDR adjusted p-value cutoff of 0.05 are shown. AveExpr, averaged expression of microarray genes; t, moderated t-statistic; B, B-statistic.
List of 526 DEGs with FDR adjusted p-value cutoffs of 0.05 corresponding to genes presented in heatmap format in Fig. 1b according to their log2L/R values.
Comparison of averaged expression values of DEGs between the two mPFC hemispheres, presented in dot plot format in Fig. 1c.
The 526 genes exhibiting laterality listed, according to their log2L/R values. Laterality genes with positive values are expressed more highly in the left mPFC, and those with negative values are expressed more highly in the right mPFC. The genes used in rank order analysis in Fig. 1e.
Differences in the expression of the housekeeping genes, Actb(b-actin), Gapdh(glyceraldehyde-3-phosphate dehydrogenase), B2m ((b2microglobinbetween control and resilient/susceptible mice, summarized as log2FC values and p-values for each gene.
GO analysis using GeneMANIA. Genes with leftward or rightward laterality were used as input to GeneMANIA, and the associated GO terms were reported.
Summary of statistics
About this article
Cite this article
Chae, S., Hong, J., Kang, K. et al. Molecular laterality encodes stress susceptibility in the medial prefrontal cortex. Mol Brain 14, 92 (2021). https://doi.org/10.1186/s13041-021-00802-w
- Molecular laterality
- Social defeat stress
- Connective tissue growth factor (CTGF)
- Medial prefrontal cortex