Next Article in Journal
The Metabolic Profile of Anchusa officinalis L. Differs According to Its Associated Arbuscular Mycorrhizal Fungi
Previous Article in Journal
Protective Effects of Bee Pollen on Multiple Propionic Acid-Induced Biochemical Autistic Features in a Rat Model
Previous Article in Special Issue
Metabolomic Associations of Asthma in the Hispanic Community Health Study/Study of Latinos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trans- and Multigenerational Maternal Social Isolation Stress Programs the Blood Plasma Metabolome in the F3 Generation

by
Joshua P. Heynen
1,2,†,
Eric J. Paxman
1,2,†,
Prachi Sanghavi
1,2,3,
J. Keiko McCreary
1,
Tony Montina
2,3,* and
Gerlinde A. S. Metz
1,2,*
1
Canadian Centre for Behavioural Neuroscience, Department of Neuroscience, University of Lethbridge, Lethbridge, AB T1K 3M4, Canada
2
Southern Alberta Genome Sciences Centre, University of Lethbridge, Lethbridge, AB T1K 3M4, Canada
3
Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, AB T1K 3M4, Canada
*
Authors to whom correspondence should be addressed.
Indicates co-first authors.
Submission received: 3 May 2022 / Revised: 12 June 2022 / Accepted: 15 June 2022 / Published: 22 June 2022
(This article belongs to the Special Issue Metabolomics of Complex Traits II)

Abstract

:
Metabolic risk factors are among the most common causes of noncommunicable diseases, and stress critically contributes to metabolic risk. In particular, social isolation during pregnancy may represent a salient stressor that affects offspring metabolic health, with potentially adverse consequences for future generations. Here, we used proton nuclear magnetic resonance (1H NMR) spectroscopy to analyze the blood plasma metabolomes of the third filial (F3) generation of rats born to lineages that experienced either transgenerational or multigenerational maternal social isolation stress. We show that maternal social isolation induces distinct and robust metabolic profiles in the blood plasma of adult F3 offspring, which are characterized by critical switches in energy metabolism, such as upregulated formate and creatine phosphate metabolisms and downregulated glucose metabolism. Both trans- and multigenerational stress altered plasma metabolomic profiles in adult offspring when compared to controls. Social isolation stress increasingly affected pathways involved in energy metabolism and protein biosynthesis, particularly in branched-chain amino acid synthesis, the tricarboxylic acid cycle (lactate, citrate), muscle performance (alanine, creatine phosphate), and immunoregulation (serine, threonine). Levels of creatine phosphate, leucine, and isoleucine were associated with changes in anxiety-like behaviours in open field exploration. The findings reveal the metabolic underpinnings of epigenetically heritable diseases and suggest that even remote maternal social stress may become a risk factor for metabolic diseases, such as diabetes, and adverse mental health outcomes. Metabolomic signatures of transgenerational stress may aid in the risk prediction and early diagnosis of non-communicable diseases in precision medicine approaches.

1. Introduction

The global rates of metabolic diseases have increased sharply over the past few decades [1]. Metabolic diseases in the mother influence the risk of adverse health outcomes in her children, including abnormal brain development, contributing to the high prevalence of impaired emotional and behavioural function and the risk of mental illness [2,3]. Mental illnesses, such as anxiety disorders or major depressive disorders, are the leading causes of disability worldwide [4]. The prevalence of neurodevelopmental disorders and mental illnesses is rapidly growing, increasing in parallel to the rise in metabolic diseases. Currently, the causal relationship between metabolic state in mother, child, and comorbid mental health issues remains unclear [5,6]; however, both conditions are linked to maternal stress and an adverse prenatal environment [7].
Maternal stress affects offspring physiology and neurodevelopment by programming the fetal hypothalamic–pituitary–adrenal (HPA) axis [7,8,9,10]. Fetal HPA axis programming can raise the lifetime risk for altered glucose and insulin metabolism, adiposity [11,12], anxiety disorders, and depression [13]. In addition, maternal anxiety and prenatal stress have been linked to low birth weight [14,15] and the development of metabolic syndromes such as type 2 diabetes, cardiovascular disease, and obesity later in life [16]. Notably, a lack of social support and social isolation represent significant stressors for pregnant women [17]. During the COVID-19 pandemic, perceived social isolation had the largest effects on the substantially elevated symptoms of anxiety (59%) and depression (37%) experienced by pregnant women in Canada [18].
Prenatal stress has been linked to altered DNA methylation marks [19] and microRNA signatures [20,21,22] in both blood and brain tissues. These epigenetic alterations enable the inheritance of metabolic and behavioural phenotypes to subsequent generations [23,24]. Thus, epigenetic markers linked to experience-dependent insulin resistance and diabetes [25,26] and altered affective states [27,28,29] may propagate across generations independent of changes in DNA sequence and affect the aging phenotype [30,31,32]. In fact, the transgenerational phenotype of high anxiety frequently coincides with altered glucose metabolism and other metabolic variations [28,29,30,33,34,35]. Furthermore, work in human cohorts of prenatal maternal stress caused by a natural disaster has shown metabolic alterations that indicate a higher risk of metabolic and cardiovascular disease [36]. Moreover, metabolic risk factors are among the most common causes of noncommunicable diseases [37], and stress critically contributes to metabolic risk [38,39]. Nevertheless, it has been difficult to disentangle the link between stress and metabolic disease risk.
Because upstream epigenetic regulation driven by prenatal or remote ancestral stress is reflected in downstream cellular metabolic functions, the metabolome, i.e., the sum of all metabolites in an organism, arguably represents a direct result of transgenerational programming [34,40,41] and provides insight into the metabolic underpinnings of epigenetically heritable diseases. Here, we used high-resolution proton nuclear magnetic resonance (1H NMR) spectroscopy and supervised and unsupervised machine learning approaches to probe robust metabolic signatures in rat blood plasma generated by transgenerational and multigenerational maternal stress in offspring from the third filial generation (F3). The models allow differentiation between the metabolomic impact of transgenerational epigenetic inheritance versus cumulative generational prenatal stress on adverse health outcomes. Both stress models were previously linked to preterm birth and metabolic disease [24], anxiety-like and depression-like symptoms [29,32,35], impaired neurodevelopmental trajectories [42,43], and accelerated biological aging [30,31]. This study reveals clinically accessible peripheral markers that may provide insight into metabolic pathways linked to the programming of adverse health outcomes and altered behaviour based on trans- and multigenerational maternal stress in rats.

2. Results

2.1. Trans- and Multigenerational Stress Generate Unique Metabolic Profiles

PCA was utilized to identify patterns in a subset of metabolites that were deemed significant using either a Mann–Whitney U test (p < 0.05) or VIAVC analysis. Metabolite profiles of the transgenerational prenatal stress (TPS) and multigenerational prenatal stress (MPS) groups were clearly separated from CONT (Figure 1A,B). Despite some overlap between the TPS and MPS profiles compared to CONT, data in Figure 1C demonstrate significant separation between the TPS and MPS groups and minimal confidence interval overlap.
There were 144 total bins created. Of these bins, TPS vs. CONT, MPS vs. CONT, and TPS vs. MPS revealed eight, eight, and thirteen significantly altered bins, respectively. Compared to CONT, both the TPS and MPS groups exhibited upregulated formate and creatine phosphate, as well as downregulated glucose (Table 1). The TPS comparison with the CONT group showed that leucine, isoleucine, alanine, 2-oxoisocaproate, and 2-hydroxyisovalerate were downregulated. Unique to the MPS group, 3-methylxanthine, threonine, and tyramine were upregulated compared to CONT, whereas betaine was downregulated (Table 1). Lastly, a comparison of TPS to MPS showed upregulation of succinate, creatine, and tyramine and downregulation of citrate, lactate, choline, alanine, and serine (Table 1).

2.2. Trans- and Multigenerational Stress Differentially Program Amino Acid Metabolism Pathways

