SLC13A5/sodium-citrate co-transporter overexpression causes disrupted white matter integrity and an autistic-like phenotype

Endoplasmic reticulum-based N ɛ -lysine acetylation serves as an important protein quality control system for the secretory pathway. Dysfunctional endoplasmic reticulum-based acetylation, as caused by overexpression of the acetyl coenzyme A transporter AT-1 in the mouse, results in altered glycoprotein ﬂ ux through the secretory pathway and an autistic-like phenotype. AT-1 works in concert with SLC25A1, the citrate/malate antiporter in the mitochondria, SLC13A5, the plasma membrane sodium/citrate symporter and ATP citrate lyase, the cytosolic enzyme that converts citrate into acetyl coenzyme A. Here, we report that mice with neuron-speci ﬁ c overexpression of SLC13A5 exhibit autistic-like behaviours with a jumping stereotypy. The mice displayed disrupted white matter integrity and altered synaptic structure and function. Analysis of both the proteome and acetyl-proteome revealed unique adaptations in the hippocampus and cortex, highlighting a metabolic response that likely plays an important role in the SLC13A5 neuron transgenic phenotype. Overall, our results support a mechanistic link between aberrant intracellular citrate/acetyl coenzyme A ﬂ ux and the development of an autistic-like phenotype. Abbreviations: ACLY = ATP citrate lyase; ASD = autism spectrum disorder; DiLeu = dimethylated leucine; DIV = days in vitro DTI = diffusion tensor imaging; DWI = diffusion-weighted imaging; ER = endoplasmic reticulum; EtOH = ethanol; FA = fractional anisotropy; FC = fear conditioning; FDR = false discovery rate; fEPSP = ﬁ eld excitatory post-synaptic potentials; GO = gene ontology; KEGG = kyoto encyclopedia of genes and genomes; LD = light – dark; LTD = long-term depression; LTP = long-term potentiation; MB = marble burying; MD = mean diffusivity; MEA = multi-electrode array; NA = numerical aperture; NODDI = neurite orientation dispersion and density imaging; NOR = novel object recognition; ODI = orientation dispersion index; OF = open ﬁ eld; nTg = neuron transgenic; PFA = paraformaldehyde; PO = propylene oxide; ROI = region of interest; SD = standard deviation; SEM = standard error of the mean; SI = social interaction; TEM = transmission electron microscopy; TBSS = tract-based spatial statistics; TCA = tricarboxylic acid cycle; WT = wild-type


Introduction
Generation and maintenance of neuronal synapses require an efficient and properly functioning secretory pathway in order to deliver the necessary proteins such as neurotransmitter receptors and ion channels. Therefore, quality control measures are vitally important in order to ensure proper function of the proteins delivered to the synapse. N ɛ -lysine acetylation within the endoplasmic reticulum (ER) serves as an important protein quality control system within the secretory pathway that (i) positively selects properly folded glycoproteins and (ii) regulates the induction of reticulophagy (ER-specific autophagy) to dispose of misfolded aggregates.  Imbalanced ER-based N ɛ -lysine acetylation profoundly alters glycoprotein flux through the secretory pathway, which impacts multiple organelles beyond the ER including the nucleus, mitochondria and lysosomal network. 11,14,[18][19][20] Aberrant ER-based N ɛ -lysine acetylation can result in multiple disease states in the nervous system. 3,7,10,11,[15][16][17]22,23 For example, loss-of-function mutations in the ER acetyl-CoA transporter AT-1/SLC33A1 are causative of hereditary spastic paraplegia 42 and are also associated with developmental disability and premature death. 22,23 On the other hand, gene duplication events of 3q25.31 containing AT-1/SLC33A1 are associated with autism spectrum disorder (ASD), dysmorphism and intellectual disability. 24 Mouse models with these genetic aberrations recapitulate the associated human diseases. 10,11,15 We recently reported that overexpression of AT-1/ SLC33A1 in mouse forebrain neurons, termed the AT-1 neuron transgenic (nTg), resulted in an autistic-like phenotype with increased dendritic branching and spine formation, altered synaptic plasticity and significant upregulation of proteins involved in synaptic generation and maintenance. 11 This phenotype was attributed to changes in glycoprotein flux through the secretory pathway due to the increased flux of acetyl-CoA from the cytosol into the ER, which altered N ɛ -lysine acetylation. In order to replenish the cytosolic acetyl-CoA, AT-1 overexpressing neurons upregulated SLC25A1/CIC, the mitochondrial citrate/malate exchanger and ATP citrate lyase (ACLY), the cytosolic enzyme that converts citrate into acetyl-CoA. 11 Thus, the neuron adapted in order to maintain the flux of cytosolic acetyl-CoA into the ER, revealing a coordinated metabolic pathway. In addition to AT-1/SLC33A1, gene duplication events of 22q11.21 (containing SLC25A1), 17q21.2 (containing ACLY) and 17p13.1 (containing SLC13A5) are all associated with ASD. 24 SLC13A5, also known as the sodium/citrate co-transporter or I'm Not Dead Yet, is a plasma membrane symporter that couples transport of extracellular dicarboxylate and tricarboxylates, mostly citrate, with four Na + cations. 25 Therefore, there appears to be a strong association between the development of ASD and aberrant intracellular citrate/acetyl-CoA flux within the cell, of which changes in ER-based N ɛ -lysine acetylation likely plays an important role. As such, we hypothesized that modulating the intracellular flux of citrate through overexpression of SLC13A5 would result in an autistic-like phenotype in the mouse, reminiscent of the AT-1 nTg mouse.
Here, we describe the generation and phenotypic characterization of a new mouse model that overexpresses SLC13A5 within forebrain neurons, using the same transgenic system as the previously reported AT-1 nTg autisticlike model. This SLC13A5 nTg mouse exhibits autistic-like behaviours with a jumping stereotypy, with evidence of disrupted white matter integrity, altered synaptic plasticity and widespread changes in the proteome and acetyl-proteome. Therefore, we conclude that the aberrant intracellular flux of citrate/acetyl-CoA within the neuron is a mechanistic driver for the development of an autistic-like phenotype.

