Marylène Bertrand

Martine Decoville

Hervé Meudal

Serge Birman

Céline Landon

↑ Ethanolamine

Ch

↑ Pyruvate

↓ Acetate

Metabolomic NMR studies at presymptomatic and symptomatic stages of Huntington's disease on a Drosophila model

Keywords: Huntington's disease, neurodegenerative diseases, Drosophila, metabolomics, NMR, presymptomatic 16-day-old 1

Huntington's disease (HD) is an inherited neurodegenerative disorder, for which diagnostic development and discovery of new therapeutic targets are urgently required.

In this study, a model of HD in Drosophila melanogaster has been used to identify metabolic biomarkers at presymptomatic and symptomatic stages of the disease. The pan-neuronal expression of a pathogenic fragment of the human Huntingtin (HTT) protein containing a 93-repeat polyglutamine expansion (Httex1p Q93) in transgenic flies induces a neuropathology with several characteristics of the human disease. The discriminant 2 metabolites between the diseased flies and their controls were identified by 1 H-nuclear magnetic resonance (NMR) and orthogonal partial least squares discriminant multivariate analysis (OPLS-DA). The experiments carried out with 10-day-old flies allowed us to identify a set of 10 biomarkers of the presymptomatic stage: NAD + , AMP, fumarate, asparagine, dimethylamine, -alanine, glutamine, succinate, glutamate and ethanol.

Remarkably, the experiments conducted with 16-day-old flies, when the symptoms of the disease were present, highlighted a different set of 6 biomarkers: phosphocholine, ethanolamine, 2-oxoglutarate, succinate, pyruvate and acetate. Our results provide a better understanding of the metabolic impairments in a widely used HD model and demonstrate that metabolism perturbations change dramatically during the development of the disease.

INTRODUCTION

Huntington's disease (HD) is an autosomal dominant inherited neurodegenerative disorder caused by a mutation in the huntingtin (HTT) gene located at band 4p16.3 on chromosome 4. [START_REF] Macdonald | A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington's disease chromosones[END_REF] Patients developing HD have an expanded CAG repeat in the first exon of the gene that includes from 36 to 121 repeats, [START_REF] Andrew | The relationship between trinucleotide (CAG) repeat length and clinical-features of Huntingtons-Disease[END_REF] leading to polyglutamine expansion in the Huntingtin protein (HTT). Depending on the length of the CAG repeat and of other genetic and environmental conditions, symptoms usually appear between the ages of 30 to 50. These symptoms include motor dysfunction, psychiatric disturbances and cognitive decline. The patients also develop other symptoms such as muscle atrophy, impaired glucose tolerance and weight loss. There is currently no treatment to stop or reverse the course of the disease, which generally leads to death within 10 to 30 years after emergence. [START_REF] Saudou | The Biology of Huntingtin[END_REF] However, some medications can limit the clinical defects. [START_REF] Carroll | Treating the whole body in Huntington's disease[END_REF][START_REF] Bates | Huntington disease[END_REF] At the molecular level, the mutant Huntingtin protein (mHTT) accumulates abnormally.

Misfolded fragments of this protein containing the polyglutamine expansion form cytoplasmic aggregates and nuclear inclusions throughout the brain, [START_REF] Labbadia | Huntington's disease: underlying molecular mechanisms and emerging concepts[END_REF] leading to progressive neuronal toxicity. The disruption of several cellular mechanisms has been reported, including alterations in gene expression and transcription, protein folding and degradation, synaptic signaling, and energy metabolism. [START_REF] Mccolgan | Huntington's disease: a clinical review[END_REF][START_REF] Caboche | Handbook of Behavioral Neuroscience[END_REF][START_REF] Illarioshkin | Molecular Pathogenesis in Huntington's Disease[END_REF] To better understand the causes and consequences of the disease, different studies have been carried out on humans and animal models to discover metabolic biomarkers of HD, [START_REF] Dumas | Metabolic Profiling and Phenotyping of Central Nervous System Diseases: Metabolites Bring Insights into Brain Dysfunctions[END_REF] most of them focusing on the brain, either post mortem [START_REF] Graham | Metabolic signatures of Huntington's disease (HD): H-1 NMR analysis of the polar metabolome in postmortem human brain[END_REF][START_REF] Patassini | Metabolite mapping reveals severe widespread perturbation of multiple metabolic processes in Huntington's disease human brain[END_REF] or in vivo. [START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF][START_REF] Unschuld | Brain metabolite alterations and cognitive dysfunction in early Huntington's disease[END_REF] Some metabolomic analyses have also been conducted in biological fluids such as plasma [START_REF] Mastrokolias | Integration of targeted metabolomics and transcriptomics identifies deregulation of phosphatidylcholine metabolism in Huntington's disease peripheral blood samples[END_REF][START_REF] Cheng | Metabolic disturbances in plasma as biomarkers for Huntington's disease[END_REF] or cerebrospinal fluid. [START_REF] Verwaest | H-1 NMR based metabolomics of CSF and blood serum: A metabolic profile for a transgenic rat model of Huntington disease[END_REF] HD animal models have been developed in various organisms such as mouse, [START_REF] Pont | Metabolic profiling for the identification of Huntington biomarkers by on-line solid-phase extraction capillary electrophoresis mass spectrometry combined with advanced data analysis tools[END_REF][START_REF] Ismailoglu | Huntingtin protein is essential for mitochondrial metabolism, bioenergetics and structure in murine embryonic stem cells[END_REF][START_REF] Zacharoff | Cortical metabolites as biomarkers in the R6/2 model of Huntington's disease[END_REF] rat, [START_REF] Verwaest | H-1 NMR based metabolomics of CSF and blood serum: A metabolic profile for a transgenic rat model of Huntington disease[END_REF][START_REF] Chang | Metabolic Profiling of 3-Nitropropionic Acid Early-Stage Huntington's Disease Rat Model Using Gas Chromatography Time-of-Flight Mass Spectrometry[END_REF][START_REF] Tsang | Metabonomic Characterization of the 3-Nitropropionic Acid Rat Model of Huntington's Disease[END_REF] sheep [START_REF] Skene | Metabolic profiling of presymptomatic Huntington's disease sheep reveals novel biomarkers[END_REF] or invertebrates such as the fruit fly Drosophila. [START_REF] Singh | NMR Spectroscopy-based Metabolomics of Drosophila Model of Huntington's Disease Suggests Altered Cell Energetics[END_REF][START_REF] Rosas-Arellano | The Tiny Drosophila Melanogaster for the Biggest Answers in Huntington's Disease[END_REF][START_REF] Marsh | Drosophila in the study of neurodegenerative disease[END_REF] Here we selected Drosophila as this model organism is well-established to study human diseases. Drosophila are easy to handle and many approved genetic tools are available for transgenesis and gene expression / inactivation in specific cell and tissue types. The Drosophila genome is 25 times smaller than the human genome, but about 75% of the genes involved in human diseases have orthologues in Drosophila. Although its brain appears to be simple compared to mammals, Drosophila display complex behaviors such as learning and memory, sleep or courtship. [START_REF] Rosas-Arellano | The Tiny Drosophila Melanogaster for the Biggest Answers in Huntington's Disease[END_REF] In addition, its nervous system is composed of various types of neurons (dopaminergic, serotonergic, glutamatergic, GABAergic …) and glia as in vertebrates. Numerous models of neurodegenerative diseases, including Parkinson's or Alzheimer's diseases, have been developed in Drosophila. [START_REF] Mcgurk | Drosophila as an In Vivo Model for Human Neurodegenerative Disease[END_REF] HD models have been described by overexpressing either in the eyes, neurons or glial cells, amino terminal fragments of human HTT containing various tracts of polyQ, or full-length human mutant HTT (128Q). [START_REF] Rosas-Arellano | The Tiny Drosophila Melanogaster for the Biggest Answers in Huntington's Disease[END_REF][START_REF] Marsh | Drosophila in the study of neurodegenerative disease[END_REF][START_REF] Mcgurk | Drosophila as an In Vivo Model for Human Neurodegenerative Disease[END_REF][START_REF] Lewis | Using Drosophila models of Huntington's disease as a translatable tool[END_REF] These models summarize the characteristic symptoms of HD, namely protein aggregation, neurodegeneration, behavioral defects and reduced lifespan.