Metabolite set enrichment analysis (MSEA) of the TPS compared to the CONT group (Figure 2A) indicated significant differences in valine, leucine, and isoleucine (branch chain amino acids) biosynthesis (p < 0.01) and degradation (p < 0.01). In addition, the glucose-alanine cycle (p < 0.01) and alanine metabolism were altered in the MSEA. Pathway topology supported the effects on branched-chain amino acid synthesis (p < 0.01) and degradation (p < 0.01) pathways but did not reveal significant pathway hits in the glucose-alanine cycle or alanine metabolic pathways (Figure 2B). The aminoacyl tRNA pathway (p < 0.01), however, was also shown to be affected by TPS.
MSEA revealed a significant impact on glycine, serine, and threonine metabolism (p = 0.013) in the MPS group when compared to CONT (Figure 3A). The pathway topology analysis reflected this significance, showing changes to the glycine, serine, and threonine metabolic pathway (p < 0.01; Figure 3B). No other pathways were found to be altered in this group, although betaine metabolism was somewhat significant in the metabolite set enrichment (p = 0.07).
When the impacts of trans- versus multigenerational stress were compared (Figure 4), the metabolite set enrichment analysis found the tricarboxylic acid (TCA) cycle (p < 0.01), methionine metabolism (p < 0.01), glycine, serine, threonine metabolism (p < 0.01), and alanine metabolism (p = 0.05) to be significantly affected. The TPS versus MPS topology analysis (Figure 4B) showed differential activation of the glycine, serine, and threonine metabolic pathway (p < 0.01) to be most significantly affected, followed by the TCA cycle (p < 0.01), the alanine, aspartate, and glutamate metabolic pathway (p < 0.01), the cyanoamino acid pathway (p = 0.03), the methane metabolic pathway (p = 0.05), and the aminoacyl tRNA pathway (p = 0.05).

2.3. Trans- and Multigenerational Stress-Induced Shifts in Energy Metabolism Are Associated with Altered Exploratory Behaviours

Table 2 presents the Spearman’s rank correlations of creatine phosphate and leucine/isoleucine to open field task metrics. Open field task scores can be found in the Supplementary Materials (Table S1). These data revealed a significant relationship between vertical exploratory time in the open field task and the concentration of creatine phosphate in the blood plasma metabolome in both TPS (Rho = −0.730, p = 0.007) and MPS rats (Rho = −0.691, p = 0.013). In the MPS group, creatine phosphate was also significantly correlated with the total distance traveled (Rho = −0.712, p = 0.009), the number of vertical moves such as rears (Rho = −0.636, p = 0.03), and the distance traveled in the center (Rho = −0.635, p = 0.003) of the open field arena (Table 2). In general, less exploratory behaviour was linked to higher creatine phosphate levels.
TPS, when compared to CONT, also displayed significant correlations between leucine/isoleucine and total distance (Rho = 0.646, p = 0.02), vertical time (Rho = 0.66, p = 0.02), and central distance (Rho = 0.596, p = 0.02) in the open field task (Table 2). MPS, when compared to CONT, showed significant correlations between leucine/isoleucine in the blood plasma and the distance traveled in the center of the open field task (Rho = 0.656, p = 0.02; Table 2). In general, lower exploratory activity was linked to reduced leucine/isoleucine levels.

3. Discussion

Metabolic risk factors are among the most common causes of noncommunicable diseases, and stress critically contributes to metabolic risk. This study shows that ancestral maternal social isolation stress induces distinct and robust metabolic profiles in adulthood, which are characterized by significant changes in energy metabolism, such as upregulated formate and creatine phosphate, and downregulated glucose. The plasma metabolome induced by transgenerational stress differs from recurrent maternal stress in terms of energy and amino acid metabolisms, cellular signaling, and immunoregulation. The findings indicate that maternal stress may lead to long-term metabolic adjustments and risks across generations of offspring, leading to altered tissue function and potentially contributing to an altered behavioural phenotype. The findings suggest that ancestral maternal stress may be a risk factor in idiopathic insulin resistance, diabetes, obesity, and other chronic and non-communicable diseases.
Substantial metabolic changes may explain the significant consequences of maternal stress on neuroendocrine and cardiometabolic functions, brain development, and emotional and cognitive impairments observed in the exposed offspring [10,22,44,45]. Moreover, it was shown that paternal stress can result in transgenerational programming of a metabolic phenotype in the F4 generation of the male lineage [34,41]. Here, we also show that trans- and multigenerational maternal stress causes distinct metabolic profiles in blood plasma, providing potentially predictive signatures of adverse health outcomes in the offspring.
Both TPS and MPS altered 21 metabolites known to be present in blood plasma. Among them, upregulation of creatine phosphate and formate represented a primary marker of ancestral stress. Creatine phosphate plays a vital role in ATP regeneration in skeletal muscle and is a rate-limiting pathway in muscle performance enabling brief, high-power activity [46] or meeting higher energy demand in response to stress [47]. On the other hand, elevated formate may trigger a metabolic switch from low to high adenine nucleotide levels, increasing the rate of glycolysis and intracellular lactate levels [48]. The involvement of these metabolites in both TPS and MPS suggests that ancestral stress may program a permanently upregulated energy metabolism, which supports a link between programmed stress vulnerability and documented low birth weights, lifetime diabetes and obesity risk, and other metabolic abnormalities in offspring exposed to stress in utero [10,49,50].
Metabolic abnormalities at birth are also considered a risk factor for mental disorders such as schizophrenia [51], depression [52], and attention deficit hyperactivity disorder [14,53]. Here we show that ancestral stress alters critical energy metabolism pathways, such as branched-chain amino acid (BCAA) biosynthesis and the tricarboxylic acid cycle, which mediate neurological and behavioural outcomes [54]. The proteinogenic BCCAs leucine, isoleucine, and valine are essential amino acids involved in glucose metabolism, neurotransmitter synthesis, and neuronal function. Thus, disruption in BCAA metabolism has been linked to abnormal brain development [54] and lifetime risk of neurodegenerative disorders, such as Alzheimer’s disease [55]. These mechanisms also implicate BCAAs in the regulation of immune functions, as nutritional interventions of both supplemented and reduced BCAA intake have been shown to affect metabolic and inflammatory functions [56,57].
When interfering and non-significant metabolites were excluded, the unsupervised PCA showed clear class separation between TPS and MPS, indicating metabolic profiles unique to each stress lineage. Interestingly, citrate and lactate were found to be significantly downregulated in the TPS group as opposed to MPS, but neither had a significantly altered citrate level compared to non-stressed controls. Both citrate and lactate are involved in energy generation, particularly in the production of ATP in the tricarboxylic acid (TCA) cycle, and in muscular ATP regeneration. If lactate levels are low in TPS animals, then metabolic acidosis and muscular fatigue may occur more quickly [58]. Similarly, low citrate levels in MPS animals may indicate energy deficiency [59]. By contrast, a relative increase in lactate in MPS animals, indicative of stress hyperlactataemia, is also secondary to anaerobic glycolysis induced by tissue hypoperfusion or hypoxia [60]. While here it may serve adaptive functions, in disease states stress hyperlactataemia is a reliable predictor of mortality [61]. This reflects the phenotype of MPS animals, whose metabolic and mental health impairments [24,31] and higher morbidity [30] exceed those of TPS animals.
TPS rats also revealed significantly downregulated choline and serine levels compared to the MPS group. Choline and its metabolites are required for several physiological functions, such as cell signaling and cholinergic transmission [62], and may be involved in blood lipid regulation and the development of cardiovascular disease and cognitive decline in aging [63]. Serine plays a central role in cellular proliferation and in nervous system development and functioning [64]. Altered levels of serine in patients with psychiatric disorders and neurological abnormalities underscore its importance in brain development and function [65]. Abnormal levels of these metabolites may help to explain anxiety-like behaviours in ancestrally stressed offspring, but it remains unclear why TPS is linked to lower levels than MPS. It has been hypothesized that cumulative stress across generations is mediated by an adaptive response to the constant presence of a stressor [66], which at times may facilitate partial physiological adaptation and resilience to recurrent stress in MPS animals [43,67].
Overall, ancestral stress had the largest impact on protein biosynthesis, energy generation, and immunity. Thus, the formation and breakdown of amino acids, as well as energy metabolism, are likely involved in the physiological and behavioural response to both TPS and MPS. These results are reinforced by the correlational analysis demonstrating associations between particular metabolites with anxiety-like behaviours. Previous studies have associated reduced exploratory behaviour, such as vertical rearing, with higher levels of anxiety [68]. Furthermore, trans- and multigenerationally stressed rats showed higher stress sensitivity and precocious onset of motor hyperactivity and risk assessment behaviours [69]. Here, anxiety-like behaviours, including reduction in the amount of vertical rearing and total and central distance traveled in the open field task, were associated with higher circulating creatine phosphate levels in the MPS group. Notably, the association between increased creatine phosphate and reduced rearing activity (vertical time), which is a behaviour linked to increased anxiety and psychomotor inhibition [70], was observed in both the MPS and TPS groups. However, many of the psychophysiological links to anxiety-like behaviours in the open field task are tenuous [71]. Nevertheless, it remains to be further determined if creatine phosphate concentrations may serve as a marker of anxiety risk in populations directly exposed to or with a family history of prenatal stress.
BCAAs such as leucine and isoleucine provided further insight into the metabolic underpinning of anxiety-like behaviours. More time spent rearing was linked to greater leucine and isoleucine concentrations in rats from the TPS group. Greater levels of leucine and isoleucine were associated with more distance travelled in the center of the open field in the MPS group. This change was mirrored in the TPS lineage with greater levels of leucine and isoleucine translating to greater total and central distance travelled in the open field task, linking remote ancestral stress to altered cellular energy metabolism and behavioural change. Thus, higher levels of leucine and isoleucine may be linked to higher stress resilience, whereas lower levels of these BCAAs may serve as markers of anxiety-like behaviours. Lower levels of BCAAs have also been linked to adverse metabolic and inflammatory functions [56,57] in addition to poorer stroke recovery, increased risk of renal disease, and type 2 diabetes [72,73,74,75]. Causal inferences of BCAA metabolism in brain function and behaviour require further investigation.
The present downstream cellular metabolic changes may reflect differential underlying epigenetic changes that are potentially heritable. The TPS generation of the maternal lineage in particular excludes direct exposure effects and unambiguously suggests heritable epigenetic transgenerational mechanisms for offspring three generations removed from maternal stress [76,77] to regulate cellular metabolism. Epigenetics may mediate the relationship between genotype and internal and external environments [33] and cellular metabolic changes may directly reflect changes in these upstream regulatory pathways. Furthermore, gestational stress may impact the developing child through the placenta or through direct or indirect effects on the fetal brain, potentially contributing to the differences observed between TPS and MPS cohorts. However, the current findings must be interpreted carefully, as this study only involved males and only included a small number of control animals.
The present study used a novel approach to quantify outcomes of ancestral stress and HPA axis dysregulation by taking a metabolomic perspective and providing possible mechanisms for health outcomes in an animal model. The metabolomic biomarkers identified in this study may provide potential biomarkers for the risk prediction and early diagnosis of non-genetically heritable stress-related diseases. The identification of predictive and diagnostic metabolic biomarkers that identify adverse health trajectories linked to maternal stress is critical for precision medicine approaches, enabling the identification of the most vulnerable individuals who are at the highest risk. Because the orchestrated regulation of metabolic pathways is essential for stress regulation, development, and successful aging, its failure may lead to cell and tissue dysfunction and result in chronic and non-communicable diseases. On the other hand, metabolic profiles are often cell-type specific [78], thus providing an opportunity for the discovery of causal pathogenic mechanisms and early diagnosis of a disease. Due to the potentially significant health impact of social isolation, understanding the metabolic consequences of ancestral stress may provide new personalized therapeutic strategies to mitigate stress-associated adverse health outcomes.

