# The Metabolomic Signature of Opa1 Deficiency in Rat Primary Cortical Neurons Shows Aspartate/Glutamate Depletion and Phospholipids Remodeling

Juan Manuel Chao de La Barca, Macarena Arrázola, Cinzia Bocca, Laetitia Arnauné-Pelloquin, Olga Iuliano, Guillaume Tcherkez, Guy Lenaers, Gilles Simard, Pascale Belenguer, Pascal Reynier

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Juan Manuel Chao de La Barca, Macarena Arrázola, Cinzia Bocca, Laetitia Arnauné-Pelloquin, Olga Iuliano, et al. The Metabolomic Signature of Opa1 Deficiency in Rat Primary Cortical Neurons Shows Aspartate/Glutamate Depletion and Phospholipids Remodeling. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-42554-7. http://02388218

HAL Id: hal-02388218

https://univ-angers.hal.science/hal-02388218

Submitted on 5 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# The Metabolomic Signature of Opa1 Deficiency in Rat Primary Cortical Neurons Shows Aspartate/ Glutamate Depletion and Phospholipids Remodeling

Juan Manuel Chao de la Barca<fnref n="1" />,<fnref n="2" />, Macarena S. Arrázola³,⁴, Cinzia Bocca<fnref n="1" />, Laetitia Arnauné-Pelloquin³, Olga Iuliano³, Guillaume Tcherkez<fnref n="5" />, Guy Lenaers<fnref n="1" />, Gilles Simard<fnref n="2" />, Pascale Belenguer³ & Pascal Reynier<fnref n="1" />,<fnref n="2" />

Pathogenic variants of *OPA1*, which encodes a dynamin GTPase involved in mitochondrial fusion, are responsible for a spectrum of neurological disorders sharing optic nerve atrophy and visual impairment. To gain insight on OPA1 neuronal specificity, we performed targeted metabolomics on rat cortical neurons with OPA1 expression inhibited by RNA interference. Of the 103 metabolites accurately measured, univariate analysis including the Benjamini-Hochberg correction revealed 6 significantly different metabolites in OPA1 down-regulated neurons, with aspartate being the most significant (*p* < 0.001). Supervised multivariate analysis by OPLS-DA yielded a model with good predictive capability (Q<fnref n="2" /> cum = 0.65) and a low risk of over-fitting (permQ2 = −0.16, CV-ANOVA *p*-value 0.036). Amongst the 46 metabolites contributing the most to the metabolic signature were aspartate, glutamate and threonine, which all decreased in OPA1 down-regulated neurons, and lysine, 4 sphingomyelins, 4 lysophosphatidylcholines and 32 phosphatidylcholines which were increased. The phospholipid signature may reflect intracellular membrane remodeling due to loss of mitochondrial fusion and/or lipid droplet accumulation. Aspartate and glutamate deficiency, also found in the plasma of OPA1 patients, is likely the consequence of respiratory chain deficiency, whereas the glutamate decrease could contribute to the synaptic dysfunction that we previously identified in this model.

Received: 19 April 2018

Accepted: 26 March 2019

Published online: 15 April 2019