In this report, we describe the changing metabolic profile of a pan-neuronal Drosophila model of HD. The purpose of our study was to identify biomarkers of HD using 1 H-nuclear magnetic resonance (NMR) at two stages of the disease, presymptomatic and symptomatic, using 10-day-old and 16-day-old HD flies, respectively. We show that the metabolic signature of the disease is significantly different at these two stages.

EXPERIMENTAL PROCEDURES

Drosophila

Culture

The fly stocks were maintained at 22°C on a standard medium (per liter: 90.25 g cornmeal, 82.5 g dry yeast, 10.75 g agar and 37.5 mL of a 10% solution of methyl-4hydroxybenzoate in ethanol). Crosses were performed at 25°C on standard medium. The F1 progeny was collected each day and maintained at 26°C on Nutri-Fly medium (Dominique Dutscher, Brumath, France). Drosophila lines were: w 1118 used as wild-type controls, elav-Gal4 (elav C155 ) [START_REF] Lin | Ectopic and increased expression of fasciclin-II alters motoneuron growth cone guidance[END_REF] and UAS-Httex1p Q93. [START_REF] Steffan | Histone deacetylase inhibitors arrest polyglutamine-dependent neurodegeneration in Drosophila[END_REF] To express Httex1p Q93 in neurons, elav-Gal4 virgin females were crossed with UAS-Httex1p-Q93 males (HD flies) (Figure 1). As controls, we crossed elav-Gal4 virgins with w 1118 males (Gal4 controls) and w 1118 virgins with UAS-Httex1p Q93 males (UAS controls). Only F1 progeny females were used for metabolite analysis in this study. The analysis was carried out in females specifically because they are larger than males and heterozygous for the X-linked elav-Gal4 insertion.

Sample Preparation

Eight to twelve samples were prepared for each batch of HD flies or controls. In order to have a good signal / noise ratio in NMR spectra, five female Drosophila were used per sample. They were flash frozen in a 1-mL microtube by immersion in liquid nitrogen. Flies were then manually crushed with a suitable conical stick in the microtube and 400 µL of ice-cold 50% acetonitrile was added as described in previous studies. [START_REF] Ott | Metabolic changes may precede proteostatic dysfunction in a Drosophila model of amyloid beta peptide toxicity[END_REF][START_REF] Sarup | The Metabolic Profile of Long-Lived Drosophila melanogaster[END_REF] The mixture was shaken with a vortex, sonicated during 5 min in an ice-cold bath and then centrifuged for 10 min at 15,300 g at 5°C. The supernatants (350 µL) were transferred to new tubes.

Acetonitrile was evaporated in a SpeedVac for 30 minutes. The tube contents were freezedried overnight and stored at -80 °C.

For NMR analyses, 300 µL of an ice-cold 50 mM phosphate buffer (pH 7.4) in D2O containing 7.5 µg/L 3-(trimethylsilyl)-propionic-2,2,3,3-d4 acid sodium salt (TSP-D4) was added to the lyophilized samples. After homogenization with vortex, centrifugation for 10 min at 15,300 g at 5°C, 250 µL of the supernatant were transferred to a 3-mm NMR tube, and kept at 5°C.

