. 2021 Jun 24;4:783. doi: 10.1038/s42003-021-02255-2

Intergenerational trauma transmission is associated with brain metabotranscriptome remodeling and mitochondrial dysfunction

Sammy Alhassen ^1,^#, Siwei Chen ^2,^3,^#, Lamees Alhassen ¹, Alvin Phan ¹, Mohammad Khoudari ¹, Angele De Silva ⁴, Huda Barhoosh ¹, Zitong Wang ¹, Chelsea Parrocha ¹, Emily Shapiro ¹, Charity Henrich ¹, Zicheng Wang ¹, Leon Mutesa ⁵, Pierre Baldi ^2,³, Geoffrey W Abbott ⁴, Amal Alachkar ^1,^3,^✉

PMCID: PMC8225861 PMID: 34168265

Abstract

Intergenerational trauma increases lifetime susceptibility to depression and other psychiatric disorders. Whether intergenerational trauma transmission is a consequence of in-utero neurodevelopmental disruptions versus early-life mother–infant interaction is unknown. Here, we demonstrate that trauma exposure during pregnancy induces in mouse offspring social deficits and depressive-like behavior. Normal pups raised by traumatized mothers exhibited similar behavioral deficits to those induced in pups raised by their biological traumatized mothers. Good caregiving by normal mothers did not reverse prenatal trauma-induced behaviors, indicating a two-hit stress mechanism comprising both in-utero abnormalities and early-life poor parenting. The behavioral deficits were associated with profound changes in the brain metabotranscriptome. Striking increases in the mitochondrial hypoxia marker and epigenetic modifier 2-hydroxyglutaric acid in the brains of neonates and adults exposed prenatally to trauma indicated mitochondrial dysfunction and epigenetic mechanisms. Bioinformatic analyses revealed stress- and hypoxia-response metabolic pathways in the neonates, which produced long-lasting alterations in mitochondrial energy metabolism and epigenetic processes (DNA and chromatin modifications). Most strikingly, early pharmacological interventions with acetyl-L-carnitine (ALCAR) supplementation produced long-lasting protection against intergenerational trauma-induced depression.

Subject terms: Molecular neuroscience, Emotion, Molecular medicine

Sammy Alhassen, Siwei Chen, et al. use mouse models to examine the effects of prenatal and postnatal stress on metabolomic and transcriptomic pathways in the brain. Their results suggest that altered mitochondrial metabolism may underlie trauma-induced behavioral deficits, and that correcting metabolism with ALCAR supplementation may protect against intergenerational transmission of traumatic stress.

Introduction

Intergenerational trauma increases lifetime susceptibility to depression and is a major risk factor for developing multiple neuropsychiatric disorders such as post-traumatic stress disorder (PTSD), autism spectrum disorder (ASD), and schizophrenia^1–4. Major depressive disorder (MDD) affects millions of people globally who have been subject to stress and trauma at some times of their lives. While stress deleteriously impacts humans worldwide and attacks all segments of societies at this moment in human history, it is disproportionately devastating to individuals in vulnerable situations such as pregnant women. Unlike other vulnerable conditions, the deleterious impacts of stress during pregnancy are doubled as they affect the mental health of the mothers and the unborn offspring. Human natural experiments provide evidence for the devastating health consequences in offspring as a result of exposure during pregnancy to existential and acute trauma such as war and natural disasters^5–7.

Whether intergenerational trauma transmission and its negative outcomes are a consequence of in utero fetal neurodevelopment disruptions or from poor maternal care by traumatized mothers is still largely ambiguous. This ambiguity is mainly due to the complexity, superposition, and inseparability of prenatal and postnatal mechanisms that are not mutually exclusive.

Converging evidence from human and animal studies supports that exposure to stress during pregnancy negatively affects maternal behavior towards, and the care of, offspring^8–10. On the other hand, maternal stress can affect offspring mental outcomes through prenatal programming, particularly epigenetic mechanisms causing permanent neuronal impairments and poor offspring outcomes that persist throughout life^11–13. Yet, it is difficult to test the mechanistic involvement of prenatal stress independent of stress effects on the mother care. The limitations in our understanding of the neurobiological and pathophysiological mechanisms of stress-induced depression, particularly the lack of measurable biochemical and molecular alterations, hamper effective pharmacological interventions to prevent or ameliorate the behavioral deficits.

Therefore, the goal of this work is to we examine the distinctive behavioral impacts and neurobiological mechanisms of in utero exposure alone (prenatal one-hit stress) versus impaired maternal care alone (postnatal one-hit stress) or combined exposure (prenatal and postnatal two-hits stress). We hypothesized that prenatal stress acts as a priming factor that synergistically interacts with subsequent early-life challenges (inadequate caregiving by traumatized mothers) to exacerbate the tendency for offspring of traumatized mothers to develop psychiatric disorders later in life.

To recapitulate intergenerational trauma experience, we modeled stress by exposing mouse dams in the last four days of their pregnancy to a predator odor. We applied cross-fostering procedures to evaluate the distinct contributions to adult behavioral deficits of prenatal trauma and subsequent postnatal inadequate caregiving. We used a comprehensive approach combining metabolomic, transcriptomic, and bioinformatic analyses to determine the molecular processes that are associated with stress-induced behavioral deficits. Based on our findings of the mechanisms involved, which revealed a role for mitochondria metabolism dysfunctions, we examined whether the behavioral deficits can be reversed by early and/or late supplementation with acetyl-L-carnitine (ALCAR).

Results

s → S pups display depressive-like behavior, social impairment, but normal cognitive functions

Compared to the offspring of control dams (c → C), the offspring of dams exposed to the predator scent (s → S) exhibited normal locomotor activities (Fig. 1a–c). The c → C and s → S mice travelled similar distances in the center and peripheral area in the open field assay (P > 0.05, unpaired t-test, Fig. 1d), and displayed similar times and number of entries to the open arm of the elevated maze (P > 0.05, unpaired t-test, Fig. 1e–h), indicating normal anxiety behavior in the s → S mice. In the social interaction assay, s → S mice displayed less interaction with the unfamiliar mice than did the c → C group (P < 0.001, two-way ANOVA, followed by Bonferroni posttest, Fig. 1i), and their D.I. was significantly lower (P < 0.001, unpaired t-test, Fig. 1j), indicating impaired sociability in the s → S mice. The s → S mice displayed depressive-like behavior, reflected by their higher immobility time in the forced swim test compared to the c → C mice (P < 0.05 unpaired t-test, Fig. 1k).

We investigated whether the prenatal stress caused cognitive dysfunction, by testing object memory, working memory, and contextual fear memory. In the novel object recognition assay, both groups spent more time exploring the new object than the old object (P < 0.01, two-way ANOVA, followed by Bonferroni posttest, Fig. 1l), the D.I. in the s → S mice was not different compared to the c → C mice (P = 0.07, unpaired t-test, Fig. 1m). In the spontaneous T-maze alternation assay, s → S and c → C mice showed a similar percentage of arm choice alternation and decision latency (P > 0.05, unpaired t-test, Fig. 1n, o). In the contextual fear-conditioning assay, the s → S group exhibited higher freezing behavior during the stimulus session (P < 0.05), but not in the retention session (P > 0.05, two-way ANOVA, followed by Bonferroni posttest, Fig. 1p). Together, these results suggest that s → S mice display normal cognitive function.

Mothers exposed to predator scent exhibited deficits in maternal behavior, and depressive-like behavior, but not locomotor activity or anxiety behavior

The exposure to the predator scent triggered fear-like responses in the mother mice, displayed by escape behavior (avoidance) and heart palpitations. Pup retrieval is a key indicator of maternal care. At postpartum day 2 (PPD2), stressed mothers (S) retrieved pups with significantly higher latencies (time to retrieve the first pup) and longer retrieval duration (retrieval of three pups) compared to the control mothers (C) (P < 0.05 and P < 0.001, for latency and duration respectively, Supplementary Fig. 1a, b). In the forced swim test on PPD5, the S mothers exhibited higher immobility time than C mothers (P < 0.001, Fig. 1c), indicating increased depressive-like behavior in these mice. Other behaviors of mothers including locomotor activity and anxiety behavior were normal (Supplementary Fig. 1d–g).

Prenatal stress results in alterations in brain metabolomics

Global brain metabolomics analysis of neonatal pups revealed that prenatal exposure to stress-induced alterations in 50 brain molecules (Fig. 2a), of which the mitochondrial metabolite 2-hydroxyglutaric acid (2-HG) displayed the highest increase in the brains of neonatal s → S pups (Fig. 2b, Supplementary Data 1). 2-HG is a hypoxia and mitochondrial dysfunction marker, and an epigenetic modifier^14–17. Therefore, its striking increase in the brains of the s → S pups within the first 24 h after birth indicates disruptions in mitochondrial respiratory functions and epigenetic processes. Alongside 2-HG elevated levels, the levels of tricarboxylic acid (TCA) cycle metabolites succinic acid and γ-aminobutyric acid (GABA) significantly increased in the brains of the neonatal s → S pups (Fig. 2b, c, Supplementary Data 1). Beside its role as a neurotransmitter, GABA regulates the TCA cycle through the GABA shunt¹⁸, in which 2-HG is formed as a by-product during conversion of GABA to succinic acid (Fig. 2c)^19–23. The increases in 2-HG, GABA, and succinate levels further support a mitochondrial metabolism dysfunction, particularly in the TCA cycle and GABA shunt, in the neonatal s → S pups.

Fig. 2 — a Heatmap of metabolites significantly changed in the brains of neonatal s → S mice (P < 0.05). b Volcano plot of metabolites significantly changed in the brains of neonatal s → S mice (P < 0.05). c Effect of prenatal exposure to stress on the tricarboxylic acid (TCA) and GABA shunt metabolism. Increased metabolites in green. d Changes of metabolites related to histidine metabolism, decreased metabolites in red. e Changes of metabolites related to lipid synthesis/metabolism in the neonatal s → S pups. Increased metabolites in green, and decreased metabolites in red. Phosphatidylcholines (PC), lysophosphatidylcholines (LPCs), phatidylethanolamines (PE), lysophosphatidylethanolamines (LPE), phosphatidylserine (PS), fatty acids (FA), arachidonic acid, glucoceramides (GlcCer), phosphtatidylglycerols (PG), and the decrease in sphingomyelins (SMs), diacylglycerols (DAG), and triglycerides (TAGs). f Pathway analysis of metabolites significantly changed in the neonatal MET pups (Supplementary Data 2 shows details of number of compounds, hits and p values). g Heatmap of metabolites significantly changed in the brains of adult s → S mice (P < 0.05). h Volcano plot showing differentially expressed metabolites in the brains from adult s → S mice. i Changes of metabolites related to lipid synthesis/metabolism in adult s → S pups. Increased metabolites in green, and decreased metabolites in red. j Pathway of metabolites significantly changed in the adult s → S pups (Supplementary Data 2 shows details of number of compounds, hits and p values). k Metabolites differentially expressed in neonatal and adult s → S mice with log₂(fold change) and −log10 (P value).

