# Metabolomic Profiling of Plasma and Erythrocytes in Sickle Mice Points to Altered Nociceptive Pathways

Klétigui Casimir Dembélé, Thomas Mintz, Charlotte Veyrat-Durebex, Floris Chabrun, Stéphanie Chupin, Lydie Tessier, Gilles Simard, Daniel Henrion, Delphine Mirebeau-Prunier, Juan Manuel Chao de La Barca, et al.

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

To cite this version:

Klétigui Casimir Dembélé, Thomas Mintz, Charlotte Veyrat-Durebex, Floris Chabrun, Stéphanie Chupin, et al.. *Metabolomic Profiling of Plasma and Erythrocytes in Sickle Mice Points to Altered Nociceptive Pathways*. Cells, 2020, 9 (6), pp.1334. 10.3390/cells9061334. [hal-02974118](https://hal.inistre.fr/hal-02974118)

HAL Id: hal-02974118

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

Submitted on 22 Dec 2023

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
<figure><img src="image_3.png" /></figure>

Article

# Metabolomic Profiling of Plasma and Erythrocytes in Sickle Mice Points to Altered Nociceptive Pathways

Klétigui Casimir Dembélé¹,²,³,†, Thomas Mintz ⁴,†, Charlotte Veyrat-Durebex ⁵,⁶,† ID,
Floris Chabrun ²,³ ID, Stéphanie Chupin ², Lydie Tessier ², Gilles Simard ², Daniel Henrion ³ ID,
Delphine Mirebeau-Prunier ²,³, Juan Manuel Chao de la Barca ²,³, Pierre-Louis Tharaux ⁴,⁷,⁸,* ID
and Pascal Reynier ²,³,* ID

1. Faculté de Pharmacie, Université des Sciences, des Techniques et des Technologies de Bamako BP, Bamako 1805, Mali; kletigui@outlook.fr

2. Département de Biochimie et Génétique, Centre Hospitalier Universitaire, 49933 Angers, France; Floris.Chabrun@chu-angers.fr (F.C.); Stephanie.Chupin@chu-angers.fr (S.C.); LyTessier@chu-angers.fr (L.T.); GiSimard@chu-angers.fr (G.S.); DePrunier@chu-angers.fr (D.M.-P.); JMChaoDeLaBarca@chu-angers.fr (J.M.C.d.l.B.)

3. Unité Mixte de Recherche (UMR) MITOVASC, Centre National de la Recherche Scientifique (CNRS) 6015, Institut National de la Santé et de la Recherche Médicale (INSERM) U1083, Université d'Angers, 49933 Angers, France; daniel.henrion@univ-angers.fr

4. Paris Cardiovascular Centre (PARCC), Institut National de la Santé et de la Recherche Médicale (INSERM), 75015 Paris, France; carmintz@gmail.com

5. Institut National de la Santé et de la Recherche Médicale (INSERM) U1253, Université François Rabelais de Tours, 37000 Tours, France; C.VEYRATDUREBEX@chu-tours.fr

6. Laboratoire de Biochimie et Biologie Moléculaire, Centre Hospitalier Universitaire de Tours, 37000 Tours, France

7. Université Paris Descartes, Sorbonne Paris Cité, 75006 Paris, France

8. Nephrology Division, Georges Pompidou European Hospital, Assistance Publique-Hôpitaux de Paris, 75015 Paris, France

* Correspondence: pierre-louis.tharaux@inserm.fr (P.-L.T.); pareynier@chu-angers.fr (P.R.)

† These authors contributed equally to this work.

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

check for
updates

Received: 7 May 2020; Accepted: 23 May 2020; Published: 26 May 2020

**Abstract:** Few data-driven metabolomic approaches have been reported in sickle cell disease (SCD) to date. We performed a metabo-lipidomic study on the plasma and red blood cells of a steady-state mouse model carrying the homozygous human hemoglobin SS, compared with AS and AA genotypes. Among the 188 metabolites analyzed by a targeted quantitative metabolomic approach, 153 and 129 metabolites were accurately measured in the plasma and red blood cells, respectively. Unsupervised PCAs (principal component analyses) gave good spontaneous discrimination between HbSS and controls, and supervised OPLS-DAs (orthogonal partial least squares-discriminant analyses) provided highly discriminant models. These models confirmed the well-known deregulation of nitric oxide synthesis in the HbSS genotype, involving arginine deficiency and increased levels of dimethylarginines, ornithine, and polyamines. Other discriminant metabolites were newly evidenced, such as hexoses, alpha-amino adipate, serotonin, kynurenine, and amino acids, pointing to a glycolytic shift and to the alteration of metabolites known to be involved in nociceptive pathways. Sharp remodeling of lysophosphatidylcholines, phosphatidylcholines, and sphingomyelins was evidenced in red blood cells. Our metabolomic study provides an overview of the metabolic remodeling induced by the sickle genotype in the plasma and red blood cells, revealing a biological fingerprint of altered nitric oxide, bioenergetics and nociceptive pathways.

**Keywords:** metabolomics; lipidomics; nociception; sickle cell disease
# 1. Introduction

Sickle cell disease (OMIM #603903) is the most common autosomal recessive disorder worldwide, caused by the homozygous pathogenic variant p.(Val6Glu) in the gene encoding the β-globin chain of hemoglobin (HBB) [1]. The increased tendency of hemoglobin S to polymerize and aggregate, enhanced in circumstances favoring hypoxia and hemoglobin deoxidation, increases the erythrocyte rigidity, producing the “sickled” red blood cells (RBC) with altered rheological properties. This leads to hemolysis, anemia and painful life-threatening vaso-occlusive crises with deleterious consequences in the organs, such as pulmonary hypertension, ischemic stroke, cutaneous leg ulceration, priapism, and renal insufficiency.

Metabolomics, the next generation of metabolic biochemistry, is a powerful hypothesis-free approach highlighting the biological imprinting of the diseases. Surprisingly, very few metabolomic studies have been performed in sickle cell disease so far [2,3]. The most documented study, performed by Darghouth et al. [4], compared the RBC from 28 adult patients with the HbSS genotype with those of 24 healthy adults (HbAA), using a chromatography-mass spectrometry based untargeted metabolomic method. Among the 89 metabolites identified, 31 exhibited significantly modified concentrations, revealing an involvement of glycolysis, pentose phosphate, glutathione, ascorbate, amino acids, polyamines, carnitine, and creatine metabolisms in the pathogenesis. Using metabolomic profiling in the whole blood and plasma of a mouse model of the disease, the Xia group progressively deciphered a central pathogenic mechanism, also found in humans, promoting the polymerization of the hemoglobin S [5,6]. They showed that increased concentrations of adenosine resulted, through the activation of the adenosine A2B receptor (ADORA2BR), in elevated concentrations of 2,3-biphosphoglycerate (2,3-DBG) and sphingosine-1-P (S1P). This activation was, however, pharmacologically reversible. The same group later showed that elevated concentrations of S1P mediated multi-tissue chronic inflammation through IL6 [7]. Lastly, a proton nuclear magnetic resonance (¹H-NMR)-based metabolomic approach was applied in the plasma of patients (n = 23) with SCD according to their level of albuminuria [8]. The concentration of six metabolites was modified according to albuminuria: betaine, proline, dimethylarginine, glutamate, leucine, and lysine.

In this study, we used a standardized and quantitative data-driven targeted metabolomics procedure to explore the metabolic fingerprints of the plasma and RBC in a SCD mouse model carrying the human mutated gene (HbSS), compared with wild type (HbAA) and heterozygous (HbAS) genotypes.

# 2. Materials and Methods

## 2.1. *Animals and Sampling Procedures*

All experiments were performed in accordance with the European Community Guiding Principles for the care and use of animals (Directive 2010/63/UE; Décret n°2013-118). Studies used humanized SCD knock-in mice, with notation B6;129-Hba<sup>tm1(HBA)</sup>Tow/Hbb<sup>tm2(HBG1,HBB*)</sup>Tow/Hbb<sup>tm3(HBG1,HBB)</sup>Tow/J, developed by Wu et al. [9]. Experimental animals (HbSS) were homozygous for mutant β globin (E to V at position 6) and expressed human HbS. Control animals (HbAS or HbAA), derived from the same colony, were heterozygous or homozygous for wild-type β globin and expressed human HbA. Our metabolomic analysis was carried out on a total of 56 mice, all studied in steady state, with an age comprised between 9-12-week-old. These mice were divided in two comparative cohorts, one comparing HbSS (n = 23; 8 males and 15 females) with HbAS (n = 14; 9 males and 5 females) and the other comparing HbSS (n = 9; 4 males and 5 females) with HbAA (n = 10; 6 males and 4 females) genotypes.

All animals were maintained in a room at constant temperature (21 ± 2 °C) with a 12 h light/dark cycle. They were fasting since the day before blood collection with free access to drinking water. 130 to 150 µL of blood were collected in the morning by the same trained operator from the mandibular vein into a heparinized tube.
Samples were centrifuged for 10 min at 5000× g and 4 °C. The supernatant (plasma) was recovered and immediately stored in 50 µL aliquots at -80 °C until analysis. None of the plasma studied showed apparent hemolysis. Then, 50 µL of the RBC pellet were collected and washed with 150 µL of NaCl 0.9%. RBCs cell pellets were immediately stored in 50 µL aliquots at -80 °C until analysis.

## 2.2. Red Blood Cell Extraction

RBC samples (50 µL) were mixed with a cold solution of PBS (15 µL) and methanol (85 µL) and transferred in pre-cooled 2.0 mL homogenization Precellys tubes prefilled with 1.4 mm diameter ceramic beads. Homogenization was performed using a Precellys homogenizer (Bertin Technologies, Montigny-le-Bretonneux, France) kept in a room at +4 °C. The supernatant (cell extract) was recovered after centrifuging the homogenate (10,000× g, 5 min at +4 °C) and kept at -80 °C until metabolomic analysis.

## 2.3. Metabolomic Analyses

Targeted quantitative metabolomic analysis was carried out with a Biocrates® Absolute IDQ p180 kit (Biocrates Life Sciences AG, Innsbruck, Austria) and a QTRAP 5500 mass spectrometer (SCIEX, Villebon-sur-Yvette, France) to quantify 188 metabolites. The full list of individual metabolites is available at http://www.biocrates.com/products/research-products/absoluteidq-p180-kit. Sample preparation and analyses were performed following the Kit User Manual.

## 2.4. Statistical Analyses

Metabolites with more than 20% of concentration values below the lower limit of quantitation (LLOQ) or above the upper limit of quantitation (ULOQ) were excluded from the raw metabolomic data before statistical analyses. Each metabolite concentration was normalized by the sum of the concentrations of all the metabolites, similar to normalization by the TIC (total ion current). Multivariate analysis was first performed using unsupervised principal component analysis (PCA) for the detection of sample grouping and outliers. Supervised orthogonal partial least squares-discriminant analysis (OPLS-DA) was then applied to maximize variation between the disease (HbSS) and control (HbAA and HbAS) groups, and to determine the metabolites contributing to this variation. Multivariate analyses for the PCA and OPLS-DA models were performed on mean-centered and unit variance-scaled (MC-UV) data. The quality of the OPLS-DA model was validated by two parameters, namely goodness of fit (R²) and goodness of prediction indicated by the cumulated Q² value (Q²cum). A threshold of 0.5 was used to determine whether an OPLS-DA model could be estimated as having a good (Q²cum ≥ 0.5) or a poor (Q²cum < 0.5) predictive capability [10]. VIP (variable importance for the projection) values summarize the importance of each variable for the OPLS-DA model, whereas the loading values are indicators of the relationship between the y vector containing the class information (i.e., HbSS or controls) and variables in the X matrix (i.e., metabolites). In the reduced model, plotting key variables based on VIP values versus loading values scaled as correlation coefficients ($p_{corr}$) yields a V-shaped graph known as the “volcano” plot, on which the most important variables are easily recognized. Multivariate data analysis was conducted using SIMCA-P v.14.0 (Umetrics, Umeå, Sweden).

# 3. Results

## 3.1. HbSS/HbAS Plasma Signature

Among the 188 metabolites analyzed by the kit, 153 were accurately measured in the plasma (the 5661 metabolites concentrations are given in Supplementary Materials, Table S1), including the hexoses (n = 1, sum of), amino acids (n = 21), biogenic amines (n = 15), carnitine, acylcarnitines (n = 12), sphingomyelins (n = 15), phosphatidylcholines (n = 75) and lysophosphatidylcholines (n = 13). Unsupervised PCA gave a good spontaneous separation of the two genotypes (R2X = 0.821; Q2cum = 0.637; Figure 1A). Supervised OPLS-DA provided a highly discriminant model for the
separation of the two groups (R2X = 0.71, R2Y = 0.945, Q₂ = 0.797, Figure 1B) with a low risk of over-fitting (CV-ANOVA *p*-value = 7.16 × 10⁻⁸).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** (plasma). PCA (A) and OPLS-DA (B) scatter plots obtained from the matrix of metabolites measured in the plasma from HbSS (blue circles, *n* = 23) and HbAS (green circles, *n* = 14) mice.</figcaption></figure>

This multivariate model included 29% (n = 44) of the accurately measured metabolites with a VIP > 1 (Figure 2). The values of VIPs and loadings are given in Table S2. These most discriminant metabolites comprised the hexoses (n = 1, sum of), nine amino acids (glycine, ornithine, proline, arginine, histidine, methionine, aspartate, alanine, and glutamine), seven biogenic amines (ADMA, total DMA, serotonin, putrescine, SDMA, kynurenine, and spermidine), six acylcarnitines (C3, C14, C16, C16:1, C18, and C18:1), five sphingomyelins, and 16 phosphatidylcholines. The sense of the modification of the concentrations of these metabolites are given in the volcano plot presented in Figure 2.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** (plasma). Volcano plot ($p_{corr}$ vs. VIP) from the OPLS-DA model shown in Figure 1B (plasma from HbSS vs. HbAS). Only the most discriminating metabolites with high VIP values ≥ 1 are listed in the boxes. Negative $p_{corr}$ values (left) indicate diminished metabolite concentrations in HbSS versus HbAS genotypes, whereas positive $p_{corr}$ values (right) indicate increased metabolite concentrations in HbSS versus HbAS genotypes.</figcaption></figure>

## 3.2. HbSS/HbAS Red Blood Cells Signature

Among the 188 metabolites analyzed by the kit, 129 were accurately measured in the RBC extracts (the 4773 metabolites concentrations are given in Table S3), including the hexoses ($n = 1$, sum of), amino acids ($n = 19$), biogenic amines ($n = 13$), carnitine, acylcarnitines ($n = 8$), sphingomyelins ($n = 14$), phosphatidylcholines ($n = 61$), and lysophosphatidylcholines ($n = 12$). Unsupervised PCA gave a good spontaneous separation of the two genotypes (R2X = 0.859; $Q_2$ = 0.678; Figure 3A). Supervised OPLS-DA provided a highly discriminant model separation of the two groups (R2X = 0.685; R2Y = 0.898; $Q_2$ = 0.787; Figure 3B) with a low risk of over-fitting (CV-ANOVA $p$-value = 7.42 × 10⁻⁹).

This multivariate model included 48% ($n = 62$) of the accurately measured metabolites with a VIP > 1 (Figure 4). The values of VIPs and Loadings are given in Table S2. The most discriminant metabolites comprised the hexoses ($n = 1$, sum of), six amino acids (glycine, glutamine, aspartate, proline, asparagine, and glutamate), four biogenic amines (alpha-AAA, putrescine, taurine, ADMA), carnitine, C3 acylcarnitine, 13 sphingomyelins, 33 phosphatidylcholines, and three lysophosphatidylcholines. The sense of the modification of the concentrations of these metabolites is given in the volcano plot presented in Figure 4.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** (RBC). PCA (A) and OPLS-DA (B) scatter plots obtained from the matrix of metabolites measured in the red blood cells from HbSS (blue circles, *n* = 23) and HbAS (green circles, *n* = 14) mice.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** (RBC). Volcano plot (*p*_corr vs. VIP) from the OPLS-DA model shown in Figure 3B (red blood cells from HbSS vs. HbAS). Only the most discriminating metabolites with high VIP values ≥ 1 are listed in the boxes. Negative *p*_corr values (left) indicate diminished metabolite concentrations in HbSS versus HbAS genotypes, whereas positive *p*_corr values (right) indicate increased metabolite concentrations in HbSS versus HbAS genotypes.</figcaption></figure>
### 3.3. Summary of the HbSS/HbAS Signature

Whole data visualization is provided in Figure 5 using the word clouds summarizing the plasma and red blood cell signatures of the HbSS/HbAS comparison.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** HbSS/HbAS signature summarization. Word clouds were calculated using metabolite VIPs to determine the size of the words, and loadings to determine the color scale, from dark green for the highest negative loadings (decreased in HbSS) and from deep red for the highest positive loadings (increased in HbSS).</figcaption></figure>

### 3.4. HbSS/HbAA Plasma Signature

Among the 188 metabolites analyzed by the kit, 131 were accurately measured in the plasma (the 2489 metabolites concentrations are given in Table S4). Unsupervised PCA gave a good spontaneous separation of the two genotypes (R2X = 0.777; Q₂ = 0.587; Figure 6A). Supervised OPLS-DA provided a highly discriminant model for separating the two groups (R2X = 0.766; R2Y = 0.917; Q₂ = 0.776; Figure 6B) with a low risk of over-fitting (CV-ANOVA p-value = 0.0023).

<figure><img src="image_10.png" /><figcaption>**Figure 6.** (plasma). PCA (A), OPLS-DA (B) scatter plots obtained from the matrix of metabolites measured in the plasmas from HbSS (blue circles, *n* = 9) and HbAA (green circles, *n* = 10) mice.</figcaption></figure>
This multivariate model included 47% (n = 62) of the accurately measured metabolites with a VIP > 1 (Figure 7). The value of VIPs and loadings are given in Table S2. The most discriminant metabolites comprised two amino acids (arginine and threonine), three biogenic amines (ADMA, putrescine and kynurenine), four sphingomyelins, 45 phosphatidylcholines, and eight lysophosphatidylcholines. The sense of the modification of the concentrations of these metabolites is given in the volcano plot presented in Figure 7.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** (plasma). Volcano scatter plots obtained from the matrix of metabolites measured in the plasmas from HbSS (blue circles, *n* = 9) and HbAA (green circles, *n* = 10) mice. Volcano plot (*p*<sub>corr</sub> vs. VIP) from the OPLS-DA model shown in B (plasma from HbSS vs. HbAA). Only the most discriminating metabolites with high VIP values ≥ 1 are listed in the boxes. Negative *p*<sub>corr</sub> values (left) indicate diminished metabolite concentrations in HbSS versus HbAA genotypes, whereas positive *p*<sub>corr</sub> values (right) indicate increased metabolite concentrations in HbSS versus HbAA genotypes.</figcaption></figure>

3.5. HbSS/HbAA Red Blood Cell Signature

Among the 188 metabolites analyzed by the kit, 113 were accurately measured in the RBC (the 2147 metabolites concentrations are given in Table S5). Unsupervised PCA gave a good spontaneous separation of the two genotypes (R2X = 0.86; Q2 = 0.749; Figure 8A). Supervised OPLS-DA provided a highly discriminant model separation of the two groups (R2X = 0.844; R2Y = 0.908; Q2 = 0.728; Figure 8B) with a low risk of over-fitting (CV-ANOVA *p*-value = 0.0066).
<figure><img src="image_12.png" /><figcaption>**Figure 8.** (RBC). PCA (A), OPLS-DA (B) scatter plots obtained from the matrix of metabolites measured in the red blood cells from HbSS (blue circles, *n* = 9) and HbAA (green circles, *n* = 10) mice.</figcaption></figure>

This multivariate model included 64% (n = 72) of the accurately measured metabolites with a VIP > 1 (Figure 9). The values of VIPs and loadings are given in Table S2. The most discriminant metabolites comprised 3 amino acids (arginine, histidine, and aspartate), five biogenic amines (ADMA, putrescine, alpha-AAA, spermine, spermidine), carnitine and acylcarnitine C3, 13 sphingomyelins, 44 phosphatidylcholines, and six lysophosphatidylcholines. The sense of the modification of the concentrations of these metabolites is given in the volcano plot presented in Figure 9.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** (RBC). Volcano scatter plots obtained from the matrix of metabolites measured in the red blood cells from HbSS (blue circles, *n* = 9) and HbAA (green circles, *n* = 10) mice. Volcano plot (*p*corr vs. VIP) from the OPLS-DA model shown in B (red blood cells from HbSS vs. HbAA). Only the most discriminating metabolites with high VIP values ≥ 1 are listed in the boxes. Negative *p*corr values (left) indicate diminished metabolite concentrations in HbSS versus HbAA genotypes, whereas positive *p*corr values (right) indicate increased metabolite concentrations in HbSS versus HbAA genotypes.</figcaption></figure>

**4. Discussion**

The contribution of nitric oxide (NO) to the vaso-occlusive process in SCD is well established [11,12]. NO is involved in the regulation of vaso-motor tone as well as in platelet aggregation and the adhesion of sickle blood cells to the activated endothelium [12]. Arginine is the precursor of NO, through NO synthases (NOS), but it can also be metabolized to ornithine and urea, through arginase activity, as the ornithine itself is a substrate for polyamines (putrescine, spermine and spermidine), proline and glutamate synthesis. SCD is characterized by a metabolic shift in favor of arginase, to the detriment of NO synthesis, due to the release of arginase from erythrocytes during hemolysis, as well as to inflammatory cytokines which reduce NO availability [12]. Our metabolomic signature is in full agreement with this metabolic shift, showing decreased concentration in the HbSS genotype of arginine in the plasma and red blood cells, whereas the direct or indirect products of arginase, such as ornithine, proline, glutamate, putrescine, spermidine, and spermine, were found to be increased in the plasma. In the HbSS red blood cells, although the putrescine was also increased, spermidine and spermine were diminished. Polyamines are known to be involved in SCD pathophysiology through their association with spectrin, which stabilizes the red blood cell membranes [13].

ADMA (asymmetric NG, NG-dimethylarginine), SDMA (symmetric dimethylarginine) and total
DMA (dimethylarginine) were globally also increased in the HbSS signatures as it has been previously
reported in the blood of patients [14]. Both ADMA and SDMA derive from intranuclear methylation of
l-arginine residuals and are released into the cytoplasm after proteolysis. Presumably, SDMA is a weak
inhibitor of eNOS, most likely weaker than ADMA. SDMA is also a potent competitor of l-arginine
transport [15], and it impairs l-arginine uptake from the perfused loop of Henle [16]. This suggests that
SDMA may compete with arginine transport and thereby still could have an indirect inhibitory effect
of NO synthesis by limiting arginine availability to NOS. In the general population, SDMA, but not
ADMA, was associated with cardiovascular mortality and all causes mortality [17]. Circulating ADMA
was found to be elevated in several conditions of endothelial dysfunctions [18], as a biomarker of poor prognostic. ADMA being an endogenous NOS inhibitor, in SCD, its increased concentration reinforces the NO deficiency [19]. These results could be due to the balance between increased production through hemolysis and haemoglobin proteolysis as well as to impaired clearance.

NO, arginine, ornithine and ADMA also play a role in pain pathophysiology and the restoration of arginine bioavailability through exogenous supplementation was shown to contribute to the treatment of the disease by increasing the plasma NO concentration and reducing the pain [20,21]. Polyamines are also involved in pain since their increased synthesis related to inflammation was shown to be involved in the development of nociception, a subcutaneous injection of putrescine, spermidine or spermine into the rat paw showing mechanical allodynia and edema [22]. In contrast, a polyamine diet has been shown to relieve pain hypersensitivity in rats [23].

In addition to this metabolomic signature centered on arginine metabolism, we also found that many other amino acids disclosed increased concentrations in the HbSS genotype, mostly in the plasma, such as aspartate, asparagine, alanine, glycine, glutamate, methionine, threonine, glutamine, and taurine. RBC of SCD patients exhibit an abnormally high abundance of N-methyl D-aspartate receptors (NMDAR), mediating excessive calcium uptake which could contribute to the crisis [24]. Interestingly, these authors proposed that glutamate may be found to be increased in this context. We here found that the glutamate, as well as the other NMDAR agonists, aspartate and glycine were all increased in HbSS. In addition to its direct role in sickle erythrocytes, NMDARs, especially those located in the dorsal horn of the spinal cord, are critically involved in nociceptive transmission and central sensitization of pain [25], which could also be impaired in the disease. Taurine, which contributes to the sensitivity to nociceptive chemical stimulation [26], is another amino acid derivative found in our sickle signatures.

Two mediators derived from the essential amino acid tryptophan showed altered concentrations, the kynurenine that was increased in the two plasma signatures and the serotonin that was reduced in the plasma of the HbSS/HbAS signature. These modifications of kynurenine and serotonin have never been reported before in our knowledge in SCD. Kynurenine metabolites such as quinolinic acid, can contribute to the development of depression via NMDA glutamatergic stimulation. The kynurenine is also an immune metabolite involved in inflammation and vasodilatation [27,28]. Interestingly, it was shown that the kynurenine metabolic pathway is an important mediator of neuropathic pain in rats that constitute a therapeutic target [29].

Serotonin or 5-hydroxytryptamine (5HT) is an intermediate product of tryptophan metabolism located primarily in the enterochromaffin cells of the intestines, brain serotonergic neurons, and blood platelets. In tissues such as the brain, inflammation increases tryptophan metabolism, thereby depleting its availability for serotonin production [30]. Serotonin is an important endogenous mediator of the descending anti-nociceptive system in the central nervous system. Indeed, modulation of painful stimuli includes the release of serotonin to inhibit pain transmission. Serotonin modulates the descending fibers from the midbrain to the dorsal horn that inhibit the transmission of the painful stimuli. In mast cells from mouse models expressing human sickle hemoglobin, there was evidence of decreased serotonin released under chronic morphine exposure [31]. Sickle hyperalgia was shown to be modulated by the descending serotonergic anti-nociceptive system in HbSS-BERK (sickle) mice, while the serotonin concentration in the spinal cord dorsal is influenced by diet and correlated with hyperalgia [32]. Serotonin is also known to enhance the interaction of platelets with circulating macrovesicles, to increase platelet activation, and to increase their overall pro-coagulant activity [33]. Taking into account the importance of coagulation and nociception mechanisms in the pathophysiology of the disease, this decreased concentration of serotonin observed in the basal state of the mouse model could either counteract the pro-coagulant state or attest to the nociception deregulation.

The concentration of the sum of hexoses was found to be reduced in the plasma and red blood cells HbSS genotype compared with HbAS. This may reflect a preferential shift to glycolytic metabolism, secondary to the low oxygen levels in tissues due to the reduced oxygen-carrying capacity of HbS. Indeed, in addition, six long-chain acylcarnitines, that are also substrates for mitochondrial oxidation,
were also found to be increased in the plasma of HbSS compared with HbAS genotypes. Such accumulation of acylcarnitines has been reported in the two previously reported metabolomic studies performed in RBC in patients [4] and transgenic mice [6]. This glycolytic shift could be associated with slightly increased insulin sensitivity in this mouse model. Indeed, we observed a reduced concentration of alpha-AAA (alpha-amino adipate) in the two RBC signatures, which was recently identified as a biomarker of insulin resistance [34]. Alpha-AAA is involved in the lysine biosynthetic pathway and has been poorly investigated in pathophysiological conditions. Its reduced concentration, combined with the low levels of the hexoses and increased long chain acylcarnitines, argues in favor of higher insulin sensitivity in this mouse model. Such increased insulin sensitivity has recently been reported in sickle patients [35]. In addition, alpha-amino adipate is known to reduce nerve injury-induced neuropathic pain in mice through glial signaling [36], indicating its potential additional role in nociception.

Carnitine, known for its antioxidant properties and propionylcarnitine (C3) were increased in the two RBC signatures, as previously reported in RBC from patients [4] and mice [6]. Carnitine and acetyl-carnitine are known to be used as a buffer and reservoir of activated acyl group for the turnover and the repair of RBC membranes during repeated cycled of erythrocyte sickling and unsickling and mechanical fragmentation under shear stress participating in HbS cell membrane fragility [37].

Lastly, the concentrations of many lysophosphatidylcholines, phosphatidylcholines, and sphingomyelins were sharply reduced by the HbSS genotype, both in the plasma and RBC. Reduced concentrations of sphingomyelin and ceramide have already been reported in the serum of patients with SCD [38]. Phosphatidylcholines were also found to be reduced in the erythrocytes from the mouse model [8]. However, in contrast with our results lysophosphatidylcholines were found to be increased in this study. This lipidomic signature could therefore reflect the remodeling of sickle RBC. It could also be involved in pain since an alteration of the phosphatidylcholines, sphingomyelins, and ceramides metabolism was also found by a metabolomic approach in the dorsal horn during chronic neuropathic pain in a rat model [39].

It is interesting to note that this murine model reveals metabolomic features closely related to the signatures observed in patients with sickle cell anemia. We recently reported the metabolomic signatures observed in the plasma and red blood cells of patients during the vaso-occlusive crisis, in comparison with the steady state [40]. Globally, the same groups of metabolites are affected in the patients and mouse model. In the plasma, the signature involving the metabolism of NO and arginine is also in the foreground in patients during the vaso-occlusive crisis, as is the energetic metabolism, the decrease concentration of serotonin and the metabolites associated with pain. In red blood cells, the phospholipid remodeling is important both in humans and mice. However, during the vaso-occlusive crisis in humans, the phosphatidylcholines and lysophosphatidylcholines are mostly increased, while they are drastically reduced in the mouse model in steady state. It is likely that this reflects the sickling of red blood cells occurring during the crisis." It would be also interesting to explore the sexual dimorphism of these sickle mice metabolomic signatures, which may explain the sex-related specificity of the clinical expression that has been reported [41].

**5. Conclusions**

In conclusion, this quantitative targeted metabo-lipidomic approach performed in plasma and RBC of a mouse model carrying the human HbSS variant revealed new deregulated pathways affecting amino acids, kynurenine, serotonin, alpha-amino adipate, hexoses, and complex lipids. These metabolic alterations associated with SCD point towards a glycolytic shift and increased insulin sensitivity in this model and to a complex set of metabolites contributing to nociception (NO, arginine, ornithine, ADMA, polyamines, glutamate, aspartate, glycine, taurine, alpha-amino adipate, kynurenine and serotonin), while the lysophosphatidylcholines, phosphatidylcholines, and sphingomyelins are potentially involved in both pain pathophysiology and RBC remodeling.
cells HbSS compared with HbAS, Table S4: Plasma HbSS compared with HbAA, Table S5: Red Blood Cells HbSS compared with HbAA.

**Author Contributions:** Designed the research: P.-L.T., P.R.; performed the research: K.C.D., T.M., C.V.-D., F.C., S.C., L.T., J.M.C.d.I.B.; analyzed the data and wrote the manuscript: K.C.D., C.V.-D., J.M.C.d.I.B., G.S., D.H., D.M.-P., P.-L.T., P.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** We acknowledge support from the School of Medicine, Pharmacy and Odontology of Bamako, the French Embassy of Bamako and Campus France (for the thesis grant accorded to Klétigui Casimir Dembélé), the Institut National de la Santé et de la Recherche Médicale (INSERM), the Centre National de la Recherche Scientifique (CNRS), the University of Angers, and the University Hospital of Angers.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Ingram, V.M. Gene mutations in human haemoglobin: The chemical difference between normal and sickle cell haemoglobin. *Nature* **1957**, *180*, 326-328. [CrossRef] [PubMed]

2. Goodman, S.R.; Pace, B.S.; Hansen, K.C.; D'alessandro, A.; Xia, Y.; Daescu, O.; Glatt, S.J. Minireview: Multiomic candidate biomarkers for clinical manifestations of sickle cell severity: Early steps to precision medicine. *Exp. Biol. Med. (Maywood)* **2016**, *241*, 772-781. [CrossRef] [PubMed]

3. Adebiyi, M.G.; Manalo, J.M.; Xia, Y. Metabolomic and molecular insights into sickle cell disease and innovative therapies. *Blood Adv.* **2019**, *3*, 1347-1355. [CrossRef] [PubMed]

4. Darghouth, D.; Koehl, B.; Madalinski, G.; Heilier, J.F.; Bovee, P.; Xu, Y.; Olivier, M.F.; Bartolucci, P.; Benkerrou, M.; Pissard, S.; et al. Pathophysiology of sickle cell disease is mirrored by the red blood cell metabolome. *Blood* **2011**, *117*, e57-e66. [CrossRef] [PubMed]

5. Zhang, Y.; Dai, Y.; Wen, J.; Zhang, W.; Grenz, A.; Sun, H.; Tao, L.; Lu, G.; Alexander, D.C.; Milburn, M.V.; et al. Detrimental effects of adenosine signaling in sickle cell disease. *Nat. Med.* **2011**, *17*, 79-86. [CrossRef]

6. Zhang, Y.; Berka, V.; Song, A.; Sun, K.; Wang, W.; Zhang, W.; Ning, C.; Li, C.; Zhang, Q.; Bogdanov, M.; et al. Elevated sphingosine-1-phosphate promotes sickling and sickle cell disease progression. *J. Clin. Invest.* **2014**, *124*, 2750-2761. [CrossRef]

7. Zhao, S.; Adebiyi, M.G.; Zhang, Y.; Couturier, J.P.; Fan, X.; Zhang, H.; Kellems, R.E.; Lewis, D.E.; Xia, Y. Sphingosine-1-phosphate receptor 1 mediates elevated IL-6 signaling to promote chronic inflammation and multitissue damage in sickle cell disease. *FASEB J.* **2018**, *32*, 2855-2865. [CrossRef]

8. Elsherif, L.; Pathmasiri, W.; McRitchie, S.; Archer, D.R.; Ataga, K.I.; Townes, T.M. Plasma metabolomics analysis in sickle cell disease patients with albuminuria-an exploratory study. *Br. J. Haematol.* **2019**, *185*, 620-623.

9. Wu, L.C.; Sun, C.W.; Ryan, T.M.; Pawlik, K.M.; Ren, J.; Townes, T.M. Correction of sickle cell disease by homologous recombination in embryonic stem cells. *Blood* **2006**, *108*, 1183-1188. [CrossRef]

10. Eriksson, L.; Johansson, E.; Kettaneh-Wold, N.; Trygg, J.; Wikström, C.; Wold, S. Part I: *Basic Principles and Applications PLS. Multi- and Magavariate Data Analysis*, 2nd ed.; Umetrics Ed: Umea, Sweden, 2006; p. 85.

11. Reiter, C.D.; Gladwin, M.T. An emerging role for nitric oxide in sickle cell disease vascular homeostasis and therapy. *Curr. Opin. Hematol.* **2003**, *10*, 99-107. [CrossRef]

12. Schnog, J.J.; Jager, E.H.; van der Dijs, F.P.; Duits, A.J.; Moshage, H.; Muskiet, F.D.; Muskiet, F.A. Evidence for a metabolic shift of arginine metabolism in sickle cell disease. *Ann. Hematol.* **2004**, *83*, 371-375. [CrossRef]

13. Natta, C.L.; Kremzner, L.T. Polyamines and membrane proteins in sickle cell disease. *Blood Cells* **1982**, *8*, 273-280. [PubMed]

14. Kato, G.J.; Wang, Z.; Machado, R.F.; Blackwelder, W.C.; Taylor, J.G., 6th; Hazen, S.L. Endogenous nitric oxide synthase inhibitors in sickle cell disease: Abnormal levels and correlations with pulmonary hypertension, desaturation, haemolysis, organ dysfunction and death. *Br. J. Haematol.* **2009**, *145*, 506-513. [CrossRef] [PubMed]

15. Tsikas, D. Does the inhibitory action of asymmetric dimethylarginine (ADMA) on the endothelial nitric oxide synthase activity explain its importance in the cardiovascular system? The ADMA paradox. *J Controv. Biomed. Res.* **2017**, *3*, 16-22. [CrossRef]
16. Closs, E.I.; Basha, F.Z.; Habermeier, A.; Forstermann, U. Interference of L-arginine analogues with L-arginine transport mediated by the y+ carrier hCAT-2B. *Nitric Oxide* **1997**, *1*, 65–73. [CrossRef]

17. Gore, M.O.; Lüneburg, N.; Schwedhelm, E.; Ayers, C.R.; Anderssohn, M.; Khera, A.; Atzler, D.; de Lemos, J.A.; Grant, P.J.; McGuire, D.K.; et al. Symmetrical dimethylarginine predicts mortality in the general population: Observations from the Dallas heart study. *Arterioscler. Thromb. Vasc. Biol.* **2013**, *33*, 2682–2688. [CrossRef]

18. Tripepi, G.; Mattace, R.F.; Sijbrands, E.; Seck, M.S.; Maas, R.; Boger, R.; Witteman, J.; Rapisarda, F.; Malatino, L.; Mallamaci, F.; et al. Inflammation and asymmetric dimethylarginine for predicting death and cardiovascular events in ESRD patients. *Clin. J. Am. Soc. Nephrol.* **2011**, *6*, 1714–1721. [CrossRef]

19. Schnog, J.B.; Teerlink, T.; van der Dijs, F.P.; Duits, A.J.; Muskiet, F.A.; CURAMA Study Group. Plasma levels of asymmetric dimethylarginine (ADMA), an endogenous nitric oxide synthase inhibitor, are elevated in sickle cell disease. *Ann. Hematol.* **2005**, *84*, 282–286. [CrossRef]

20. Bakshi, N.; Morris, C.R. The role of the arginine metabolome in pain: Implications for sickle cell disease. *J. Pain Res.* **2016**, *9*, 167–175.

21. Eleutério, R.M.N.; Nascimento, F.O.; Araújo, T.G.; Castro, M.F.; Filho, T.P.A.; Filho, P.A.M.; Eleutério, J., Jr.; Elias, D.B.D.; Lemes, R.P.G. Double-Blind Clinical Trial of Arginine Supplementation in the Treatment of Adult Patients with Sickle Cell Anaemia. *Adv. Hematol.* **2019**, *2019*, 4397150. [CrossRef]

22. Silva, M.A.; Klafke, J.Z.; Rossato, M.F.; Gewehr, C.; Guerra, G.P.; Rubin, M.A.; Ferreira, J. Role of peripheral polyamines in the development of inflammatory pain. *Biochem. Pharmacol.* **2011**, *82*, 269–277. [CrossRef]

23. Rivat, C.; Richebé, P.; Laboureyras, E.; Laulin, J.P.; Havouis, R.; Noble, F.; Moulinoux, J.P.; Simonnet, G. Polyamine deficient diet to relieve pain hypersensitivity. *Pain* **2008**, *137*, 125–137. [CrossRef] [PubMed]

24. Hänggi, P.; Makhro, A.; Gassmann, M.; Schmugge, M.; Goede, J.S.; Speer, O.; Bogdanova, A. Red blood cells of sickle cell disease patients exhibit abnormally high abundance of N-methyl D-aspartate receptors mediating excessive calcium uptake. *Br. J. Haematol.* **2014**, *167*, 252–264. [CrossRef] [PubMed]

25. Zhou, H.Y.; Chen, S.R.; Pan, H.L. Targeting N-methyl-D-aspartate receptors for treatment of neuropathic pain. *Expert Rev. Clin. Pharmacol.* **2011**, *4*, 379–388. [CrossRef] [PubMed]

26. Lötsch, J.; Hummel, T.; Warskulat, U.; Coste, O.; Häussinger, D.; Geisslinger, G.; Tegeder, I. Congenital taurine deficiency in mice is associated with reduced sensitivity to nociceptive chemical stimulation. *Neuroscience* **2014**, *259*, 63–70. [CrossRef] [PubMed]

27. Van der Leek, A.P.; Yanishevsky, Y.; Kozyrskyj, A.L. The Kynurenine Pathway as a Novel Link between Allergy and the Gut Microbiome. *Front. Immunol.* **2017**, *8*, 1374. [CrossRef]

28. Wang, Y.; Liu, H.; McKenzie, G.; Witting, P.K.; Stasch, J.P.; Hahn, M.; Changsirivathanathamrong, D.; Wu, B.J.; Ball, H.J.; Thomas, S.R.; et al. Kynurenine is an endothelium-derived relaxing factor produced during inflammation. *Nat. Med.* **2010**, *16*, 279–285.

29. Rojewska, E.; Ciapała, K.; Piotrowska, A.; Makuch, W.; Mika, J. Pharmacological Inhibition of Indoleamine 2,3-Dioxygenase-2 and Kynurenine 3-Monooxygenase, Enzymes of the Kynurenine Pathway, Significantly Diminishes Neuropathic Pain in a Rat Model. *Front. Pharmacol.* **2018**, *9*, 724. [CrossRef]

30. Ballas, S.K. Current issues in sickle cell pain and its management. *Hematol. Am. Soc. Hematol. Educ. Program.* **2007**, *1*, 97–105. [CrossRef]

31. Manning, B.M.; Hebbel, R.P.; Gupta, K.; Haynes, C.L. Carbon-fiber microelectrode amperometry reveals sickle-cell-induced inflammation and chronic morphine effects on single mast cells. *ACS Chem. Biol.* **2012**, *7*, 543–551. [CrossRef]

32. Tran, H.; Jarrett, S.; Nguyen, J.; Badgaiyan, R.; Gupta, K. Activation of the Central Serotonergic System Reduces Hyperalgesia in Sickle Mice. *Blood* **2016**, *128*, 266. [CrossRef]

33. Lopez-Vilchez, I.; Diaz-Ricart, M.; White, J.G.; Escolar, G.; Galan, A.M. Serotonin enhances platelet procoagulant properties antheir activation induced during platelet tissue factor upake. *Cardiovasc. Res.* **2009**, *84*, 309–316. [CrossRef]

34. Wang, T.J.; Ngo, D.; Psychogios, N.; Dejam, A.; Larson, M.G.; Vasan, R.S.; Ghorbani, A.; O'Sullivan, J.; Cheng, S.; Rhee, E.P.; et al. 2-Aminoadipic acid is a biomarker for diabetes risk. *J. Clin. Invest.* **2013**, *123*, 4309–4317. [CrossRef] [PubMed]

35. Akinlade, K.S.; Kumuyi, A.S.; Rahamon, S.K.; Olaniyi, J.A. Insulin Sensitivity, Inflammation, and Basal Metabolic Rate in Adults with Sickle Cell Anemia. *Int. J. Appl. Basic. Med. Res.* **2018**, *8*, 106–110. [CrossRef] [PubMed]
36. Chen, G.; Luo, X.; Qadri, M.Y.; Berta, T.; Ji, R.R. Sex-Dependent Glial Signaling in Pathological Pain: Distinct Roles of Spinal Microglia and Astrocytes. *Neurosci. Bull.* **2018**, *34*, 98–108. [CrossRef]

37. Arduini, A.; Mancinelli, G.; Radatti, G.L.; Dottori, S.; Molajoni, F.; Ramsay, R.R. Role of carnitine and carnitine palmitoyltransferase as integral components of the pathway for membrane phospholipid fatty acid turnover in intact human erythrocytes. *J. Biol. Chem.* **1992**, *267*, 12673–12681.

38. Aslan, M.; Kıraç, E.; Kaya, S.; Özcan, F.; Salim, O.; Küpesiz, O.A. Decreased Serum Levels of Sphingomyelins and Ceramides in Sickle Cell Disease Patients. *Lipids* **2018**, *53*, 313–322. [CrossRef]

39. Patti, G.J.; Yanes, O.; Shriver, L.P.; Courade, J.P.; Tautenhahn, R.; Manchester, M.; Siuzdak, G. Metabolomics implicates altered sphingolipids in chronic pain of neuropathic origin. *Nat. Chem. Biol.* **2012**, *8*, 232–234. [CrossRef]

40. Dembélé, K.C.; Veyrat-Durebex, C.; Aldiouma, G.; Chupin, S.; Goïta, Y.; Homedan, C.; Mirebeau-Prunier, D.; Cissé, B.; Simard, G.; Diallo, D.; et al. Sickle cell disease: Metabolomic profiles of vaso-occlusive crisis in plasma and erythrocytes. *J. Clin. Med.* **2020**, *9*, 1092. [CrossRef]

41. Ceglie, G.; Di Mauro, M.; Tarissi De Jacobis, I.; de Gennaro, F.; Quaranta, M.; Baronci, C.; Villani, A.; Palumbo, G. Gender-Related Differences in Sickle Cell Disease in a Pediatric Cohort: A Single-Center Retrospective Study. *Front Mol. Biosci.* **2019**, *6*, 140. [CrossRef]

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

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).