4. Materials and Methods

4.1. Experimental Design

4.1.1. Animal Model

The study involved 24 adult male Long-Evans rats obtained from the F3 generation of three lineages: yoked controls (CONT, n = 8), transgenerational prenatal stress (TPS; n = 8), and multigenerational prenatal stress (MPS; n = 8). Four samples from the control group did not yield enough plasma for metabolomics analysis, resulting in a reduced control group size (CONT, n = 4). TPS rats were the third filial (F3) generation of a lineage in which only the first filial (F1) generation was stressed prenatally; subsequent generations remained unstressed (Figure 5). MPS rats were the F3 generation of a lineage in which each consecutive generation (F1, F2 and F3) was prenatally stressed.

4.1.2. Gestational Stress

Pregnant dams were stressed by social isolation, which has been shown to result in mild psychosocial stress in rats [79] and profound epigenetic and behavioural phenotypes [69,80]. Rats were housed in pairs from weaning until postnatal (P) day 90. Then, dams were separated and housed alone, one per cage, throughout pregnancy until weaning of their offspring [69,80]. Pairing for mating began at P90 and the maximum time spent in social isolation prior to conception was 30 days. Control rats remained housed in pairs, except for the period from gestational day 21 to lactational day 21 (time of weaning), to allow for undisturbed rearing of their offspring.

4.1.3. Breeding Colony

Animals were bred in-house for at least five generations prior to the beginning of the experiment. Each generation F0-F4 was outcrossed to avoid inbreeding. Distinct lineages were monitored through the JAX Colony Management System (JCMS; Jackson Laboratory, Bar Harbour, ME, USA). One male offspring per litter was randomly selected for this metabolomics study. Bystander effects of stress were avoided by using designated testing and housing spaces. All housing, handling, testing, and tissue sampling conditions were harmonized and carefully controlled across generations.
The rats were housed in polycarbonate shoebox cages on corn cob bedding under a 12 h light/dark cycle with lights on at 7:30 AM. The room temperature was maintained at 20 °C with relative humidity at 30%. All procedures were performed in accordance with the guidelines of the Canadian Council on Animal Care and approved by the University of Lethbridge Animal Welfare Committee.

4.2. Behavioural Testing

Anxiety-like behaviours towards a novel environment and risk taking were assessed in an open field test at P90. Rats were individually placed into an Accuscan activity monitoring box (Accuscan Instruments Inc., Columbus, OH, USA; 42 × 42 × 30 cm) and recorded for 10 min. VersaMax™ software (Accuscan Instruments Inc., OH, USA) monitored the animals’ horizontal and vertical activity to be analyzed using VersaDat™ software (Accuscan Instruments Inc., Columbus, OH, USA; [29]). Vertical (rearing activity) and horizontal (distance travelled) exploration were considered.

4.3. Sample Collection and Preparation

At the age of 140 days (P140), 600 μL of blood was collected from the lateral tail vein using a 23-gauge butterfly needle coated with lithium heparin, while rats were anaesthetized using 4% isoflurane. All blood samples were collected at the same time of day, between 8:00 and 9:00 AM in the morning. EDTA and citrate collection tubes were avoided, as they give additional spectral signals in 1H NMR spectroscopy [81]. Blood was transferred to centrifuge tubes and plasma was obtained by centrifugation at 1600× g for 15 min at 4 °C. Samples were stored at −80 °C until further processing.
Water soluble metabolites were extracted from plasma using a methanol precipitation protocol [82]. Methanol was added in a 2:1 ratio (550 μL buffered plasma + 1.10 mL methanol) in 2.0 mL Eppendorf tubes. These were vortexed for 5 s, incubated at −20 °C for 20 min, and centrifuged at 12,000× g for 30 min at 4 °C. Supernatant was decanted to fresh tubes and allowed to dry under a gentle stream of nitrogen gas until all liquid had evaporated. Dried pellets were then re-suspended in 600 μL of phosphate buffer calibrated to pH 7.4 and vortexed for 10 s, or until the pellets were completely dissolved. The phosphate buffer was prepared as a 4:1 ratio of KH2PO4:K2HPO4 in a 4:1 H2O:D2O solution to obtain a final concentration of 0.5 M. The D2O contained 0.03% w/v trimethylsilyl propanoic acid (TSP) as a chemical shift reference for 1H NMR spectroscopy. To maintain metabolite integrity, 0.03% w/v of sodium azide (NaN3) was also added to the buffer as an antimicrobial agent. Finally, samples were centrifuged again at 12,000× g for 5 min at 4 °C and 550 μL of re-suspended supernatant was added to 5 mm NMR tubes for data acquisition.

4.4. NMR Data Acquisition and Processing

NMR spectra were collected on a 700 MHz Bruker Avance III HD spectrometer (Bruker, Milton, ON, Canada; Figure S1). The 1-D NOESY gradient water suppression pulse sequence ‘noesygpr1d’ was used. Each sample was run for 512 scans to a total acquisition size of 256 k. The spectra were zero filled to 512 k, automatically phased and baseline corrected, and line-broadened by 0.3 Hz. The processed spectra were then exported to MATLAB (The MathWorks, Natick, MA, USA) for statistical analysis. Spectra were binned using Dynamic Adaptive Binning [83]. Each spectrum was normalized to the total unit area of all spectral bins and the spectral region corresponding to water was removed from data normalization. The data set was then Pareto scaled to reduce the influence of intense peaks while emphasizing weaker ones and log transformed [84]. All peaks were referenced to TSP (0.00δ).

4.5. Statistical Analysis