Other metabolites that were altered in the brains of the neonatal s → S pups include molecules involved in glycolysis and energy metabolism such as myo-inositol (increased) and palmitoyl-carnitine (decreased) (Fig. 2b Supplementary Data 1). The decrease in palmitoyl-carnitine reflects reduced oxidative capacity and mitochondrial ATP production^24–27. The mitochondrial dysfunction in neonatal s → S pups is further indicated by the marked decrease of carnosine (>3-fold)^28,29, which was associated with decreases of other histidine metabolism components (Fig. 2b, d, Supplementary Data 1,2). Carnosine (β‐alanyl‐L‐histidine) is a dipeptide that has mitochondrial pH-buffering capacity and enhancing activity on the respiratory chain complexes, and mitochondrial energy production^28–31. The depletion of these energy-related molecules, thus, suggests that the oxygen–glucose deprivation in the neonatal s → S pups’ brains creates a condition of excessive need for carnosine and its derivatives, to restore glycolysis and ATP production, and that the need exceeds the resources available of these molecules³². On the other hand, carnosine depletion might result in a reduction in the buffering of mitochondrial pH, thus contributing to the cellular and particularly mitochondrial acidity, associated with the production of 2-HG¹⁵.

Metabolites of monoamine neurotransmitters’ pathways displayed decreases in the brains of neonatal s → S pups including noradrenaline and dopamine precursor (tyrosine), serotonin precursor (tryptophan), and histamine-metabolism components (Fig. 2b).

The major lipidomic changes included increases in triacylglycerols (TAGs), cholesterol, ceramides (Cer), and decreases in acyl-carnitine (palmitoyl-carnitine), and membrane phospholipids (Fig. 2a, e). Bioinformatic analysis of metabolites altered in the brains of neonatal s → S pups substantiated enrichment of pathways associated with epigenetic processes, energy metabolism, mitochondrial functions, fatty acid oxidation, and complex lipid metabolism (FDR < 0.05, Fig. 2f, Supplementary Data 2).

In the adult brains, 2-HG again exhibited the highest change (>6-fold increase) in the brains of the s → S pups, indicating mitochondrial dysfunction (Fig. 2g, h, Supplementary Data 1). Other metabolites whose levels increased are molecules involved in energy and lipid metabolism/storage (myo-inositol, glycerol, TAGs, and cholesterol). There were profound abnormalities in the membrane, mitochondria, and signaling lipidomics. Except for cholesterol, glycerol, TAGs and diacylglycerol (DGs), which were increased in the S adult brains, all other lipids were decreased (Fig. 2h–i). Among the notable metabolites whose levels decreased in the s → S mice were aspartic acid, palmitoyl-carnitine, and S-Adenosyl-methionine (SAM).

Bioinformatic analysis revealed similar biochemical pathways enrichments in the brain’ metabolites of adult s → S mice to those significantly enriched in the neonatal s → S pups, indicating early but long-lasting metabolomic signatures produced by prenatal stress (FDR < 0.05, Fig. 2j, Supplementary Data 2). This is further supported by the finding that seven metabolites exhibited changes in the same directions in both neonatal and adult s → S mice: 2-HG, myo-inositol, N8-Acetylspermidine, palmitoyl-carnitine, cholesterol, LPC 20:4, and methylhistidine (Fig. 2k). Hence, metabolic pathways enriched in the brains’ metabolites of both neonatal and adult s → S mice included epigenetic/oncogenic action of 2-HG and succinate, mitochondrial GABA metabolism, carnitine synthesis, mitochondrial beta-oxidation of long-chain fatty acids, catecholamine biosynthesis, glutathione metabolism, inositol metabolism, histidine metabolism, methionine metabolism, and phosphatidylinositol-phosphate metabolism (Fig. 2f, j).

Prenatal stress results in distinct alterations in brain transcriptomics

Global brain transcriptomic analysis of differentially expressed genes (DEGs) revealed a subset of genes that exhibited a ≥ 1.5-fold change in the neonatal s → S mice (P < 0.05); Avp, Egr1, C1ql1, Fos, Crhbp (downregulated) and Baiap2l1, mt-Tqm, Irak1bp1, and 11 microRNAs (Mirs) (upregulated) (Fig. 3a, Supplementary Data 3). To verify the microarray results, we conducted RT-qPCR on six genes (Egr1, Fos, Avp, Vgf, Mcm6, and C1ql1) that were shown by microarray to display significant changes in the neonatal s → S mice (P < 0.05); the RT-qPCR results confirmed the microarray findings (Fig. 3b).

Fig. 3 — a Volcano plot showing differentially expressed genes in brains from neonatal s → S versus c → C mice. Baiap2l1 exhibited the largest expression change among the upregulated genes (green symbols, log(fold) > 0.5); Thy1 exhibited the largest expression change among the downregulated genes (red symbols log(fold) < −0.5). b Real-time qPCR validation of microarray analysis in the brains of neonatal mice. Relative global mRNA levels of Egr1, Fos, Avp, Vgf, Mcm6, and C1ql1 genes in the mice brain assessed by quantitative real-time PCR (n = 5). Results were normalized to Gapdh. Unpaired student test (Egr1: t = 6.346, P = 0.0002; Fos: t = 4.027, P = 0.0038; Avp: t = 6.807, P = 0.0001; Vgf: t = 5.172, P = 0.0009, Mcm6: t = 5.950, P = 0.0003; C1q11: t = 2.563, P = 0.0335): neonatal c → C vs neonatal s → S, *P < 0.05, **P < 0.01, ***P < 0.01. Data are presented as means ± S.E.M. c, d Significantly enriched biological pathways in the upregulated (c) and downregulated (d) genes in the neonatal mice with number of genes’ hits. Hypergeometric test, P < 0.05. e, f Differential (e) TFs and (f) RBPs, whose target gens are DEGs in the neonatal mice. STRING network visualization, using Markov clustering algorithm (MCL), with the number of their differential target genes. g Volcano plot showing differentially expressed genes (DEGs) in brains from adult s → S versus c → C mice. h, i Significantly enriched biological pathways in the upregulated (h) and downregulated (i) genes in the adult mice with number of genes’ hits. Hypergeometric test, P < 0.05. j, k Differential (j) TFs and (k) RBPs whose target gens are DEGs in adult s → S mice. STRING network visualization, using Markov clustering algorithm (MCL), with the number of their differential target genes. l Volcano plot of genes showing similar and opposite directions of change in neonatal and adult s → S pups. m–p Immunohistochemistry analysis of Arc in the brains of adult s → S and c → C mice; m, o Representative images of Arc-immunoreactivity in the (m) PCX, (o) Amyg, n, p quantification of Arc immuno-positive cells in the (n) PCX, (p) Amyg. Unpaired t-test, PCX: t = 2.83, P = 0.029, n = 3 C 5 s; Amyg: t = 6.165, P = 0.0008, n = 3 C 5 s; Scale bar = 50 µM. Values represent mean ± SEM. piriform cortex (PCX), and amygdala (Amyg). q, r Network graphs visualizing the metabolite-transcriptome interaction in (q) neonatal s → S pups and (r) adult s → S pups, using “spring” layout from the NetworkX package with Python. Each blue node presents a metabolite, and each red node presents the enzyme-coding gene. If a gene was involved in the metabolic reaction of a metabolite, an edge was drawn to connect these two nodes. The size of the node is proportional to its connectedness in the graph: the more connected the node is, the bigger its size.

Pathway analysis using hypergeometric test and Pathway Common Database and the ConsensusPathDB in neonatal s → S mice revealed the enrichment of mRNA processing, extracellular matrix organization, mitochondrial transport, GPCR signaling, glutamate transmission, glycolytic processes, ATP production, sympathetic system activation, and circadian regulation of gene expression (Fig. 3c, d Supplementary Data 4). After correction for multiple testing (FDR q < 0.05), mRNA processing remained enriched (Supplementary Data 4).

Aligned with the gene regulation enriched pathways, 57 transcription factors (TFs) and 13 RNA-binding proteins (RBPs) were differentially expressed in the s → S mice (Supplementary Data 5). Bioinformatic analyses using MotifMap and ChiPseq databases reveal that many of the DEGs in the neonatal s → S mice are indeed targets (differential target genes (DTGs)) for more than one of the differential TFs and RBPs (Fig. 3e, f, Supplementary Fig. 2a–h, Supplementary Data 5).

In the adult s → S mice, the top downregulated genes were Fos, Cacna2d3, Zfp874b, mt-Ts1, and several small nucleolar RNAs and ribosomal proteins (Fig. 3g, Supplementary Data 3). Among the top upregulated genes were mitochondrial genes and genes known to regulate energy metabolism: Glp1r, Tfap2b/d, Gpx3, Calb2, Alox8, Chrna3, Galr1, Ebf1, Paqr6, Gdpd2, Akr1b10, Syn3, Aqp4, Il6ra, and Tacr3 (Fig. 3g, Supplementary Data 3).

The hypergeometric test identified enrichment of pathways that are associated with cellular responses to stress, oxidative phosphorylation, electron transport chain, glycolysis, TCA cycle and respiratory electron transport, acyl-carnitine pathways, HIF-1-alpha transcription factor network, oxidative stress-induced senescence, circadian gene expression, DNA and histone methylation, sympathetic signaling in cardiovascular system, and epigenetic regulation of gene expression (P < 0.05, Fig. 3h, i, Supplementary Data 4). After correction for multiple testing (FDR q < 0.05), mitochondrial Electron Transport Chain remained significantly enriched (Supplementary Data 4).

Aligned with the gene regulation enriched pathways, 52 TFs (40 upregulated and 12 downregulated), and nine RBPs were differentially expressed in the brains of the adult s → S mice (Fig. 3j, k, Supplementary Data 4). Several differential genes were identified as targets for these differential TFs and RBPs (Fig. 3j, k, S3a–h, Supplementary Data 4).

Eighty-three genes followed the same change direction in neonatal and adult s → S mice, whereas 62 genes exhibited opposite alterations (Fig. 3l). The genes exhibiting the greatest degree of differential expression that changed in the same direction in neonatal and adult s → S mice were Baiap2l1 and Fos (Fig. 3l).

Given the alteration in many genes involved in synaptic plasticity and glutamate transmission in both neonatal and adult s → S mice, we analyzed relative expression of activity-regulated cytoskeleton-associated protein (Arc/Arg3.1), a critical gene for proper glutamatergic synaptic function and plasticity^33,34. We found that in the s → S adult pups, there were lower numbers of Arc positive cells in the amygdala, hippocampus, striatum, and piriform cortex, with no change in the frontal cortex and nucleus accumbens, compared to c → C adult pups (Fig. 3m–p, Supplementary Fig. 4a–h), though the total numbers of cells in these regions were not altered in s → S group, evidenced by the similar number of DAPI-stained nuclei.

Bioinformatic integration of transcriptomic and metabolomic data substantiates mitochondrial and epigenetic mechanisms

Our metabolic profiling was in good agreement with gene expression data. The activities of the metabolites and their corresponding genes involved in the metabolic reactions were linked together through known enzymatic reactions from HMDB. The linkage analysis of the changes in gene expression levels with the changes in metabolite levels in the neonatal pups confirm the master biological pathways including mRNA modification, response to stress and hypoxia, oxidative-reduction process, phospholipid metabolic process, and DNA demethylation. The molecular functions enrichment analysis identified enzymatic activities of 2-ketoglutarate-dependent dioxygenase, phospholipase, and cytidine deaminase (Fig. 3q, r, Supplementary Data 6). The linkage analysis of the DEGs with the changes in metabolite levels in the adult pups revealed master pathways including chromatin modification, lipid metabolism, and glutamate pathways. Enriched chromatin modification (histone methylation (H4-R3 and H3-K27)) was associated with enrichment in histone methyltransferase activity. Enriched lipid metabolism was linked with phospholipase activity, and phospholipid transporter activity (transferring acyl group) (Fig. 3q, r Supplementary Data 6).