NMR Acquisition and Processing

The NMR measurements were performed at 25°C on a Bruker Avance-III HD 700 spectrometer, equipped with a 5-mm HCN triple resonance cryoprobe. Bruker's TopSpin software 3.2 package was used for the acquisition and processing of NMR spectra. [START_REF] Macdonald | A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington's disease chromosones[END_REF] H-NMR spectra were acquired using a cpmgpr1d pulse program, which attenuates broad signals from high-molecular-weight components (D1 = 4s; 128 scans; spectral width: 11 ppm; experimental time: 14 min). The water signal was suppressed by presaturation of the water peak. For assignment purposes, bidimensional COSY, TOCSY and 13 C-HSQC spectra were acquired (see Figure S1 in supporting information). The first spectrum was manually phased. The following spectra were automatically phased and referenced to TSP resonance at 0.0 ppm and their baselines were adjusted. A moderate local baseline correction [START_REF] Zhang | Baseline correction using adaptive iteratively reweighted penalized least squares[END_REF] and a peak alignment of the NMR spectra [START_REF] Bloemberg | Improved parametric time warping for proteomics[END_REF][START_REF] Wehrens | Chemometrics with R : Multivariate Data Analysis in the Natural Sciences and Life Sciences[END_REF] was performed with NMRProcFlow open source software (https://www.nmrprocflow.org/ [START_REF] Jacob | NMRProcFlow: a graphical and interactive tool dedicated to 1D spectra processing for NMR-based metabolomics[END_REF] ). The NMR spectra (2 zones excluding the water peak's region from 4.86 to 4.75 ppm) were segmented into variable buckets using the "intelligent bucketing" method [START_REF] De Meyer | NMR-based characterization of metabolic alterations in hypertension using an adaptive, intelligent binning algorithm[END_REF] of the NMRProcFlow software.

The spectral regions from 1.38 -1.25 ppm and 0.93-0.85 ppm (signals of residual lipids) and the spectral regions having a Signal / Noise Ratio (SNR) lower than 3 were excluded.

All NMR spectra were processed with nmrprocflow in exactly the same way so that they could be compared reliably, and were exported into a unique data matrix.

Statistical Analysis

This data matrix was further split into three sets of data, allowing three analyses: i) the presymptomatic analysis compares the 10-day-old fly samples, carrying the disease (HD flies) or not (controls); ii) the symptomatic analysis compares the 16-day-old fly samples, carrying the disease or not and; iii) the change analysis compares the HD flies at the two stages of the disease.

Each set of data was imported in Workflow4Metabolomics 3.0 Galaxy online infrastructure for the statistical analysis achievement (https://galaxy.workflow4metabolomics.org/ [START_REF] Giacomoni | Workflow4Metabolomics: a collaborative research infrastructure for computational metabolomics[END_REF] ).

Before the statistical analyses, the buckets in each NMR spectrum were normalized to the total sum of the spectral intensity. For each set of data, multivariate as well as univariate statistical analyses were performed to check the quality and the dispersion of samples, the quality of variables, and to find the discriminant variables between the groups.

Multivariate Statistical Analyses

OPLS-DA (orthogonal partial least-squares discriminant analysis) was used as multivariate analysis in order to find discriminant variables of two groups of samples. [START_REF] Trygg | Chemometrics in metabonomics[END_REF] Three parameters were followed to evaluate the quality of OPLS-DA models: R2X and R2Y to quantify the goodness-of-fit and Q2Y to quantify predictive accuracy of the model.

Permutation tests (100 cycles) were conducted to assess the robustness. Discriminant variables were defined by a threshold of variable importance of projection (VIP) coefficients in the prediction greater than 1, derived from the OPLS-DA model. This method offers the advantage to remove disturbing variation.

Univariate Statistical Analyses

Through the Workflow4Metabolomics online infrastructure [START_REF] Giacomoni | Workflow4Metabolomics: a collaborative research infrastructure for computational metabolomics[END_REF] univariate statistical analysis were performed, in order to check the significance of the variables. The factor of interest being qualitative with 2 levels, t-test with p-value significance threshold of 0.05 was performed and p and q-values were calculated. [START_REF] Colquhoun | An investigation of the false discovery rate and the misinterpretation of p-values[END_REF] The q-value gives the expected positive false discovery rate obtained by rejecting the null hypothesis for any result with an equal or smaller q-value. For metabolites with several peaks, common or average p and qvalues of the different peaks are given.

Metabolite Identification

The NMR peaks highlighted by VIPs upper than 1 were then assigned. The corresponding metabolites were identified from their 1 H and [START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF] C chemical shifts, with databases such as HMDB, [START_REF] Wishart | HMDB 4.0: the human metabolome database for 2018[END_REF] BRMB 42 and a home-made database. For unambiguous identifications, when needed, some reference compounds were directly added to the samples to identify precisely the affected peaks of interest. It is worth mentioning that metabolites displaying several peaks in the 1 H-NMR spectrum were retained only when all the peaks were present and when all or most of the peaks were discriminant.

Discrimination performance of biomarkers

Using the open source MetaboanalystR software [START_REF] Xia | MetaboAnalyst 3.0-making metabolomics more meaningful[END_REF] (https://www.metaboanalyst.ca) prediction models and ROC (Receiver Operating Characteristic) curves were performed to identify a set of metabolites capable of classifying unknown samples as HD or control samples with high sensitivity and specificity. A matrix composed of all samples (10 and 16-day-old HD and control flies) and of the 46 variables / peaks of the discriminating metabolites was created and ROC curves were constructed to select the best metabolite combination. A model with 15 metabolites (one variable by metabolite) was created.

Finally, the performance of the model was evaluated with ROC curve and logistic regression analyses were performed with selected variables and samples.

RESULTS

Nervous system expression of the first exon of human HTT with 93 CAG repeats (Httex1p Q93) has been extensively characterized in previous works as a reliable HD model in Drosophila. [START_REF] Steffan | Histone deacetylase inhibitors arrest polyglutamine-dependent neurodegeneration in Drosophila[END_REF][START_REF] Lievens | Expanded polyglutamine peptides disrupt EGF receptor signaling and glutamate transporter expression in Drosophila[END_REF][START_REF] Lievens | AKT-sensitive or insensitive pathways of toxicity in glial cells and neurons in Drosophila models of Huntington's disease[END_REF][START_REF] Pallos | Inhibition of specific HDACs and sirtuins suppresses pathogenesis in a Drosophila model of Huntington's disease[END_REF][START_REF] Zhang | Inactivation of Drosophila Huntingtin affects long-term adult functioning and the pathogenesis of a Huntington's disease model[END_REF][START_REF] Besson | Enhanced Neuronal Glucose Transporter Expression Reveals Metabolic Choice in a HD Drosophila Model[END_REF][START_REF] Besson | Increased energy metabolism rescues glia-induced pathology in a Drosophila model of Huntington's disease[END_REF] Here we used elav-Gal4 as a driver, allowing expression of the polyglutaminated protein in all neurons in adult flies which causes reduced viability and neuronal degeneration in adult flies. Because Httex1p Q93 was not only expressed in the brain but also in the ventral nerve cord and peripheral neurons with this driver, metabolites were extracted from the entire Drosophila and not only from the head. In this study, we performed NMR experiments on 10-and 16-day-old HD flies and on their two controls: 10-day-old flies are at a presymptomatic stage, which means that they do not show a significantly higher death rate than the controls, while the survival of 16-day-old HD flies is reduced to about 50% of that of the controls, indicating that they have reached a fully symptomatic stage. [START_REF] Lievens | AKT-sensitive or insensitive pathways of toxicity in glial cells and neurons in Drosophila models of Huntington's disease[END_REF] A typical NMR spectrum obtained with 16-day-old HD flies is shown in Figure 2.

To eliminate the contribution of each transgene to modifications of the metabolite patterns, two controls were processed in parallel. The first one corresponds to flies carrying only the elav-Gal4 transgene (Gal4 control) and the second one to flies carrying only the UAS-Httex1p Q93 transgene (UAS control). For each experiment, statistical analysis was performed by comparing HD flies with Gal4 controls then with UAS controls, and finally HD flies with the associated Gal4 and UAS controls. Only VIPs common to the three comparisons were retained for further analysis. This method excludes the variables which are discriminant for one type of controls only and keeps the most reliable variables.

Presymptomatic stage biomarkers

An OPLS-DA model was used to find the discriminant metabolites between the 10-dayold HD flies and the controls (Figure 3). Only discriminant variables (VIPs>1) from OPLS-DA and showing confidence interval at the 95% level were retained: i) 157 variables were identified between HD flies and the associated two controls, ii) 228 variables between HD flies and Gal4 controls, and finally iii) 189 variables between HD flies and UAS controls.