The overall structure of scaled spectra was visualized and compared across control and experimental groups using principal component analysis (PCA). A Shapiro–Wilk test for data normality was applied and all data were found to be non-parametric. A Mann–Whitney U test was subsequently used to determine which spectral bins were significantly altered in comparison groups [85]. Bonferroni–Holm correction was applied to all univariate statistical tests to correct for multiple comparisons. Bins with p-values less than or equal to 0.05 were deemed significant and considered as potential markers defining substantial alterations in metabolite concentration between groups.
PCA and Mann–Whitney U tests were carried out using the web-based Metaboanalyst software [86] and MATLAB® (MathWorks, Natick, MA, USA), respectively. The variable importance analysis based on random variable combination (VIAVC) algorithm [87] was utilized to determine which set of metabolite bins led to the best discrimination between the comparison groups (referred to as the best subset). This machine learning, multivariate statistical algorithm makes use of random data combinations and iteratively tests each combination using both PLS-DA and the area under the curve of the receiver operator characteristic curve (ROC) to determine the best subset of metabolites leading to class separation [88]. PCA score plots (Figure 1) were generated using only those bins identified as significant by Mann–Whitey U and VIAVC testing (Table 1). Pathway topology analysis and metabolite set enrichment analysis were carried out using Metaboanalyst [89,90]. The Rattus Norvegicus pathway library was utilized for both analyses, the hypergeometric test was selected for the over-representation analysis, and relative-betweenness centrality was chosen for the pathway topology analysis.
Spearman’s rank correlations were carried out using SPSS 26 for Windows (IBM Corporation, Armonk, NY, USA) to determine the relationship between normalized metabolite concentrations and behavioural data. This analysis concerned only metabolites determined to be significantly altered between groups, based on the multivariate and univariate statistical tests outlined above. Samples with p-values equal to or less than 0.05 were considered significant.

4.6. Metabolite Identification

A spectral database of pure metabolites was generated and used to identify most metabolites found in the NMR spectra. The Human Metabolome Database [91,92,93] was used to supplement and aid the creation of our spectral library, and to identify some substances that were not obtained for the creation of this database. Furthermore, tools of the Chenomx 8.2 NMR Suite (Chenomx Inc., Edmonton, AB, Canada) were used for identifying and quantifying NMR metabolites which allowed for the spectral deconvolution of biofluid samples into individual components.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo12070572/s1, Figure S1: Labelled 1H NMR spectra of rat plasma metabolome; Table S1: Open field task scores.

Author Contributions

Conceptualization, J.K.M. and G.A.S.M.; methodology, J.P.H., E.J.P., J.K.M., P.S. and T.M.; software, T.M; formal analysis, J.P.H., E.J.P., P.S., J.K.M. and T.M.; writing—original draft preparation, J.P.H., E.J.P., P.S., T.M. and G.A.S.M.; writing—review and editing, J.P.H., T.M. and G.A.S.M.; supervision, T.M. and G.A.S.M.; funding acquisition, T.M. and G.A.S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Sciences and Engineering Research Council of Canada (NSERC), grant #363195 (G.A.S.M.), and the University of Lethbridge Health Research Accelerator Fund (HRAF; T.M. and G.A.S.M.). E.J.P. was supported by a CIHR CGS-M award.

Institutional Review Board Statement