Distinct behavioral impacts of in utero exposure vs. impaired maternal care

Since mothers exposed to stress during pregnancy exhibited impaired maternal care and increased depressive-like behavior, we applied cross-fostering procedures to examine the distinctive effects of the in utero exposure vs. impaired maternal care upon the behavioral phenotype in adulthood (Fig. 4a). In the cross-fostering behavioral studies, mice were tested only in social interaction and forced swim assays, since these were the only two behaviors that were impaired in s → S. We found that pups born to control mothers and reared by stressed mothers (c → S mice) displayed similar behavioral phenotypes to those born to and raised by their stressed biological mothers: the c → S mice displayed social deficits (P < 0.001, two-way ANOVA followed by Bonferroni posttest, Fig. 4b), and increased immobility time in the forced swim test (P < 0.01, one-way ANOVA, followed by Tukey posttest, Fig. 4c). These results support a crucial role for the early life environment in the observed phenotypes. However, the caregiving by control females of mice born to stressed females (s → C) did not reverse the behavioral deficits. Thus, although s → C mice spent more time with the unfamiliar mice than the empty cup, the difference in time spent was not statistically significant (Fig. 4b). Similarly, caregiving by control mothers slightly reduced immobility time of s → C mice in the forced swim test; however, the levels were not significantly different from those of s → S mice (P > 0.05, when compared to c → C or s → S mice, Fig. 4c).

Fig. 4 — a Schematic of the Cross-fostering experiment design. b Time mice spent interacting with empty cup and unfamiliar mouse in the social interaction assay (n = 8 c → C, 8 s → S, 8 c → S, 7 s → C). Two-way ANOVA: object effect (F_1,54 = 7.255, P = 0.0094), exposure to stress effect (F_3,54 = 0.1260, P = 0.9443), stress x object interaction (F_3,54 = 3.422, P = 0.0235) followed by Bonferroni post hoc test: empty cup vs unfamiliar mouse, **P < 0.01, ns not significant (P > 0.05). Data are presented as means ± S.E.M. c Time mice spent immobile in the forced swim assay (n = 10 c → C, 10 s → S, 10 c → S, 7 s → C). One way ANOVA revealed a significant stress effect (F_3,33 = 7.531, P = 0.0006) followed by Tukey post hoc test: c → C vs s → S vs c → S, **P < 0.01, ns not significant (P > 0.05). d Volcano plot showing DEGs in brains from adult c → S versus c → C mice. e Volcano plot showing DEGs in brains from adult s → C versus c → C mice. f, g Volcano plot of genes showing similar direction of change in c → S and s → S: (f) upregulated and (g) downregulated genes. h, i Volcano plot of genes showing similar direction of change in s → C and s → S: (h) upregulated and (i) downregulated genes. j Volcano plot showing DEGs in brains from adult s → C versus s → S mice. k, l Volcano plot of genes showing similar direction of change in c → S and s → C mice: (k) downregulated and (l) upregulated genes. m, n Enriched biological pathways in the DEGs of (g) c → S and (h) s → C mice with number of genes’ hits. Hypergeometric test, P < 0.05. o Venn diagram showing the overlap of DEGs in all mouse groups. p Overlap of DEGs in adult mouse groups that received one hit of stress, prenatally only (s → C), postnatally only (c → S) or two hits of stress (s → S). q Overlap of DEGs in neonatal and adult mice that received one hit of stress prenatally with or without postnatal hit of stress.

Distinct Transcriptomic impacts of in utero exposure vs. impaired maternal care

Transcriptomic analysis in adult c → S and s → C mice revealed 640 and 973 DEGs, respectively, compared to c → C mice, of which 108 and 69 genes in adult c → S and s → C mice, respectively, were altered in the same direction of the s → S DEGs (Fig. 4d, e, Supplementary Data 3).

Among the top DEGs that displayed the same direction of change in both s → S and c → S mice were Il6ra, Aqp4, Atp13a5, Akr1b10, Tspan18, and Paqr6 (upregulated) (Fig. 4f), and Scarna6, Snora34, Mcm6, Akr1c20, Mir1291, and Crh, (downregulated) (Fig. 4g). On the other hand, among the top DEGs that displayed similar changes in s → S and s → C mice were Il6ra, Baiap2l1, Akr1b10, Gstp1, and Lxn (upregulated) and Fos, Mcm6, Fbl, Stk-ps2, and Esf1 (Fig. 4h, i). Most interestingly, caregiving by control mothers reversed the alterations in 559 genes (Fig. 4j). Among the top DEGs in s → C group compared to s → S group are Gpx3, Galr1, Chrna3, Alox8, Calb2, Glp1r, Tfap2d, and Mab21l2, which were downregulated, and Atp5e, Crhbp, Cck, and Pop5, which were upregulated (Fig. 4j).

Mice exposed to one hit of stress either prenatally (s → C mice) or postnatally (c → S mice) shared 103 DEGs, which exhibited changes in the same directions in both s → C and c → S mice (Fig. 4k, l).

Hypergeometric distribution and Fisher’s exact test revealed common significant enriched pathways between s → C and c → S mice, including mRNA translation (initiation, elongation, and termination) and mRNA decay (FDR q < 0.05, Fig. 4m, n, Supplementary Data 4).

Further, 20 DEGs overlapped across the three adult groups that were exposed to one hit (s → C and c → S mice) or two hits (s → S) of stress (Fig. 4o, p). When considering the neonatal transcriptomic changes, we found only four genes that exhibited early- and long-lasting expression changes in the animals that were exposed at least prenatally to stress (neonatal s → S, s → S, and s → C) (Fig. 4q). Only one gene (Mcm6) displayed the same change in all animals that were exposed to stress prenatally and/or postnatally and was changed as early as 24-hr after birth.

Transcriptomic-overlapping between human MDD and the four stress mouse groups

The transcriptomics linkage analysis of the DEGs in the four groups of S mice with gene expression changes in human major depressive disorder (MDD) revealed 203, 137, 105, and 74 DEGs in the neonatal s → S, adult s → S, s → C, and c → S mice respectively that overlap with human MDD (Fig. 5, Supplementary Data 7). Fisher exact test revealed significant overlaps between human MDD and all the stress mouse groups, with P values of 1.33E–31, 3.34E–21, 1.09E–19, and 4.58E–11 for intersections between MDD and neonatal s → S, adult s → S, s → C, and c → S respectively.

Fig. 5 — a Top common DEGs and their log(fold change) in MDD patients and 24-h s → S mice with same direction changes. b Top common DEGs and their log(fold change) in MDD patients and adult s → S mice with same direction changes. c Top common DEGs and their log(fold change) in MDD patients and c → S mice with same direction changes. d Top common DEGs and their log(fold change) in MDD patients and s → C mice with same direction changes.

Acetyl-L-carnitine (ALCAR) protected against and reversed depressive-like behavior induced by prenatal trauma

Given that the behavioral deficits in the s → S mice were associated with a decrease in palmitoyl-carnitine, increase in 2-HG, and impairments in the mitochondrial functions, lipid β-oxidation and phospholipid metabolism, we tested whether the behavioral deficits could be rescued by an agent, acetyl-L-carnitine (ALCAR), that is known to restore these functions. ALCAR is known to elicit its antidepressant effects through a number of mechanisms including rescuing mitochondrial functions and lipid β-oxidation metabolism, improving glutamate transmission, and enhancing chromatin acetylation by providing an acetyl group to chromatin^35,36. ALCAR was supplemented in the drinking water of s → S mice either from weaning to adulthood (3–8 weeks), or for one week in adulthood (7–8 weeks) (Fig. 6a). ALCAR supplementation for 1 week during adulthood rescued the depressive-like behavior (P < 0.001, Fig. 6b) in s → S mice. One week after ALCAR cessation, however, the anti-depressant effect of ALCAR was diminished (Fig. 6c). ALCAR supplementation from weaning rendered s → S mice resistant to developing depressive-like behavior (P < 0.01, Fig. 6b). Strikingly, ALCAR antidepressant effects lasted 1 week after its cessation (P < 0.001, Fig. 6c).

Fig. 6 — a Schematic of the two ALCAR administration regimes, and the chemical structure of ALCAR. b Time mice spent immobile in the forced swim assay following the treatment of ALCAR for 5 weeks from weaning, and 1 week during the adult stage, (n = 9 c → C, 8 s → S, 7 s → S + ALCAR from Weaning, 10 s → S + ALCAR in Adulthood). One-way ANOVA: treatment effect (F_3,30 = 9.72, P = 0.0001) followed by Tukey post hoc test, **P < 0.01, ***P < 0.001, ns not significant (P > 0.05). c Time mice spent immobile in the forced swim assay 1 week after the cessation of ALCAR treatment, One-way ANOVA revealed a significant treatment effect (F_3,30 = 32, P = 0.0001) followed by Tukey post hoc test, ***P < 0.001, ns not significant (P > 0.05). Data are presented as means ± S.E.M.

Discussion

The long-term devastating impacts of prenatal traumatic stress are well established in humans and animals^1,3,4,37. Yet, the mechanistic involvement of prenatal stress independent of stress effects on maternal care is largely unknown.

Here, we present results demonstrating that exposure to trauma during pregnancy induces in the offspring long-lasting depressive-like behavior and social deficits, without affecting cognitive functions. Using cross-fostering experiments, we demonstrate that these behavioral deficits are associated with divergent and convergent mechanisms of both in-utero and early-life parenting environments. Thus, the exposure to postnatal early-life hit of stress, through raising of normal pups by traumatized mothers, produces a behavioral phenotype that is similar to that induced by a double-hit of stress (in pups born and raised by their biological traumatized mothers). Good caregiving by normal mothers of pups that were exposed to prenatal hit of stress, however, did not completely prevent the manifestation of trauma- induced behavioral deficits. Thus, our data demonstrate distinctive behavioral impacts of one-hit exposure to stress (either prenatally through or postnatally) or two-hits stress (prenatal and postnatal exposures).

Associated with the distinct behavioral deficits we found, through metabolomic, transcriptomic and bioinformatic analyses, mechanisms that involve stress- and hypoxia-response energy metabolic pathways, especially mitochondrial ATP production. These acute responses seem to have resulted in long-lasting adaptations in glycolysis, homeostasis of energy lipids, and epigenetic processes pertaining to DNA and chromatin modifications, as evidenced by the disruptions of these pathways in the brains of adult mice. We therefore propose a model through which stress exposure during pregnancy induces, in progeny, early and long-lasting mitochondrial metabolism dysfunctions and epigenetic changes, associated with social behavioral deficits (Fig. 7a). According to our model, the striking increase in the mitochondrial metabolite and epigenetic modifier 2-hydroxyglutarate (2-HG) in the brains of neonatal mice, whose mothers were exposed to extreme stress during pregnancy, likely forms the first step of the consequent stress-response events.

Fig. 7 — a Prenatal exposure to stress elicits depression through a number of mechanisms that may begin with hypoxia, which results in compromised glycolysis, energy deprivation, mitochondrial dysfunction (oxidative phosphorylation and ATP production), and altered epigenetic processes. b Possible mechanisms through which ALCAR may produce its antidepressant effects.

Exposure to extreme stress during pregnancy is known to activate the maternal fight-or-flight sympathetic system, and produce fetal reduction of oxygen (hypoxia), which might be detected in the infants several months postnatally^38–42. The mitochondrial metabolite 2-HG is an established marker for the hypoxia, defective mitochondria, decreased electron transport chain activity, and acidic cellular conditions^14,15,43,44. Accumulation of 2-HG in hypoxia occurs as a result of the TCA cycle dysfunction and, in turn, 2-HG inhibits electron transport and glycolysis to counterbalance the deleterious consequences of hypoxia⁴³. Therefore, 2-HG accumulation in the neonatal s → S pups, accompanied by succinate accumulation and activated GABA shunt, indicates a hypoxic state and a disruption of the mitochondrial metabolism. The hypoxic state in newborn s → S mice seems to trigger adaptation mechanisms, particularly in metabolic and energy demand^45–47. Accordingly, acute stress-hypoxia responses in the neonatal s → S mice affect glycolysis and shift mitochondrial-energy metabolism (disruption in ATP production, accumulation of TAGs and depletion of membrane phospholipids). In addition to its role in hypoxia, 2-HG is an epigenetic modifier, which regulate epigenetic programming through its competitive inhibition of 2-ketoglutarate (2-KG)-dependent enzymes⁴⁸. 2-KG-dependent enzymes known to be inhibited by 2-HG include TET(1-3) DNA demethylase enzymes, Jumonji family of histone lysine demethylases (JmjC-KDMs), and prolyl hydroxylases (PHDs 1–3), which hydroxylate proline in hypoxia-inducible factors (HIFs), as a post-translational modification^49–53. Through these mechanisms, 2-HG alters DNA and histone methylation^49,50,54,55, and increases HIFs accumulation, resulting in selective upregulation of HIF-target genes^56,57. In support of HIFs accumulation, we found lower levels in the neonatal s → S pups of 4-hydroxyproline, a marker for the degradation of the proline-containing HIFs⁵⁸. The reduced levels of 4-hydroxyproline in the neonatal s → S pups indicate decreased hydroxylation and, hence, lower degradation of HIFs, causing HIFs accumulation.

HIFs target genes include practically all genes of the glycolytic pathway^59–61. While these metabolic and epigenetic changes aim to protect the brain from hypoxia during development, their effects may persist and confer the susceptibility for psychiatric disorders including depression later in life. First, the perturbation of the glycolysis pathways during stress triggers the activation of GABA shunt as an alternative pathway to enhance energy production^23,62–64. Ultimately, the stress-associated depletion of the brain’s energy resources leads to alteration in GABA and glutamate transmission and results in impaired neuronal plasticity underlying depression^65,66, reflected by the decreased Arc expression in brain regions involved in emotion and stress. Second, the disruption of glycerolipids homeostasis exacerbate neural dysfunction and has acute and long-lasting impacts through directing the lipid flux toward energy storage rather than membrane extension or signaling. Lastly, while the early life 2-HG-triggered epigenetic modifications serve to offset mitochondrial dysfunction in response to stress⁶⁶, they may produce long-lasting adaptations by altering methylation states and the expression of key genes involved in neurodevelopment, neuron maturation and differentiation, axon-genesis, and synaptic plasticity^67,68.

Our findings indicate a two-hits mechanistic of convergence as well as divergence of some key metabolic, epigenetic and neuronal pathways induced by prenatal and early life trauma, which are involved in the pathophysiology of depression. This is supported by the overlap in specific sets of genes and pathway networks among s → S, s → C and c → S mice. Our finding of common genes in human depression and the four experimental animal groups (neonatal s → S, adult s → S, s → C and c → S) provides a translational aspect of our model and support for our convergent pathways hypothesis of the emergence of stress-induced depression. Strikingly, among the different stress mouse groups, the neonatal s → S mice exhibited the highest overlap in DEGs with human MDD. This finding strongly indicates the importance of our model in identifying biomarkers that serve as predictive markers in early life for the risk of the development of depression in adulthood, and in providing a potential target for early preventive approaches and therapeutic strategies.

Thus, a pharmacological intervention that restores the perturbed pathways at the appropriate time of life would be predicted to protect against depression. ALCAR possesses unique features that make it an ideal candidate, providing a rare opportunity for preventive rather than therapeutic intervention for depression. First, ALCAR has long been known to enhance mitochondrial function and facilitate ATP production^35,36. Second, ALCAR promotes transportation of fatty acids into the mitochondria for subsequent oxidation and ATP generation and, thus, corrects the lipid profile, and directs the flow of lipids toward membrane production and energy production and TAGs reduction^69,70. Third, ALCAR produces epigenetic modifications that involve histone acetylation (through providing acetyl group from ALCAR to chromatin)^35,36. Through its epigenetic actions, ALCAR regulates the expression of key genes crucial for synaptic plasticity^35,36. Fourth, ALCAR has been shown to correct glutamate transmission in animal models of stress^35,36. Lastly, and most importantly, several randomized clinical studies demonstrated the effectiveness of ALCAR supplementation to decrease depressive symptoms^71,72. Indeed, the antidepressant effects of ALCAR have been speculated to occur through more than one of these mechanisms including enhancing mitochondrial function and facilitate ATP production^35,36.

We found that ALCAR antidepressant effects outlasted its treatment end when it was administered early in life (at weaning time), but not when it was administered in adulthood. The timing and duration of ALCAR administration might be deterministic factors in the long-lasting anti-depressant effect of ALCAR, given that animals administered ALCAR from weaning time received the drug for 5 weeks, whereas those administered ALCAR in postnatal week 7 received it for 1 week. However, considering that 1 week cessation of ALCAR is sufficient as a washout period, the long-lasting effect of ALCAR seems to be more likely due to the administration timing (weaning time vs adulthood) rather than duration (1 week vs 5 weeks).

This long-lasting anti-depressant effect of ALCART after its cessation suggests epigenetic mechanisms, through which ALCART likely increases histone acetylation causing chromatin structure remodeling. Since prenatal stress induced in s → S mice changes in all pathways that can be improved by ALCAR (mitochondrial functions and lipid β-oxidation metabolism, glutamate transmission, and epigenetic modifications), and given the distinctive effects of 1 week vs 5 weeks of ALCAR administration, it is likely that ALCAR elicited its antidepressant effects on s → S mice via more than one mechanism (Fig. 7b). These findings provide a unique opportunity for utilizing ALCAR as a prophylactic supplementation to young children and adolescents to protect against depression in high-risk populations. Remarkably, ALCAR levels are low in depression patients, and the declines are greater in patients with a history of childhood trauma and emotional neglect⁷³.

Given the unique features of ALCAR, this natural supplement can represent an innovative and unique prophylactic and therapeutic strategy, should it be administered at the appropriate time of life. This has the potential to change the lives of millions of people who suffer from major depression or have the risk of developing this disabling disorder, particularly those in which the depression arose from prenatal traumatic stress.

Our study is, however, not without limitations. First, ideally, metabolomics and transcriptomics analyses should be conducted on specific brain regions. However, we started with the whole brain since several regions might be involved in stress-induced behavioral deficits including prefrontal cortex, hippocampus, amygdala, ventral tegmental area, nucleus accumbens, septum, and piriform cortex. Second, ideally, we would have wanted to test the effects of ALCAR on both depressive-like behavior and social interaction. However, the need to test the animals in both of these behavioral assays after administering the supplement for a specific time proved impossible to run the two assays on the same animals. Region-specific metabolomics and transcriptomic analyses as well as the effects of ALCAR on social interaction are the focus of our ongoing work.

In conclusion, we demonstrate that intergenerational trauma induces social deficits and depressive-like behavior through divergent and convergent mechanisms of both in utero and early-life parenting environments. We establish 2-HG as an early predictive biomarker for trauma-induced behavioral deficits and demonstrate that early pharmacological correction of mitochondria metabolism dysfunction by ALCAR can permanently reverse the behavioral deficits.

Methods

Animals and breeding procedures

Eight-week-old Swiss Webster mice were obtained from Charles River Laboratories (Wilmington, MA). Mice were maintained on a 12-h light/dark cycle (light on at 7:00 am). One male and two female mice were group-housed, and pregnant mice were individually housed from the fourteenth day of pregnancy (gestational day 14) until delivery. After weaning (postnatal day 21), male mice were group-housed separately (Fig. 1a). All experimental procedures were approved by the Institutional Animal Care and Use Committee of University of California, Irvine and were performed in compliance with national and institutional guidelines for the care and use of laboratory animals.

Predator scent (PS) exposure and cross-fostering

On gestational day 17, the pregnant mice were exposed to a predator scent (2,5-dihydro-2,4,5-trimethylthiazoline, or TMT, a constituent of fox urine) for 1 h on each of four consecutive days in a separate cage under the fume hood and were then returned to their home cage. Bedding mixed with a few drops of distilled water (Control, C) or 0.02% TMT (Stress, S) was placed in a 60 × 15 mm polystyrene plate placed in each cage to achieve scent exposure^74–77.

Maternal behavior: pups’ retrieval assay with postpartum mice

Maternal behavior was examined in predator-odor-exposed dams and control dams using pups’ retrieval assay⁷⁸. The assay was conducted on PPD2, in which, a maximum of 5 min of video recording time to calculate the mother’s latency and duration. The mother was temporarily removed from the cage, and pups were removed except for three pups, which were placed on each of the three corners of the cage (not the nest corner). The mother was then returned to her nest and the latency and duration of pups’ retrieval were recorded.

Cross-fostering mouse model and animal care

Pups (control, c or stress, s) of control (C) and stress (S) mothers were raised by their biological mother (c → C and s → S) or were cross-fostered to the opposite mother treatment (s → C and c → S, for control mothers with stress pups or stress mothers with control pups, respectively) within 24-h of birth. The pups remained with either their biological or foster mother until the weaning day (postnatal day 21). Subsequently, the male mice were then selected for the study and group-housed in groups ranging from 3 to 5.

Behavioral analyses

Male mice were tested from postnatal week 8 to week 13 with a battery of behavioral paradigms in the following order: locomotion and stereotypy/open field, social interaction, spontaneous T-maze alternation, novel object recognition, elevated plus maze, forced swim, contextual fear conditioning (Fig. 1a). The sequence of specific assays spaced by 3-6 days inter-assay interval was adapted from previous reports^79,80. In subsequent cross-fostering behavioral studies, mice were tested only in social interaction and forced swim assays. In ALCAR behavioral studies, mice were tested only in one behavioral assay (forced swim).

Locomotor activity

Mice were placed into a locomotion test chamber (Med Associates, Inc.) for 90 min as we previously described⁸¹. The first 30 min was allotted for the animal to habituate to the chamber before the test. The horizontal, vertical, and stereotypic activities were then recorded for the remaining 60 min and were then analyzed by the Activity Monitor 5 software (Med Associates, Inc.).

Open field assay

Prior to testing, the open field test chamber (40 × 40 cm, Med Associates, Inc.) was sectioned into a central zone, a 24 × 24 cm square in the middle of the test chamber, and a peripheral zone, the remaining outer area⁸². Mice were placed in the chamber and the time the animal spent within 10 min in the central or peripheral zones was recorded. Using the predetermined zone areas, the center-to-periphery exploration time ratio was assessed by the Activity Monitor 5 software (Med Associates, Inc.).

Elevated plus maze

A standard elevated plus-maze, made of grey Plexiglas, was placed in a sound-proof observation room with controlled light (200 Lux) on the central platform of the maze⁸³. Animals were given a 30-min period to habituate to the room before being tested. During the test, the mice were placed in the center of the plus facing an open arm and were given 5 min to explore. The behavior was recorded and scored by two independent observers blind to the animal treatments. The animals were scored based on time spent in the closed and open arms and the number of entries to the closed and open arms.

Social interaction

The social interaction test was performed using the three-chambered apparatus⁸². Two wire mesh cups were placed in the middle of both the right and left chamber (one on each side). In one of the cups on either side, a control mouse (an unfamiliar mouse of the same strain, gender, and age with no prior contact with the subject mouse) was placed. The subject mouse was placed in the middle chamber and was given 5 min to explore with the dividing doors closed. The dividing doors were then removed to allow the mouse to freely travel between the three chambers for 10 min. The duration of direct contact between the subject mouse and both the empty cup and the cup with the mouse was measured. The relative exploration time was recorded and expressed by a discrimination index: [D.I. = (T_mouse − T_empty)/(T_mouse + T_empty) × 100%]. The test was video recorded and analyzed by the ANY-MAZE software (Stoelting Co.).

Spontaneous T-maze alternation assay

The T-maze has three designated parts: the main stem, the side arms, and the starting area (AccuScan Instruments, Inc.)⁸². The mice were placed in the starting area behind a sliding door, blocking the main stem and side arms and were given a 30 s acclimation period before the start of each trial. After the acclimation period, the sliding door was removed, and the mice can choose either the left or right-side arm. After a choice (all four paws in the chosen arm), the sliding door was closed behind the animal, allowing the mouse to explore the chosen arm for 30 s, before being manually returned to the start area for the next trial. A total of eight trials were completed, allowing for a total of seven total possible alternations. The alternation percentage was calculated as 100 x (number of alternations/7). The time for the animal to decide on a side arm (decision latency) was also recorded.

Novel object recognition assay

The novel object recognition (NOR) assay consists of training and testing phases⁸¹. Prior to training, the mice were handled for 1–2 min a day for 3 days and were given 10 min a day for 3 days inside the empty experimental apparatus to habituate to the environment. During the training phase, mice were placed in the experimental apparatus, and were given 10 min to explore and examine two identical objects. Twenty-four hours later, one of the familiar objects was replaced with a novel object and the animals were given 5 min to explore the chamber. During both the training and testing phases, the duration and number of times the mice interact with the familiar and novel objects were recorded individually. The relative exploration time was recorded and expressed by a D.I.: Equation (1) [D.I. = (T_novel − T_familiar)/(T_novel + T_familiar) × 100%]. Tests were video recorded and analyzed by the ANY-MAZE software (Stoelting Co.).

Forced swim

The forced swim assay was performed as previously described⁸⁴. During the 6-min assay, the first 2 min of activity was discarded. For the last 4 min, mice were videotaped and the immobility time was recorded. The ANY-MAZE software was used to record and analyze immobility (Stoelting Co.).

Contextual fear conditioning assay

The Contextual Fear Conditioning Assay was performed as previously described⁸¹. The assay was split into a training and retention session. Prior to the training session, mice were handled for 1–2 min per day for 3 days. During the training session, mice were placed in the conditioning chamber (TSE Systems, Inc.), with a striped pattern, for 2.5 min, before receiving a 0.7 mA foot shock for 2 s. After an additional 30 s, the mice were returned to their home cages. The freezing behavior (lack of movement for at least 3 s in a 5 s interval) pre-shock and post-shock were measured. Twenty-four hours after training, the mice were placed in the chamber again, with the same striped pattern, for 5 min and the freezing behavior, without shock, and was assessed. Freezing behavior was scored as freezing (1) or not (0) for every 5 s interval and the percentage of freezing behavior was calculated as 100 x (number of freezing intervals/total intervals).

Acetyl-L-carnitine (ALCAR) treatment

The ALCAR treatment protocol was adapted from a previously reported protocol⁷⁰. The s → S mice received ALCAR (0.3%, Sigma Aldrich) in the drinking water either on postnatal day 21 (PPD 21 weaning day) and continued to week 8, or on week 7 and continued to week 8. Control groups received normal tap water. Forced swim tests were carried out on the last day of the treatment and 7 days after discontinuing ALCAR.

Immunohistochemistry

Immunohistochemistry methods to examine the Arc positive (Arc + ) neurons were conducted as we previously described⁸⁴. Mice were perfused under isoflurane anesthesia with 4% paraformaldehyde (PFA), and brains were removed, post-fixed in 4% PFA for overnight at 4 °C, and stored in 30% sucrose. Coronal sections were cut at 20 µm thickness, and three sections were used from each region of interest. Brain regions were defined according to their anatomy using the mouse brain atlas⁸⁵. Sections were blocked with 4% normal goat serum in PBS with 0.3% Triton X-100 for 1 h, and were then incubated in the blocking buffer that contains the primary antibodies (Arc 1:500, Sigma). The sections were then washed with PBS, incubated with the secondary antibody (1:500) for 1 h, then with 4′,6-diamidino-2-phenylindole (DAPI) for 5 min, and were then mounted on slides. Imaging was carried out using Leica Sp8 TCS confocal microscope (UCI optical biology core facility). Arc positive (Arc⁺) neurons were counted in the bilateral areas of each section, and the mean counts of three non-consecutive sections per brain of 3–5 brains were calculated. Cell counts were performed using ImageJ⁸⁶, and confirmed manually by two persons blind to the experiment conditions.

Brain metabolite analyses

Tissues of one brain hemisphere (cortical-subcortical/mouse) were collected from neonatal and adult mice, and extracted following the protocols first published in⁸⁷. Metabolic profiling was carried out by the West Coast Metabolomics Center (University of California Davis). Three metabolic platforms were profiled: (1) primary metabolites including hydroxyl acids, purines, pyrimidines, carbohydrates and sugar phosphates, amino acids, and aromatics, (2) lipids, and (3) biogenic amines and methylated and acetylated amines⁸⁸.

mRNA microarray analysis

Microarray experiments and analysis were performed as we previously described^81,89. Total RNA was extracted (Qiagen RNA extraction kit) from one brain hemisphere of the 24-h and 13 weeks’ old pups according to the manufacturer’s protocol. RNA samples with A260/A280 absorbance ratios between 2.00 and 2.20 were reverse-transcribed into cDNA and analyzed by “whole-transcript transcriptomics” using the GeneAtlas microarray system (Affymetrix) and manufacturer’s protocols. MoGene 2.1 ST array strips (Affymetrix) were used to hybridize to newly synthesized ss-cDNA. Each array comprised 770,317 distinct 25-mer probes to probe an estimated 28,853 transcripts, with a median 27 probes per gene. Only annotated genes were included in the differential analysis.

Real time quantitative polymerase chain reaction (RT-qPCR)

RNA was reverse-transcribed into cDNA using the Transcriptor First Strand cDNA synthesis kit (Roche Molecular Systems, Inc) following manufacturer’s instructions. RT-qPCR was performed using the LightCycler480 SYBR Green I Master kit (Roche Molecular Systems, Inc), and analyzed using the amplification reactions were carried out on a LightCycler® 480 instrument with gene specific primers. Samples from 5 animals were run in duplicates and averaged. The comparative CT method was used to obtain the Ct values from each RNA sample and normalized to a housekeeping gene (GAPDH)⁹⁰. The difference between the Ct values (ΔCt) of the gene of interest and the housekeeping gene was calculated for each sample and the difference in the ΔCt values between the experimental and control samples ΔΔCt was calculated. Fold change was obtained as 2^(−ΔΔCt)⁸¹. The following primers were used (5′–3′): MCM6: Forward CGACAGCTTGAGAGCATGATCC, Reverse TGACATCAGGCGTCTCTACACG; AVP Forward GCTACTTCCAGAACTGCCCAAG, Reverse CAGCAGATGCTTGGTCCGAAGC; VGF Forward CGGACTGGAACGAAGTACTC, Reverse GTTTGGGAAGGGGTGATAGG; C1QL1: Forward AG TATGTGGGCAGACCTCTGCA, CCAGCTTGATGAAGACCTCGTC; FOS: Forward GATGTTCT CGGGTTTCAACG, Reverse GGAGAAGGAGTCGGCTGG; EGR1: Forward GTCACCTCCGCC TCCCCCAACC, Reverse GGCGGCGATAAGAATGCTGAAGGA; GAPDH: Forward CATCCCA GAGCTGAAC, Reverse TCAGATGCCTGCTTCAC.

Bioinformatic analysis of metabolomics and transcriptomics, and metabolites-transcriptomic Integration

The bioinformatic analysis of the differential metabolites differential genes was conducted as we described previously⁸⁸. A differential analysis was performed between the s → S and c → C groups for both 24-h old pups and 13-week old adult mice groups using the Cyber-T program^91,92 to identify the top up- or downregulated metabolites, using p-value cutoff at 0.05. The metabolic pathway enrichment analysis was performed on the differentially expressed metabolites using Fisher exact test, which is based on a hypergeometric distribution, to identify enriched pathways from the Small Molecule Pathway Database^93,94. The inputs of the 2 × 2 contingency table used for the Fisher exact test were: number of metabolite hits, number of metabolite non-hits, number of non-hits in metabolites associated with the pathway, and number of background metabolites. Pathways with p-value less than 0.05 were considered enriched in the given list of metabolites. Differentially expressed genes in s → S and c → C groups for both neonatal and adult mice were identified using Transcriptome Analysis Console (TAC) with p-value cutoff at 0.05, and the transcriptomic data underwent further analyses. These analyses were: (1) pathway enrichment analysis using fisher exact test on pathways from the Pathway Common Database⁹³ and the ConsensusPathDB⁹⁴, and STRING analysis (https://string-db.org/); (2) upstream transcription factor binding site (TFBS) enrichment analysis of promoters using the promoters and their target genes identified by the MotifMap database^95,96 and the CHiPSeq database from ENCODE;⁹⁷ (3) identification of transcription factors (TFs) and RNA-binding proteins (RBP) and analyses of their differentially expressed downstream targets using MotifMap and CHiPSeq databases.

For Integrated metabolite-transcriptomic analysis, metabolites were identified by InChI Key or name in Human Metabolome Database (HMDB)⁹⁸. Enzymes and enzyme-coding genes involved in metabolic reactions of the metabolites were queried from the biological properties section of the database. The transcriptomic expression of those enzyme-coding genes in s → S and c → C groups was extracted from the microarray data. In this way, the activity of the metabolites and their enzyme-coding genes in the s → S and c → C groups were linked together for holistic analysis. Network graphs were constructed with nodes being a metabolite or a gene, and edges being a known reaction between the metabolite and enzyme coded by the gene.

Integrated mouse-human transcriptomic analysis

The DEGs in both adult and newborn mouse groups from the transcriptomic analysis were compared with DEGs in human depression⁹⁹. The overlaps between DEGs common to stress mice and human depression were computed for adult and newborn mouse groups using Fisher exact test.

Statistics and reproducibility

Statistical analyses of behavioral data were carried out using GraphPad Prism (GraphPad Software, Inc.). Data were presented as means ± S.E.M. Results were analyzed by unpaired student t-test, one-way, and two-way ANOVA followed by the appropriate post hoc comparisons, and P < 0.05 was considered statistically significant.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Supplementary information

Supplementary Information^{(2.2MB, pdf)}

42003_2021_2255_MOESM2_ESM.pdf^{(173.7KB, pdf)}

Description of Additional��Supplementary Files

Supplementary Data 1^{(80.1KB, xlsx)}

Supplementary Data 2^{(18.1KB, xlsx)}

Supplementary Data 3^{(1MB, xlsx)}

Supplementary Data 4^{(77.1KB, xlsx)}

Supplementary Data 5^{(30.8KB, xlsx)}

Supplementary Data 6^{(21KB, xlsx)}

Supplementary Data 7^{(48KB, xlsx)}

Supplementary Data 8^{(16.8KB, xlsx)}

Supplementary Data 9^{(27.1KB, xlsx)}

Supplementary Data 10^{(36.1KB, xlsx)}

Reporting Summary^{(317.8KB, pdf)}

Acknowledgements

The authors would like to thank Dr. Olivier Civelli for helpful discussions. We would like to thank the IIE-SRF for their support to Amal Alachkar. The work of AA was supported by the department of pharmaceutical sciences-UCI. The work of SC and PB was supported in part by NIH grant GM123558. The work of GWA was supported by NIGMS grant GM130377.

Author contributions

A.A. designed the experiments and wrote the manuscript; L.M. helped with designing the experiments; S.A., A.P., Z.W., C.P., C.H., M.K., and Z.W. conducted behavioral experiments on the pups; L.A. conducted behavioral experiments on the mothers; E.C. and H.B. conducted the immunostaining experiments, A.S. and G.A. conducted the transcriptomics analysis, S.C. and P.B. conducted the bioinformatic analysis.

Data availability

Our microarray transcriptomic data have been deposited into GEO public repository: https://www.ncbi.nlm.nih.gov/geo/ (accession number GSE171275). The raw data of transcriptomic and metabolomic studies are available in the Supplementary Data 1 and 3, and source data underlying figures are available in Supplementary Data 8. Other datasets generated during the current study are available from the corresponding author upon reasonable request.

Code availability

The source code is freely downloadable from: http://cybert.ics.uci.edu/.

Competing interests

The authors declare no competing interests.

Footnotes

Peer review information Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Sammy Alhassen, Siwei Chen.

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-021-02255-2.

References

1.Van den Bergh BR, et al. Antenatal maternal anxiety is related to HPA-axis dysregulation and self-reported depressive symptoms in adolescence: a prospective study on the fetal origins of depressed mood. Neuropsychopharmacology. 2008;33:536–545. doi: 10.1038/sj.npp.1301450. [DOI] [PubMed] [Google Scholar]
2.Yehuda R, et al. Maternal, not paternal, PTSD is related to increased risk for PTSD in offspring of Holocaust survivors. J. Psychiatr. Res. 2008;42:1104–1111. doi: 10.1016/j.jpsychires.2008.01.002. [DOI] [PMC free article] [PubMed] [Google Scholar]
3.Kinney DK, et al. Autism prevalence following prenatal exposure to hurricanes and tropical storms in Louisiana. J. Autism Dev. Disord. 2008;38:481–488. doi: 10.1007/s10803-007-0414-0. [DOI] [PubMed] [Google Scholar]
4.Meyer U, et al. Relative prenatal and postnatal maternal contributions to schizophrenia-related neurochemical dysfunction after in utero immune challenge. Neuropsychopharmacology. 2008;33:441–456. doi: 10.1038/sj.npp.1301413. [DOI] [PubMed] [Google Scholar]
5.Musanabaganwa C, et al. Burden of post-traumatic stress disorder in postgenocide Rwandan population following exposure to 1994 genocide against the Tutsi: A meta-analysis. J. Affect Disord. 2020;275:7–13. doi: 10.1016/j.jad.2020.06.017. [DOI] [PMC free article] [PubMed] [Google Scholar]
6.Fossion P, et al. Transgenerational transmission of trauma in families of Holocaust survivors: the consequences of extreme family functioning on resilience, sense of coherence, anxiety and depression. J. Affect Disord. 2015;171:48–53. doi: 10.1016/j.jad.2014.08.054. [DOI] [PubMed] [Google Scholar]
7.Stein AD, et al. Maternal exposure to the Dutch famine before conception and during pregnancy: quality of life and depressive symptoms in adult offspring. Epidemiology. 2009;20:909–915. doi: 10.1097/EDE.0b013e3181b5f227. [DOI] [PMC free article] [PubMed] [Google Scholar]
8.Zoubovsky SP, et al. Chronic psychosocial stress during pregnancy affects maternal behavior and neuroendocrine function and modulates hypothalamic CRH and nuclear steroid receptor expression. Transl. Psychiatry. 2020;10:6. doi: 10.1038/s41398-020-0704-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
9.Van den Bergh, B. R. H. et al. Prenatal developmental origins of behavior and mental health: The influence of maternal stress in pregnancy. Neurosci. Biobehav. Rev. 117, 26–64 (2020). [DOI] [PubMed]
10.Meek LR, et al. Effects of stress during pregnancy on maternal behavior in mice. Physiol. Behav. 2001;72:473–479. doi: 10.1016/S0031-9384(00)00431-5. [DOI] [PubMed] [Google Scholar]
11.Kingsbury M, et al. Stressful life events during pregnancy and offspring depression: evidence from a prospective cohort study. J. Am. Acad. Child Adolesc. Psychiatry. 2016;55:709–716 e2. doi: 10.1016/j.jaac.2016.05.014. [DOI] [PubMed] [Google Scholar]
12.Lehrner A, et al. Maternal PTSD associates with greater glucocorticoid sensitivity in offspring of Holocaust survivors. Psychoneuroendocrinology. 2014;40:213–220. doi: 10.1016/j.psyneuen.2013.11.019. [DOI] [PMC free article] [PubMed] [Google Scholar]
13.Darnaudery M, Maccari S. Epigenetic programming of the stress response in male and female rats by prenatal restraint stress. Brain Res Rev. 2008;57:571–585. doi: 10.1016/j.brainresrev.2007.11.004. [DOI] [PubMed] [Google Scholar]
14.Intlekofer AM, et al. Hypoxia Induces Production of L-2-Hydroxyglutarate. Cell Metab. 2015;22:304–311. doi: 10.1016/j.cmet.2015.06.023. [DOI] [PMC free article] [PubMed] [Google Scholar]
15.Intlekofer AM, et al. L-2-Hydroxyglutarate production arises from noncanonical enzyme function at acidic pH. Nat. Chem. Biol. 2017;13:494–500. doi: 10.1038/nchembio.2307. [DOI] [PMC free article] [PubMed] [Google Scholar]
16.Hunt RJ, et al. Mitochondrial stress causes neuronal dysfunction via an ATF4-dependent increase in L-2-hydroxyglutarate. J. Cell Biol. 2019;218:4007–4016. doi: 10.1083/jcb.201904148. [DOI] [PMC free article] [PubMed] [Google Scholar]
17.Laukka T, et al. Fumarate and succinate regulate expression of hypoxia-inducible genes via TET enzymes. J. Biol. Chem. 2016;291:4256–4265. doi: 10.1074/jbc.M115.688762. [DOI] [PMC free article] [PubMed] [Google Scholar]
18.Schousboe A, Sonnewald U, Waagepetersen HS. Differential roles of alanine in GABAergic and glutamatergic neurons. Neurochem Int. 2003;43:311–315. doi: 10.1016/S0197-0186(03)00017-2. [DOI] [PubMed] [Google Scholar]
19.Struys EA, et al. Metabolism of gamma-hydroxybutyrate to d-2-hydroxyglutarate in mammals: further evidence for d-2-hydroxyglutarate transhydrogenase. Metabolism. 2006;55:353–358. doi: 10.1016/j.metabol.2005.09.009. [DOI] [PubMed] [Google Scholar]
20.Ravasz D, et al. Catabolism of GABA, succinic semialdehyde or gamma-hydroxybutyrate through the GABA shunt impair mitochondrial substrate-level phosphorylation. Neurochem Int. 2017;109:41–53. doi: 10.1016/j.neuint.2017.03.008. [DOI] [PubMed] [Google Scholar]
21.Salminen A, et al. Hypoxia and GABA shunt activation in the pathogenesis of Alzheimer’s disease. Neurochem Int. 2016;92:13–24. doi: 10.1016/j.neuint.2015.11.005. [DOI] [PubMed] [Google Scholar]
22.Chowdhury GM, et al. Altered cerebral glucose and acetate metabolism in succinic semialdehyde dehydrogenase-deficient mice: evidence for glial dysfunction and reduced glutamate/glutamine cycling. J. Neurochem. 2007;103:2077–2091. doi: 10.1111/j.1471-4159.2007.04887.x. [DOI] [PubMed] [Google Scholar]
23.Kang TC, et al. Spatial and temporal alterations in the GABA shunt in the gerbil hippocampus following transient ischemia. Brain Res. 2002;944:10–18. doi: 10.1016/S0006-8993(02)02596-9. [DOI] [PubMed] [Google Scholar]
24.Kerner J, Hoppel C. Fatty acid import into mitochondria. Biochim Biophys. Acta. 2000;1486:1–17. doi: 10.1016/S1388-1981(00)00044-5. [DOI] [PubMed] [Google Scholar]
25.Longo N, Amat di San Filippo C, Pasquali M. Disorders of carnitine transport and the carnitine cycle. Am. J. Med Genet C. Semin Med Genet. 2006;142C:77–85. doi: 10.1002/ajmg.c.30087. [DOI] [PMC free article] [PubMed] [Google Scholar]
26.Furuno T, et al. Roles of long chain fatty acids and carnitine in mitochondrial membrane permeability transition. Biochem Pharm. 2001;62:1037–1046. doi: 10.1016/S0006-2952(01)00745-6. [DOI] [PubMed] [Google Scholar]
27.Branca D, et al. Stabilizing action of L-carnitine on the energy-linked processes of mitochondria isolated from perfused rat liver. Biochem Pharm. 1986;35:2839–2841. doi: 10.1016/0006-2952(86)90201-7. [DOI] [PubMed] [Google Scholar]
28.Corona C, et al. Effects of dietary supplementation of carnosine on mitochondrial dysfunction, amyloid pathology, and cognitive deficits in 3xTg-AD mice. PLoS One. 2011;6:e17971. doi: 10.1371/journal.pone.0017971. [DOI] [PMC free article] [PubMed] [Google Scholar]
29.Macedo LW, et al. Acute carnosine administration increases respiratory chain complexes and citric acid cycle enzyme activities in cerebral cortex of young rats. Mol. Neurobiol. 2016;53:5582–5590. doi: 10.1007/s12035-015-9475-9. [DOI] [PubMed] [Google Scholar]
30.Smith EC. The buffering of muscle in rigor; protein, phosphate and carnosine. J. Physiol. 1938;92:336–343. doi: 10.1113/jphysiol.1938.sp003605. [DOI] [PMC free article] [PubMed] [Google Scholar]
31.Harding JW, O’Fallon JV. The subcellular distribution of carnosine, carnosine synthetase, and carnosinase in mouse olfactory tissues. Brain Res. 1979;173:99–109. doi: 10.1016/0006-8993(79)91099-0. [DOI] [PubMed] [Google Scholar]
32.Baek SH, et al. Modulation of mitochondrial function and autophagy mediates carnosine neuroprotection against ischemic brain damage. Stroke. 2014;45:2438–2443. doi: 10.1161/STROKEAHA.114.005183. [DOI] [PMC free article] [PubMed] [Google Scholar]
33.Alberi L, et al. Activity-induced Notch signaling in neurons requires Arc/Arg3.1 and is essential for synaptic plasticity in hippocampal networks. Neuron. 2011;69:437–444. doi: 10.1016/j.neuron.2011.01.004. [DOI] [PMC free article] [PubMed] [Google Scholar]
34.Shepherd JD, et al. Arc/Arg3.1 mediates homeostatic synaptic scaling of AMPA receptors. Neuron. 2006;52:475–484. doi: 10.1016/j.neuron.2006.08.034. [DOI] [PMC free article] [PubMed] [Google Scholar]
35.Nasca C, et al. Role of the astroglial glutamate exchanger xCT in ventral Hippocampus in resilience to stress. Neuron. 2017;96:402–413 e5. doi: 10.1016/j.neuron.2017.09.020. [DOI] [PubMed] [Google Scholar]
36.Nasca C, et al. L-acetylcarnitine causes rapid antidepressant effects through the epigenetic induction of mGlu2 receptors. Proc. Natl Acad. Sci. USA. 2013;110:4804–4809. doi: 10.1073/pnas.1216100110. [DOI] [PMC free article] [PubMed] [Google Scholar]
37.Perroud N, et al. The Tutsi genocide and transgenerational transmission of maternal stress: epigenetics and biology of the HPA axis. World J. Biol. Psychiatry. 2014;15:334–345. doi: 10.3109/15622975.2013.866693. [DOI] [PubMed] [Google Scholar]
38.Morishima HO, Yeh MN, James LS. Reduced uterine blood flow and fetal hypoxemia with acute maternal stress: experimental observation in the pregnant baboon. Am. J. Obstet. Gynecol. 1979;134:270–275. doi: 10.1016/S0002-9378(16)33032-0. [DOI] [PubMed] [Google Scholar]
39.Vlassaks E, et al. The effects of fetal and perinatal asphyxia on neuronal cytokine levels and ceramide metabolism in adulthood. J. Neuroimmunol. 2013;255:97–101. doi: 10.1016/j.jneuroim.2012.09.011. [DOI] [PubMed] [Google Scholar]
40.Capani F, et al. Protein ubiquitination in postsynaptic densities after hypoxia in rat neostriatum is blocked by hypothermia. Exp. Neurol. 2009;219:404–413. doi: 10.1016/j.expneurol.2009.06.007. [DOI] [PubMed] [Google Scholar]
41.Grimaldi M, et al. Early changes in the synapses of the neostriatum induced by perinatal asphyxia. Nutr. Neurosci. 2012;15:103–110. doi: 10.1179/1476830511Y.0000000026. [DOI] [PubMed] [Google Scholar]
42.Morishima HO, Pedersen H, Finster M. The influence of maternal psychological stress on the fetus. Am. J. Obstet. Gynecol. 1978;131:286–290. doi: 10.1016/0002-9378(78)90602-6. [DOI] [PubMed] [Google Scholar]
43.Oldham WM, et al. Hypoxia-Mediated Increases in L-2-hydroxyglutarate Coordinate the Metabolic Response to Reductive Stress. Cell Metab. 2015;22:291–303. doi: 10.1016/j.cmet.2015.06.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
44.Nadtochiy SM, et al. Acidic pH Is a metabolic switch for 2-hydroxyglutarate generation and signaling. J. Biol. Chem. 2016;291:20188–20197. doi: 10.1074/jbc.M116.738799. [DOI] [PMC free article] [PubMed] [Google Scholar]
45.Wheaton WW, Chandel NS. Hypoxia. 2. Hypoxia regulates cellular metabolism. Am. J. Physiol. Cell Physiol. 2011;300:C385–C393. doi: 10.1152/ajpcell.00485.2010. [DOI] [PMC free article] [PubMed] [Google Scholar]
46.Murray AJ. Oxygen delivery and fetal-placental growth: beyond a question of supply and demand? Placenta. 2012;33:e16–e22. doi: 10.1016/j.placenta.2012.06.006. [DOI] [PubMed] [Google Scholar]
47.Smith TF, et al. Pre- and perinatal ischemia-hypoxia, the ischemia-hypoxia response pathway, and ADHD risk. Behav. Genet. 2016;46:467–477. doi: 10.1007/s10519-016-9784-4. [DOI] [PubMed] [Google Scholar]
48.Xu W, et al. Oncometabolite 2-hydroxyglutarate is a competitive inhibitor of alpha-ketoglutarate-dependent dioxygenases. Cancer Cell. 2011;19:17–30. doi: 10.1016/j.ccr.2010.12.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
49.Janke, R., Iavarone, A. T. & Rine, J. Oncometabolite D-2-Hydroxyglutarate enhances gene silencing through inhibition of specific H3K36 histone demethylases. Elife, 6, e22451 (2017). [DOI] [PMC free article] [PubMed]
50.Shim EH, et al. L-2-Hydroxyglutarate: an epigenetic modifier and putative oncometabolite in renal cancer. Cancer Disco. 2014;4:1290–1298. doi: 10.1158/2159-8290.CD-13-0696. [DOI] [PMC free article] [PubMed] [Google Scholar]
51.Bailey, P. S. J. & Nathan, J. A. Metabolic Regulation of Hypoxia-Inducible Transcription Factors: The Role of Small Molecule Metabolites and Iron. Biomedicines, 6, (2018). [DOI] [PMC free article] [PubMed]
52.Losman JA, et al. (R)-2-hydroxyglutarate is sufficient to promote leukemogenesis and its effects are reversible. Science. 2013;339:1621–1625. doi: 10.1126/science.1231677. [DOI] [PMC free article] [PubMed] [Google Scholar]
53.Turcan S, et al. IDH1 mutation is sufficient to establish the glioma hypermethylator phenotype. Nature. 2012;483:479–483. doi: 10.1038/nature10866. [DOI] [PMC free article] [PubMed] [Google Scholar]
54.Schvartzman JM, et al. 2-hydroxyglutarate inhibits MyoD-mediated differentiation by preventing H3K9 demethylation. Proc. Natl Acad. Sci. USA. 2019;116:12851–12856. doi: 10.1073/pnas.1817662116. [DOI] [PMC free article] [PubMed] [Google Scholar]
55.Chowdhury R, et al. The oncometabolite 2-hydroxyglutarate inhibits histone lysine demethylases. EMBO Rep. 2011;12:463–469. doi: 10.1038/embor.2011.43. [DOI] [PMC free article] [PubMed] [Google Scholar]
56.Koivunen P, et al. Transformation by the (R)-enantiomer of 2-hydroxyglutarate linked to EGLN activation. Nature. 2012;483:484–488. doi: 10.1038/nature10898. [DOI] [PMC free article] [PubMed] [Google Scholar]
57.Komotar RJ, et al. IDH1 and IDH2 mutations in gliomas and the associated induction of hypoxia-inducible factor and production of 2-hydroxyglutarate. Neurosurgery. 2010;66:N20–N21. doi: 10.1227/01.neu.0000369899.41915.67. [DOI] [PubMed] [Google Scholar]
58.Jaakkola P, et al. Targeting of HIF-alpha to the von Hippel-Lindau ubiquitylation complex by O2-regulated prolyl hydroxylation. Science. 2001;292:468–472. doi: 10.1126/science.1059796. [DOI] [PubMed] [Google Scholar]
59.Miska J, et al. HIF-1alpha is a metabolic switch between glycolytic-driven migration and oxidative phosphorylation-driven immunosuppression of tregs in glioblastoma. Cell Rep. 2019;27:226–237 e4. doi: 10.1016/j.celrep.2019.03.029. [DOI] [PMC free article] [PubMed] [Google Scholar]
60.Del Rey MJ, et al. Hif-1alpha knockdown reduces glycolytic metabolism and induces cell death of human synovial fibroblasts under normoxic conditions. Sci. Rep. 2017;7:3644. doi: 10.1038/s41598-017-03921-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
61.Agrawal A, et al. Normoxic stabilization of HIF-1alpha drives glycolytic metabolism and regulates aggrecan gene expression in nucleus pulposus cells of the rat intervertebral disk. Am. J. Physiol. Cell Physiol. 2007;293:C621–C631. doi: 10.1152/ajpcell.00538.2006. [DOI] [PubMed] [Google Scholar]
62.Seo JY, et al. Neuroprotection of ebselen against ischemia/reperfusion injury involves GABA shunt enzymes. J. Neurol. Sci. 2009;285:88–94. doi: 10.1016/j.jns.2009.05.029. [DOI] [PubMed] [Google Scholar]
63.Yogeeswari P, Sriram D, Vaigundaragavendran J. The GABA shunt: an attractive and potential therapeutic target in the treatment of epileptic disorders. Curr. Drug Metab. 2005;6:127–139. doi: 10.2174/1389200053586073. [DOI] [PubMed] [Google Scholar]
64.Kang TC, et al. The altered expression of GABA shunt enzymes in the gerbil hippocampus before and after seizure generation. Neurochem Int. 2003;42:239–249. doi: 10.1016/S0197-0186(02)00079-7. [DOI] [PubMed] [Google Scholar]
65.Morava E, Kozicz T. Mitochondria and the economy of stress (mal)adaptation. Neurosci. Biobehav Rev. 2013;37:668–680. doi: 10.1016/j.neubiorev.2013.02.005. [DOI] [PubMed] [Google Scholar]
66.Picard M, et al. An energetic view of stress: Focus on mitochondria. Front Neuroendocrinol. 2018;49:72–85. doi: 10.1016/j.yfrne.2018.01.001. [DOI] [PMC free article] [PubMed] [Google Scholar]
67.Trollmann R, Gassmann M. The role of hypoxia-inducible transcription factors in the hypoxic neonatal brain. Brain Dev. 2009;31:503–509. doi: 10.1016/j.braindev.2009.03.007. [DOI] [PubMed] [Google Scholar]
68.Cohen SS, et al. Impaired cerebral angiogenesis in the fetal lamb model of persistent pulmonary hypertension. Int J. Dev. Neurosci. 2014;38:113–118. doi: 10.1016/j.ijdevneu.2014.08.003. [DOI] [PMC free article] [PubMed] [Google Scholar]
69.Ferreira GC, McKenna MC. L-Carnitine and Acetyl-L-carnitine Roles and Neuroprotection in Developing Brain. Neurochem Res, 2017;42:1661–1675. doi: 10.1007/s11064-017-2288-7. [DOI] [PMC free article] [PubMed] [Google Scholar]
70.Cherix, A. et al., Metabolic signature in nucleus accumbens for anti-depressant-like effects of acetyl-L-carnitine. Elife, 9, (2020). [DOI] [PMC free article] [PubMed]
71.Wang SM, et al. A review of current evidence for acetyl-l-carnitine in the treatment of depression. J. Psychiatr. Res. 2014;53:30–37. doi: 10.1016/j.jpsychires.2014.02.005. [DOI] [PubMed] [Google Scholar]
72.Veronese N, et al. Acetyl-L-Carnitine Supplementation and the Treatment of Depressive Symptoms: A Systematic Review and Meta-Analysis. Psychosom. Med. 2018;80:154–159. doi: 10.1097/PSY.0000000000000537. [DOI] [PubMed] [Google Scholar]
73.Nasca C, et al. Acetyl-l-carnitine deficiency in patients with major depressive disorder. Proc. Natl Acad. Sci. USA. 2018;115:8627–8632. doi: 10.1073/pnas.1801609115. [DOI] [PMC free article] [PubMed] [Google Scholar]
74.Belzung C, et al. Behavioral and neurochemical changes following predatory stress in mice. Neuropharmacology. 2001;41:400–408. doi: 10.1016/S0028-3908(01)00072-7. [DOI] [PubMed] [Google Scholar]
75.Mashoodh R, Sinal CJ, Perrot-Sinal TS. Predation threat exerts specific effects on rat maternal behaviour and anxiety-related behaviour of male and female offspring. Physiol. Behav. 2009;96:693–702. doi: 10.1016/j.physbeh.2009.01.001. [DOI] [PubMed] [Google Scholar]
76.Mueller BR, Bale TL. Early prenatal stress impact on coping strategies and learning performance is sex dependent. Physiol. Behav. 2007;91:55–65. doi: 10.1016/j.physbeh.2007.01.017. [DOI] [PubMed] [Google Scholar]
77.St-Cyr S, et al. Maternal predator odour exposure programs metabolic responses in adult offspring. Sci. Rep. 2018;8:8077. doi: 10.1038/s41598-018-26462-w. [DOI] [PMC free article] [PubMed] [Google Scholar]
78.Alachkar A, et al. Inactivation of the melanin concentrating hormone system impairs maternal behavior. Eur. Neuropsychopharmacol. 2016;26:1826–1835. doi: 10.1016/j.euroneuro.2016.08.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
79.Paylor R, et al. The use of behavioral test batteries, II: effect of test interval. Physiol. Behav. 2006;87:95–102. doi: 10.1016/j.physbeh.2005.09.002. [DOI] [PubMed] [Google Scholar]
80.McIlwain KL, et al. The use of behavioral test batteries: effects of training history. Physiol. Behav. 2001;73:705–717. doi: 10.1016/S0031-9384(01)00528-5. [DOI] [PubMed] [Google Scholar]
81.Alachkar A, et al. Prenatal one-carbon metabolism dysregulation programs schizophrenia-like deficits. Mol. Psychiatry. 2018;23:282–294. doi: 10.1038/mp.2017.164. [DOI] [PMC free article] [PubMed] [Google Scholar]
82.Vawter MP, et al. Melanin concentrating hormone signaling deficits in schizophrenia: association with memory and social impairments and abnormal sensorimotor gating. Int J. Neuropsychopharmacol. 2020;23:53–65. doi: 10.1093/ijnp/pyz051. [DOI] [PMC free article] [PubMed] [Google Scholar]
83.Sanathara NM, et al. Melanin concentrating hormone modulates oxytocin-mediated marble burying. Neuropharmacology. 2018;128:22–32. doi: 10.1016/j.neuropharm.2017.09.008. [DOI] [PMC free article] [PubMed] [Google Scholar]
84.Phan J, et al. Mating and parenting experiences sculpture mood-modulating effects of oxytocin-MCH signaling. Sci. Rep. 2020;10:13611. doi: 10.1038/s41598-020-70667-x. [DOI] [PMC free article] [PubMed] [Google Scholar]
85.G. Paxinos, K. F. The Mouse Brain in Stereotaxic Coordinates (Academic Press, 2001).
86.Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods. 2012;9:671–675. doi: 10.1038/nmeth.2089. [DOI] [PMC free article] [PubMed] [Google Scholar]
87.Matyash V, et al. Lipid extraction by methyl-tert-butyl ether for high-throughput lipidomics. J. Lipid Res. 2008;49:1137–1146. doi: 10.1194/jlr.D700041-JLR200. [DOI] [PMC free article] [PubMed] [Google Scholar]
88.Chen S, et al. Metabolomic and transcriptomic signatures of prenatal excessive methionine support nature rather than nurture in schizophrenia pathogenesis. Commun. Biol. 2020;3:409. doi: 10.1038/s42003-020-01124-8. [DOI] [PMC free article] [PubMed] [Google Scholar]
89.Heerdt PM, et al. Transcriptomic analysis reveals atrial KCNE1 down-regulation following lung lobectomy. J. Mol. Cell Cardiol. 2012;53:350–353. doi: 10.1016/j.yjmcc.2012.05.010. [DOI] [PMC free article] [PubMed] [Google Scholar]
90.Nagasaki H, et al. Differential expression of the thyrostimulin subunits, glycoprotein alpha2 and beta5 in the rat pituitary. J. Mol. Endocrinol. 2006;37:39–50. doi: 10.1677/jme.1.01932. [DOI] [PubMed] [Google Scholar]
91.Kayala MA, Baldi P. Cyber-T web server: differential analysis of high-throughput data. Nucleic Acids Res. 2012;40:W553–W559. doi: 10.1093/nar/gks420. [DOI] [PMC free article] [PubMed] [Google Scholar]
92.Baldi P, Long AD. A Bayesian framework for the analysis of microarray expression data: regularized t -test and statistical inferences of gene changes. Bioinformatics. 2001;17:509–519. doi: 10.1093/bioinformatics/17.6.509. [DOI] [PubMed] [Google Scholar]
93.Cerami EG, et al. Pathway Commons, a web resource for biological pathway data. Nucleic Acids Res. 2011;39:D685–D690. doi: 10.1093/nar/gkq1039. [DOI] [PMC free article] [PubMed] [Google Scholar]
94.Kamburov A, et al. ConsensusPathDB–a database for integrating human functional interaction networks. Nucleic Acids Res. 2009;37:D623–D628. doi: 10.1093/nar/gkn698. [DOI] [PMC free article] [PubMed] [Google Scholar]
95.Liu Y, et al. MotifMap-RNA: a genome-wide map of RBP binding sites. Bioinformatics. 2017;33:2029–2031. doi: 10.1093/bioinformatics/btx087. [DOI] [PubMed] [Google Scholar]
96.Daily K, et al. MotifMap: integrative genome-wide maps of regulatory motif sites for model species. BMC Bioinforma. 2011;12:495. doi: 10.1186/1471-2105-12-495. [DOI] [PMC free article] [PubMed] [Google Scholar]
97.Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74. doi: 10.1038/nature11247. [DOI] [PMC free article] [PubMed] [Google Scholar]
98.Wishart DS, et al. HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. 2018;46:D608–D617. doi: 10.1093/nar/gkx1089. [DOI] [PMC free article] [PubMed] [Google Scholar]
99.Gandal MJ, et al. Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap. Science. 2018;359:693–697. doi: 10.1126/science.aad6469. [DOI] [PMC free article] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplementary Information^{(2.2MB, pdf)}

42003_2021_2255_MOESM2_ESM.pdf^{(173.7KB, pdf)}

Description of Additional Supplementary Files

Supplementary Data 1^{(80.1KB, xlsx)}

Supplementary Data 2^{(18.1KB, xlsx)}

Supplementary Data 3^{(1MB, xlsx)}

Supplementary Data 4^{(77.1KB, xlsx)}

Supplementary Data 5^{(30.8KB, xlsx)}

Supplementary Data 6^{(21KB, xlsx)}

Supplementary Data 7^{(48KB, xlsx)}

Supplementary Data 8^{(16.8KB, xlsx)}

Supplementary Data 9^{(27.1KB, xlsx)}

Supplementary Data 10^{(36.1KB, xlsx)}

Reporting Summary^{(317.8KB, pdf)}

Data Availability Statement

The source code is freely downloadable from: http://cybert.ics.uci.edu/.

PERMALINK

Intergenerational trauma transmission is associated with brain metabotranscriptome remodeling and mitochondrial dysfunction

Sammy Alhassen

Siwei Chen

Lamees Alhassen

Alvin Phan

Mohammad Khoudari

Angele De Silva

Huda Barhoosh

Zitong Wang

Chelsea Parrocha

Emily Shapiro

Charity Henrich

Zicheng Wang

Leon Mutesa

Pierre Baldi

Geoffrey W Abbott

Amal Alachkar

Abstract

Introduction

Results

s → S pups display depressive-like behavior, social impairment, but normal cognitive functions

Fig. 1. The s → S pups display depressive-like behavior, social impairment, but normal cognitive functions.

Mothers exposed to predator scent exhibited deficits in maternal behavior, and depressive-like behavior, but not locomotor activity or anxiety behavior

Prenatal stress results in alterations in brain metabolomics

Fig. 2. Prenatal exposure to stress causes changes in brain metabolites.

Prenatal stress results in distinct alterations in brain transcriptomics

Fig. 3. Prenatal exposure to stress causes alteration in brain transcriptomic signature.

Bioinformatic integration of transcriptomic and metabolomic data substantiates mitochondrial and epigenetic mechanisms

Distinct behavioral impacts of in utero exposure vs. impaired maternal care

Fig. 4. Cross-fostering effects on behavioral phenotype and transcriptomic profile.

Distinct Transcriptomic impacts of in utero exposure vs. impaired maternal care

Transcriptomic-overlapping between human MDD and the four stress mouse groups

Fig. 5. Integrated mouse-human transcriptomic analysis reveals overlap between genes whose expressions changed in mice exposed to stress and MDD patients.

Acetyl-L-carnitine (ALCAR) protected against and reversed depressive-like behavior induced by prenatal trauma

Fig. 6. Acetyl-L-carnitine (ALCAR) reverses depressive-like behavior after its cessation.

Discussion

Fig. 7. Model.

Methods

Animals and breeding procedures

Predator scent (PS) exposure and cross-fostering

Maternal behavior: pups’ retrieval assay with postpartum mice

Cross-fostering mouse model and animal care

Behavioral analyses

Locomotor activity

Open field assay

Elevated plus maze

Social interaction

Spontaneous T-maze alternation assay

Novel object recognition assay

Forced swim

Contextual fear conditioning assay

Acetyl-L-carnitine (ALCAR) treatment

Immunohistochemistry

Brain metabolite analyses

mRNA microarray analysis

Real time quantitative polymerase chain reaction (RT-qPCR)

Bioinformatic analysis of metabolomics and transcriptomics, and metabolites-transcriptomic Integration

Integrated mouse-human transcriptomic analysis

Statistics and reproducibility

Reporting summary

Supplementary information

Acknowledgements

Author contributions

Data availability

Code availability

Competing interests

Footnotes

Supplementary information

References

Associated Data

Supplementary Materials

Data Availability Statement

ACTIONS

PERMALINK

RESOURCES

Similar articles

Cited by other articles

Links to NCBI Databases