All the OPLS-DA analyses performed led to a good or excellent separation of the groups and had large R2X (0.80; 0.77 and 0.72, respectively) and R2Y (0.99; 0.98 and 0.96, respectively) which indicated high reproducibility and a good fitting of the model to the data. The Q2Y parameters, calculated by cross validation were close to 1 (0.93; 0.93 and 0.94, respectively) which are indicative of the high predictivity of the models. The discriminant variables of the three OPLS-DA were compared to each other using a Venn diagram and the 79 variables common to the three groups were retained for further analysis. The metabolites identified for the 10-days-old flies from the 79 discriminant variables are presented in Table 1 (VIP score of each of the discriminant metabolites is displayed in the Figure S2 as supporting information). All the q-values are very low, which means that the identified metabolites are statistically highly significant. At this stage of the disease, all the identified metabolic biomarkers are lower in HD flies than in controls.

Symptomatic stage biomarkers

The same procedure as for the 10-day-old flies was conducted for the 16-day-old flies (Figure 4). The discriminant variables (VIPs) >1 from OPLS-DA and showing confidence interval at the 95% level were retained: i) 156 variables were found between HD flies and Only 6 metabolites were unambiguously identified from the 46 discriminant variables, and are presented in Table 2 (the VIP score of each of the discriminant metabolites is displayed in the Figure S3 as supporting information). Note that at this stage of the disease, four metabolic biomarkers have increased in HD flies, whereas two have decreased.

Table 2. Metabolites identified from the 46 common discriminant variables with VIP >1, for the 16-day-old flies. For metabolites with several peaks, the q-value is common to the different peaks. Downward arrows (↓) and upward arrows (↑) indicate that the proportion of metabolites was lower or, respectively, higher for HD flies than for controls. s-singlet, t-triplet, m-multiplet. ov-overlapped. * 0.01< Value < 0.05; ** 0.001< Value < 0.01; *** 0.0001< Value < 0.001, Value < 0.0001.

Analysis of biomarker change

In order to follow the change of the identified metabolites with the progression of the disease, the comparison of 10-day-old and 16-day-old HD flies was carried out. The OPLS-DA analysis highlighted 177 discriminant variables (Figure 5), making it possible to identify 15 metabolites (Table 3). Note that all the metabolites previously identified (Tables 1 and2) were discriminant between 10-and 16-day-old flies, with the exception of pyruvate. Taurine, which was not identified in the previous analysis, is now discriminant between 10-and 16-day-old HD flies.

To make sure that the change of biomarkers is well correlated with the progression of the disease and not with the age of flies, the same analysis was performed between 10-and 16-day-old control flies. Results are presented in Table 3. We observed that three metabolites (glutamate, acetate and ethanol) changed in the same way in diseased and control flies according to age.