Animals
Mice were housed in standard cages provided by the University Laboratory Animal Resources and grouped with 1-5 littermates per cage. Animals were supplied standard chow and water ad libitum. All animal experiments were performed in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals approved by the Institutional Animal Care and Use Committee of the University of Wisconsin-Madison (protocol #M005120; #M005327; #M019011). SLC13A5 nTg mice used throughout the study were heterozygous for both the TRE-SLC13A5 and the Camk2a-tTA transgenes. Non-transgenic, wild-type (WT) littermates were used as controls throughout the study. The specific age and sex of animals used for experiments are noted in the figure legends.

Behaviour testing
All behavioural assays were conducted within the Waisman Center Behavioral Testing Service (Madison, WI, USA) with the experimenter blind to the genotype of the mice during testing. A 30-min acclimation period in the testing room was performed prior to each behaviour assay. The following behavioural assays have been previously described 11 : marble burying (MB) assay, novel object recognition (NOR), social interaction (SI) and fear conditioning (FC) paradigm.

Open field exploration
Each mouse was removed from its home cage and placed in the centre of the arena and allowed to explore for a duration of 30 min. Each mouse received one open field (OF) session. The Omnitech Fusion system with photobeams continuously monitored and recorded the animal's placement during exploration. Quantified parameters included total distance travelled (cm), vertical activity episode count, total ambulation time (s) and distance travelled in the centroid (cm). Data were recorded using the Omnitech Fusion system with a centre ratio zone map.

Light/dark exploration
Each mouse was placed into a split arena for a 10-min assay. Each mouse received one session. Time spent (s) and number of entries into the light and dark portions of the arena were recorded.

Jumping activity
Videos were taken during the dark cycle with mice in their home cages using a red light for visualization. During each bout of jumping, the number of jumps per mouse was quantified in a 50 s time window.

Primary neuron culturing
Primary neurons were harvested, cultured and stained as previously described 11 on micro-dissected cortex and hippocampus; cultures were maintained for a maximum duration of 28 days. For immunocytofluorescence analysis of synapse formation, the following primary antibodies were used: Syn-1 (Cell Signal Technologies; 5927S; 1:200), Homer1 (Synaptic Systems; 160-011; 1:100), Psd-95 (Thermo Fisher; MAI-045; 1:500), Vglut2 (Synaptic Systems; 135-043; 1:500) and NF-200 (Thermo Fisher; PA3-16753; 1:1000). Images were collected on a Nikon A1 inverted confocal microscope using NIS-Elements AR version 5.11.01 software with 405 nm (blue channel), 488 nm (green channel), 561 nm (red channel) and 640 nm (far red) laser wavelengths with the Galvano scan head. Multi-z-stack images (1024 × 1024 pixels with 25 z-steps every 0.15 μm) were acquired using a 60× oil objective [numerical aperture (NA) = 1.4; 0.21 μm/pixel] at a pinhole size of 39.59 μm. Neurons and neurites were identified by positive NeuN and NF-200 staining, respectively. The .nd2 images were imported into Imaris (Bitplane; Version 9.5) and converted to native .ims format. For neuron morphology, the Filament Tracer module with the Autopath method was used to trace dendrites starting from the soma with a thinnest diameter of 1.5 μm; dendrite seed points were removed within 30 μm of the soma. Dendritic spines were subsequently detected with a seed point diameter of 0.7 μm and maximum length of 7 μm. A Sholl analysis with 1 μm-spaced spheres was used to quantify dendritic branching. For the assessment of synapse formation, signal from the nucleus was first removed by masking the pre-and post-synaptic marker channels using a surface reconstruction of the DAPI signal (setting voxels within the surface to zero). Two-micrometre diameter spots were fit to the pre-and post-synaptic marker signals and spots were considered colocalized if within 1 μm of each other. The spot counts per image were normalized to neurite volume, which was estimated from surface reconstruction of the NF-200 signal (neurons and soma) less the surface reconstruction of the DAPI signal (nucleus).
For multi-electrode array (MEA), 48-well MEA plates (Axion Biosystems; M768-tMEA-48B-5) were pre-coated in filter-sterilized 0.1% polyethyleneimine (Sigma-Aldrich; 181978-5G) diluted in borate buffer (Thermo Fisher; 28341) for 1 h. After rinsing 4 times with sterile water, the plates were dried overnight in the biosafety cabinet. The following day, embryonic hippocampal and cortical cell suspensions were generated as described above and supplemented with mouse laminin (Thermo Fisher; 23017015; 1 μg/ml). A total of 50 000 cells/well in a volume between 5 and 10 μl were plated in the centre of the well. After a 1-h incubation at 37°C and 5% CO 2 , 200 μl neuron culture media was added per well. The cells were maintained at 37°C and 5% CO 2 with a half media change every 3-4 days and spontaneous activity recordings were taken every 7 days in vitro (DIV) using the Maestro Pro Multiwell Multielectrode Array and Impedance System (Axion Biosystems). Media changes were avoided on the day of recording. Before recording, the plate equilibrated in the MEA recording chamber for 5 min at 37°C and 5% CO 2 . Data were collected for a 10-min duration using the Axis Navigator software (Axion Biosystems; Version 2.0.2.5) with the Neural Real Time configuration for continuous spontaneous activity. A band-pass filter of 3000 Hz (lowpass) to 200 Hz (high-pass) was applied with a variable threshold spike detector at +6 standard deviations (SD) of the root mean squared of the background noise. A minimum spike rate of 5 spikes/min was used to determine an active electrode and a mature network was considered to be a well with at least 8 out of 16 active electrodes. Bursts were detected in mature networks using the inter-spike interval threshold with a maximum inter-spike interval of 100 ms and minimum number of five spikes. Network bursts were detected in mature networks with a maximum inter-spike interval of 100 ms, minimum number of 50 spikes and minimum of 35% participating electrodes. A synchronicity window of 20 ms was applied. Mean firing rate (Hz) and burst rate (Hz) represent neural excitability while network burst rate (Hz) and network synchronicity (index value between 0 and 1) represent network synchronization.

Diffusion-weighted imaging
On P120, mice were anaesthetized with isoflurane then transcardially perfused with phosphate-buffered solution Raw data files were converted to NIfTI format. FMRIB Software Library (FSL) was used to correct for eddy current artefacts. A diffusion-weighted imaging (DWI)-based mouse brain atlas was used as a template and to define region of interests (ROIs), including the left and right hippocampus, amygdala and corpus callosum. 26 Multishell diffusion data were fit with the Microstructure Diffusion Toolbox to models of neurite orientation dispersion and density imaging (NODDI) and ActiveAx. 27 An additional compartment of isotropic restriction was included to account for potential fixative effects as recommended. 28 Whole-brain voxel-wise differences between WT and SLC13A5 nTg animals for all MC-DWI indices were determined as previously described with tract-based spatial statistics (TBSS). 29 An fractional anisotropy (FA) threshold of 0.2 was applied for the creation of the skeleton and permutation test results for multiple comparisons and threshold-free cluster enhancement was implemented with FSL's randomize to compare the SLC13A5 nTg group to the WT group, with P , 0.05 as the threshold for significance. Region-based differences between WT and SLC13A5 nTg animals were determined by calculating the indices of DWI in each diffusion model [FA, mean diffusivity (MD), NDI, ODI, IC, EC and stat] for each ROI and comparing the values between WT and SLC13A5 nTg animals via an unpaired t-test (significance P , 0.05); statistically significant differences were then determined after controlling for multiple comparisons with the Benjamini-Hochberg procedure with the false discovery rate (FDR) set to 0.05.

Electrophysiology
Extracellular recordings of field excitatory post-synaptic potentials (fEPSPs), long-term potentiation (LTP) and longterm depression (LTD) were conducted as previously described, 11 with the following modifications. Slice preparation solution contained the following: 124 mM NaCl, 1.25 mM NaH 2 PO 4 , 3 mM KCl, 25 mM NaHCO 3 , 10 mM glucose, 1 mM sodium ascorbate, 3 mM kynurenic acid, 3.6 mM MgSO 4 and 0.8 mM CaCl 2 . Recording artificial CSF contained the following: 124 mM NaCl, 1.25 mM NaH 2 PO 4 , 3 mM KCl, 25 mM NaHCO 3 , 15 mM glucose, 0.8 mM sodium ascorbate, 1.3 mM MgSO 4 and 2.5 mM CaCl 2 . All solutions were buffered to pH 7.3 when saturated with carbogen and had a confirmed osmolality between 294 and 297 mOsm. Recordings were taken from coronal slices using fire-polished borosilicate glass recording pipettes filled with 1 M NaCl (3-5 MΩ) and Pt/Ir concentric bipolar stimulating electrodes. For LTD induction, 15 μM N-methyl-D-aspartate (NMDA) was perfused over the slice for 5 min at a rate of 3 ml/min. Potentiation and depression were defined as the mean fEPSP slope during the last 10 min of the recording divided by the average of the last 10 min of baseline immediately preceding induction of LTP or LTD.

Synaptosome preparation
Crude synaptosomes (pre-and post-synaptic) were prepared from the fresh or snap-frozen hippocampus and cortex using Syn-PER synaptic protein extraction reagent (Thermo Fisher; 87793) supplemented with protease inhibitor cocktail (Millipore; 11836170001) following manufacturer's instructions.
Golgi staining images were pre-processed in ImageJ (Version 1.52) by first inverting then using the Background Subtract tool three consecutive times with a pixel size of 100, 50 and 25. Images were saved as .tiff, imported into Imaris (Bitplane; Version 9.5) and converted into native .ims format. Secondary dendritic branches were semimanually traced using the Autopath method of the Filament Tracer module with a diameter of 0.25 μm; automatic dendrite volume detection was enabled. To reduce computation time for spine detection, a mask was first created to eliminate signal distant from the reconstructed dendrites. Spines were identified with a minimum diameter of 0.25 μm, maximum length of 5 μm and enabling for detection of branch spines. Dendrite spine density (in spines per 10 μm dendrite length) and spine volume (μm 3 ) were extracted for external analysis.

Transmission electron microscopy
Following CO 2 euthanasia, the brain was extracted and the left hippocampus micro-dissected in PBS. The tissue was subsequently fixed in 2.5% glutaraldehyde in 0.1 M phosphate buffer (PB) overnight at 4°C. The fixed samples were rinsed 5 × 5 min in PB and post-fixed in 1% osmium tetroxide, 1% potassium ferrocyanide in 0.1 M PB for 1 h at room temperature, then rinsed in PB as before. Dehydration was performed in a graded EtOH (ethanol) series (35%, 50%, 70%, 80%, 90% for 10 min each step, 95% for 20 min, 100% for 2 × 10 min) at room temperature and 100% EtOH at 4°C overnight then transitioned in propylene oxide (PO) 2 × 7 min at room temperature. Fully dehydrated samples were infiltrated in increasing concentrations of PolyBed 812 (Polysciences Inc.) and PO mixtures. Embedding and polymerization took place in fresh PolyBed 812 for 24 h at 60°C. The samples were sectioned on a Leica EM UC6 ultramicrotome at 100 nm, collected on formvar coated 2 × 1 mm slot Cu grids (EMS Hatfield, PA) and post-stained with uranyl acetate and lead citrate. The sectioned samples were viewed at 80 kV on a Philips CM120 transmission electron microscope equipped with AMT BioSprint12 digital camera (AMT Imaging Systems).
Images were analysed on ImageJ (Version 1.52) using the GRatio plugin. 31 The G-ratio reported was calculated using the traced perimeter values (GRatioPerimeter). The axon diameter was calculated by dividing the traced axon perimeter by pi.

Proteomics
Cortical and hippocampal tissues were lysed in 8 M urea buffer containing protease and phosphatase inhibitors (Thermo Fisher) followed by a minute of sonication. The lysate was reduced with 10 mM dithiothreitol at 56°C for 30 min, then alkylated with 20 mM iodoacetamide at room temperature in the dark for an additional 15 min. Approximately 100 μg of protein was then digested with trypsin with a trypsin/protein ratio of 1:100 overnight at 37°C. Peptides were extracted and dried (SpeedVac; Eppendorf) and concentrations of the peptide mixture were measured. Samples were lyophilized and stored at −80°C until further processing.
Dimethylated leucine (DiLeu) tags were synthesized as previously described. 32 DiLeu tags were activated in anhydrous DMF combined with DMTMM and NMM at 0.7× molar ratio and vortexed at room temperature for 45 min. After centrifugation, the supernatant was used immediately for peptide labelling. DiLeu labelling was performed by addition of labelling solution at a 20:1 tag to digested peptide ratio by weight and vortexed at room temperature for 2 h. The labelling reactions were quenched by addition of hydroxylamine to a concentration of 0.25% and the labelled peptide samples were dried in vacuo. The samples were combined and cleaned with SCX SpinTips (Protea Biosciences) and desalted with Omix C18 pipet tips (Agilent).
The liquid chromatography-tandem mass spectrometry detection system consisted of a nanoflow high-performance liquid chromatograph instrument (Dionex UltiMate 3000 UPLC system; Thermo Fisher) coupled to a Q Exactive HF Orbitrap mass spectrometer (Thermo Fisher) with a nanoelectrospray ion source (Thermo Fisher). In brief, 0.5 μg of peptide mixture dissolved in buffer A [0.1% formic acid (FA)] was loaded onto a 75 μm × 15 cm fabricated column filled with 1.7 μm Bridged Ethylene Hybrid packing materials (130 Å; Waters) over a 126 min linear gradient of 3-45% Mobile Phase B (buffer A, 0.1% FA in water; buffer B, 0.1% FA in ACN) with a flow rate of 300 nl/min. The MS analysis was performed in a data-dependent manner using an Orbitrap mass analyser. For a full mass spectrometry survey scan, the target value was 1 × 10 5 and the scan ranged from 300 to 1500 m/z at a resolution of 60 000, with a maximum injection time of 100 ms. For the MS2 scan, up to 15 of the most intense precursor ions from a survey scan were selected for MS/MS and detected by the Orbitrap at a mass resolution of 15 000 at m/z 400. Only precursor ions with charge states of 2-6 were selected for fragmentation by higherenergy collision dissociation with a normalized collision energy of 30%. The automatic gain control for MS/MS was set to 8e3, with maximum ion injection times of 100 ms. Dynamic exclusion time was 45 s and the window for isolating the precursors was 1.4m/z.
Protein and peptide identification and quantification were conducted through MaxQuant version 1.5.3.8. 33 Raw files were searched against the Uniprot mouse reviewed database (August 2019) using the integrated Andromeda search engine with FDR , 1% at the peptide and protein level. Trypsin was selected as the enzyme with at most two missed cleavages. A reverse database for the decoy search was generated automatically in MaxQuant. Enzyme specificity was set to 'Trypsin', with a minimum number of seven amino acids were required for peptide identification. Static modifications included carbamidomethylation of cysteine residues (+57.02146 Da) and DiLeu labelling on N-terminus and lysine residues (+145.12801 Da). Dynamic modifications included oxidation of methionine residues (+15.99492 Da) and deamidation of asparagine and glutamine residues (+0.98402 Da). The first search mass tolerance was 20 ppm and the main search peptide tolerance was 4.5 ppm. The FDRs of the peptide-spectrum matches and proteins were set to ,1%. Quantitation was performed using Perseus software. 34 Briefly, the raw reporter ion intensity in each DiLeu channel was corrected for isotope impurities 32 and normalized for mixing differences by equalizing the total signal in each channel. In cases where no signal was detected in a channel, the missing value was assigned with the noise level of the original spectrum (noise-band capping of missing channels), and the resultant intensity was not corrected for impurities or normalized for uneven mixing.

Stoichiometry of protein acetylation
Cortical and hippocampal tissue was prepared and analysed as previously described. 18 Acetyl stoichiometry values are reported as a range from 0 to 1 representing 0-100% of a detected lysine site being endogenously acetylated, respectively. Changes in acetyl stoichiometry values (Δ acetyl stoichiometry) are consistently reported as the nTg value less the WT value. Pathway analysis and network plot construction were conducted using the R package enrichplot with an overrepresentation analysis using the Mus musculus organism database. 35 Relevant parameters included a minimum and a maximum number of genes for a category of 5 and 2000, respectively, with the Benjamini-Hochberg method for multiple test adjustment to a FDR of 0.05.

Statistical analysis
Data analysis was performed using GraphPad Prism version 9.0.1. Data are expressed as mean + SD unless otherwise specified. Comparison of the means was performed using an unpaired t-test for two groups and ordinary one-way or two-way ANOVA for ≥3 groups followed by either Tukey-Kramer (comparison between all groups) or Dunnett's (comparison to one control group) multiple comparisons test. If no sex differences were found via ANOVA testing, data were combined for ease of visualization. The details of each statistical test are described in the figure legends. Grubb's test was used to remove outliers, which are determined at P , 0.05. Differences in the mean were declared statistically significant if P , 0.05 and the following statistical significance indicators are used throughout the article: *P , 0.05; **P , 0.005; ***P , 0.0005.

Data availability
The proteomics data that support the findings of this study have been deposited to the ProteomeXchange Consortium (ID number PXD026624). The acetyl-proteomics data that support the findings of this study have been deposited to the ProteomeXchange Consortium (ID number PXD025424) and the MassIVE partner repository (ID number MSV000087209). The authors declare that all other data supporting the findings of this study are available within the paper and its Supplementary Data/Tables.

SLC13A5 nTg mice exhibit an autistic-like phenotype with jumping stereotypy
To evaluate our hypothesis that increased citrate flux into the cytosol can cause an autistic-like phenotype, we generated transgenic C57BL/6J mice that specifically overexpress human SLC13A5 in forebrain neurons, henceforth referred to as the nTg, using the same expression system as the previously reported AT-1 nTg model. 11 To ensure that our transgenic system was working as expected, we performed RT-qPCR on cDNA from isolated adult neurons using primers specific for human SLC13A5; as expected, we observed a PCR product only with SLC13A5 nTg cDNA and not in WT controls (Fig. 1A). Additionally, immunofluorescence staining of hippocampal brain slices revealed a robust increase in SLC13A5 protein in the SLC13A5 nTg slices that co-localized with the neuronal marker NeuN (Fig. 1B). These data support the successful generation of a transgenic mouse line that exhibits forebrain neuron-specific overexpression of SLC13A5.
Upon generation of the transgenic mouse line, it became apparent that the SLC13A5 nTg mice exhibited repetitive bouts of repetitive jumping behaviour in their home cage seen both during the day and night cycle while their WT littermates did not (Supplementary Movie 1). This behaviour was formally quantified with the mice undisturbed in their home cage by counting the number of jumps in a 50 s time period ( Fig. 2A). Interestingly, in an OF assay, the SLC13A5 nTg mice did not exhibit a change in vertical activity compared to WT controls, suggesting this jumping stereotypy was suppressed in an unfamiliar environment (Fig. 2B). The OF assay also revealed a decrease in the overall total distance travelled, distance travelled in the centre of the arena and ambulatory time (Fig. 2B). In a light-dark (LD) assay, the SLC13A5 nTg mice entered the light side less frequently but did not spend significantly more time on either side compared to WT controls (Fig. 2C). In sum, these data show that SLC13A5 nTg mice exhibit repetitive jumping behaviour while in their home cage but have reduced ambulatory activity in a novel environment.
We also performed several behavioural assays to gauge both learning and memory formation as well as autistic-like features like we did in the AT-1 nTg model. 11 In a MB assay, the SLC13A5 nTg mice buried less marbles compared to WT controls (Fig. 2D), which was also observed in the AT-1 nTg mice. 11 Furthermore, the SLC13A5 nTg mice spent more time investigating both novel and familiar objects in a NOR paradigm; however, they did not display a change in the percentage of time investigating the novel object (Fig. 2E). In a FC paradigm, the SLC13A5 nTg mice compared to WT controls displayed reduced freezing behaviour after conditioning (Day 2), but not during training (Day 1), revealing reduced learning and memory formation in the transgenic mice (Fig. 2F). Moreover, the SLC13A5 nTg mice preferred interacting with a mouse over an empty cup like WT mice in a SI assay, but when presented with a familiar or novel mouse, the SLC13A5 nTg mice did not significantly distinguish the two like WT controls (Fig. 2G). Overall, the behavioural abnormalities observed in the SLC13A5 nTg mice are similar to the AT-1 nTg in the MB, FC and SI assays, 11 leading us to conclude the SLC13A5 nTg mouse exhibits an autistic-like phenotype.

White matter integrity is disrupted in SLC13A5 nTg mice
Changes in white matter microstructure are commonly observed in individuals with ASD. 36,37 Therefore, we desired to assess the white matter integrity of our autistic-like SLC13A5 nTg model by performing ex vivo DWI using both ROI and TBSS analyses with three different diffusion models: diffusion tensor imaging (DTI), NODDI and ActiveAx. We selected the right/left hippocampus and right/left amygdala ROIs based upon the behavioural aberrations observed and included the right/left corpus callosum as comparison ROIs where we did not expect to observe differences. The DTI analysis revealed a decrease in FA within the left hippocampus and right/left amygdala (Fig. 3A); likewise, our TBSS analysis revealed multiple clusters of voxels with decreased FA (Fig. 3B). We also observed an increase in DTI MD in the TBSS analysis, but no changes in the ROI analysis survived corrections for multiple comparisons ( Fig. 3A and B). NODDI modelling did not reveal any changes in WT versus SLC13A5 nTg mice ( Supplementary  Fig. 1). ActiveAx revealed an increase in the 'extraneuronal/extracellular' EC parameter in the left hippocampus as well as a decrease in the 'intra-neuronal' IC parameter in the right/left amygdala and increase in the right corpus callosum (Fig. 3A). No changes in the ActiveAx stat parameter were observed (Supplementary Fig. 1). Taken together, these results suggest disruptions in the integrity of white matter that can be from a number of causes, including changes within the axonal tracts or myelin structure or composition. 38 To address potential changes in myelination, we first performed the basic histologic assessment of our SLC13A5 nTg mice and found no obvious anatomic differences ( Supplementary Fig. 2A). Additionally, we performed immunofluorescence staining for oligodendrocytes, which also did not reveal any changes from WT mice (Supplementary Fig. 2B and C). Finally, we processed the left hippocampus from WT and SLC13A5 nTg mice for transmission electron microscopy (TEM) for high-resolution  Open field (OF) assay. WT, n = 10; Tg, n = 10; equal males/ females. Total distance: P = 0.0158 via two-way ANOVA (genotype × sex; F statistic = 7.291 for genotype factor). Distance in centroid: P = 0.0005 via two-way ANOVA (genotype × sex; F statistic = 18.73 for genotype factor). Ambulatory time: P = 0.0023 via two-way ANOVA (genotype × sex; F statistic = 13.16 for genotype factor). (C) Light-dark (LD) exploration. WT, n = 10; Tg, n = 10; equal males/females. Light entries: P = 0.0302 via two-way ANOVA (genotype × sex; F statistic = 5.654 for genotype factor). (D) Marble burying (MB) assay. WT, n = 10; Tg, n = 10; equal males/females. Total: P = 0.0053 via two-way ANOVA (genotype × sex; F statistic = 10.420 for genotype factor). 50-99%: P = 0.0042 via two-way ANOVA (genotype × sex; F statistic = 11.160 for genotype factor). (E) Novel object recognition (NOR). WT, n = 10; Tg, n = 10; equal males/females. Investigation time (novel object): P = 0.0114 via two-way ANOVA (genotype × sex; F statistic = 8.172 for genotype factor). Investigation time (all objects): P = 0.0230 via two-way ANOVA (genotype × sex; F statistic = 6.325 for genotype factor). (F) Fear conditioning (FC). WT, n = 10; Tg, n = 10; equal males/females. Day 2 (challenge): P = 0.0147 via two-way ANOVA (genotype × sex; F statistic = 7.479 for genotype factor). (G) Social interaction (SI) assay. WT, n = 8 with three females and five males; Tg, n = 10, equal males/females. WT sociability: P = 0.0007 via two-way ANOVA (mouse/empty × sex; F statistic = 20.520 for mouse/empty factor). SLC13A5 nTg sociability: P , 0.0001 via two-way ANOVA (mouse/empty × sex; F statistic = 79.49 for mouse/empty factor). WT recognition: P = 0.0003 via two-way ANOVA (familiar/novel × sex; F statistic = 24.61 for familiar/novel factor). All data are from mice at 4-7 months of age.  imaging to visualize myelinated axons. We specifically examined for differences in myelination of the perforant path axons within the dentate gyrus stratum moleculare. There was a trend for smaller axon diameter in the SLC13A5 nTg, but this did not achieve statistical significance (Fig. 3C). However, by plotting the G-ratio versus axon diameter, we observed significantly more small-diameter axons with higher G-ratios in the SLC13A5 nTg mice, supporting reduced thickness of myelin in the transgenic mice (Fig. 3C). Overall, these data reveal widespread disruptions in white matter integrity in the SLC13A5 nTg mice that are potentially a result of aberrant myelination.

SLC13A5 nTg mice have abnormal synaptic plasticity and spontaneous electrical activity
In addition to alterations in white matter microstructure, we assessed both synaptic structure and function in our SLC13A5 nTg mice, as the AT-1 nTg model displayed changes in synaptic density and biochemical composition. 11 We first performed hippocampal brain slice electrophysiology to assess synaptic plasticity, specifically probing the Schaffer collateral synapse from CA3 onto CA1 pyramidal neurons. We did not observe any difference in LTP (Fig. 3D), but there was a robust increase in LTD in the SLC13A5 nTg mice compared to WT controls (Fig. 3E). Furthermore, we used the Golgi staining method to visualize neuron morphology in vivo, like for the AT-1 nTg, to evaluate for potential explanations of enhanced LTD in the SLC13A5 nTg mice. We did not observe obvious changes in neuron morphology, namely the density of dendritic branching, at low magnification ( Supplementary Fig. 3). We then took high-magnification images of the CA1 pyramidal neurons to quantify dendritic spine density and morphology of their apical dendrites. We did not observe differences in dendritic spine density, but there was an increase in spine volume observed in the SLC13A5 nTg mice compared to WT controls (Fig. 3F). Finally, we evaluated protein expression levels in crude synaptosomes focussing on the proteins that significantly increased in the AT-1 nTg model. 11 To our surprise, we found minimal differences in both the cortex and hippocampus in male and female SLC13A5 nTg mice versus WT controls ( Supplementary  Fig. 4). Therefore, in vivo, our SLC13A5 nTg mice exhibited altered Schaffer collateral-to-CA1 synaptic plasticity with a change in CA1 pyramidal neuron apical dendritic spine morphology. We also uncovered several phenotypic differences from the AT-1 nTg model.
To better characterize neuron morphology and synapse formation, we cultured primary embryonic neurons from WT and SLC13A5 nTg littermates in vitro, performing parallel experiments on both cortical and hippocampal neurons. We expected to see similar changes in neuron morphology like was observed for the AT-1 nTg model, namely increased dendritic branching and spine formation. 11 After 15 DIV, we were surprised to find that SLC13A5 nTg cortical neurons did not exhibit changes from WT in dendritic branching, dendritic spine density or spine morphology (Fig. 4A). In addition to quantifying dendritic spine density, we also assessed spontaneous synapse formation 15 DIV by immunofluorescence staining for the pre-synaptic marker Syn-1 and post-synaptic marker Psd-95. When normalized to neurite volume, we observed a decrease in puncta density for both Syn-1 and Psd-95 as well as a reduction in co-localization between the puncta, suggesting reduced synaptic density in the SLC13A5 nTg cortical neurons compared to WT (Fig. 4B). Finally, to assess the activity of our neurons, we plated them on MEAs and monitored their spontaneous activity every 7 days for a total of 28 DIV. We monitored the spontaneous activity of our neurons by examining the number of active electrodes, mean firing rate and burst frequency; the network activity, implying functional neuron-to-neuron communication, was quantified by network burst frequency and synchronicity index. It is important to note that the activity parameters (mean firing rate, burst frequency, network burst frequency and network synchronicity) required a network to have at least 8 out of 16 active electrodes. At every time point, the SLC13A5 nTg cortical neurons had significantly more active electrodes compared to WT and the mean firing rate was significantly higher at 14 and 21 DIV (Fig. 4C). Note that activity parameter comparisons could not be made at seven DIV since there were no WT networks that exhibited at least eight active electrodes. Moreover, the SLC13A5 nTg cortical networks exhibited a higher burst frequency at 14 and 21 DIV; higher network burst frequency at 14 and 28 DIV and higher synchronicity index at 14 DIV (Fig. 4C). Our hippocampal primary neurons exhibited similar findings to the cortical neurons in morphology (Fig. 5A), but the hippocampal neurons had almost no changes in synaptic density or spontaneous electrical activity ( Fig. 5B and C). Overall, our data show that the SLC13A5 nTg mice exhibit changes in synaptic plasticity and dendritic spine morphology in vivo as well as synaptic density and spontaneous electrical activity in vitro, with notable differences from the previously characterized AT-1 nTg autisticlike model and between cortical and hippocampal-derived neurons.

SLC13A5 overexpression results in differential proteomic and acetyl-proteomic adaptations
To evaluate the molecular underpinning of the SLC13A5 nTg phenotype, we conducted quantitative proteomics on both hippocampal and cortical tissue. The total numbers of proteins detected were 1592 and 1580 in the hippocampus and cortex, respectively; of these proteins detected, 607 and 107 in the hippocampus and cortex, respectively, exhibited an expression level that was significantly different   plots displaying all quantified acetylpeptides in the SLC13A5 nTg versus WT hippocampus (top) and cortex (bottom) across the three subcellular fractions (cytosolic, nuclear/membrane and chromatin-associated). Change in acetylation stoichiometry (Δ acetyl stoichiometry) is calculated as the SLC13A5 nTg stoichiometry value less the WT stoichiometry value. Statistically significant acetylpeptides (P , 0.05) are shown in green with all other acetylpeptides in grey. n = 4 male mice per genotype at 4 months of age. (B) Overlap of the statistically significant acetylpeptides across the three subcellular fractions in both the hippocampus and cortex. (C) Overlap between the hippocampus and cortex of all the statistically significant acetylpeptides (combined from the three subcellular fractions). A heat map of the nine overlapping acetylpeptides between the two tissue types is shown. (D) The fold enrichment of KEGG pathways determined from proteins harbouring the acetylation sites that were significantly changed from WT in the hippocampus and cortex. The top 10 categories sorted by enrichment score are shown that have been filtered by a FDR score of 0.05. (E) Gene-network plots of proteins harbouring the acetylation sites that were significantly changed from WT in the hippocampus and cortex. Plots constructed using an overrepresentation analysis using the gene ontology cellular component function database. The dot size of each network category is scaled by the number of overlapping proteins within the category. The top 10 categories sorted by enrichment score are shown that have been filtered by a FDR score of 0.05. Since proteins could be found in multiple fractions and potentially acetylated on multiple lysine sites, the Δ acetyl stoichiometry value displayed in the network plot is as follows: (i) if the acetylpeptide was detected in multiple subcellular fractions, the fraction selected for illustration was prioritized as follows: cytoplasmic, nuclear/membrane, then chromatin-associated and (ii) if multiple sites per protein exhibited significantly changed acetyl stoichiometry, the lowest number lysine residue was selected for illustration. SLC13A5/NaCT overexpression and autism (P , 0.05) from WT controls (Fig. 6A). Interestingly, the overall fold change in the detected proteome was significantly different between the hippocampus and cortex (Fig. 6B) and only 37 of the significantly changed proteins overlapped between the two tissue types (Fig. 6C). Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of the proteins that exhibited statistically significant changes from WT resulted in multiple overlapping categories related to glucose metabolism, including oxidative phosphorylation, carbon metabolism and citrate cycle [tricarboxylic acid cycle (TCA) cycle] (Fig. 6D). These same proteins were also analysed in the gene ontology (GO) cellular component function database to construct a gene-network plot, which revealed several categories that directly apply to the SLC13A5 nTg phenotype including myelin sheath, synaptic vesicle and neuron-to-neuron synapse (Fig. 6E).
To complement the proteome, we also characterized the acetyl-proteome due to hypothesized changes in protein acetylation from increased transport of citrate into the cytosol. Acetylation can dramatically impact protein function independent of protein expression level, 39 which may play an important role in the SLC13A5 nTg phenotype. As such, we performed subcellular fractionation (cytoplasmic, nuclear/ membrane and chromatin-associated) on hippocampal and cortical tissue in order to increase our detection resolution. 40 In all fractions combined, we detected 2506 and 3886 protein acetylation sites (hereafter referred to as acetylpeptides) in the hippocampus and cortex, respectively. The following number of acetylpeptides exhibited a statistically significant change in acetylation stoichiometry from WT controls (P , 0.05): cytoplasmic, 26 in the hippocampus and 39 in the cortex; nuclear/membrane, 18 in the hippocampus and 45 in the cortex; chromatin-associated, 52 in the hippocampus and 27 in the cortex (Fig. 7A and B). Three acetylpeptides were detected in multiple fractions (Fig. 7B). When the acetylpeptides from all three fractions were combined, nine of these overlapped between the hippocampus and cortex, with many exhibiting opposite changes in acetyl stoichiometry (Fig. 7C). To assess the biological impact of these acetylation changes, a KEGG pathway analysis was conducted on the proteins harbouring the significantly changed acetylation sites in both the hippocampus and cortex. Again, there was a predominance of glucose metabolism-related categories including carbon metabolism and citrate cycle (TCA cycle) (Fig. 7D). Finally, the GO cellular component gene-network plots were strikingly similar between the hippocampus and cortex despite the minimal overlap between the acetylation sites, containing several relevant categories such as myelin sheath, neuron-to-neuron synapse and post-synaptic density (Fig. 7E). Similar to the proteomic analysis, the SLC13A5 nTg acetyl-proteome revealed differential adaptations between the hippocampus and cortex with changes consistent with the autistic-like phenotype. Finally, both the proteome and acetyl-proteome might also highlight a compensatory response from glia.
Global profiling of both the proteome and acetyl-proteome in the SLC13A5 nTg mouse allowed us to assess the overlap between the significantly affected proteins in both analyses. There was a larger overlap between the proteome and acetylproteome in the hippocampus versus the cortex with 29 and 3 commonly shared proteins, respectively (Fig. 8A). To better illustrate this, we constructed a manual cluster plot using the proteins found within the single overlapping GO cellular component myelin sheath category between the proteome and acetyl-proteome (see gene-network plots in Figs 6E and 7E). This plot demonstrates that a majority of the proteins were regulated at either the level of the proteome or the acetyl-proteome, and a handful of proteins exhibited complex changes both in the protein expression level and acetylation stoichiometry, sometimes in multiple lysine sites (Fig. 8B). It is interesting to note that, proportionally, changes in the hippocampus were driven more by changes in protein abundance, whereas changes in the cortex were driven more by changes in acetylation. Finally, we sorted the proteins in both the acetyl-proteome and proteome analyses by their subcellular localization, revealing that 40% of the proteins impacted by SLC13A5 overexpression flux through the secretory pathway (i.e. localize to the ER, Golgi apparatus, membrane or are secreted; Supplementary Fig. 5A).
GO cellular component gene-network plots using these secretory pathway-related proteins revealed categories similar to before, such as myelin sheath, neuron-to-neuron synapse and post-synaptic density, emphasizing the importance of this subset of proteins in the SLC13A5 nTg phenotype (Supplementary Fig. 5B). Taken together, the SLC13A5 nTg hippocampus and cortex displayed differential adaptations in the proteome and acetyl-proteome, with changes in metabolism being a potentially important mechanism underlying the mouse phenotype.

Discussion
Here we reported that overexpression of SLC13A5 restricted to forebrain neurons results in an autistic-like phenotype with jumping stereotypy. The mice displayed disruptions in white matter integrity as well as aberrant synaptic structure and function both in vivo and in vitro. Finally, the proteome and acetyl-proteome uniquely adapted to SLC13A5 overexpression within the hippocampus and cortex, which implicated the potential importance of metabolic changes in the phenotype observed.
Prior study of the autistic-like AT-1 nTg mouse led us to the conclusion that increased flux of acetyl-CoA from the cytosol to the ER may be a mechanistic driver of ASD. 11 This is supported by the observation that gene duplication events in AT-1/SLC33A1, SLC25A1, SLC13A5 and ACLY are all associated with ASD. 24 Importantly, AT-1 overexpressing neurons exhibited cellular adaptations to maintain the cytosolic-to-ER flux of acetyl-CoA by upregulation of SLC25A1 and ACLY, which act in concert to deliver additional citrate to the cytosol for conversion to acetyl-CoA. 11 The results presented here affirm that SLC13A5 overexpression also results in an autistic-like phenotype, likely from increased citrate/acetyl-CoA availability in the cytosol and flux of acetyl-CoA into the ER. Indeed, 40% of the proteins affected in the SLC13A5 nTg acetyl-proteome and proteome localized to the secretory pathway and were relevant to the phenotypic features of the mouse model. Targeting the ER acetylation machinery, such as by inhibition of the N ɛ -lysine acetyltransferases ATase1 and ATase2, 2,10,12,16 is an active area of investigation to assess the direct contribution of changes in ER-based acetylation in the development of the autistic-like phenotype in the SLC13A5 nTg mouse.
Alterations in membrane physiology or organelle adaptations outside of the ER as a result of SLC13A5 overexpression are also likely to play a role in the phenotype that we observed. First, SLC13A5 functions as a sodium cotransporter, carrying in four Na + for every dicarboxylate or tricarboxylate molecule imported into the cell, 25 which results in membrane depolarization. This impact on membrane potential may have substantial implications for neuron physiology-especially when SLC13A5 is expressed at high levels-that likely triggers changes within the neuron to account for this continuous source of membrane depolarization. Our proteomic dataset supports this as we saw changes in expression levels of Na + /K + -ATPase subunits, namely Atp1a1, Atp1a2, Atp1a3, Atp1b1 and Atp1b2, as well as other proteins involved in membrane excitability including Slc8a2, Scn2a and Scn2b. It is also important to point out that human SLC13A5 is low-affinity/high-capacity transporter for citrate, which is in contrast to mouse Slc13a5 that is a high-affinity/low-capacity transporter 41,42 ; this results in significantly more citrate import into the cell over the same amount of time for the human transporter compared to the mouse transporter. Furthermore, within the neuron, citrate participates in the synthesis of the neurotransmitters acetylcholine, glutamate and gammaaminobutyric acid (GABA), 42 and the altered availability of citrate within the cytosol may impact one or all of these neurotransmitter systems. Interestingly, we did see expression level changes in several genes related to neurotransmission including the glutamate receptors Grm2, Grm5 and Grin2a; GABA receptor subunits Gabra1 and Gabbr1; the diazepam binding inhibitor Dbi; and vesicle-associated proteins Snap25 and Sygr3.
The critical role of human SLC13A5 in neuron physiology is emphasized by loss-of-function mutations associated with early infantile epileptic encephalopathy-25, 43,44 which may in part be due to imbalances in GABA synthesis. Likewise, deletion of mouse Slc13a5 results in a propensity for epileptic seizures and quantitative proteomic analysis revealed disruptions in GABA and serotonin synthesis as well as adaptations in oxidative phosphorylation and lipid metabolism. 45 Our proteome and acetyl-proteome analyses strongly suggest changes in glucose metabolism within the SLC13A5-overexpressing neurons, which are not surprising given the central role of citrate and acetyl-CoA in cellular energetics as well as the aforementioned mouse Slc13a5 knockout study. We propose that the neuron is likely accommodating for membrane depolarization by excessive Na + /citrate co-transport, which would require additional ATP to run the Na + /K + -ATPase. It would not be surprising if the SLC13A5 nTg mouse exhibits an increased capacity for oxidative phosphorylation to generate additional ATP. There may be uncharacterized adaptations in the SLC13A5 nTg mitochondria that are a result of altered histone acetylation and methylation signatures, which was observed in the AT-1 nTg model. 11 Dissecting the phenotype of the SLC13A5 nTg mouse revealed unexpected but intriguing differences from the AT-1 nTg mouse. First, the SLC13A5 nTg mice exhibited a jumping stereotypy while in their home cage, which was not observed in the AT-1 nTg model but is a feature of other rodent ASD models including the C58/J inbred mouse strain and Jakmip1 knockout mouse. 11,46,47 Second, unlike the AT-1 nTg, the SLC13A5 nTg mice displayed reduced locomotion as demonstrated by the OF and LD assays, but oddly this behaviour was not observed in the NOR or FC assays as demonstrated by overall increased time exploring objects and similar freezing behaviour on the training day, respectively. It appears that the SLC13A5 nTg behaviour is variable depending on the context, such as repetitively jumping when undisturbed in the home cage versus overly engaged in exploration during the NOR assay. It is also important to mention that the cognitive deficits we observed in the FC and SI assays are in contrast to improved cognitive function observed in a mouse Slc13a5 knockout model, 48 suggesting the importance of balanced Slc13a5 function in learning and memory function. We have not excluded the possibility that the unusual behaviour observed in the SLC13A5 nTg model could fit into other categories of neurodevelopmental diseases, such as attention deficit/hyperactivity disorder, which has considerable overlap with ASD. [49][50][51] Further behavioural testing, such as motor impulsiveness, could better categorize the behavioural phenotype of the SLC13A5 nTg mouse. 52 Finally, SLC13A5 nTg mice did not exhibit changes in hippocampal LTP, dendritic branching and spine formation or synaptosome protein expression like the AT-1 nTg model. 11 It is important to point out that we used crude synaptosomes unlike in the AT-1 nTg study where purified post-synaptic densities were used for Western blotting; therefore, our synaptosomes may have not been pure enough to detect subtle changes in protein expression. The SLC13A5 nTg mice did, however, have increased dendritic spine volume in vivo. Therefore, even with the similar behavioural abnormalities between the two autistic-like models, the underlying cellular and molecular changes are rather distinct. This likely highlights the importance of other mechanisms beyond ER-based N ɛ -lysine acetylation in the SLC13A5 nTg mice that are contributing to the development of the phenotype, such as the differential effect of the AT-1 and SLC13A5 manipulations on the cytosolic pools of citrate and acetyl-CoA, and consequent functional responses.
Via ex vivo DWI and TEM, we found disrupted white matter integrity in the SLC13A5 nTg mice with evident changes in myelination in the hippocampus. What remains to be elucidated is whether these changes are due to an exclusive problem in the neurons or a defect in neuron-tooligodendrocyte communication. Our proteome and acetylproteome suggest that SLC13A5 overexpression results in widespread changes in both protein expression level and acetylation status of multiple genes involved in myelination, but the specific cell type exhibiting those changes was not resolved. Alterations in neuronal secretory proteins or metabolites like N-acetylaspartate could impact oligodendrocyte function at multiple points during neurodevelopment, 53,54 and more detailed exploration could reveal a new function of ER-based N ɛ -lysine acetylation in this neuron-to-glia communication.
Both our primary neuron and proteome/acetyl-proteome analyses revealed unique adaptations in the hippocampus versus cortex to SLC13A5 overexpression. This was evident by changes observed in synaptic density and spontaneous electrical activity in primary neurons derived from the cortex but not hippocampus, as well as different sets of changes in protein expression and acetyl stoichiometry between the cortex and hippocampus. It is important to note that our hippocampal primary neurons are more active with networks maturing more quickly compared to our cortical neurons, potentially from culture conditions that were originally optimized for hippocampal primary neuron cultures. Additionally, the population of neuron subtypes within the hippocampal and cortical cultures may be relevant, especially if membrane potential and neurotransmitter synthesis are impacted by SLC13A5 overexpression.
In conclusion, by overexpressing SLC13A5 in mouse forebrain neurons, we have reinforced the conclusion that aberrant citrate/acetyl-CoA flux within the neuron is a mechanistic driver of an autistic-like phenotype, resulting in disrupted white matter integrity, changes in synaptic structure and plasticity and widespread adaptations in the proteome and acetyl-proteome. The metabolic versus ER-based N ɛ -lysine acetylation contributions to the SLC13A5 nTg phenotype remain to be elucidated.

Supplementary material
Supplementary material is available at Brain Communications online.