Metabolomic Analysis of Plasma from GABAB(1) Knock-Out Mice Reveals Decreased Levels of Elaidic Trans-Fatty Acid

Mice lacking the GABAB(1) subunit of gamma-aminobutyric acid (GABA) type B receptors exhibit spontaneous seizures, hyperalgesia, hyperlocomotor activity, and memory impairment. Although mice lacking the GABAB(1) subunit are viable, they are sterile, and to generate knockout (KO) mice, it is necessary to cross heterozygous (HZ) mice. The aim of our study was to detect the metabolic differences between the three genotypes of GABAB(1) KO mice in order to further characterize this experimental animal model. Plasma samples were collected from wild-type (WT), HZ, and KO mice. Samples were analyzed by means of a gas chromatography-mass spectrometry (GC-MS) platform. Univariate t-test, and partial least square discriminant analysis (PLS-DA) were performed to compare the metabolic pattern of different genotypes. The metabolomic analysis highlighted differences between the three genotypes and identified some metabolites less abundant in KO mice, namely elaidic acid and other fatty acids, and chiro-inositol.


Introduction
Gamma-aminobutyric acid (GABA) regulates inhibitory neurotransmission in the adult mammalian brain acting on GABA A and GABA B receptors, and both types of receptors are involved in the fine and precise regulation of neuronal excitability. GABA B receptors are G protein-coupled receptors and exist as heterodimers composed by the R1 and R2 subunits encoded by different genes [1]. Accumulating evidences suggest that GABA receptors may promote cell migration, differentiation, and synaptogenesis, therefore affecting Central Nervous System (CNS) development. Accordingly, the knocking down of GABA B receptors in mice, by eliminating the gene coding for the R1 subunit (GABA B(1) knock-out), causes behavioral alterations and abnormalities including epilepsy, impaired memory, hyperalgesia, hyperthermia, and hyperactivity [2][3][4][5][6]. Moreover, recent findings in our laboratory demonstrate that both GABA B(1) knockout (KO) and heterozygous (HZ) mice consumed high amounts of ethanol and prefer to drink alcohol rather than water [7] and could be used as an experimental model to study alcohol addiction.
All these observations suggest a rearrangement in the neuronal networks and altered neuronal excitability in the CNS of these mice [7] accompanied by changes in expression of a number of genes. Accordingly, GABA B(1) KO mice are characterized by altered levels of mRNA encoding for the different GABA A receptor subunits [7] as well as growth factors, transcription factors, and enzymes encoding Metabolites 2020, 10, 484 2 of 12 genes (Follesa et al. unpublished observations), and these alterations in gene expression may be regulated by changes in DNA methylation [8,9]. Nevertheless, little is known about the functional consequences that these changes can produce, and to the best of our knowledge, no information is available on metabolomic studies performed in GABA B(1) KO mice. To further characterize this mouse strain, the aims of our research were to investigate the metabolic differences between plasma from three genotypes of mice obtained by crossing GABA B(1) HZ mice. The ultimate goal of this study was to characterize GABA B(1) KO mice and identify metabolic routes than could be exploited to further understand the molecular, cellular, and functional alterations that could play a pivotal role in epilepsy, impaired memory, hyperalgesia, hyperthermia, hyperactivity [2][3][4][5][6], changes in gene expression, and drinking behavior [7] observed in these mice and give some new insights for intervention to treat and/or restore these CNS alterations. Biological samples, such as blood, can provide comprehensive information being minimally invasive to obtain, rendering it a good medium for clinical diagnostics.