Discrimination performance of biomarkers

We have performed ROC analysis to evaluate the predictive value of the discriminant metabolites. ROC curves were constructed with different combination (2 to 46 variables).

The results of multi-ROC analyses showed that the models with 10 or more variables have an area under the curve (AUC) greater than 95% (see Figure S4). The AUC value close to 1 shows high sensitivity and specificity of the model. Then, a set of 15 discriminant metabolites were manually selected based on their LASSO frequencies and of their mean importance measures (phosphocholine, glutamine, pyruvate, taurine, -alanine, succinate, acetate, 2-oxoglutarate, ethanol, ethanolamine, NAD + , AMP, fumarate, glutamate and dimethylamide). The significance of this biomarker model was evaluated using ROC curve and logistic regression (Figure S5). Then different sets of samples were tested to see if the set of the 15 variables was able to sort the HD samples from the controls. The predicted class probabilities for each sample showed a good classification accuracy where 29 out of 30 samples were correctly classified (example on Figure S6).

DISCUSSION

In this study, a transgenic model of HD in Drosophila melanogaster has been used to identify metabolic biomarkers at presymptomatic (10-day-old flies) and symptomatic (16day-old flies) stages of the disease. Our study has been performed on the whole body of the flies because in this model the mutant Htt has been expressed in all fly neurons, which are present in the head as components of the brain, but also in other parts of the fly body, such as the ventral nerve cord in the thorax or the enteric nervous system in the abdomen.

A complete analysis of HD-associate defects in our model therefore required to study metabolites in the entire fly. Moreover, brain activity and state are known to influence the physiology of the whole organism, both in health and disease conditions. The discriminant metabolites between the disease carrier flies and their controls were identified by 1 H-NMR and OPLS-DA multivariate analyses. This allowed a better understanding of the metabolic impairments associated with the disease and to follow the metabolite disturbance over its progression.

The three statistical analyses reported here enabled us to: i) identify 10 biomarkers of the presymptomatic stage; ii) identify 6 biomarkers of the symptomatic stage and finally, iii) monitor the change of these biomarkers during HD-like pathogenesis. It is worth mentioning that some metabolites could not be unambiguously identified from discriminant variables. Therefore, the final lists of metabolites only report the compounds that we could identified with certainty. Moreover, the regions corresponding to residual lipid peaks in NMR spectra were removed from the matrix before the statistical analyses, in order not to affect the results. As a consequence, the lactate peaks, at 1.33 and 1.35 ppm, could not be considered in our work although lactate is a known HD biomarker which has been highlighted in several previous studies. [START_REF] Verwaest | H-1 NMR based metabolomics of CSF and blood serum: A metabolic profile for a transgenic rat model of Huntington disease[END_REF][START_REF] Singh | NMR Spectroscopy-based Metabolomics of Drosophila Model of Huntington's Disease Suggests Altered Cell Energetics[END_REF][START_REF] Tsang | Metabolic characterization of the R6/2 transgenic mouse model of Huntington's disease by high-resolution MAS (1)H NMR spectroscopy[END_REF] The discriminating performance of these metabolites was tested on the whole samples by an exploratory analysis using multivariate ROC curves. A minimum of ten metabolites is required to have a good discriminant ability. A set of 15 metabolites allows to class between HD and control samples with a good accuracy and to predict them.