OPA1 is a dynamin GTPase controlling the mitochondrial network and cristae organization through its role in the inner mitochondrial membrane structuration<fnref n="1" />. More than 300 *OPA1* gene variants have been reported since its first description as the main contributor to dominant optic atrophies (DOA, MIM#165500)<fnref n="2" />⁻⁴. Clinical studies have progressively shown that OPA1 dysfunction was responsible for a large spectrum of neurological disorders involving optic, auditory and peripheral nerves as well as the brain<fnref n="5" />,⁶. This has highlighted the importance of OPA1, and more generally of mitochondrial dynamics, in neuronal plasticity⁷.

OPA1 inactivation affects many functions, such as mitochondrial fusion, mitochondrial cristae organization and cell death⁸,⁹, oxidative phosphorylation (OXPHOS) and energy production¹⁰, mitochondrial DNA maintenance¹¹⁻¹³, calcium fluxes¹⁴, reactive oxygen species production¹⁵,¹⁶, inflammation¹⁷, ageing¹⁷,¹⁸, mitophagy and mitochondrial renewal¹⁸⁻²⁰.

To explore further the neuronal specificity characterizing OPA1 disorders, we have developed a rat primary neuron culture model, with a siRNA-induced OPA1 haploinsufficiency²¹, mimicking the common pathological

<footnote n="1">¹Equipe Mitolab, Institut MITOVASC, CNRS 6015, INSERM U1083, Université d'Angers, Angers, France.</footnote>

<footnote n="2">²Département de Biochimie et Génétique, Centre Hospitalier Universitaire, Angers, France. ³Centre de Recherches sur la Cognition animale (CRCA), Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, Toulouse, France. ⁴Present address: Center for Integrative Biology, Facultad de Ciencias, Universidad Mayor, Santiago, Chile.</footnote>

<footnote n="5">⁵Research School of Biology, The Australian National University, Canberra, 2601, ACT, Australia. Pascale Belenguer and Pascal Reynier contributed equally. Correspondence and requests for materials should be addressed to J.M.C.d.l.B. (email: JMChaoDeLaBarca@chu-angers.fr)</footnote>
mechanism of DOA. This model allowed us disclosing dendritic and synaptic impairments associated with altered mitochondrial respiration and reduced abundance of mitochondria along dendrites¹⁶,²¹. Accordingly, synaptic defects and dendropathy were also found in vivo in retinal ganglion cells (RGC) from a DOA mouse model²²,²³. Furthermore, OPA1 loss induces pro-oxidative cellular conditions, leading to increased cell death¹⁶ and reduced mitophagy²⁴ in rat cortical neurons.

Metabolomics is particularly useful to explore neuro-metabolic disorders. We recently used a non-targeted metabolomics approach on blood samples from OPA1 patients and showed a specific metabolic signature that includes the alteration of purine metabolism as well as that of amino acids and fatty acids²⁵. We also applied targeted metabolomics on 9 tissues of 3 month old *Opa1delTTAG/+* mice¹⁸, revealing a pre-symptomatic discriminant metabolomics signature in optic nerves, while other tissues were not affected at this age, thereby demonstrating the high sensitivity of optic nerve metabolism to OPA1 dysfunction²⁶. In the later study, the metabolomics signature was characterized by a decrease in sphingomyelins and lysophosphatidylcholines, and an alteration of metabolites involved in neuro-protection and neuronal metabolism, suggesting myelin sheath alterations and axonal dysfunction, respectively.

In order to explore the neuronal specificity of the OPA1 clinical expression, we use here targeted metabolomics to characterize the metabolism of cultured rat cortical neurons in which OPA1 expression was inhibited by siRNA interference. In contrast to skin fibroblasts isolated from patients with similar OPA1 haploinsufficiency, the cortical neurons were considerably affected.

## Material and Methods

**Cell culture of rat embryonic cortical neurons.** All animal procedures were approved by the "Animal Experimentation Ethics Committee" of the *Centre National de la Recherche Scientifique/Fédération de Recherche de Biologie de Toulouse* (C2EA-01) under the protocol number 01024-01.7. All experiments were performed in accordance with the relevant guidelines and regulations. Pregnant Wistar rats were delivered 48 h before sacrifice and kept with food and water *ad libitum*. Embryos were removed at Day 17 (D17) from pentobarbital anaesthetized females (Ceva Santé Animale, Libourne, France) and cortices were dissected as described²¹. Cortical cells (5 10⁶) obtained from D17 embryonic cortices were electroporated using the Rat Neuron Nucleofector® Kit (Amaxa®, Lonza, Basel, Switzerland) with either control luciferase-targeting (D-001210-02, Dharmacon, Lafayette, Colorado, USA) or OPA1-targeting (D-099769-01, Dharmacon) small interfering RNA (3 µg), leading to siLUC (controls) and siOPA1 cells. Transfected cells were plated on poly-D-lysine coated Petri dishes (1.25 10⁶ cells/35 mm) grown in Neurobasal A-25 medium supplemented with B27, L-glutamine, penicillin/streptomycin, amphotericin and lactic acid as described²¹. Cell pellets (3 10⁶) obtained after 9 days of culture were conserved at -80 °C until Western-blot and metabolomics analyses were performed.

**Western-Immuno-blotting.** Proteins were extracted from primary cultured cortical neurons as described¹⁶. Protein extracts (10 µg) were subjected to SDS-PAGE and blotted onto nitrocellulose membranes, which were treated first with anti-OPA1 (1/500, BD Biosciences, Le Pont de Claix, France) and anti-GAPDH (1/500, Merck Millipore MAB374, Molsheim, France) primary antibodies, and then with horseradish peroxidase-conjugated secondary antibodies (1/5000, Abcam ab6789 and ab6721, Paris, France) as described²¹. After enhanced chemiluminescent detection, signals were analyzed using ChemiDoc™ MP Imaging System and quantified by Image LabTM Software (Biorad®, Marnes la Coquette, France).

**Metabolite extraction.** Each cell pellet (3 10⁶ cells, n = 10 for each siOPA1 and siLUC pairs) was dissolved in 100 µl of cold methanol/PBS (85:15, v/v). The mixture was transferred to a 0.5 ml homogenizer tube prefilled with ceramic beads. Cell lysis was achieved in a Precellys® 24 homogenizer (Bertin instruments, Montigny-Le-Bretonneux, France) by two cycles of grinding (2 × 20 s at 6,500 rpm with a break of 20 s, followed by 30 s at 6,000 rpm) at 4 °C. The homogenate was centrifuged at 10,000 g for 5 min at 4 °C, the supernatant was then evaporated, and the resulting pellet stored at -80 °C until further analysis.

**Metabolomics analysis.** Targeted quantitative metabolomics analyses were carried out as described²⁷, using the Biocrates® Absolute IDQ p180 kit (Biocrates Life sciences AG, Innsbruck, Austria). This kit uses mass spectrometry (QTRAP 5500, SCIEX, Villebon-sur-Yvette, France) to quantify up to 188 different endogenous molecules (Supplementary Table S1). Flow injection analysis coupled with tandem mass spectrometry (FIA-MS/MS) was used to analyze carnitine, acylcarnitines, lipids and hexoses, whereas liquid chromatography (LC) was used to separate amino acids and biogenic amines before mass spectrometry quantitation.

All reagents used in this analysis were of LC-MS grade and purchased from VWR (Fontenay-sous-Bois, France) and Merck (Molsheim, France). Sample preparation and analysis were performed following the Kit User Manual. Briefly, after thawing on ice, the pellet was dissolved in 40 µL methanol. After vigorous mixing, 10 µL of each sample were added to a well of a 96-well plate. Metabolites were extracted and derivatized for the quantitation of amino acids and biogenic amines. Extracts were finally diluted with MS running solvent before FIA and LC-MS/MS analysis. Three quality controls (QCs) composed of three concentrations of a human plasma samples: low (QC1), medium (QC2) and high (QC3), were used to evaluate the performance of the analytical assay. A seven-point serial dilution of calibrators was added to the plate of the kit to generate calibration curves for the quantification of amino acids and biogenic amines.

**Statistical analyses.** Raw data were examined before statistical analyses to exclude metabolites having more than 20% of their concentration values below the lower limit of quantitation (LLOQ) or above the upper limit of quantitation (ULOQ). To account for differences in metabolite concentration amongst samples merely due to differences in the number of cells used, each metabolite concentration belonging to a particular sample was divided by the sum of all the metabolites in that sample (row normalization).
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Box plot showing the OPA1 protein content relative to GAPDH protein in neurons transfected with siLUC or siOPA1, and cultured 9 days *in vitro*. Data represent the mean ± SEM of 10 independent experiments, and were statistically treated with Mann-Whitney **p** < 0.01. AU = Arbitrary Units. One representative immunoblot is shown on the right panel.</figcaption></figure>

Univariate analysis was performed using non-parametric Mann-Whitney-Wilcoxon test (thereafter simply referred to as 'Wilcoxon test') for quantitative variables and the Fischer exact test for qualitative variables. The Benjamini-Hochberg correction was used to compare multiple metabolite concentrations to keep the false discovery rate below 5%. Otherwise differences were considered significant when $p \le 0.05$.

Multivariate analysis was performed using SIMCA-P v.14.1 (Umetrics, Umeå, Sweden). Unsupervised principal component analysis (PCA) was employed to detect strong outliers and visualize sample clusters. A supervised method called orthogonal projection to latent structures by means of partial least square-discriminant analysis (OPLS-DA) was used to find the linear combination of metabolite concentration in the X (matrix of metabolites) space that best correlates to the Y response vector (binary classifier: siOPA1/siLUC). The quality of models was appreciated using two parameters: $R^2$ (goodness of fit) and $Q^{2}_{cum}$ (goodness of prediction). In the model with the best predictive capabilities (i.e., $Q^{2}_{cum} > 0.5$), the risk of overfitting was measured by the predictive capabilities ($Q^{2}_{perm}$) of the model obtained when elements of the response vector Y were randomly permuted (permutation test) (forming a new Y vector, $Y_{perm}$). A non-over-fitted model is characterized by poor predictive capabilities (i.e., negative $Q^{2}_{cum-perm}$) when the permutation test yields models with no correlation between Y and $Y_{perm}$. In the model retained, variables (here, metabolites) were selected according to their variable importance for the projection (VIP) and their loading rescaled as a correlation coefficient ($p_{cor}$) between the original variable in the X matrix and the predictive component of interest. For a given metabolite $m$ and a unique predictive component $tp$, VIP$_m$ summarizes its importance for the OPLS-DA model, whereas $p_{cor,m}$ indicates the magnitude of the correlation between $m$ and $tp$. Thus, important variables (i.e., VIP > 1) having a high absolute $p_{cor}$ value were retained since they were considered to be critical for group discrimination in predictive models. Plotting the VIP value versus $p_{cor}$ for all variables (the “volcano” plot) enables visualization of the importance of each variable for group discrimination in the OPLS model.

In addition to the information obtained from individual metabolites, 25 sums and ratios of metabolites with biological significance have been calculated according to the RatioExplorer module provided by Biocrates® (Supplementary Table S2).

## Results

**OPA1 inactivation.** To gain information on the impact of OPA1 inactivation on the metabolome of neuronal cells, we used embryonic rat cortical neurons grown in primary culture transfected with small interfering RNA against OPA1 (siOPA1; n = 10) or control small interfering RNA (siLUC; n = 10). After 9 days of culture *in vitro*, a Western-blot analysis using anti-OPA1 antibodies revealed that OPA1 expression, normalized to GAPDH, was reduced by 41.3% on average in siOPA1 compared with siLUC ($p < 0.0029$) (Fig. 1 and Supplementary Fig. S3).

**Univariate analysis.** After validation of the kit plate based on QC samples, 103 (54.8%) metabolites were found to be accurately measured and consequently used for statistical analyses (Supplementary Table S4).

Univariate analysis using the Wilcoxon test showed that 36 metabolites were significantly different ($p < 0.05$) in siOPA1 compared to siLUC control neurons, while only 6 of them (decrease of aspartate and increase of lysoPC 20:3, lysoPC 20:4, SM 24:1, SM (OH) 14:1 and PC ae 36:4) remained significantly different after having applied the Benjamini-Hochberg correction with an overall false discovery rate (FDR) of 0.046 (Fig. 2).

**Multivariate analysis.** Multivariate analysis by unsupervised principal component analysis (PCA) showed neither spontaneous grouping nor any strong outliers in the first principal plot (PC1 vs. PC2) (Fig. 3A). By contrast, supervised multivariate analysis by OPLS-DA showed a clear differentiation between siOPA1 and siLUC neurons, with high predictive capabilities ($Q^{2}_{cumulated} = 0.65$) and good performances in the permutation test ($R^2X = 0.86$ and $R^2Y = 0.95$; $Q^2: (-0.16)$), and the CV-ANOVA tests ($p$-value: 0.036) (Fig. 3B).

The representation of VIP (variable importance for the projection) values against correlation coefficient ($p_{corr}$) (volcano plot, Fig. 4) showed 46 important metabolites with a VIP value $\ge 1$ and as such, considered as being crucial for group discrimination. Amongst these, aspartate, glutamate and threonine were decreased in siOPA1 compared to siLUC neurons, whereas 43 metabolites were increased. Lysine was the only increased amino acid
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Box plot showing the 6 metabolites significantly different in siOPA1 (n = 10) compared to siLUC controls (n = 10) neurons using the univariate Wilcoxon test after Benjamini-Hochberg correction. Values in the y-axis have no dimension as they represent relative concentrations. ***p < 0.001; **p < 0.01. PC ae: Alkyl-acyl phosphatidylcholine; lysoPC: lysophosphatidylcholine; SM: sphingomyelin; SM(OH): hydroxy-sphingomyelin.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** PCA (A) and OPLS-DA (B) scatter plots obtained from the matrix of metabolites for the 10 samples from siLUC (blue circles) and the 10 samples from siOPA1 (green circles). (A) PCA shows neither clear grouping nor outlier in the first principal plan, the green point appearing outside the ellipse being not a strong outlier according to Hotelling's T2 range. (B) There is a clear between-group discrimination in the OPLS-DA plot along the predictive latent variable (p LV). Legend: PC1,2: Principal Components 1 and 2; o LV: first orthogonal latent variable; p LV: predictive latent variable.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Volcano plot ($p_{corr}$ vs. VIP) from the OPLS-DA model. Only the most discriminating metabolites having high VIP values ≥ 1 (indicated by the horizontal red line) have been labelled. Negative $p_{corr}$ values (left) indicate diminished metabolite concentrations in siOPA1 neurons versus siLuc neurons, whereas positive $p_{corr}$ values (right) indicate increased metabolite concentrations in siOPA1 neurons compared to the control group. The metabolomic signature accompanying the loss of OPA1 expression is associated to decreased levels of the amino acids (green bubbles) aspartate (Asp), glutamate (Glu) and threonine (Thr) and to increased levels of 34 phosphatidylcholines (PC, orange bubbles), 4 lysophosphatidylcholines (lysoPC, blue bubbles), 4 sphingomyelins (SM and SM(OH)), yellow bubbles) and the amino acid lysine (Lys, green bubble). Red-rimmed bubbles indicate metabolites significantly different between siOPA1 and siLUC neurons in univariate analysis (Wilcoxon test) after Benjamini-Hochberg correction.</figcaption></figure>

<table><caption>**Table 1.** Metabolite sums and ratios calculated with amounts accurately determined and listed according to their P-value (from the smallest to the largest). In order to determine significance of each sum or ratio, observed p-values should be lower than the corrected alpha threshold. Legend: SFA: saturated fatty acids; MUFA: mono-unsaturated fatty acids; PUFA: poly-unsaturated fatty acids; UFA: unsaturated fatty acid, equal to the sum of MUFA and PUFA; Methionine-SO: sulfoxidized methionine; PC: phosphatidylcholine; aa: diacyl; ae: acyl-alkyl.</caption><thead><tr><th>Ratio or sum</th><th>Fold Change (SiOPA1/ SiLUC)</th><th>Observed p-value</th><th>Corrected alpha threshold</th></tr></thead><tbody><tr><td>UFA/SFA aa</td><td>0.68</td><td><strong>0.00021</strong></td><td><strong>0.0089</strong></td></tr><tr><td>UFA/SFA ae</td><td>0.83</td><td><strong>0.00209</strong></td><td><strong>0.0132</strong></td></tr><tr><td>SFA ae</td><td>1.66</td><td>0.02323</td><td>0.0170</td></tr><tr><td>SFA aa</td><td>1.66</td><td>0.02881</td><td>0.0204</td></tr><tr><td>Putrescine/Spermidine</td><td>1.23</td><td>0.05381</td><td>0.0236</td></tr><tr><td>Spermine/Spermidine</td><td>1.51</td><td>0.15641</td><td>0.0266</td></tr><tr><td>Methionine-SO/Methionine</td><td>1.14</td><td>0.19032</td><td>0.0294</td></tr><tr><td>Citrulline/Arginine</td><td>0.27</td><td>0.19280</td><td>0.0321</td></tr><tr><td>MUFA ae</td><td>1.27</td><td>0.24745</td><td>0.0348</td></tr><tr><td>PUFA ae</td><td>1.20</td><td>0.24745</td><td>0.0374</td></tr><tr><td>Total lysoPC/ Total PC</td><td>0.75</td><td>0.24745</td><td>0.0399</td></tr><tr><td>MUFA aa</td><td>1.17</td><td>0.31499</td><td>0.0423</td></tr><tr><td>PUFA aa</td><td>1.03</td><td>0.63053</td><td>0.0447</td></tr><tr><td>Tyrosine/Phenylalanine</td><td>0.97</td><td>0.63053</td><td>0.0470</td></tr></tbody></table>

while 4 sphingomyelins, 4 lysophosphatidylcholines and 34 phosphatidylcholines were all increased in siOPA1 neurons. Importantly, the 6 significant metabolites identified by the univariate analysis were the most important contributors of the 46 discriminating metabolites identified in the multivariate model.

Fourteen metabolites sums and/or ratios were calculated with available metabolite contents. Only unsaturated to saturated diacyl- and acyl-alkyl phosphatidylcholine (PC aa and PC ae, respectively) were significantly different between siOPA1 and siLUC neurons after Benjamini-Hochberg correction, with lower values for both ratios in siOPA1 neurons (Table 1).

Discussion

Our results show that the metabolome of embryonic rat cortical neurons is significantly impacted by siRNA-mediated decreased OPA1 expression as compared to controls. Taken as a whole, the OPA1 metabolomics signature is characterized by a decrease in aspartate, glutamate and threonine, and an increase in lysine and 42 glycerophospholipids.
The decrease in aspartate concentration was the most significant alteration identified by both univariate (*p* < 0.001) and multivariate (VIP > 1.9) analyses. This important metabolic feature associated to OPA1 silencing has also been found in our recent study of blood from OPA1 patients²⁵ and in *Opa1*⁻/⁻ mouse embryonic fibroblasts²⁷, thus in three samples from three different species. In contrast, the patient fibroblasts carrying heterozygous OPA1 pathogenic variants did not revealed a discriminant metabolomic signature²⁸. Glutamate, the metabolism of which is tightly coupled to aspartate metabolism via the Krebs cycle and aminotransferases, was also reduced in siOPA1 neurons. This feature was also observed in patients and mouse embryonic fibroblasts. In the mouse model of OPA1 haploinsufficiency, aspartate concentration was invariant while glutamate concentration was respectively lower and higher in the optic nerve of symptomatic 11 month-old mice or asymptomatic 3 month-old mice²⁶.

Two mechanisms could explain the alteration of aspartate homeostasis upon OPA1 deficiency. First, it has been reported that OPA1 directly interacts with the two mitochondrial aspartate/glutamate carriers (SLC25A12 and SLC25A13) involved in the exchange of aspartate and glutamate across the inner mitochondrial membrane²⁹. Therefore, the decrease in OPA1 expression may alter the cytosolic-mitochondrial partitioning of glutamate and aspartate, and consequently their homeostasis. Second, aspartate deficiency has been shown to be a consequence of mitochondrial respiratory inhibition and impaired oxidative phosphorylation³⁰,³¹. This would agree with the reduced mitochondrial membrane potential and the reduced expression of OXPHOS proteins that were previously reported in siOPA1 rat cortical neurons¹⁶,²¹ and more generally in OPA1 mutated conditions¹⁰. It should be noted that these two mechanisms are not mutually exclusive, because mitochondrial solute carriers use the inner membrane potential provided by OXPHOS as an energy supply for amino acids transport across this membrane.

Since aspartate is used as an amino-nitrogen donor during *de novo* glutamate synthesis in glial cells, the decrease in aspartate could lead to depletion in glutamate which in turn may impair synaptic transmission. Glutamate is indeed considered to be the main neurotransmitter of most excitatory synapses in the brain, through its interaction with the NMDA receptor (N-methyl-D-aspartate receptor). Interestingly, we have shown that synaptic plasticity is altered in siOPA1 neurons²¹ and others have shown that OPA1 expression influences NMDA receptor expression in the retina of a DOA mouse model³². Thus, the altered synaptogenesis plasticity related to OPA1 may be due to decreased NMDA receptor excitation as a consequence of aspartate/glutamate reduced concentrations. In turn, the lysine being a glutamate precursor through the saccharopine pathway in the mammalian central nervous system³³, the impairment of glutamate biosynthesis may be responsible for the accumulation of lysine observed in the signature.

The increase in four sphingomyelins, four lysophosphatidylcholines and 34 phosphatidylcholines shows that OPA1 deficiency impacts considerably on phospholipidic composition of cortical neurons. These phospholipids represented a large proportion (42/105 = 39%) of the phospholipid species measured in the present targeted metabolomic analysis. Validating our results, such increased phosphatidylcholines was also confirmed in the same model using another mass spectrometer based lipidomic approach (data not shown). We previously reported that 12 phosphatidylcholines and one lysophosphatidylcholines also increased in the optic nerves of 11-month old *Opa1*⁻/⁻ female mice²⁶. Also, in siOPA1 neurons identical to those used here, we previously showed that OPA1 silencing led to an alteration of the mitochondrial network that became highly fragmented²¹. The reconfiguration of phospholipid composition observed here may therefore reflect this fragmentation. In fact, the involvement of OPA1 in membrane fusion is mediated by a direct interaction with phospholipids in the inner mitochondrial membrane. Accordingly in yeast, Mgm1p, the OPA1 ortholog, promotes changes in phospholipid composition in liposomes to trigger membrane fusion³⁴.

The phospholipid remodeling may also be linked to lipid droplets accumulating when mitochondrial fusion or OXPHOS activity is inhibited. Both the levels of triacylglycerols and the number and volume of lipid droplets increase upon OPA1 down-regulation in adipocytes and mouse embryonic fibroblasts²⁵,³⁶. Lipid droplets accumulation has also been found in the skeletal muscle of a mouse model with a conditional *Opa1*⁺/⁻ knock-out¹⁷. In addition, mitochondrial fusion is required in starved cells to redistribute lipid droplets-derived fatty acids across the mitochondrial network, to ensure optimal oxidation, and OPA1 dysfunction redirects fatty acids back to lipid droplets in both starved and non-starved cells³⁶. In addition, it has been shown that a fraction of OPA1 protein ectopically localizes at the surface of lipid droplets in adipocytes, playing a direct role in the adrenergic control of lipolysis by anchoring the protein kinase A at the surface of lipid droplets³⁷. Thus, the phospholipid signature we found here in siOPA1 neurons may simply result from lipid droplets that could accumulate, impacting fatty acids redistribution³⁶ and/or lipolysis³⁷.

The unsaturated to saturated diacyl- and acyl-alkyl phosphatidylcholine ratios also significantly decreased in siOPA1 compared to siLUC neurons, suggesting a general increase in the saturation of fatty acid chains of phosphatidylcholines and thus a decrease in the desaturases activities. This effect is likely related to oxidative stress that occurs in siOPA1 neurons¹⁶. In fact, both activity and expression of desaturases are regulated by the intracellular redox state³⁸ and a low degree of fatty acid unsaturation may be a mechanism to protect cellular membranes from lipid peroxidation³⁹.

Taken as a whole, our study shows that OPA1 deficiency in rat cortical neurons defines a highly specific metabolomic signature. It includes a general reconfiguration of glycerophospholipids, likely reflecting the alteration of mitochondrial fusion and the accumulation of lipid droplets. The second feature of the OPA1 signature in neurons was a sharp reduction in aspartate and glutamate concentrations, also found in other models of OPA1 dysfunction. This aspartate/glutamate signature is an attractive candidate to explain the neuronal specificity of the OPA1 clinical expression, through the perturbation of the NMDA receptor signaling.
## References

1. Olichon, A. et al. Loss of OPA1 perturbates the mitochondrial inner membrane structure and integrity, leading to cytochrome c release and apoptosis. *J. Biol. Chem.* **278**, 7743–7746 (2003).

2. Delettre, C. et al. Nuclear gene Opa1, encoding a mitochondrial dynamin-related protein, is mutated in dominant optic atrophy. *Nat. Genet.* **26**, 207–210 (2000).

3. Alexander, C. et al. Opa1, encoding a dynamin-related GTPase, is mutated in autosomal dominant optic atrophy linked to chromosome 3q28. *Nat. Genet.* **26**, 211–215 (2000).

4. Ferre, M. et al. Improved locus-specific database for Opa1 mutations allows inclusion of advanced clinical data. *Hum. Mutat.* **36**, 20–25 (2015).

5. Chao de la Barca, J. M. et al. Opa1-related disorders: Diversity of clinical expression, modes of inheritance and pathophysiology. *Neurobiol. Dis.* **9**, 20–26 (2016).

6. Yu-Wai-Man, P. et al. A neurodegenerative perspective on mitochondrial optic neuropathies. *Acta Neuropathol.* **132**, 789–806 (2016).

7. Bertholet, A. M. et al. Mitochondrial fusion/fission dynamics in neurodegeneration and neuronal plasticity. *Neurobiol. Dis.* **9**, 3–19 (2016).

8. Frezza, C. et al. Opa1 controls apoptotic cristae remodeling independently from mitochondrial fusion. *Cell* **126**, 177–189 (2006).

9. Olichon, A. et al. Opa1 alternate splicing uncouples an evolutionary conserved function in mitochondrial fusion from a vertebrate restricted function in apoptosis. *Cell Death Differ.* **14**, 682–692 (2007).

10. Chevrollier, A. et al. Hereditary optic neuropathies share a common mitochondrial coupling defect. *Ann. Neurol.* **63**, 794–798 (2008).

11. Amati-Bonneau, P. et al. Opa1 mutations induce mitochondrial DNA instability and optic atrophy 'plus' phenotypes. *Brain* **131**, 338–351 (2008).

12. Hudson, G. et al. Mutation of Opa1 causes dominant optic atrophy with external ophthalmoplegia, ataxia, deafness and multiple mitochondrial DNA deletions: a novel disorder of mtDNA maintenance. *Brain* **131**, 329–337 (2008).

13. Elachouri, G. et al. OPA1 links human mitochondrial genome maintenance to mtDNA replication and distribution. *Genome Res.* **21**, 12–20 (2011).

14. Dayanithi, G. et al. Characterization of Ca2+ signalling in postnatal mouse retinal ganglion cells: involvement of Opa1 in Ca2+ clearance. *Ophthalmic Genet.* **31**, 53–65 (2010).

15. Tang, S., Le, P. K., Tse, S., Wallace, D. C. & Huang, T. Heterozygous mutation of Opa1 in Drosophila shortens lifespan mediated through increased reactive oxygen species production. *PLoS One* **4**, e4492 (2009).

16. Millet, A. M. et al. Loss of functional Opa1 unbalances redox state: implications in dominant optic atrophy pathogenesis. *Ann. Clin. Transl. Neurol.* **3**, 408–421 (2016).

17. Tezze, C. et al. Age-Associated Loss of Opa1 in Muscle Impacts Muscle Mass, Metabolic Homeostasis, Systemic Inflammation, and Epithelial Senescence. *Cell Metab.* **25**, 1374–1389.e6 (2017).

18. Sarzi, E. et al. The human Opa1 delTTAG mutation induces premature age-related systemic neurodegeneration in mouse. *Brain* **135**, 3599–3613 (2012).

19. White, K. E. et al. Opa1 deficiency associated with increased autophagy in retinal ganglion cells in a murine model of dominant optic atrophy. *Invest. Ophthalmol. Vis. Sci.* **50**, 2567–2571 (2009).

20. Kane, M. S. et al. Autophagy controls the pathogenicity of Opa1 mutations in dominant optic atrophy. *J. Cell. Mol. Med.* **21**, 2284–2297 (2017).

21. Bertholet, A. M. et al. OPA1 loss of function affects *in vitro* neuronal maturation. *Brain* **136**, 1518–1533 (2013).

22. Williams, P. A., Morgan, J. E. & Votruba, M. Opa1 deficiency in a mouse model of dominant optic atrophy leads to retinal ganglion cell dendropathy. *Brain* **133**, 2942–2951 (2010).

23. Williams, P. et al. Opa1 is essential for retinal ganglion cell synaptic architecture and connectivity. *Brain* **135**, 493–505 (2012).

24. Moulis, M. F. et al. OPA1 haploinsufficiency induces a BNIP3-dependent decrease in mitophagy in neurons: relevance to Dominant Optic Atrophy. *J. Neurochem.* **140**, 485–494 (2017).

25. Bocca, C. et al. Plasma Metabolomic Signature Involving Purine Metabolism in Human Optic Atrophy 1 (OPA1)-Related Disorders. *Invest. Ophthalmol. Vis. Sci.* **59**, 185–195 (2018).

26. Chao de la Barca, J. M. et al. Targeted Metabolomics Reveals Early Dominant Optic Atrophy Signature in Optic Nerves of Opa1 delTTAG/+ Mice. *Invest. Ophthalmol. Vis. Sci.* **58**, 812–820 (2017).

27. Bocca, C. et al. The Metabolomic Bioenergetic Signature of Opa1-Disrupted Mouse Embryonic Fibroblasts Highlights Aspartate Deficiency. *Sci. Rep.* **8**, 11528 (2018).

28. Chao de la Barca, J. M. et al. The metabolomic signature of Leber's hereditary optic neuropathy reveals endoplasmic reticulum stress. *Brain* **139**, 2864–2876 (2016).

29. Patten, D. A. et al. OPA1-dependent cristae modulation is essential for cellular adaptation to metabolic demand. *EMBO J.* **33**, 2676–2691 (2014).

30. Birsoy, K. et al. An Essential Role of the Mitochondrial Electron Transport Chain in Cell Proliferation Is to Enable Aspartate Synthesis. *Cell* **162**, 540–551 (2015).

31. Sullivan, L. B. et al. Supporting Aspartate Biosynthesis Is an Essential Function of Respiration in Proliferating Cells. *Cell* **162**, 552–563 (2015).

32. Nguyen, D. et al. A new vicious cycle involving glutamate excitotoxicity, oxidative stress and mitochondrial dynamics. *Cell Death Dis.* **2**, e240 (2011).

33. Papes, F., Surpili, M. J., Langone, F., Trigo, J. R. & Arruda, P. The essential amino acid lysine acts as precursor of glutamate in the mammalian central nervous system. *FEBS Lett.* **488**, 34–8 (2001).

34. Rujiviphat, J., Meglei, G., Rubinstein, J. L. & McQuibban, G. A. Phospholipid association is essential for dynamin-related protein Mgm1 to function in mitochondrial membrane fusion. *J. Biol. Chem.* **284**, 28682–28686 (2009).

35. Kita, T. et al. Possible role of mitochondrial remodelling on cellular triacylglycerol accumulation. *J. Biochem.* **146**, 787–796 (2009).

36. Rambold, A. S., Cohen, S. & Lippincott-Schwartz, J. Fatty acid trafficking in starved cells: regulation by lipid droplet lipolysis, autophagy, and mitochondrial fusion dynamics. *Dev Cell* **32**, 678–692 (2015).

37. Pidoux, G. et al. Optic atrophy 1 is an A-kinase anchoring protein on lipid droplets that mediates adrenergic control of lipolysis. *EMBO J.* **30**, 4371–4386 (2011).

38. Meesapyodsuk, D. & Qiu, X. The front-end desaturase: structure, function, evolution and biotechnological use. *Lipids* **47**, 227–37 (2012).

39. Naudí, A. et al. Membrane lipid unsaturation as physiological adaptation to animal longevity. *Front. Physiol.* **4**, 372 (2013).

## Acknowledgements

This work was supported by Institut National de la Santé et de la Recherche Médicale (INSERM), Centre National de la Recherche Scientifique (CNRS), University of Angers, University Hospital of Angers, University of Toulouse Paul Sabatier. This work was also supported by the following patients' foundations: "Kjer France", "Fondation VISIO", "Ouvrir les Yeux", "Union Nationale des Aveugles et Déficients Visuels", "Association contre les Maladies Mitochondriales", "Retina France", and "Association Point de Mire".
**Author Contributions**

J.M.C.B., P.B. and P.R. designed the study. J.M.C.B., M.S.A., C.B., L.A.P. and I.O. performed the experiments. All authors analyzed and interpreted data. All authors contributed to writing the manuscript and approved the final version. P.B. and P.R. contributed equally to this work.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-42554-7.

**Competing Interests:** The authors declare no competing interests.

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

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019