Results
Plasma samples from 82 mice (32 WT, 41 HZ, and 9 KO) underwent metabolomic analysis. The three possible genotypes were compared in different combinations: WT vs. HZ, WT vs. KO, and HZ vs. KO.
WT vs. HZ: After a preliminary outlier removal, the groups resulted formed by 29 WT and 31 HZ. Univariate analysis did not reveal any significant metabolite. Partial least square discriminant analysis (PLS-DA) analysis allowed to build a statistically validated model (p < 0.01): the relative two-dimensional (2D) score plot and the corresponding VIP (Variable Importance in Projection) score were plotted and are reported in Figure 1a,b, respectively.
Metabolites 2020, 10, x FOR PEER REVIEW 2 of 13 different GABAA receptor subunits [7] as well as growth factors, transcription factors, and enzymes encoding genes (Follesa et al. unpublished observations), and these alterations in gene expression may be regulated by changes in DNA methylation [8,9]. Nevertheless, little is known about the functional consequences that these changes can produce, and to the best of our knowledge, no information is available on metabolomic studies performed in GABAB(1) KO mice. To further characterize this mouse strain, the aims of our research were to investigate the metabolic differences between plasma from three genotypes of mice obtained by crossing GABAB(1) HZ mice. The ultimate goal of this study was to characterize GABAB(1) KO mice and identify metabolic routes than could be exploited to further understand the molecular, cellular, and functional alterations that could play a pivotal role in epilepsy, impaired memory, hyperalgesia, hyperthermia, hyperactivity [2][3][4][5][6], changes in gene expression, and drinking behavior [7] observed in these mice and give some new insights for intervention to treat and/or restore these CNS alterations. Biological samples, such as blood, can provide comprehensive information being minimally invasive to obtain, rendering it a good medium for clinical diagnostics.

Results
Plasma samples from 82 mice (32 WT, 41 HZ, and 9 KO) underwent metabolomic analysis. The three possible genotypes were compared in different combinations: WT vs. HZ, WT vs. KO, and HZ vs. KO.
WT vs. HZ: After a preliminary outlier removal, the groups resulted formed by 29 WT and 31 HZ. Univariate analysis did not reveal any significant metabolite. Partial least square discriminant analysis (PLS-DA) analysis allowed to build a statistically validated model (p < 0.01): the relative twodimensional (2D) score plot and the corresponding VIP (Variable Importance in Projection) score were plotted and are reported in Figure 1a and 1b, respectively.
WT vs. KO: Two groups formed by 29 (WT) and 9 (KO) were submitted to statistical analysis. The t-test revealed the elaidic acid as statistically different between groups (p = 4.1417 × 10 −4 , FDR (False Discovery Rate) = 0.0207). The unequivocal identification of elaidic acid, in comparison with the more abundant cis-isomer oleic acid, is reported as Supplementary Figures S1-S3.
The PLS-DA analysis allowed to build the model whose plots are reported in Figure 2. As evident from the VIP score plot, all metabolites resulted more abundant in WT, except for mannitol ( Figure 2b). The PLS-DA analysis allowed to build the model whose plots are reported in Figure 2. As evident from the VIP score plot, all metabolites resulted more abundant in WT, except for mannitol ( Figure 2b). HZ vs. KO: Two groups formed by 34 (HZ) and 9 (KO) were submitted to statistical analysis. Again, the t-test revealed the elaidic acid as statistically different between groups (p = 6.9305 × 10 −4 , FDR = 0.034652). The PLS-DA analysis allowed to build the model whose plots are reported in Figure 3. HZ vs. KO: Two groups formed by 34 (HZ) and 9 (KO) were submitted to statistical analysis. Again, the t-test revealed the elaidic acid as statistically different between groups (p = 6.9305 × 10 −4 ,   The statistical parameters of the three PLS-DA models are reported in Table 1, and the corresponding metabolite variations are summarized in Table 2. Fatty acids, namely elaidic, linoleic, myristic, and palmitoleic acids, appeared to be less abundant in KO mice in comparison with both WT and HZ. Analogous behavior was shown by chiro-inositol and an undefined inositol isomer (other than myo-and scyllo-inositol). Due to the apparent similarity of WT and HZ, the combined group WT + HZ was compared with KO (69 vs. 9 samples), obtaining the PLS-DA model reported in Figure 4. This model was best described by the first three components, and the cross validation provided accuracy = 0.88732, R 2 = 0.6844, Q 2 = 0.23364. The permutation test on 100 permutations gave p = 0.4.
As expected, all metabolites resulted more abundant in the group formed by WT and HZ, except for trimethanolmethylamine, 2-monostearine, and ethanolamine (Figure 4b).   Elaidic acid resulted as the statistically most significant metabolite in all comparisons involving KO mice, being found in lower amounts in plasma when compared either with HZ or WT. Elaidic acid, the trans isomer of oleic acid, mainly derives from diet, being a component of ruminant fat, meat, and dairy products. To investigate the source of elaidic acid in our samples, the mice food was analyzed. The GC-MS analysis of the lipophilic phase revealed a fatty acid composition in agreement with what was declared by the food manufacturer, with linoleic, oleic, and palmitic acid being the most abundant. We found a small amount of elaidic acid, around 1% of total fatty acid content, most likely derived from the whey powder declared by the manufacturer in the food formula.

Discussion
The comparison of the metabolite content of plasma samples from different genotype mice revealed several interesting features in the animals lacking GABA B receptors. KO mice were characterized by lower amounts of fatty acids when compared with HZ or WT (alone or grouped together). In particular, the most important fatty acids were elaidic, linoleic, myristic, and palmitoleic. The greatest impact can be ascribed to elaidic acid, its different concentration between KO and either WT or HZ is clearly revealed also by univariate (t-test) analysis. Elaidic acid is the trans-isomer of C18 mono-unsaturated oleic acid, its occurrence in biologic samples is mainly ascribed to exogenous sources. Accordingly, elaidic acid is usually found in ruminant-derived food [10], and is degraded by β-oxidation in mitocondria, as reported by Yu [11], although with a less efficient mechanism than saturated or cis-unsaturated fatty acids. Fluctuation in the concentrations of these compounds may affect cellular functions that could be involved in the pathophysiology and behavioral phenotype observed in these KO mice. It has been reported that elaidic acid supplementation on pregnant mice induced a pro-inflammatory and adipogenic transcriptional profile compared with the oleic cis-isomer, leading to global adipose tissue DNA hypermethylation in the progeny [8]. Moreover, elaidic acid and linolelaidic acid enhanced extracellular ATP-induced apoptosis, with elevated activation of the apoptosis signal-regulating kinase 1 ASK1 pathway in a macrophage-like cell line. Trans fatty acids promote extracellular ATP-induced apoptosis by targeting ASK1 and suggest new TFA-associated pathways leading to inflammatory signal transduction and cell death leading to onset and progression of atherosclerosis [12]. We cannot rule out the possibility that trans fatty acids may affect apoptosis also in other cell types, including neurons, and that the reduction of trans fatty acid concentrations we observed in GABA B(1) KO mice could play a role in regulating the increased hippocampal neurogenesis measured in GABA B(1) KO mice (Follesa et al., unpublished observation) by reducing hippocampal apoptosis. These alterations in hippocampal cell proliferation might be part of a neuronal circuitry rearrangement in the CNS of GABA B(1) KO mice that involves also the gene expression of GABA A receptors and related behavior [7] and elaidic acid may alter the function of GABA A receptors. Accordingly, it was found that a mixture of eight fatty acids, comprising elaidic acid, produce anxiolytic-like effects comparable to diazepam in Wistar rats, exerting these effects through GABA A receptors [13]. In agreement with these observations, GABA B(1) KO mice showed an anxiety behavior [14] as well as depressive [15,16] and the lower levels of elaidic acid we observed may exacerbate these condition in GABA B(1) KO mice.
The lower amounts of fatty acids in KO mice could be ascribed to a higher degradation in GABA B(1) KO mice as these mice have altered metabolic activity. These metabolic alterations may be necessary to counteract the higher energy expenditure [17]. In our investigation, we did not measure food intake, but we did not observe significant changes in body weight between the three genotypes. Other studies demonstrated that GABA B(1) KO mice consume about 15% more food compared to control WT mice and this was accompanied by a strong reduction in plasma leptin levels [18], but at the same time, they show a 50% reduction of white adipose tissue (WAT), and increase in the weight of organs such as heart and brain [18]. This higher energy expenditure need may be the result of hyperlocomotor activity, wild running, as well as frequent spontaneous seizures. Accordingly, GABA B(1) KO mice are characterized by spontaneous seizures, hyperalgesia, hyperlocomotor [2][3][4][5][6] activity, and memory impairment [19]. On the other hand trans fatty acids may act on DNA methylation that could regulate the changes in gene expression observed in these mice [7]. All these allostatic adaptive mechanisms, adopted to survive the lack of GABA B receptors, may also be responsible for the large variability observed in some of the parameters we measured in these germinal knock-out mice.
Two other metabolites found in lower amounts in KO mice (and also in HZ when compared with WT) were chiro-inositol and a not precisely identified inositol isomer (other than myo-and scyllo-inositol). Chiro-inositol is one of the nine stereoisomers of 1,2,3,4,5,6-hexahydroxycyclohexane, and myo-inositol is the most commonly found in nature. Chiro-inositol in vivo is converted in myo-inositol [20], and is a component of an inositolphosphoglycan insulin mediator found in rat liver, stimulating pyruvate dehydrogenase phosphatase, kinetically similar to insulin. Its diet supplementation was found to increase the effect of insulin in patients with the polycystic ovary syndrome, improving ovulatory function and decreasing serum androgen concentrations, blood pressure, and plasma triglyceride concentrations [21]. In view of the above considerations, we cannot exclude that GABA B(1) KO mice could have a dysregulation in insulin and insulin-mediated pathways that needs to be further investigated in future studies.
The limitation of this study essentially lies in the low number of samples, especially for the KO group. This is a pilot study and the data reported, although very preliminary, allow to extract some interesting data leading the way to supplementary studies on bigger cohorts.

Materials and Methods
All procedures were performed in accordance with the European legislation EU Directive 2010/63, the "Guide for Care and Use of Laboratory Animals" adopted by NIH, USA (8th edition, 2011) and were approved by the Animal Ethics Committee of the

Animals and Sample Study Population
Balb/c GABA B(1) germinal KO mice were generously provided by Dr. Martin Gassmann (University of Basel, Switzerland) and have been previously described [19]. In order to prevent genetic drift, after about 20 generations (F20), the Balb/c HZ mice were crossed with WT mice of the FVB (Friend Virus B NIH Jackson) strain to obtain new HZ mice (F1). The F1 HZ mice where then crossed again to keep the colony until the next F20. Therefore, the strain of animal used in our laboratories is different from the original [7] due to the crossing with the FVB strain. Over the years, our colony of mice was maintained by crossing HZ mice since GABA B(R1) knock-out mice are sterile and therefore, to obtain KO mice, we need to cross HZ ± GABA B(R1) mice [22]. The mice used in the present study derived from independent mating of the F12-F16 generation. Litters were genotyped by polymerase chain reaction (PCR) using specific primers as described previously [7,23] and only males were selected to be utilized at 8 to 10 weeks of age with a mean ± SEM (standard error mean) of 9.474 ± 0.1254 weeks with a range from 8 to 10.71 weeks, and body weights ranging from 21 to 31 g, and the means ± SEM for the different genotypes were WT 26.75 ± 0.50 (n = 32), HZ 25.75 ± 0.41 (n = 41), and KO 26.56 ± 0.91 (n = 9). In agreement with other studies [2,18], we did not observe any statistical differences in body weight between the three different genotypes (One-way analysis of variance (ANOVA), p-value = 0.2872).

Samples Preparation and GC-MS Analysis
A total number of 82 mice were used: 32 WT, 41 HZ, and 9 KO. Blood was collected from the trunk into heparinized tubes and centrifuged at 900× g for 10 min at 4 • C, as previously described [24]. The resulting plasma supernatant was frozen at −80 • C until GC-MS analysis. Plasma samples were analyzed as previously reported for human samples [25]. The method required slight adjustments for the low amount of plasma available. In brief, 50 µL of plasma were deproteinized with cold methanol (200 µL) and centrifuged for 10 min at 16,900× g. The upper phase (150 µL) was transferred in glass vials (1.5 mL) with glass insert (200 µL) and evaporated to dryness. Methoxyamine hydrochloride (30 µL, 0.24 M solution in pyridine) was added to each sample and left for 17 h at room temperature. MSTFA (N-Methyl-N-trimethylsilyltrifluoroacetamide, 30 µL) was added and reacted for 1 h at room temperature. The analytical platform was an Agilent 5977B GC/MS interfaced to the GC 7890B (Agilent Technologies, Palo Alto, CA, USA), with a DB-5 ms column (Agilent J&W Scientific, Folsom, CA, USA). Instrumental parameters were previously reported [25]. One µL of each sample was injected in split mode (1:2), sample depth 6 mm. Mass spectra were acquired in full scan mode using 2.28 scans/s with a mass range of 50-700 Amu. Each chromatogram was analyzed by the free software AMDIS (Automated Mass spectral Deconvolution and Identification System) (http: //chemdata.nist.gov/mass-spc/amdis), allowing the automatic identification of each chromatographic peak by comparison of the relative mass spectra and retention times with those stored in an in-house library of 255 metabolites. Other metabolites, not included in the in-house library, were identified using NIST14 (National Institute of Standards and Technology's mass spectral database) and the Golm Metabolome Database (GMD) (http://gmd.mpimp-golm.mpg.de/). In this way, 67 compounds were detected: following the identification levels defined by the Metabolomics Standards Initiative (MSI) [26], 53 were "confidently identified compounds" (level 1), 3 "putatively annotated compounds" (level 2), 9 "putatively annotated compound class" (level 3), and 2 "unknown compounds". The Excel data sheet obtained (67 metabolites × 82 samples) was then submitted to statistical analysis.

Statistical Analysis
The data matrix was processed with the web-based tool MetaboAnalyst 4.0 (http://www. metaboanalyst.ca/) [27]. Missing values were replaced with the half of the minimum positive values in the original data, data were normalized by sum, log transformed, and Pareto scaled. Univariate analysis (t-test), principal component analysis (PCA), partial least square discriminant analysis (PLS-DA), and its associated variable importance in projection (VIP) score were obtained. PCA allowed outlier detection: outlier samples were removed before building the PLS-DA models. PLS-DA models were tested with the leave-one-out cross validation (LOOCV) method for the evaluation of statistical parameters (accuracy, correlation coefficient R 2 , cross validation coefficient Q 2 ), determining the optimal number of components for each model description. Due to the tendency of PLS-DA to over fit data, all models were validated in order to understand whether the separation is statistically significant using the B/W permutation test [28]. Statistical differences in body weight by one-way ANOVA followed by Dunnett's multiple comparisons test were performed using GraphPad Prism version 8.0.0 for Windows, GraphPad Software, San Diego, California USA, www.graphpad.com.

Mice Food Analysis
A piece of food (0.45 g) was ground in a glass mortar and stirred in a round bottom flask with water (10 mL), chloroform (10 mL), and methanol (10 mL) for 2 h. The two phases were separated in a separatory funnel (upper layer hydrophilic phase, lower layer lipophilic phase), dried, and derivatized. The hydro-phase was methoxymated/silylated and the lipo-phase was methylated, and both were analyzed by GC-MS, as previously reported [29].

Conclusions
Plasma samples from GABA B(1) KO mice were compared with HZ and WT from a metabolomic point of view. The KO phenotype resulted characterized by lower amounts of plasmatic fatty acids, the most important being elaidic acid, and two inositol isomers. The different metabolite content may regulate cellular and molecular modifications that could be related with some of the peculiar behaviors of KO mice, mainly related to CNS impairment.