The metabolites identified in our study can be classified in three groups: i) the first group includes the metabolites directly involved in energy metabolism (fumarate, succinate, 2oxoglutarate (or α-ketoglutarate), pyruvate, acetate, NAD + and AMP), ii) the second group corresponds to amino acids and derivatives (glutamate, glutamine, asparagine, β-alanine and taurine), and iii) the third group corresponds to miscellaneous metabolites (dimethylamine, acetate, phosphocholine, ethanol and ethanolamine). Most of these discriminant metabolites have been reported in the literature as biomarkers of HD, in particular glutamate, [START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF][START_REF] Unschuld | Brain metabolite alterations and cognitive dysfunction in early Huntington's disease[END_REF][START_REF] Graham | Novel Metabolite Biomarkers of Huntington's Disease As Detected by High-Resolution Mass Spectrometry[END_REF][START_REF] Behrens | Impaired glutamate transport and glutamate-glutamine cycling: downstream effects of the Huntington mutation[END_REF] glutamine, [START_REF] Patassini | Metabolite mapping reveals severe widespread perturbation of multiple metabolic processes in Huntington's disease human brain[END_REF][START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF][START_REF] Verwaest | H-1 NMR based metabolomics of CSF and blood serum: A metabolic profile for a transgenic rat model of Huntington disease[END_REF][START_REF] Zacharoff | Cortical metabolites as biomarkers in the R6/2 model of Huntington's disease[END_REF][START_REF] Singh | NMR Spectroscopy-based Metabolomics of Drosophila Model of Huntington's Disease Suggests Altered Cell Energetics[END_REF][START_REF] Tsang | Metabolic characterization of the R6/2 transgenic mouse model of Huntington's disease by high-resolution MAS (1)H NMR spectroscopy[END_REF][START_REF] Behrens | Impaired glutamate transport and glutamate-glutamine cycling: downstream effects of the Huntington mutation[END_REF][START_REF] Joyner | Probing the Metabolic Aberrations Underlying Mutant Huntingtin Toxicity in Yeast and Assessing Their Degree of Preservation in Humans and Mice[END_REF] phosphocholine, [START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF][START_REF] Mastrokolias | Integration of targeted metabolomics and transcriptomics identifies deregulation of phosphatidylcholine metabolism in Huntington's disease peripheral blood samples[END_REF][START_REF] Cheng | Metabolic disturbances in plasma as biomarkers for Huntington's disease[END_REF][START_REF] Skene | Metabolic profiling of presymptomatic Huntington's disease sheep reveals novel biomarkers[END_REF] acetate, 50,53 succinate 17,24,53 and pyruvate. [START_REF] Singh | NMR Spectroscopy-based Metabolomics of Drosophila Model of Huntington's Disease Suggests Altered Cell Energetics[END_REF][START_REF] Naia | Insulin and IGF-1 regularize energy metabolites in neural cells expressing full-length mutant huntingtin[END_REF][START_REF] Papsdorf | Polyglutamine toxicity in yeast induces metabolic alterations and mitochondrial defects[END_REF] Other metabolites, such as leucine, isoleucine, valine, [START_REF] Skene | Metabolic profiling of presymptomatic Huntington's disease sheep reveals novel biomarkers[END_REF][START_REF] Cheng | Metabolic disturbances in plasma as biomarkers for Huntington's disease[END_REF][START_REF] Mochel | Early Energy Deficit in Huntington Disease: Identification of a Plasma Biomarker Traceable during Disease Progression[END_REF] myo-inositol, [START_REF] Graham | Metabolic signatures of Huntington's disease (HD): H-1 NMR analysis of the polar metabolome in postmortem human brain[END_REF][START_REF] Sturrock | A Longitudinal Study of Magnetic Resonance Spectroscopy Huntington's Disease Biomarkers[END_REF][START_REF] Tkac | Neurochemical changes in Huntington R6/2 mouse striatum detected by in vivo H-1 NMR spectroscopy[END_REF] N-acetyl-aspartate (NAA), [START_REF] Unschuld | Brain metabolite alterations and cognitive dysfunction in early Huntington's disease[END_REF][START_REF] Sturrock | A Longitudinal Study of Magnetic Resonance Spectroscopy Huntington's Disease Biomarkers[END_REF][START_REF] Van Den Bogaard | Exploratory 7-Tesla magnetic resonance spectroscopy in Huntington's disease provides in vivo evidence for impaired energy metabolism[END_REF] gaminobutyric acid (GABA) [START_REF] Dumas | Metabolic Profiling and Phenotyping of Central Nervous System Diseases: Metabolites Bring Insights into Brain Dysfunctions[END_REF][START_REF] Patassini | Metabolite mapping reveals severe widespread perturbation of multiple metabolic processes in Huntington's disease human brain[END_REF][START_REF] Tsang | Metabonomic Characterization of the 3-Nitropropionic Acid Rat Model of Huntington's Disease[END_REF] were quite often pointed out, but they did not appear as discriminant in our analyses. The closest report to our study is the recent work of Singh and collaborators [START_REF] Singh | NMR Spectroscopy-based Metabolomics of Drosophila Model of Huntington's Disease Suggests Altered Cell Energetics[END_REF] on a Drosophila model of HD displaying retinal deformations. Among the metabolites they identified by NMR analysis performed on methanolic extracts of heads, four are common to ours: NAD + , pyruvate, succinate and glutamine. These authors used the GMR-Gal4 driver, which selectively expresses in the fly compound eyes, explaining why the analysis was focused on heads, unlike our model in which mutant Htt is expressed in all neurons present over the entire body. Our analysis was carried out, therefore, on a larger diversity of neuronal subtypes than that used in Singh's study. The eye photoreceptors are highly specialized cells that are different from all other neurons, so one can expect that the results obtained in the retina will not be identical to those obtained in brain or other tissues. However, it is interesting to note the similarities between our results and those of Singh et al., particularly regarding energetic pathways.

Brain energy supply is well known to be dysfunctional in HD. [START_REF] Tkac | Homeostatic adaptations in brain energy metabolism in mouse models of Huntington disease[END_REF][START_REF] Mochel | Energy deficit in Huntington disease: why it matters[END_REF] Numerous evidence suggests a dysfunction of mitochondria in HD, and more specifically of the electron transport chain. However, other energy generating pathways are also implicated. In particular, the tricarboxylic acid (TCA) cycle is disrupted, [START_REF] Verwaest | H-1 NMR based metabolomics of CSF and blood serum: A metabolic profile for a transgenic rat model of Huntington disease[END_REF][START_REF] Dubinsky | Towards an Understanding of Energy Impairment in Huntington's Disease Brain[END_REF] and compensatory pathways are used to fuel the TCA cycle, such as fatty acid oxidation or glutamate oxidation for example. It has been shown that increasing glucose transport in neurons significantly enhanced Drosophila survival in the Httex1p Q93 model. [START_REF] Besson | Enhanced Neuronal Glucose Transporter Expression Reveals Metabolic Choice in a HD Drosophila Model[END_REF] The same authors observed that overexpression of glucose-6-phosphate dehydrogenase (G6PD), the key enzyme in the pentose phosphate (PPP) pathway, also markedly extended the lifespan of HD flies.

Our results highlight metabolites involved in such bioenergetic pathways: fumarate, succinate and 2-oxoglutarate, which are intermediates of the TCA cycle, and pyruvate, which is the end product of glycolysis. Acetate is also involved in energy metabolism as it may derived from fatty acid oxidation, a pathway used when carbohydrates are low.

Polyzos et al. showed in a mouse model of HD (HdhQ(150/150)) that the level of glucose was lower in the striatum of HD mice, and that astrocytes reprogramed their metabolism to use fatty acids as a fuel. [START_REF] Polyzos | Metabolic Reprogramming in Astrocytes Distinguishes Region-Specific Neuronal Susceptibility in Huntington Mice[END_REF] NAD + and NADP + are coenzymes in enzymatic oxidationreduction (redox) reactions such as glycolysis and the TCA cycle. They are also involved in the electron transport chain. A decrease in NAD + was already noted in two studies with two different animal models of HD : Drosophila 24 and the mouse. [START_REF] Toczek | An impaired metabolism of nucleotides underpins a novel mechanism of cardiac remodeling leading to Huntington's disease related cardiomyopathy[END_REF] Among the amino acids identified at presymptomatic stage, glutamate and glutamine are non-essential amino acids, which are involved in brain functions. In neurons, glutamate acts as an excitatory neurotransmitter. In astrocytes, it is used to produce glutamine, which is transported into neurons and then converted into glutamate. However, glutamate may also participate to fuel the TCA cycle by generating 2-oxoglutarate. At 10 days, HD flies have a lower level of glutamate and glutamine than controls, which may suggest that these two amino acids are preferentially used to produce energy. We also noticed a decrease in asparagine at presymptomatic stage. This was already reported in the metabolic profiling of presymptomatic HD in sheep [START_REF] Skene | Metabolic profiling of presymptomatic Huntington's disease sheep reveals novel biomarkers[END_REF] and in the plasma of HD patients. [START_REF] Gruber | Huntington' diseaseimbalance of amino acid levels in plasma of patients and mutation carriers[END_REF] Like glutamate, asparagine can also participates in energy metabolism as an intermediate of oxaloacetate.

Finally, other metabolites such as phosphocholine and ethanolamine were found to increase at symptomatic stage in our model. They are precursors of phosphatidylcholine and phosphatidylethanolamine, which are major components of cellular and mitochondrial membranes and more specifically of the neuronal membranes. The increase in phosphocholine and ethanolamine observed in this study likely reflects therefore an altered turnover of membranes resulting from neurodegeneration. Alterations in phosphocholine levels have previously been reported in the blood of HD patients [START_REF] Mastrokolias | Integration of targeted metabolomics and transcriptomics identifies deregulation of phosphatidylcholine metabolism in Huntington's disease peripheral blood samples[END_REF] and in various mouse models. [START_REF] Pepin | In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington's disease[END_REF][START_REF] Zacharoff | Cortical metabolites as biomarkers in the R6/2 model of Huntington's disease[END_REF][START_REF] Peng | Characterization of Behavioral, Neuropathological, Brain Metabolic and Key Molecular Changes in zQ175 Knock-In Mouse Model of Huntington's Disease[END_REF] HD is a progressive disease that gets worse with age and it is tempting to think that metabolic alterations will also worsen with age. However, Dubinsky [START_REF] Dubinsky | Towards an Understanding of Energy Impairment in Huntington's Disease Brain[END_REF] suggested that the compensation phenomena that arise in HD evolve in the course of the disease. In this study, we have highlighted the discriminating metabolites in presymptomatic and symptomatic flies. Comparative analysis between the 10-day-old HD flies and 16-day-old HD flies identifies 15 discriminants metabolites. Among them, three (glutamate, acetate and ethanol) were also discriminant between 10-day-old and 16-day-old control flies suggesting that the variation in the rate of these metabolites is more a consequence of age rather than of the course of the disease. With the exception of taurine, these metabolites are those identified in 10-day-old and 16-day-old HD flies. Most of these metabolites were discriminant only in 10-day-old HD flies or only in 16-day-old HD flies when comparing to their controls. For example, fumarate has a lower level when compared with controls in 10-day-old HD flies, but was not discriminant when the 16-dayold HD flies were compared to their controls. These results suggest that discriminating metabolites at 10 days returned to a normal level at 16 days and that new discriminating metabolites appeared at this stage. Regarding succinate, its level changes in the opposite direction compared to the controls between these two stages, decreasing in 10-day-old HD flies, but increasing in 16-day-old HD flies. It is interesting to note that taurine is not discriminating when comparing 10-day-old and 16-day-old HD flies with their controls, but emerges as a discriminating metabolite when comparing 10-day-old and 16-day old HD flies. The discriminating metabolites retained in the results are those common from the three comparison pairs HD flies / Gal4 controls, HD flies / UAS controls and HD flies / sum of controls, as indicated. In this case taurine was not discriminant, but when comparing 10-day-old and 16-day-old HD flies only with the sum of the controls or only with the Gal4 control taurine was discriminant, its level being lower in 10-day-old HD flies and higher in 16-day-old HD flies. If the taurine level changes in an opposite direction in 10-day-old and 16-day-old HD flies, taurine can become a discriminating metabolite between these two stages of the disease.

Taurine is a sulfur-containing amino-acid displaying multiple functions. [START_REF] Tsang | Metabonomic Characterization of the 3-Nitropropionic Acid Rat Model of Huntington's Disease[END_REF][START_REF] Cheng | Metabolic disturbances in plasma as biomarkers for Huntington's disease[END_REF][START_REF] Tkac | Neurochemical changes in Huntington R6/2 mouse striatum detected by in vivo H-1 NMR spectroscopy[END_REF] It is involved in brain development, neuromodulation, osmoregulation, stabilization of membranes, inflammation, neuroprotection and it has also an antioxidant function. [START_REF] Menzie | Taurine and central nervous system disorders[END_REF] Recently, it has been shown that taurine has a protective effect in models of neurodegenerative disorders and in particular in a rat HD model. [START_REF] Jakaria | Taurine and its analogs in neurological disorders: Focus on therapeutic potential and molecular mechanisms[END_REF][START_REF] Tadros | Neuroprotective effect of taurine in 3-nitropropionic acid-induced experimental animal model of Huntington's disease phenotype[END_REF] The change in taurine level during the progression of the disease in our Drosophila model, could reflect this role in neuroprotection.

As a whole, our results describe that some metabolites decrease at early stage and increase at a later stage, demonstrating that the impairments in metabolite pathways change over the course of the disease. Such observations have already been noted in the literature. [START_REF] Dubinsky | Towards an Understanding of Energy Impairment in Huntington's Disease Brain[END_REF][START_REF] Polyzos | Metabolic Reprogramming in Astrocytes Distinguishes Region-Specific Neuronal Susceptibility in Huntington Mice[END_REF] A majority of the metabolites identified in this study are related to energy metabolism, they are intermediates of the TCA cycle or intermediates of compensatory energy mediating pathways such as fatty acid oxidation or oxidative glutamate pathway.

During the progression of the disease, the cells adapt their metabolism to compensate the deficit in energy, but adopt different compensatory pathways based on accumulated metabolic and transcriptional alterations. The different energy supply pathways are tightly regulated and interconnected. For example, AMP activates several enzymes of the TCA cycle (pyruvate dehydrogenase, citrate synthase and isocitrate dehydrogenase) and the ratio NAD + / NADH influences the pyruvate / lactate ratio. [START_REF] Dubinsky | Towards an Understanding of Energy Impairment in Huntington's Disease Brain[END_REF] Other studies have reported different expression rates of enzymes involved in energy biosynthesis pathways. For example, Zabel et al. [START_REF] Zabel | A Large Number of Protein Expression Changes Occur Early in Life and Precede Phenotype Onset in a Mouse Model for Huntington Disease[END_REF] have shown the up-regulation of enzymes involved in glycolysis and gluconeogenesis in 2 weeks R6/2 mice. Using the same model, Perluigi et al. [START_REF] Perluigi | Proteomic analysis of protein expression and oxidative modification in R6/2 transgenic mice[END_REF] have demonstrated an increasing expression level of dihydrolipoamide S-succinyl-transferase and aspartate aminotransferase over the course of disease (10-week-old mice) and a decrease in pyruvate dehydrogenase expression level in 10-week-old HD transgenic mice compared with young (4-week-old) mice. In addition, expression of mhtt may affect transcriptional activity and may result in alteration in expression of proteins involved in energy metabolism as demonstrated by Naia et al. [START_REF] Naia | Histone Deacetylase Inhibitors Protect Against Pyruvate Dehydrogenase Dysfunction in Huntington's Disease[END_REF] who found in striatal HdhQ 111/111 cells an increase of pyruvate dehydrogenase kinase expression leading to a decrease of pyruvate dehydrogenase activity.

CONCLUSION

In HD, changes in metabolite levels precede structural brain alterations. Therefore, the study of molecular biomarkers at a presymptomatic stage is essential to understand better the pathogenesis and to find new therapeutic targets. Our work identifies 16 discriminatory metabolites at presymptomatic and symptomatic stages of HD modeled in Drosophila melanogaster. These metabolites confirm dysfunctions in energy production leading to mitochondrial dysfunction and a disruption of glycolysis, TCA and glutamate-glutamine cycles. Our results suggest that a set of 15 metabolites characteristic of both stages of the disease could be used to distinguish HD flies from controls, indicating that they have a predictive value. As defined here on a Drosophila model, among this 15 metabolites, nine metabolites, namely NAD + , fumarate, taurine, ethanolamine, dimethylamine, -

Figure 1 .

 1 Figure 1. Fly crossing scheme used to generate the three sets of samples. Disease model

Figure 2 .

 2 Figure 2. Example of 1 H-NMR spectrum with the identified metabolites. This spectrum

Figure 3 .

 3 Figure 3. Discriminant variables at the presymptomatic stage (10-day-old flies). OPLS-

Table 1 .

 1 Metabolites identified from the 79 common discriminant variables with VIP >1, for the 10-day-old flies. Downward arrows (↓) indicate that the ratio of metabolite was lower for HD flies than for controls. s-singlet, d-doublet, dd-doublet of doublet, t-triplet, m-multiplet. ** 0.001< Value < 0.01; ***0.0001< Value < 0.001, **** Value < 0.0001.

  the associated two controls, ii) 136 variables between HD flies and Gal4 controls and finally iii) 191 variables between HD flies and UAS controls. All the OPLS-DA analyses performed lead to a good or excellent separation of the groups. Correct R2X (0.57; 0.71 and 0.50, respectively) and large R2Y (0.95; 0.99 and 0.98, respectively) indicated high reproducibility and a good fitting of the model to the data. Q2Y parameters were calculated by cross validation and were close to 1 (0.87; 0.89 and 0.97, respectively), which demonstrates the high predictivity of models. Furthermore, the discriminant variables produced by the three OPLS-DA analyses were compared to each other using a Venn diagram and the 46 common variables were retained.

Figure 4 .

 4 Figure 4. Discriminant variables at the symptomatic stage (16-day-old flies). OPLS-DA

Figure 5 .Table 3 .

 53 Figure 5. Discriminant variables over the course of the disease (10-and 16 day-old flies).

Figure S1 -

 S1 Figure S1 -13 C-HSQC spectra of 16-day-old HD flies.

Figure S2 -

 S2 Figure S2 -VIP score of the metabolites for the 10-day-old flies extracted from the

Figure S3 -

 S3 Figure S3 -VIP score of the metabolites for the 16-day-old flies extracted from the

Figure S4 -

 S4 Figure S4 -Plot of ROC curves for several biomarker models based on their average

Figure S5 -Figure S6 .

 S5S6 Figure S5 -Plot of the ROC curve for the 15 biomarker model

ACKNOWLEDGMENTS

The authors acknowledge the CNRS and the "Centre Val de Loire" region for their financial support.

alanine, succinate, phosphocholine and 2-oxoglutarate, shown to be independent of age, may be considered first to follow the progression of HD in flies. They could be considered as potential candidates to be tested in other HD models and in human patients.