The study was conducted in accordance with the recommendations specified in the Canadian Council on Animal Care Guidelines and approved by the University of Lethbridge Animal Welfare Committee (protocol # 1715, approved 21 November 2017).

Data Availability Statement

The data presented in this study are available upon request from the corresponding authors, as it has not been uploaded to an online database.

Acknowledgments

The authors would like to thank Michael Opyr for his assistance with both the data processing scripts and the spectrometer. We also thank the University of Lethbridge for use of the Magnetic Resonance Facility, which provided access to and training on the 700 MHz NMR spectrometer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saklayen, M.G. The global epidemic of the metabolic syndrome. Curr. Hypertens. Rep. 2018, 20, 12. [Google Scholar] [CrossRef] [Green Version]
  2. Casas, M.; Project, O.B.O.T.I.; Forns, J.; Martínez, D.; Guxens, M.; Fernandez-Somoano, A.; Ibarluzea, J.; Lertxundi, N.; Murcia, M.; Rebagliato, M.; et al. Maternal pre-pregnancy obesity and neuropsychological development in pre-school children: A prospective cohort study. Pediatr. Res. 2017, 82, 596–606. [Google Scholar] [CrossRef] [Green Version]
  3. Mina, T.H.; Lahti, M.; Drake, A.J.; Forbes, S.; Denison, F.C.; Räikkönen, K.; Norman, J.E.; Reynolds, R.M. Maternal lipids in pregnancy are associated with increased offspring cortisol reactivity in childhood. Psychoneuroendocrinology 2017, 83, 79–83. [Google Scholar] [CrossRef]
  4. James, S.L.; Abate, D.; Abate, K.H.; Abay, S.M.; Abbafati, C.; Abbasi, N.; Abbastabar, H.; Abd-Allah, F.; Abdela, J.; Abdelalim, A.; et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018, 392, 1789–1858. [Google Scholar] [CrossRef] [Green Version]
  5. Contu, L.; Hawkes, C.A. A Review of the impact of maternal obesity on the cognitive function and mental health of the offspring. Int. J. Mol. Sci. 2017, 18, 1093. [Google Scholar] [CrossRef] [Green Version]
  6. Godfrey, K.M.; Reynolds, R.; Prescott, S.L.; Nyirenda, M.; Jaddoe, V.W.V.; Eriksson, J.G.; Broekman, B.F.P. Influence of maternal obesity on the long-term health of offspring. Lancet Diabetes Endocrinol. 2016, 5, 53–64. [Google Scholar] [CrossRef] [Green Version]
  7. Charil, A.; Laplante, D.P.; Vaillancourt, C.; King, S. Prenatal stress and brain development. Brain Res. Rev. 2010, 65, 56–79. [Google Scholar] [CrossRef]
  8. King, S.; Laplante, D.P. The effects of prenatal maternal stress on children’s cognitive development: Project Ice Storm. Stress 2005, 8, 35–45. [Google Scholar] [CrossRef] [PubMed]
  9. Ping, E.Y.; Laplante, D.P.; Elgbeili, G.; Hillerer, K.M.; Brunet, A.; O’Hara, M.W.; King, S. Prenatal maternal stress predicts stress reactivity at 2½ years of age: The Iowa Flood Study. Psychoneuroendocrinology 2015, 56, 62–78. [Google Scholar] [CrossRef] [PubMed]
  10. Harris, A.; Seckl, J. Glucocorticoids, prenatal stress and the programming of disease. Horm. Behav. 2011, 59, 279–289. [Google Scholar] [CrossRef]
  11. Dancause, K.N.; Vilar, M.; Wilson, M.; Soloway, L.E.; Dehuff, C.; Chan, C.; Tarivonda, L.; Regenvanu, R.; Kaneko, A.; Lum, J.K.; et al. Behavioral risk factors for obesity during health transition in Vanuatu, South Pacific. Obesity 2012, 21, E98–E104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Reynolds, R.M. Corticosteroid-mediated programming and the pathogenesis of obesity and diabetes. J. Steroid Biochem. Mol. Biol. 2010, 122, 3–9. [Google Scholar] [CrossRef] [PubMed]
  13. Van den Bergh, B.R.H.; van den Heuvel, M.I.; Lahti-Pulkkinen, M.; Braeken, M.; de Rooij, S.R.; Entringer, S.; Hoyer, D.; Roseboom, T.; Räikkönen, K.; King, S.; et al. Prenatal developmental origins of behavior and mental health: The influence of maternal stress in pregnancy. Neurosci. Biobehav. Rev. 2020, 117, 26–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Maric, N.P.; Dunjic, B.; Stojiljkovic, D.J.; Britvic, D.; Jasovic-Gasic, M. Prenatal stress during the 1999 bombing associated with lower birth weight—a study of 3,815 births from Belgrade. Arch. Women’s Ment. Health 2009, 13, 83–89. [Google Scholar] [CrossRef]
  15. Lederman, S.A.; Rauh, V.; Weiss, L.; Stein, J.L.; Hoepner, L.; Becker, M.; Perera, F.P. The effects of the World Trade Center event on birth outcomes among term deliveries at three Lower Manhattan hospitals. Environ. Health Perspect. 2004, 112, 1772–1778. [Google Scholar] [CrossRef] [Green Version]
  16. Li, N.; Wang, W.; Wu, G.; Wang, J. Nutritional support for low birth weight infants: Insights from animal studies. Br. J. Nutr. 2017, 117, 1390–1402. [Google Scholar] [CrossRef] [Green Version]
  17. Reid, K.M.; Taylor, M.G. Social support, stress, and maternal postpartum depression: A comparison of supportive relationships. Soc. Sci. Res. 2015, 54, 246–262. [Google Scholar] [CrossRef]
  18. Lebel, C.; MacKinnon, A.; Bagshawe, M.; Tomfohr-Madsen, L.; Giesbrecht, G. Elevated depression and anxiety symptoms among pregnant individuals during the COVID-19 pandemic. J. Affect. Disord. 2020, 277, 5–13. [Google Scholar] [CrossRef]
  19. Lei, C.-L.; Massart, R.; Suderman, M.J.; Machnes, Z.; Elgbeili, G.; Laplante, D.P.; Szyf, M.; King, S. DNA methylation signatures triggered by prenatal maternal stress exposure to a natural disaster: Project Ice Storm. PLoS ONE 2014, 9, e107653. [Google Scholar] [CrossRef] [Green Version]
  20. Zucchi, F.C.R.; Yao, Y.; Ward, I.D.; Ilnytskyy, Y.; Olson, D.M.; Benzies, K.; Kovalchuk, I.; Kovalchuk, O.; Metz, G.A.S. Maternal stress induces epigenetic signatures of psychiatric and neurological diseases in the offspring. PLoS ONE 2013, 8, e56967. [Google Scholar] [CrossRef]
  21. Hollins, S.; Cairns, M.J. MicroRNA: Small RNA mediators of the brains genomic response to environmental stress. Prog. Neurobiol. 2016, 143, 61–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Meaney, M.J.; Szyf, M. Environmental programming of stress responses through DNA methylation: Life at the interface between a dynamic environment and a fixed genome. Dialog-Clin. Neurosci. 2005, 7, 103–123. [Google Scholar] [CrossRef]
  23. Babenko, O.; Kovalchuk, I.; Metz, G.A.S. Stress-induced perinatal and transgenerational epigenetic programming of brain development and mental health. Neurosci. Biobehav. Rev. 2015, 48, 70–91. [Google Scholar] [CrossRef] [PubMed]
  24. Yao, Y.; Robinson, A.M.; Zucchi, F.C.R.; Robbins, J.C.; Babenko, O.; Kovalchuk, O.; Kovalchuk, I.; Olson, D.M.; Metz, G.A.S. Ancestral exposure to stress epigenetically programs preterm birth risk and adverse maternal and newborn outcomes. BMC Med. 2014, 12, 121. [Google Scholar] [CrossRef] [Green Version]
  25. Aerts, L.; Van Assche, F. Animal evidence for the transgenerational development of diabetes mellitus. Int. J. Biochem. Cell Biol. 2006, 38, 894–903. [Google Scholar] [CrossRef] [PubMed]
  26. Thamotharan, M.; Garg, M.; Oak, S.; Rogers, L.M.; Pan, G.; Sangiorgi, F.; Lee, P.W.N.; Devaskar, S.U. Transgenerational inheritance of the insulin-resistant phenotype in embryo-transferred intrauterine growth-restricted adult female rat offspring. Am. J. Physiol. Metab. 2007, 292, E1270–E1279. [Google Scholar] [CrossRef] [Green Version]
  27. Skinner, M.K.; Anway, M.D.; Savenkova, M.I.; Gore, A.C.; Crews, D. Transgenerational epigenetic programming of the brain transcriptome and anxiety behavior. PLoS ONE 2008, 3, e3745. [Google Scholar] [CrossRef]
  28. Gapp, K.; Jawaid, A.; Sarkies, P.; Bohacek, J.; Pelczar, P.; Prados, J.; Farinelli, L.; Miska, E.; Mansuy, I.M. Implication of sperm RNAs in transgenerational inheritance of the effects of early trauma in mice. Nat. Neurosci. 2014, 17, 667–669. [Google Scholar] [CrossRef] [Green Version]
  29. Kiss, D.; Ambeskovic, M.; Montina, T.; Metz, G.A.S. Stress transgenerationally programs metabolic pathways linked to altered mental health. Cell. Mol. Life Sci. 2016, 73, 4547–4557. [Google Scholar] [CrossRef]
  30. Ambeskovic, M.; Ilnytskyy, Y.; Kiss, D.; Currie, C.; Montina, T.; Kovalchuk, I.; Metz, G.A. Ancestral stress programs sex-specific biological aging trajectories and non-communicable disease risk. Aging 2020, 12, 3828–3847. [Google Scholar] [CrossRef]
  31. Ambeskovic, M.; Babenko, O.; Ilnytskyy, Y.; Kovalchuk, I.; Kolb, B.; Metz, G.A.S. Ancestral stress alters lifetime mental health trajectories and cortical neuromorphology via epigenetic regulation. Sci. Rep. 2019, 9, 6389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ambeskovic, M.; Roseboom, T.J.; Metz, G.A. Transgenerational effects of early environmental insults on aging and disease incidence. Neurosci. Biobehav. Rev. 2020, 117, 297–316. [Google Scholar] [CrossRef] [PubMed]
  33. Feinberg, A.P. Phenotypic plasticity and the epigenetics of human disease. Nature 2007, 447, 433–440. [Google Scholar] [CrossRef] [PubMed]
  34. Van Steenwyk, G.; Roszkowski, M.; Manuella, F.; Franklin, T.B.; Mansuy, I.M. Transgenerational inheritance of behavioral and metabolic effects of paternal exposure to traumatic stress in early postnatal life: Evidence in the 4th generation. Environ. Epigenetics 2018, 4, dvy023. [Google Scholar] [CrossRef]
  35. McCreary, J.K.; Erickson, Z.T.; Paxman, E.; Kiss, D.; Montina, T.; Olson, D.M.; Metz, G.A.S. The rat cumulative allostatic load measure (rCALM): A new translational assessment of the burden of stress. Environ. Epigenetics 2019, 5, dvz005. [Google Scholar] [CrossRef]
  36. Paxman, E.J.; Boora, N.S.; Kiss, D.; Laplante, D.P.; King, S.; Montina, T.; Metz, G.A.S. Prenatal maternal stress from a natural disaster alters urinary metabolomic profiles in Project Ice Storm participants. Sci. Rep. 2018, 8, 12932. [Google Scholar] [CrossRef]
  37. Esmailnasab, N.; Moradi, G.; Delaveri, A. Risk factors of non-communicable diseases and metabolic syndrome. Iran. J. Public Health 2012, 41, 77–85. [Google Scholar]
  38. Hjemdahl, P. Stress and the Metabolic Syndrome: An interesting but enigmatic association. Circulation 2002, 106, 2634–2636. [Google Scholar] [CrossRef] [Green Version]
  39. Poplawski, J.; Radmilovic, A.; Montina, T.D.; Metz, G.A.S. Cardiorenal metabolic biomarkers link early life stress to risk of non-communicable diseases and adverse mental health outcomes. Sci. Rep. 2020, 10, 13295. [Google Scholar] [CrossRef]
  40. Bouatra, S.; Aziat, F.; Mandal, R.; Guo, A.C.; Wilson, M.R.; Knox, C.; Bjorndahl, T.C.; Krishnamurthy, R.; Saleem, F.; Liu, P.; et al. The human urine metabolome. PLoS ONE 2013, 8, e73076. [Google Scholar] [CrossRef] [Green Version]
  41. van Steenwyk, G.; Gapp, K.; Jawaid, A.; Germain, P.; Manuella, F.; Tanwar, D.K.; Zamboni, N.; Gaur, N.; Efimova, A.; Thumfart, K.M.; et al. Involvement of circulating factors in the transmission of paternal experiences through the germline. EMBO J. 2020, 39, e104579. [Google Scholar] [CrossRef] [PubMed]
  42. McCreary, J.K.; Truica, L.S.; Friesen, B.; Yao, Y.; Olson, D.M.; Kovalchuk, I.; Cross, A.R.; Metz, G.A. Altered brain morphology and functional connectivity reflect a vulnerable affective state after cumulative multigenerational stress in rats. Neuroscience 2016, 330, 79–89. [Google Scholar] [CrossRef] [PubMed]
  43. Skelin, I.; Needham, M.; Molina, L.; Metz, G.; Gruber, A. Multigenerational prenatal stress increases the coherence of brain signaling among cortico–striatal–limbic circuits in adult rats. Neuroscience 2015, 289, 270–278. [Google Scholar] [CrossRef] [PubMed]
  44. Seckl, J.R. Glucocorticoid “programming” and PTSD risk. Ann. N. Y. Acad. Sci. 2006, 1071, 351–378. [Google Scholar] [CrossRef]
  45. Meaney, M.J.; Szyf, M.; Seckl, J.R. Epigenetic mechanisms of perinatal programming of hypothalamic-pituitary-adrenal function and health. Trends Mol. Med. 2007, 13, 269–277. [Google Scholar] [CrossRef]
  46. Baker, J.S.; McCormick, M.C.; Robergs, R.A. Interaction among skeletal muscle metabolic energy systems during intense exercise. J. Nutr. Metab. 2010, 2010, 905612. [Google Scholar] [CrossRef] [Green Version]
  47. Allen, P.J. Creatine metabolism and psychiatric disorders: Does creatine supplementation have therapeutic value? Neurosci. Biobehav. Rev. 2012, 36, 1442–1462. [Google Scholar] [CrossRef] [Green Version]
  48. Oizel, K.; Tait-Mulder, J.; Fernandez-De-Cossio-Diaz, J.; Pietzke, M.; Brunton, H.; Lilla, S.; Dhayade, S.; Athineos, D.; Blanco, G.R.; Sumpton, D.; et al. Formate induces a metabolic switch in nucleotide and energy metabolism. Cell Death Dis. 2020, 11, 310. [Google Scholar] [CrossRef]
  49. Raaijmakers, A.; Zhang, Z.-Y.; Claessens, J.; Cauwenberghs, N.; van Tienoven, T.P.; Wei, F.-F.; Jacobs, L.; Levtchenko, E.; Pauwels, S.; Kuznetsova, T.; et al. Does extremely low birth weight predispose to low-renin hypertension? Hypertension 2017, 69, 443–449. [Google Scholar] [CrossRef]
  50. Cherak, S.J.; Giesbrecht, G.F.; Metcalfe, A.; Ronksley, P.E.; Malebranche, M.E. The effect of gestational period on the association between maternal prenatal salivary cortisol and birth weight: A systematic review and meta-analysis. Psychoneuroendocrinology 2018, 94, 49–62. [Google Scholar] [CrossRef]
  51. Gunther, N.; Drukker, M.; Feron, F.; Van Os, J. Association of mental health problems in childhood with prenatal and postnatal physical growth. Eur. Psychiatry 2005, 20, 277–286. [Google Scholar] [CrossRef] [PubMed]
  52. Gale, C.R.; Martyn, C.N. Birth weight and later risk of depression in a national birth cohort. Br. J. Psychiatry 2004, 184, 28–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Breslau, N. Psychiatric sequelae of low birth weight. Epidemiol. Rev. 1995, 17, 96–106. [Google Scholar] [CrossRef] [PubMed]
  54. Alexandre-Gouabau, M.-C.; Moyon, T.; David-Sochard, A.; Fenaille, F.; Cholet, S.; Royer, A.-L.; Guitton, Y.; Billard, H.; Darmaun, D.; Rozé, J.-C.; et al. Comprehensive preterm breast milk metabotype associated with optimal infant early growth [attern. Nutrients 2019, 11, 528. [Google Scholar] [CrossRef] [Green Version]
  55. Li, H.; Ye, D.; Xie, W.; Hua, F.; Yang, Y.; Wu, J.; Gu, A.; Ren, Y.; Mao, K. Defect of branched-chain amino acid metabolism promotes the development of Alzheimer’s disease by targeting the mTOR signaling. Biosci. Rep. 2018, 38, BSR20180127. [Google Scholar] [CrossRef] [Green Version]
  56. Fontana, L.; Cummings, N.E.; Apelo, S.I.A.; Neuman, J.C.; Kasza, I.; Schmidt, B.; Cava, E.; Spelta, F.; Tosti, V.; Syed, F.A.; et al. Decreased consumption of branched-chain amino acids improves metabolic health. Cell Rep. 2016, 16, 520–530. [Google Scholar] [CrossRef] [Green Version]
  57. Rosa, L.; Scaini, G.; Furlanetto, C.B.; Galant, L.S.; Vuolo, F.; Dall’Igna, D.M.; Schuck, P.F.; Ferreira, G.C.; Dal-Pizzol, F.; Streck, E.L. Administration of branched-chain amino acids alters the balance between pro-inflammatory and anti-inflammatory cytokines. Int. J. Dev. Neurosci. 2015, 48, 24–30. [Google Scholar] [CrossRef]
  58. Robergs, R.A.; Ghiasvand, F.; Parker, D. Biochemistry of exercise-induced metabolic acidosis. Am. J. Physiol. Integr. Comp. Physiol. 2004, 287, R502–R516. [Google Scholar] [CrossRef]
  59. Zheng, S.; Yu, M.; Lu, X.; Huo, T.; Ge, L.; Yang, J.; Wu, C.; Li, F. Urinary metabonomic study on biochemical changes in chronic unpredictable mild stress model of depression. Clin. Chim. Acta 2010, 411, 204–209. [Google Scholar] [CrossRef]
  60. Garcia-Alvarez, M.; Marik, P.; Bellomo, R. Stress hyperlactataemia: Present understanding and controversy. Lancet Diabetes Endocrinol. 2014, 2, 339–347. [Google Scholar] [CrossRef]
  61. Kushimoto, S.; Akaishi, S.; Sato, T.; Nomura, R.; Fujita, M.; Kudo, D.; Kawazoe, Y.; Yoshida, Y.; Miyagawa, N. Lactate, a useful marker for disease mortality and severity but an unreliable marker of tissue hypoxia/hypoperfusion in critically ill patients. Acute Med. Surg. 2016, 3, 293–297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Fischer, L.M.; da Costa, K.-A.; Kwock, L.; Galanko, J.; Zeisel, S.H. Dietary choline requirements of women: Effects of estrogen and genetic variation. Am. J. Clin. Nutr. 2010, 92, 1113–1119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Leermakers, E.T.; Moreira, E.M.; Jong, J.C.K.-D.; Darweesh, S.K.; Visser, T.; Voortman, T.; Bautista, P.K.; Chowdhury, R.; Gorman, D.; Bramer, W.M.; et al. Effects of choline on health across the life course: A systematic review. Nutr. Rev. 2015, 73, 500–522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Tabatabaie, L.; Klomp, L.; Berger, R.; de Koning, T. l-Serine synthesis in the central nervous system: A review on serine deficiency disorders. Mol. Genet. Metab. 2010, 99, 256–262. [Google Scholar] [CrossRef]
  65. De Koning, T.J.; Snell, K.; Duran, M.; Berger, R.; Poll-The, B.-T.; Surtees, R. l-Serine in disease and development. Biochem. J. 2003, 371, 653–661. [Google Scholar] [CrossRef]
  66. Nederhof, E.; Ormel, J.; Oldehinkel, A.J. Mismatch or Cumulative Stress: The pathway to depression is conditional on attention style. Psychol. Sci. 2014, 25, 684–692. [Google Scholar] [CrossRef]
  67. Faraji, J.; Soltanpour, N.; Ambeskovic, M.; Zucchi, F.C.R.; Beaumier, P.; Kovalchuk, I.; Metz, G.A.S. Evidence for ancestral programming of resilience in a two-hit stress model. Front. Behav. Neurosci. 2017, 11, 89. [Google Scholar] [CrossRef] [Green Version]
  68. Costall, B.; Jones, B.; Kelly, M.; Naylor, R.; Tomkins, D. Exploration of mice in a black and white test box: Validation as a model of anxiety. Pharmacol. Biochem. Behav. 1989, 32, 777–785. [Google Scholar] [CrossRef]
  69. McCreary, J.K.; Erickson, Z.T.; Hao, Y.; Ilnytskyy, Y.; Kovalchuk, I.; Metz, G.A.S. Environmental intervention as a therapy for adverse programming by ancestral stress. Sci. Rep. 2016, 6, 37814. [Google Scholar] [CrossRef] [Green Version]
  70. Faraji, J.; Soltanpour, N.; Jafari, S.Y.; Moeeini, R.; Pakdel, S.; Moharreri, A.; Metz, G.A. Stress inhibits psychomotor performance differently in simple and complex open field environments. Horm. Behav. 2014, 65, 66–75. [Google Scholar] [CrossRef]
  71. Ennaceur, A. Tests of unconditioned anxiety—Pitfalls and disappointments. Physiol. Behav. 2014, 135, 55–71. [Google Scholar] [CrossRef] [PubMed]
  72. Wang, X.; Liu, T.; Song, H.; Cui, S.; Liu, G.; Christoforou, A.; Flaherty, P.; Luo, X.; Wood, L.; Wang, Q.M. Targeted metabolomic profiling reveals association between altered amino acids and poor functional recovery after stroke. Front. Neurol. 2019, 10, 1425. [Google Scholar] [CrossRef] [PubMed]
  73. Tofte, N.; Suvitaival, T.; Trost, K.; Mattila, I.M.; Theilade, S.; Winther, S.A.; Ahluwalia, T.S.; Frimodt-Møller, M.; Legido-Quigley, C.; Rossing, P. Metabolomic assessment reveals alteration in polyols and branched chain amino acids associated with present and future renal impairment in a discovery cohort of 637 persons with type 1 diabetes. Front. Endocrinol. 2019, 10, 818. [Google Scholar] [CrossRef] [PubMed]
  74. Gonzalez-Franquesa, A.; Burkart, A.M.; Isganaitis, E.; Patti, M.-E. What have metabolomics approaches taught us about type 2 diabetes? Curr. Diabetes Rep. 2016, 16, 74. [Google Scholar] [CrossRef] [PubMed]
  75. Lee, J.; Choi, J.-Y.; Kwon, Y.-K.; Lee, D.; Jung, H.-Y.; Ryu, H.-M.; Cho, J.-H.; Ryu, D.H.; Kim, Y.-L.; Hwang, G.-S. Changes in serum metabolites with the stage of chronic kidney disease: Comparison of diabetes and non-diabetes. Clin. Chim. Acta 2016, 459, 123–131. [Google Scholar] [CrossRef]
  76. Skinner, M.K. What is an epigenetic transgenerational phenotype?: F3 or F2. Reprod. Toxicol. 2008, 25, 2–6. [Google Scholar] [CrossRef] [Green Version]
  77. Zucchi, F.C.R.; Yao, Y.; Metz, G.A. The secret language of destiny: Stress imprinting and transgenerational origins of disease. Front. Genet. 2012, 3, 96. [Google Scholar] [CrossRef] [Green Version]
  78. Cable, J.; Pourquié, O.; Wellen, K.E.; Finley, L.W.S.; Aulehla, A.; Gould, A.P.; Teleman, A.; Tu, W.B.; Garrett, W.S.; Miguel-Aliaga, I.; et al. Metabolic decisions in development and disease—a Keystone Symposia report. Ann. N. Y. Acad. Sci. 2021, 1506, 55–73. [Google Scholar] [CrossRef]
  79. Hawkley, L.C.; Gu, Y.; Luo, Y.-J.; Cacioppo, J.T. The mental representation of social connections: Generalizability extended to Beijing adults. PLoS ONE 2012, 7, e44065. [Google Scholar] [CrossRef]
  80. McCreary, J.K.; Erickson, Z.T.; Metz, G.A. Environmental enrichment mitigates the impact of ancestral stress on motor skill and corticospinal tract plasticity. Neurosci. Lett. 2016, 632, 181–186. [Google Scholar] [CrossRef]
  81. Barton, R.H.; Waterman, D.; Bonner, F.W.; Holmes, E.; Clarke, R.; Procardis, C.; Nicholson, J.K.; Lindon, J.C. The influence of EDTA and citrate anticoagulant addition to human plasma on information recovery from NMR-based metabolic profiling studies. Mol. BioSyst. 2009, 6, 215–224. [Google Scholar] [CrossRef] [PubMed]
  82. Gowda, G.A.N.; Raftery, D. Quantitating metabolites in protein precipitated serum using NMR spectroscopy. Anal. Chem. 2014, 86, 5433–5440. [Google Scholar] [CrossRef] [PubMed]
  83. Anderson, P.E.; Mahle, D.A.; Doom, T.E.; Reo, N.V.; DelRaso, N.J.; Raymer, M.L. Dynamic adaptive binning: An improved quantification technique for NMR spectroscopic data. Metabolomics 2010, 7, 179–190. [Google Scholar] [CrossRef]
  84. Craig, A.; Cloarec, O.; Holmes, E.; Nicholson, A.J.K.; Lindon, J.C. Scaling and normalization effects in NMR spectroscopic metabonomic data sets. Anal. Chem. 2006, 78, 2262–2267. [Google Scholar] [CrossRef] [PubMed]
  85. Goodpaster, A.M.; Romick-Rosendale, L.; Kennedy, M.A. Statistical significance analysis of nuclear magnetic resonance-based metabonomics data. Anal. Biochem. 2010, 401, 134–143. [Google Scholar] [CrossRef]
  86. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [Green Version]
  87. Yun, Y.-H.; Liang, F.; Deng, B.-C.; Lai, G.-B.; Gonçalves, C.M.V.; Lu, H.-M.; Yan, J.; Huang, X.; Yi, L.-Z.; Liang, Y.-Z. Informative metabolites identification by variable importance analysis based on random variable combination. Metabolomics 2015, 11, 1539–1551. [Google Scholar] [CrossRef]
  88. Fawcett, T. An Introduction to ROC analysis. Pattern Recogn. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  89. Xia, J.; Wishart, D.S. MetPA: A web-based metabolomics tool for pathway analysis and visualization. Bioinformatics 2010, 26, 2342–2344. [Google Scholar] [CrossRef] [Green Version]
  90. Xia, J.; Wishart, D.S. MSEA: A web-based tool to identify biologically meaningful patterns in quantitative metabolomic data. Nucleic Acids Res. 2010, 38, W71–W77. [Google Scholar] [CrossRef] [Green Version]
  91. Wishart, D.S.; Knox, C.; Guo, A.C.; Eisner, R.; Young, N.; Gautam, B.; Hau, D.D.; Psychogios, N.; Dong, E.; Bouatra, S.; et al. HMDB: A knowledgebase for the human metabolome. Nucleic Acids Res. 2008, 37, D603–D610. [Google Scholar] [CrossRef] [PubMed]
  92. Wishart, D.S.; Jewison, T.; Guo, A.C.; Wilson, M.; Knox, C.; Liu, Y.; Djoumbou, Y.; Mandal, R.; Aziat, F.; Dong, E.; et al. HMDB 3.0—The Human Metabolome Database in 2013. Nucleic Acids Res. 2013, 41, D801–D807. [Google Scholar] [CrossRef] [PubMed]
  93. Wishart, D.S.; Tzur, D.; Knox, C.; Eisner, R.; Guo, A.C.; Young, N.; Cheng, D.; Jewell, K.; Arndt, D.; Sawhney, S.; et al. HMDB: The Human Metabolome Database. Nucleic Acids Res. 2007, 35, D521–D526. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Principle component analysis (PCA) plots showing unsupervised separation between the groups: (A) transgenerationally stressed rats (TPS) vs. non-stressed controls (CONT); (B) multigenerationally stressed rats (MPS) vs. non-stressed controls (CONT); and (C) transgenerationally stressed rats (TPS) vs. multigenerationally stressed rats (MPS). Each PCA was carried out using only the bins that were determined to be significantly altered. Each point (triangle, cross, circle, or square) represents one individual based on the list of blood plasma metabolites found to be statistically significant via a Mann–Whitney U test and VIAVC analysis. Coloured ellipses represent 95% confidence intervals. x and y axes show principal components 1 and 2, respectively, with brackets indicating the variance explained by each principal component. Note that the metabolomes of both TPS and MPS differed substantially from CONT and from each other.
Figure 1. Principle component analysis (PCA) plots showing unsupervised separation between the groups: (A) transgenerationally stressed rats (TPS) vs. non-stressed controls (CONT); (B) multigenerationally stressed rats (MPS) vs. non-stressed controls (CONT); and (C) transgenerationally stressed rats (TPS) vs. multigenerationally stressed rats (MPS). Each PCA was carried out using only the bins that were determined to be significantly altered. Each point (triangle, cross, circle, or square) represents one individual based on the list of blood plasma metabolites found to be statistically significant via a Mann–Whitney U test and VIAVC analysis. Coloured ellipses represent 95% confidence intervals. x and y axes show principal components 1 and 2, respectively, with brackets indicating the variance explained by each principal component. Note that the metabolomes of both TPS and MPS differed substantially from CONT and from each other.
Metabolites 12 00572 g001
Figure 2. Pathway-associated metabolite set enrichment analysis (MSEA) for TPS vs. CONT. (A) Plot showing affected biological processes in transgenerationally stressed (TPS) compared to non-stressed (CONT) rats based on metabolites identified as significantly altered between these groups. The p-value for each pathway is shown using the heatmap on the right of the figure and the fold enrichment shows how many times greater than chance the process is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in TPS and non-stressed CONT rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis gives the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. TPS had the highest impact on valine, leucin, and isoleucine biosynthesis.
Figure 2. Pathway-associated metabolite set enrichment analysis (MSEA) for TPS vs. CONT. (A) Plot showing affected biological processes in transgenerationally stressed (TPS) compared to non-stressed (CONT) rats based on metabolites identified as significantly altered between these groups. The p-value for each pathway is shown using the heatmap on the right of the figure and the fold enrichment shows how many times greater than chance the process is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in TPS and non-stressed CONT rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis gives the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. TPS had the highest impact on valine, leucin, and isoleucine biosynthesis.
Metabolites 12 00572 g002
Figure 3. Pathway-associated metabolite set enrichment analysis (MSEA) for MPS vs. CONT. (A) Plot displaying affected biological processes in multigenerationally stressed (MPS) compared to non-stressed control (CONT) rats based on metabolites identified as significantly different between these groups. The p-value for each pathway is indicated in the gradient on the right side. The fold enrichment shows how many times greater than chance the pathway is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in MPS and non-stressed CONT rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis indicates the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. MPS had the highest impact on the glycine, serine, and threonine metabolism.
Figure 3. Pathway-associated metabolite set enrichment analysis (MSEA) for MPS vs. CONT. (A) Plot displaying affected biological processes in multigenerationally stressed (MPS) compared to non-stressed control (CONT) rats based on metabolites identified as significantly different between these groups. The p-value for each pathway is indicated in the gradient on the right side. The fold enrichment shows how many times greater than chance the pathway is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in MPS and non-stressed CONT rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis indicates the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. MPS had the highest impact on the glycine, serine, and threonine metabolism.
Metabolites 12 00572 g003
Figure 4. Pathway-associated metabolite set enrichment analysis (MSEA) for TPS vs. MPS. (A) Plot showing affected biological processes in transgenerationally stressed (TPS) compared to multigenerationally stressed (MPS) rats, based on metabolites identified as significantly altered between these groups. The p-value for each pathway is shown using the heatmap on the right of the figure and the fold enrichment shows how many times greater than chance the process is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in TPS vs. MPS rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis indicates the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. The largest difference in pathway activity between TPS and MPS was observed in the glycine, serine, and threonine metabolism.
Figure 4. Pathway-associated metabolite set enrichment analysis (MSEA) for TPS vs. MPS. (A) Plot showing affected biological processes in transgenerationally stressed (TPS) compared to multigenerationally stressed (MPS) rats, based on metabolites identified as significantly altered between these groups. The p-value for each pathway is shown using the heatmap on the right of the figure and the fold enrichment shows how many times greater than chance the process is involved. (B) Metabolomic pathway analysis showing all matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis in TPS vs. MPS rats. The y-axis shows the negative natural log of p, such that a higher value on the y-axis gives a lower p-value. The x-axis indicates the pathway impact, which correlates to the number of metabolite hits in a particular pathway. Only metabolic pathways with p < 0.05 are labeled. The largest difference in pathway activity between TPS and MPS was observed in the glycine, serine, and threonine metabolism.
Metabolites 12 00572 g004
Figure 5. Illustration summarizing the experimental design generating lineages of transgenerational prenatal stress (TPS) and multigenerational prenatal stress (MPS) through filial generation zero (F0) to F3.
Figure 5. Illustration summarizing the experimental design generating lineages of transgenerational prenatal stress (TPS) and multigenerational prenatal stress (MPS) through filial generation zero (F0) to F3.
Metabolites 12 00572 g005
Table 1. p-values of blood plasma metabolites found to be significant in male Long-Evans rats in either a Mann–Whitney U test, the variable importance analysis based on random variable combination (VIAVC), or both. Rats were either non-stressed (CONT), transgenerationally-stressed (TPS), or multigenerationally-stressed (MPS). Up- or downregulation of the metabolites is also indicated. Metabolites with multiple peaks are represented as metabolite.1, metabolite.2… metabolite.n.
Table 1. p-values of blood plasma metabolites found to be significant in male Long-Evans rats in either a Mann–Whitney U test, the variable importance analysis based on random variable combination (VIAVC), or both. Rats were either non-stressed (CONT), transgenerationally-stressed (TPS), or multigenerationally-stressed (MPS). Up- or downregulation of the metabolites is also indicated. Metabolites with multiple peaks are represented as metabolite.1, metabolite.2… metabolite.n.
GroupMetaboliteNMR Chemical Shift Range of Bin (ppm)Mann-Whitney U TestVIAVCVIP ScoreRegulation
F3-TPS vs. F3-CONTCreatine phosphate3.035641–3.0281.62 × 10−24.40 × 10−91.51Up
Formate8.5343705–8.4422.83 × 10−22.09 × 10−61.51Up
Glucose.13.715–3.7041.09 × 10−17.02 × 10−60.87Down
Leucine.1, isoleucine.1, 2-hydroxyisovalerate.10.9586–0.94767.27 × 10−21.12 × 10−470.84Down
Alanine1.498–1.48784.85 × 10−22.35 × 10−200.81Down
Glucose.23.526–3.5141233.68 × 10−11.20 × 10−70.80Down
Leucine.2, isoleucine.20.9812075–0.96824.85 × 10−21.85 × 10−300.72Down
2-Hydroxyisovalerate.2, 2-oxoisocaproate0.9476–0.91336155.70 × 10−16.83 × 10−220.40Down
F3-MPS vs. F3-CONTSinglet at 8.38 ppm8.442–8.3885841.09 × 10−13.46 × 10−371.54Up
Formate8.5343705–8.4427.27 × 10−23.03 × 10−541.14Up
Creatine phosphate3.035641–3.0282.83 × 10−22.39 × 10−21.09Up
3-Methylxanthine8.388584–8.06758753.68 × 10−14.02 × 10−311.02Up
Threonine3.5998105–3.594.85 × 10−21.66 × 10−200.94Up
Tyramine.17.0614405–6.056.83 × 10−11.25 × 10−200.82Up
Glucose, betaine3.28–3.26840657.27 × 10−22.81 × 10−310.64Down
Tyramine.27.2946365–7.06144059.33 × 10−13.90 × 10−270.35Up
F3-TPS vs. F3-MPSCitrate.12.68–2.67286.50 × 10−25.94 × 10−111.64Down
Citrate.22.543–2.52698.30 × 10−23.34 × 10−121.40Down
Citrate.32.5269135–2.5111.05 × 10−17.66 × 10−111.31Down
Citrate.42.6727785–2.622971.05 × 10−12.25 × 10−81.17Down
Choline3.22–3.2032.07 × 10−21.75 × 10−130.91Down
Lactate4.137413–4.126813.79 × 10−2-0.91Down
Succinate2.412–2.4021.95 × 10−11.01 × 10−100.89Up
Alanine.11.487783–1.4768.30 × 10−29.29 × 10−160.89Down
Alanine.21.498–1.48782.34 × 10−13.88 × 10−100.83Down
Serine.1, creatine3.971–3.9372.34 × 10−14.95 × 10−80.77Up
Serine.2 3.8517095–3.8494.99 × 10−25.3 × 10−30.65Down
Alanine.33.8066745–3.7985751.61 × 10−11.27 × 10−80.45Down
Tyramine7.2946365–7.06144052.34 × 10−12.99 × 10−80.36Up
Table 2. Spearman’s rank correlations of blood plasma metabolite concentrations of creatine phosphate and leucine/isoleucine in relation to open field (OF) behavioural measurements. Correlation coefficients were obtained utilizing metabolite concentrations from comparisons of non-stressed (CONT) versus multigenerationally-stressed (MPS) Long-Evans rats and non-stressed (CONT) versus transgenerationally-stressed (TPS) Long-Evans rats. *—indicates a correlation with a p-value less than 0.05; **—indicates a correlation with a p-value less than 0.01.
Table 2. Spearman’s rank correlations of blood plasma metabolite concentrations of creatine phosphate and leucine/isoleucine in relation to open field (OF) behavioural measurements. Correlation coefficients were obtained utilizing metabolite concentrations from comparisons of non-stressed (CONT) versus multigenerationally-stressed (MPS) Long-Evans rats and non-stressed (CONT) versus transgenerationally-stressed (TPS) Long-Evans rats. *—indicates a correlation with a p-value less than 0.05; **—indicates a correlation with a p-value less than 0.01.
Creatine PhosphateLeucine/Isoleucine
TestTPS vs. CONTMPS vs. CONTTPS vs. CONTMPS vs. CONT
Rhop-valueRhop-valueRhop-valueRhop-value
Total Distance −0.712 **0.0090.646 *0.023
Number of Vertical Moves −0.636 *0.026
Vertical Time−0.730 **0.007−0.691 *0.0130.660 *0.02
Central Distance −0.635 *0.0260.674 *0.0160.656 *0.02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Heynen, J.P.; Paxman, E.J.; Sanghavi, P.; McCreary, J.K.; Montina, T.; Metz, G.A.S. Trans- and Multigenerational Maternal Social Isolation Stress Programs the Blood Plasma Metabolome in the F3 Generation. Metabolites 2022, 12, 572. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12070572

AMA Style

Heynen JP, Paxman EJ, Sanghavi P, McCreary JK, Montina T, Metz GAS. Trans- and Multigenerational Maternal Social Isolation Stress Programs the Blood Plasma Metabolome in the F3 Generation. Metabolites. 2022; 12(7):572. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12070572

Chicago/Turabian Style

Heynen, Joshua P., Eric J. Paxman, Prachi Sanghavi, J. Keiko McCreary, Tony Montina, and Gerlinde A. S. Metz. 2022. "Trans- and Multigenerational Maternal Social Isolation Stress Programs the Blood Plasma Metabolome in the F3 Generation" Metabolites 12, no. 7: 572. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12070572

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop