Article

# Preventing the Increase in Lysophosphatidic Acids: A New Therapeutic Target in Pulmonary Hypertension?

Thomas Duflot <fnref n="1" />,* , Ly Tu <fnref n="2" />,<fnref n="3" /> , Matthieu Leuillier <fnref n="4" />, Hind Messaoudi <fnref n="4" />, Déborah Groussard <fnref n="4" />, Guillaume Feugray <fnref n="5" />, Saïda Azhar <fnref n="4" />, Raphaël Thuillet <fnref n="2" />,<fnref n="3" />, Fabrice Bauer <fnref n="6" />, Marc Humbert <fnref n="2" />,<fnref n="3" /> , Vincent Richard <fnref n="1" />, Christophe Guignabert <fnref n="2" />,<fnref n="3" /> and Jérémy Bellien <fnref n="1" />

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

check for
updates

**Citation:** Duflot, T.; Tu, L.; Leuillier, M.; Messaoudi, H.; Groussard, D.; Feugray, G.; Azhar, S.; Thuillet, R.; Bauer, F.; Humbert, M.; et al. Preventing the Increase in Lysophosphatidic Acids: A New Therapeutic Target in Pulmonary Hypertension? *Metabolites* **2021**, *11*, 784. https://doi.org/10.3390/metabo11110784

Academic Editors:
Hector Gallart-Ayala,
Justine Bertrand-Michel and
Olivier Berdeaux

Received: 22 October 2021
Accepted: 11 November 2021
Published: 17 November 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

Copyright: © 2021 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 (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Cardiovascular diseases (CVD) are the leading cause of premature death and disability in humans that are closely related to lipid metabolism and signaling. This study aimed to assess whether circulating lysophospholipids (LPL), lysophosphatidic acids (LPA) and monoacylglycerols (MAG) may be considered as potential therapeutic targets in CVD. For this objective, plasma levels of 22 compounds (13 LPL, 6 LPA and 3 MAG) were monitored by liquid chromatography coupled with tandem mass spectrometry (HPLC/MS<fnref n="2" />) in different rat models of CVD, i.e., angiotensin-II-induced hypertension (HTN), ischemic chronic heart failure (CHF) and sugen/hypoxia(SuHx)-induced pulmonary hypertension (PH). On one hand, there were modest changes on the monitored compounds in HTN (LPA 16:0, 18:1 and 20:4, LPC 16:1) and CHF (LPA 16:0, LPC 18:1 and LPE 16:0 and 18:0) models compared to control rats but these changes were no longer significant after multiple testing corrections. On the other hand, PH was associated with important changes in plasma LPA with a significant increase in LPA 16:0, 18:1, 18:2, 20:4 and 22:6 species. A deleterious impact of LPA was confirmed on cultured human pulmonary smooth muscle cells (PA-SMCs) with an increase in their proliferation. Finally, plasma level of LPA(16:0) was positively associated with the increase in pulmonary artery systolic pressure in patients with cardiac dysfunction. This study demonstrates that circulating LPA may contribute to the pathophysiology of PH. Additional experiments are needed to assess whether the modulation of LPA signaling in PH may be of interest.

**Keywords:** lysophospholipids; lysophosphatidic acids; cardiovascular diseases; HPLC-MS/MS; rodent models; pulmonary hypertension; chronic heart failure; hypertension

## 1. Introduction

Cardiovascular diseases (CVD) are the leading cause of mortality and a major contributor to disability. In 2019, 523 million cases and 18.6 million deaths were reported corresponding to a 100% and a 50% increase since 1990, respectively. Since CVD remain the leading cause of disease burden in the world and still rise in all countries, the need to find new therapeutic target is of utmost importance [1]. CVD strongly rely to atherosclerosis

<footnote n="1">1 UNIROUEN, INSERM U1096, CHU Rouen, Department of Pharmacology, Normandie University, F-76000 Rouen, France; vincent.richard@univ-rouen.fr (V.R.); jeremy.bellien@chu-rouen.fr (J.B.)</footnote>

<footnote n="2">2 INSERM UMR_S 999, Hôpital Marie Lannelongue, F-92350 Le Plessis-Robinson, France; lyieng@gmail.com (L.T.); raphael.thuillet@inserm.fr (R.T.); marc.humbert@aphp.fr (M.H.); christophe.guignabert@inserm.fr (C.G.)</footnote>

<footnote n="3">3 School of Medicine, Université Paris-Saclay, Le Kremlin-Bicêtre, F-92290 Châtenay-Malabry, France</footnote>

<footnote n="4">4 UNIROUEN, INSERM U1096, Normandie University, F-76000 Rouen, France; leuillier.matthieu@gmail.com (M.L.); hindiimessaoudi@gmail.com (H.M.); deborah.groussard@laposte.net (D.G.); saida.azhar1@univ-rouen.fr (S.A.)</footnote>

<footnote n="5">5 UNIROUEN, INSERM U1096, CHU Rouen, Department of General Biochemistry, Normandie University, F-76000 Rouen, France; guillaume.feugray@chu-rouen.fr</footnote>

<footnote n="6">6 UNIROUEN, INSERM U1096, CHU Rouen, Department of Cardiology, Normandie University, F-76000 Rouen, France; fabrice.bauer@chu-rouen.fr</footnote>

<footnote>* Correspondence: thomas.duflot@chu-rouen.fr; Tel.: +33-2-32-88-84-91</footnote>
and both are characterized by an increase in oxidative stress and chronic low-grade in-
flammatory condition associated to lipid homeostasis dysregulation [2]. The effectiveness
of therapeutic strategies based on circulating lipids reduction such as statins or PCSK-
9 inhibitors confirms the close link between atherosclerosis, CVD and lipid metabolism
and it is now well admitted that those strategies significantly reduce long-term mortality
and morbidity, especially in the elderly [3–5]. Most lipids circulate through bloodstream,
either free or bound to other molecules including albumin or lipoproteins. Lipoproteins
are macromolecular complexes of lipids and proteins with a central hydrophobic core
containing non-polar lipids such as triglycerides and cholesterol esters surrounded by
a hydrophilic layer composed of phospholipids, free cholesterol and apolipoproteins.
CVD trigger the production of oxidized low-density lipoprotein (ox-LDL) and oxidized
lipoprotein(a) (oxLp(a)) which become highly enriched in oxidized phospholipids (ox-PL)
displaying pro-atherogenic properties [6–8]. Ox-PL are especially prone to hydrolysis by
the lipoprotein associated phospholipase A₂ (Lp-PLA₂), also known as serum platelet-
activating factor-acetyl hydrolase (PAF-AH), an inflammatory marker of CVD, leading
to the release of both oxidized fatty acids and lysophospholipids (LPL) [9,10]. These
compounds exert deleterious biological effects due to their interaction with nucleic acids,
phospholipids and proteins promoting atherosclerosis and cardiovascular events [11]. Inter-
estingly, bioactive LPL, and in particular lysophosphatidic acids (LPA) have potent effects
on vascular cells, promoting vacoconstrictor, pro-inflammatory and proliferating effects
that are known to contribute to atherosclerosis [12,13]. However, the role of LPL in the
pathophysiology of many other CVD has been poorly evaluated [14–17]. LPA metabolism
is complex with various anabolic and catabolic pathways. There are four major enzymatic
pathways for LPA production: the extracellular LPL-autotaxin (LPL-ATX) pathway (1),
the phosphatidic acid-phospholipase pathway (PA-PLA₁/PLA₂) (2), the monoacylglycerol
kinase (MAGK) pathway (3) and the *de novo* glycerophosphate acyltransferase (GPAT)
pathway (4). Amongst these pathways responsible for LPA production, the last three
are intracellular and reversible with the intervention of lysophosphatidic acid acyltrans-
ferases (LPAAT), lipid phosphate phosphatases (LPP) and lysophospholipases (lysoPL) for
pathway (2), (3) and (4), respectively (Figure 1) [18].

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Biochemical pathways of lysophosphatidic acid metabolism. LPA, substrates/products and enzymes are depicted in green, blue and orange, respectively. ATX: autotaxin, G3P: glycero-3-phosphate, GPAT: glycerophosphate acyltransferase, LPA: lysophosphatidic acid LPAAT: lysophosphatidic acid acyltransferases, LPL: lysophospholipids, LPP: lysophosphatidic acid phosphatase, lysoPL: lysophospholipase, MAG: monoacylglycerol, MAGK: monoacylglycerol kinase, PA: phosphatidic acid, PLA: phospholipase.</figcaption></figure>
Thus, this study aims to evaluate circulating LPL, LPA and MAG levels in different rat models of CVD, i.e., hypertension (HTN), heart failure (HF) and pulmonary hypertension (PH) in order to provide new insights to the interest of targeting LPA metabolism with pharmacological compounds.

## 2. Results

### 2.1. Angiotensin-II Induced Hypertension (Ang-II HTN)

As expected, administration of Ang-II for 28 days induced an increase in aldosterone level (242 ± 108 vs. 1030 ± 633 pM, *p* = 0.004; Figure 2A) and in systolic blood pressure (SBP) (117 ± 4 vs. 176 ± 11 mm Hg, *p* < 0.001; Figure 2B) compared to control rats. Ang-II HTN also led to a decrease in left ventricular fractional shortening (LVFS) (45 ± 3 vs. 26 ± 3%, *p* < 0.001; Figure 2C) and an increase in heart weight (1.54 ± 0.13 vs. 1.76 ± 0.13 g, *p* = 0.004, Figure 2D), showing the development of cardiac dysfunction and hypertrophy.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Evaluation of cardiovascular parameters in control and hypertensive (HTN) rats. (A) Aldosterone levels (pM); (B) Systolic blood pressure (SBP, mm Hg); (C) Left ventricular fractional shortening (LVFS, %) and (D) Heart weight (g). N = 8 for each group. *p*-values computed by Student's *t*-test.</figcaption></figure>
Amongst the 22 monitored compounds, four LPA species exhibited significant increase in the HTN group compared with control group (LPA(16:0), LPA(16:1), LPA(18:1) and LPA(20:4)). However, these increases were no longer statistically significant after Benjamini-Hochberg (BH) correction to control the false discovery rate (Table 1) [19].

<table><caption>**Table 1.** Impact of Ang-II HTN on LPA, LPL and MAG levels.</caption><thead><tr><th>Analyte</th><th>Control Group (N = 8)</th><th>HTN Group (N = 8)</th><th>p-Value</th><th>Adjusted p-Value <sup>a</sup></th></tr></thead><tbody><tr><td>LPA(16:0)</td><td>0.11 [0.11–0.12]</td><td>0.15 [0.13–0.22]</td><td><strong>0.0229</strong></td><td>0.241</td></tr><tr><td>LPA(18:0)<sup>b</sup></td><td>0.09 [0.08–0.10]</td><td>0.11 [0.10–0.14]</td><td>0.0832</td><td>0.341</td></tr><tr><td>LPA(18:1)</td><td>0.09 [0.09–0.10]</td><td>0.13 [0.10–0.14]</td><td><strong>0.0438</strong></td><td>0.241</td></tr><tr><td>LPA(18:2)<sup>b</sup></td><td>1.78 [1.57–1.97]</td><td>2.20 [1.56–2.46]</td><td>0.163</td><td>0.398</td></tr><tr><td>LPA(20:4)<sup>b</sup></td><td>1.91 [1.50–2.13]</td><td>2.56 [2.14–2.92]</td><td><strong>0.0308</strong></td><td>0.241</td></tr><tr><td>LPA(22:6)<sup>b</sup></td><td>0.10 [0.08–0.11]</td><td>0.14 [0.11–0.17]</td><td>0.0931</td><td>0.341</td></tr><tr><td>LPC(16:0)</td><td>16.7 [16.2–17.4]</td><td>17.3 [15.0–18.0]</td><td>0.568</td><td>0.781</td></tr><tr><td>LPC(16:1)<sup>b</sup></td><td>1.24 [1.10–1.40]</td><td>1.58 [1.40–1.77]</td><td><strong>0.0438</strong></td><td>0.241</td></tr><tr><td>LPC(18:0)<sup>b</sup></td><td>29.7 [27.3–32.9]</td><td>30.4 [28.7–31.7]</td><td>0.914</td><td>0.914</td></tr><tr><td>LPC(18:1)</td><td>7.24 [6.36–7.53]</td><td>7.22 [6.83–8.14]</td><td>0.343</td><td>0.629</td></tr><tr><td>LPC(18:2)<sup>b</sup></td><td>31.5 [28.5–32.6]</td><td>30.9 [28.2–32.0]</td><td>0.749</td><td>0.867</td></tr><tr><td>LPC(20:4)<sup>b</sup></td><td>11.1 [10.3–12.1]</td><td>12.6 [11.6–13.7]</td><td>0.148</td><td>0.398</td></tr><tr><td>LPC(22:6)<sup>b</sup></td><td>0.34 [0.31–0.43]</td><td>0.43 [0.38–0.46]</td><td>0.269</td><td>0.538</td></tr><tr><td>LPE(16:0)</td><td>0.54 [0.47–0.61]</td><td>0.62 [0.58–0.68]</td><td>0.118</td><td>0.371</td></tr><tr><td>LPE(18:0)<sup>b</sup></td><td>0.80 [0.76–0.93]</td><td>0.95 [0.89–1.00]</td><td>0.205</td><td>0.451</td></tr><tr><td>LPE(18:1)<sup>b</sup></td><td>0.22 [0.21–0.27]</td><td>0.25 [0.22–0.27]</td><td>0.376</td><td>0.636</td></tr><tr><td>LPE(18:2)<sup>b</sup></td><td>0.45 [0.44–0.52]</td><td>0.46 [0.41–0.49]</td><td>0.439</td><td>0.690</td></tr><tr><td>LPE(20:4)<sup>b</sup></td><td>0.16 [0.15–0.17]</td><td>0.17 [0.16–0.20]</td><td>0.836</td><td>0.876</td></tr><tr><td>LPE(22:6)<sup>b</sup></td><td>0.11 [0.09–0.12]</td><td>0.12 [0.10–0.15]</td><td>0.612</td><td>0.792</td></tr><tr><td>MAG(18:1)</td><td>0.13 [0.10–0.16]</td><td>0.12 [0.11–0.15]</td><td>0.805</td><td>0.876</td></tr><tr><td>MAG(18:2)<sup>b</sup></td><td>0.67 [0.57–1.00]</td><td>0.64 [0.58–0.79]</td><td>0.673</td><td>0.823</td></tr><tr><td>MAG(20:4)<sup>b</sup></td><td>0.02 [0.02–0.03]</td><td>0.02 [0.01–0.02]</td><td>0.508</td><td>0.745</td></tr></tbody></table>

## 2.2. Chronic Heart Failure (CHF)

Definitive coronary artery ligation-induced CHF was defined by the decrease in LVFS (54 ± 4 vs. 26 ± 5% for control and CHF groups, respectively, *p* < 0.001, Figure 3A), the increase in heart weight (1.48 ± 0.12 vs. 1.68 ± 0.15 g, *p* = 0.008, Figure 3B), which concerned both the left (1.13 ± 0.09 vs. 1.27 ± 0.13 g, *p* = 0.016, Figure 3C) and right (0.24 ± 0.04 vs. 0.30 ± 0.05 g, *p* = 0.009, Figure 3D) ventricles but no change in mean pulmonary artery pressure (17.1 ± 1.1 vs. 17.0 ± 1.4 mm Hg, *p* = 0.861).

Regarding LPA metabolism, LPA(16:0), LPE(16:0) and LPE(18:0) were increase in CHF rats compared to control rats, whereas LPC(18:1) was decrease. However, these differences were no longer significant after BH correction (Table 2).

<footnote><sup>a</sup> *p*-values adjusted according to Benjamini and Hochberg multiple testing. <sup>b</sup> Analytes without analytical standard. HTN: hypertension, LPA: lysophosphatidic acid, LPC: lysophosphatidylcholine, LPE: lysophosphatidylethanolamine, MAG: monoacylglycerol. *p*-values computed by Student's *t*-test. *p*-values < 0.05 are highlighted in bold.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Evaluation of cardiac parameters in control and chronic heart failure (CHF) rats. (A) Left ventricular fractional shortening (LVFS, %), (B) Heart weight (g), (C) Left ventricle weight (g) and (D) Right ventricle weight (g). N = 8 and 9 for Control and CHF, respectively. *p*-values computed by Student's *t*-test.</figcaption></figure>
<table><caption>**Table 2.** Impact of CHF on LPA, LPL and MAG levels.</caption><thead><tr><th>Analyte</th><th>Control Group (N = 8)</th><th>CHF Group (N = 9)</th><th>p-Value</th><th>Adjusted p-Value <sup>a</sup></th></tr></thead><tbody><tr><td>LPA(16:0)</td><td>0.058 [0.046–0.060]</td><td>0.082 [0.056–0.103]</td><td><strong>0.025</strong></td><td>0.186</td></tr><tr><td>LPA(18:0)<sup>b</sup></td><td>0.050 [0.048–0.057]</td><td>0.070 [0.065–0.080]</td><td>0.336</td><td>0.615</td></tr><tr><td>LPA(18:1)</td><td>0.069 [0.063–0.076]</td><td>0.067 [0.054–0.075]</td><td>0.757</td><td>0.779</td></tr><tr><td>LPA(18:2)<sup>b</sup></td><td>0.87 [0.77–0.95]</td><td>1.00 [0.88–1.06]</td><td>0.193</td><td>0.499</td></tr><tr><td>LPA(20:4)<sup>b</sup></td><td>1.24 [1.04–1.31]</td><td>1.21 [1.14–1.56]</td><td>0.204</td><td>0.499</td></tr><tr><td>LPA(22:6)<sup>b</sup></td><td>0.066 [0.062–0.070]</td><td>0.068 [0.061–0.084]</td><td>0.447</td><td>0.615</td></tr><tr><td>LPC(16:0)</td><td>13.0 [11.3–14.2]</td><td>12.3 [11.1–13.2]</td><td>0.302</td><td>0.615</td></tr><tr><td>LPC(16:1)<sup>b</sup></td><td>1.06 [0.98–1.37]</td><td>0.93 [0.71–1.60]</td><td>0.750</td><td>0.779</td></tr><tr><td>LPC(18:0)<sup>b</sup></td><td>19.4 [15.7–20.5]</td><td>18.7 [17.1–21.4]</td><td>0.608</td><td>0.743</td></tr><tr><td>LPC(18:1)</td><td>5.02 [4.61–5.64]</td><td>3.97 [3.3–4.48]</td><td><strong>0.016</strong></td><td>0.186</td></tr><tr><td>LPC(18:2)<sup>b</sup></td><td>19.7 [18.7–21.6]</td><td>19.4 [17.1–20.5]</td><td>0.430</td><td>0.615</td></tr><tr><td>LPC(20:4)<sup>b</sup></td><td>8.97 [8.27–10.63]</td><td>8.4 [8.16–8.81]</td><td>0.350</td><td>0.615</td></tr><tr><td>LPC(22:6)<sup>b</sup></td><td>0.386 [0.331–0.393]</td><td>0.325 [0.274–0.362]</td><td>0.196</td><td>0.499</td></tr><tr><td>LPE(16:0)</td><td>0.55 [0.52–0.58]</td><td>0.70 [0.61–0.72]</td><td><strong>0.017</strong></td><td>0.186</td></tr><tr><td>LPE(18:0)<sup>b</sup></td><td>0.74 [0.63–0.79]</td><td>0.82 [0.78–0.88]</td><td><strong>0.046</strong></td><td>0.250</td></tr><tr><td>LPE(18:1)<sup>b</sup></td><td>0.25 [0.23–0.27]</td><td>0.27 [0.25–0.30]</td><td>0.425</td><td>0.615</td></tr><tr><td>LPE(18:2)<sup>b</sup></td><td>0.38 [0.37–0.42]</td><td>0.44 [0.37–0.48]</td><td>0.383</td><td>0.615</td></tr><tr><td>LPE(20:4)<sup>b</sup></td><td>0.158 [0.143–0.169]</td><td>0.168 [0.167–0.201]</td><td>0.089</td><td>0.390</td></tr><tr><td>LPE(22:6)<sup>b</sup></td><td>0.12 [0.11–0.14]</td><td>0.14 [0.12–0.15]</td><td>0.478</td><td>0.619</td></tr><tr><td>MAG(18:1)</td><td>0.07 [0.06–0.09]</td><td>0.11 [0.10–0.12]</td><td>0.126</td><td>0.462</td></tr><tr><td>MAG(18:2)<sup>b</sup></td><td>0.18 [0.15–0.21]</td><td>0.17 [0.11–0.26]</td><td>0.756</td><td>0.779</td></tr><tr><td>MAG(20:4)<sup>b</sup></td><td>0.0061 [0.0049–0.0084]</td><td>0.0068 [0.0045–0.0094]</td><td>0.779</td><td>0.779</td></tr></tbody></table>

## 2.3. *Pulmonary Hypertension (PH) Induced by Sugen/Hypoxia (SuHx)*

Administration of the VEGF-receptor antagonist sugen (SU5416) associated with chronic hypoxia was used to induce severe PH. The SuHx rat model is a well-recognized animal model of severe PH (Group 1 PH). In this model, mean arterial pulmonary pressure (mPAP) was markedly increased compared with the control group (48 ± 7 vs. 16 ± 1 mm Hg, *p* < 0.001, Figure 4A). Heart weight was also increased (1.7 ± 0.2 vs. 1.2 ± 0.1 g, *p* < 0.001, Figure 4B), which not concerned the left (0.99 ± 0.06 vs. 0.96 ± 0.10 g, *p* = 0.472, Figure 4C) but only the right (0.73 ± 0.14 vs. 0.25 ± 0.03 g, *p* < 0.001, Figure 4D) ventricle. In addition, cardiac output was decreased in the PH group compared with the control group (53 ± 9 vs. 98 ± 4 mL/min, *p* < 0.001, Figure 4E).

PH was associated with an increase in all monitored LPA that remains significant after BH correction except for LPA (18:0), and a trend tower an increase in LPE(16:0) with an adjusted *p*-value of 0.056 (Table 3 and Figure 5).

<footnote><sup>a</sup> *p*-values adjusted according to Benjamini and Hochberg multiple testing. <sup>b</sup> Analytes without analytical standard. HTN: hypertension, LPA: lysophosphatidic acid, LPC: lysophosphatidylcholine, LPE: lysophosphatidylethanolamine, MAG: monoacylglycerol. *p*-values computed from Student's *t*-test. *p*-values < 0.05 are highlighted in bold.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evaluation of pulmonary and cardiac parameters in control and PH rats. (A) Mean pulmonary arterial pressure (mPAP, mm Hg). N = 6 and 5 for Control and PH, respectively, (B) Heart weight (g). N = 7 and 5 for Control and PH, respectively, (C) Left ventricle weight (g). N = 7 and 5 for Control and PH, respectively, (D) Right ventricle weight (g). N = 7 and 5 for Control and PH, respectively, and (E) cardiac output (mL/min). N = 6 and 5 for Control and PH, respectively. PH: pulmonary hypertension. *p*-values computed by Student's *t*-test.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Impact of PH on lysophosphatidic acid (LPA) species (A) LPA(16:0), (B) LPA(18:0); (C) LPA(18:1), (D) LPA(18:2), (E) LPA(20:4), (F) LPA(22:6). *p*-values computed by Student's *t*-test. N = 10 and 7 for Control and PH, respectively.</figcaption></figure>
<table><caption>**Table 3.** Impact of PH induced by SuHx on LPA, LPL and MAG levels.</caption><thead><tr><th>Analyte</th><th>Control Group (N = 10)</th><th>PH Group (N = 7)</th><th>p-Value</th><th>Adjusted p-Value <sup>a</sup></th></tr></thead><tbody><tr><td>LPA(16:0)</td><td>0.091 [0.080-0.12]</td><td>0.21 [0.19-0.26]</td><td><strong>&lt;0.001</strong></td><td><strong>0.003</strong></td></tr><tr><td>LPA(18:0)<sup>b</sup></td><td>0.13 [0.11-0.17]</td><td>0.19 [0.17-0.25]</td><td>0.080</td><td>0.219</td></tr><tr><td>LPA(18:1)</td><td>0.062 [0.056-0.085]</td><td>0.156 [0.125-0.182]</td><td><strong>&lt;0.001</strong></td><td><strong>0.001</strong></td></tr><tr><td>LPA(18:2)<sup>b</sup></td><td>0.99 [0.80-1.29]</td><td>2.31 [2.03-4.00]</td><td><strong>&lt;0.001</strong></td><td><strong>0.001</strong></td></tr><tr><td>LPA(20:4)<sup>b</sup></td><td>1.22 [1.05-1.51]</td><td>2.18 [1.99-2.76]</td><td><strong>&lt;0.001</strong></td><td><strong>0.004</strong></td></tr><tr><td>LPA(22:6)<sup>b</sup></td><td>0.064 [0.057-0.082]</td><td>0.137 [0.119-0.160]</td><td><strong>&lt;0.001</strong></td><td><strong>&lt;0.001</strong></td></tr><tr><td>LPC(16:0)</td><td>16.8 [16.0-17.3]</td><td>16.7 [14.9-17.2]</td><td>0.357</td><td>0.604</td></tr><tr><td>LPC(16:1)<sup>b</sup></td><td>1.96 [1.625-2.07]</td><td>1.99 [1.585-2.86]</td><td>0.348</td><td>0.604</td></tr><tr><td>LPC(18:0)<sup>b</sup></td><td>32.8 [31.3-35.9]</td><td>33.7 [30.7-36.0]</td><td>0.749</td><td>0.827</td></tr><tr><td>LPC(18:1)</td><td>7.01 [6.66-7.32]</td><td>6.80 [6.11-7.31]</td><td>0.485</td><td>0.686</td></tr><tr><td>LPC(18:2)<sup>b</sup></td><td>28.7 [26.3-30.2]</td><td>28.1 [26.8-30.0]</td><td>0.669</td><td>0.827</td></tr><tr><td>LPC(20:4)<sup>b</sup></td><td>11.1 [10.3-11.8]</td><td>9.3 [7.6-10.4]</td><td>0.231</td><td>0.462</td></tr><tr><td>LPC(22:6)<sup>b</sup></td><td>0.36 [0.32-0.39]</td><td>0.36 [0.33-0.39]</td><td>0.693</td><td>0.827</td></tr><tr><td>LPE(16:0)</td><td>0.78 [0.73-0.79]</td><td>0.81 [0.80-1.03]</td><td><strong>0.015</strong></td><td>0.056</td></tr><tr><td>LPE(18:0)<sup>b</sup></td><td>1.14 [1.05-1.28]</td><td>1.39 [1.23-1.48]</td><td>0.053</td><td>0.168</td></tr><tr><td>LPE(18:1)<sup>b</sup></td><td>0.29 [0.27-0.34]</td><td>0.38 [0.335-0.42]</td><td>0.149</td><td>0.328</td></tr><tr><td>LPE(18:2)<sup>b</sup></td><td>0.53 [0.47-0.59]</td><td>0.62 [0.53-0.66]</td><td>0.457</td><td>0.686</td></tr><tr><td>LPE(20:4)<sup>b</sup></td><td>0.23 [0.21-0.25]</td><td>0.22 [0.18-0.24]</td><td>0.131</td><td>0.320</td></tr><tr><td>LPE(22:6)<sup>b</sup></td><td>0.158 [0.143-0.167]</td><td>0.176 [0.146-0.179]</td><td>0.752</td><td>0.827</td></tr><tr><td>MAG(18:1)</td><td>0.095 [0.066-0.156]</td><td>0.089 [0.058-0.113]</td><td>0.499</td><td>0.686</td></tr><tr><td>MAG(18:2)<sup>b</sup></td><td>0.70 [0.40-1.10]</td><td>0.60 [0.41-0.98]</td><td>0.860</td><td>0.860</td></tr><tr><td>MAG(20:4)<sup>b</sup></td><td>0.014 [0.011-0.018]</td><td>0.015 [0.006-0.023]</td><td>0.806</td><td>0.844</td></tr></tbody></table>

## 2.4. Impact of LPA on Pulmonary Artery Smooth Muscle Cell (PA-SMC) Proliferation

To better evaluate whether LPA species may be involved in the pathophysiology of PH, we investigated their impact on the proliferation of human pulmonary artery smooth muscle cells (PA-SMC) derived from PH and control patients using 5-bromo-2-deoxyuridine (BrdU) incorporation (Figure 6A,B, respectively). Of interest, the conducted ANOVA was significant (*p* = 0.006 and *p* < 0.001, respectively) and all tested LPA species (LPA(18:1), LPA(18:2) and LPA(20:4)) were found to increase PA-SMC proliferation, a phenomenon that was more pronounced in PH PA-SMC.

## 2.5. Relationship between LPA Levels and Index of PH in Patients with Cardiac Dysfunction

We aimed to determine if LPA levels are associated with estimated pulmonary artery systolic pressure (PAPs) determined by transthoracic echocardiography in 90 patients with cardiac dysfunction followed in the Department of Cardiology of Rouen University Hospital. The estimated PAPs was 36 ± 14 mm Hg and its increase was positively associated with LPA(16:0) levels (*p* = 0.013; Figure 7).

<footnote><sup>a</sup> *p*-values adjusted according to Benjamini and Hochberg multiple testing. <sup>b</sup> Analytes without analytical standard. HTN: hypertension, LPA: lysophosphatidic acid, LPC: lysophosphatidylcholine, LPE: lysophosphatidylethanolamine, MAG: monoacylglycerol. *p*-values computed by Student's *t*-test. *p*-values < 0.05 are highlighted in bold.</footnote>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Impact of LPA species on cultured pulmonary artery smooth muscle cell (PA-SMC) proliferation assessed by 5-bromo-2-deoxyuridine (BrdU) incorporation. (A) PA-SMC derived from control (CTRL) patients (N = 4 in triplicates). (B) PA-SMC derived from a patient with pulmonary hypertension (PH) (8 replicates of 1 patient). *p*-values computed by Dunnett's post-hoc test after ANOVA significance.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Association between estimated pulmonary artery systolic pressure (mm Hg) and relative abundance of plasma LPA(16:0). N = 90, *p*-value computed by Student's *t*-test.</figcaption></figure>

**3. Discussion**

**3.1. LPL, LPA and MAG Quantitation**

The present analytical method aimed to monitor 22 compounds. The choice has been made to investigate compounds with and without analytical standards to have the widest possible overview of LPA metabolism based on known fragmentation patterns for each class. This leads to the limit that results are expressed as compound-to-IS area under the curve ratio and not as absolute quantitation. Furthermore, since we did not have specific deuterated internal standard for each compound, we cannot be sure of the reliability of the correction of matrix effects. However, for each model, all samples were prepared and analyzed uninterruptedly resulting in the absence of inter-batch variability [20]. Further-
more, the use of these IS, even if not being the best-suited, allowed to decrease within-batch variability normalizing area under the curve of the monitored compounds. Moreover, the low dispersibility of the intra-group values allows us to be sufficiently confident about the robustness and reliability of the results obtained. Finally, coefficients of variation for all monitored compounds did not exceed 30% for both repeatability and reproducibility except for LPC(22:6), LPE(18:0), MAG(18:1) and MAG(20:4).

Precautions should be taken in the preparation and analysis of these compounds. First of all, since LPA are highly acidic compounds, it is strongly recommended to perform an acidic extraction to avoid poor recovery. This could be carried out using hydrochloric or formic acid. Here, we used formic acid to prevent production of LPA by LPL hydrolysis that could be induced by too acidic solutions. Furthermore, it is mandatory to chromatographically resolve LPA and LPL that possess the same esterified fatty acid since collision-induced dissociation may artificially produce LPA. The natural abundance of LPA being a tiny fraction of LPL, artificial production or detection of only a minor part of LPL to LPA, either during extraction or MS analysis, could dramatically increase the concentration of these compounds [21].

In our study, LPL and MAG were monitored in plasma samples mainly reflecting extracellular production by ATX and distribution through lipoproteins and albumin. Since LPA are also constitutive of oxLDL, quantitation of albumin and lipoproteins, especially oxLDL, as well as expression and/or activity of key enzymes involved in LPA metabolism could have been worth to better assess relationships within LPA metabolism

## 3.2. CVD Models

All evaluated models have been well characterized with an increase in SBP, a decrease in LVFS and an increase in mPAP for the HTN, CHF and PH model, respectively. As observed in humans, the three tested models were also associated with the development of cardiac hypertrophy and dysfunction. Of note, although the development of PH secondary to left heart disease is frequent, we did not observe an increase in mPAP three months after the surgery used to induce CHF in rats.

The development of HTN induced a significant increase of 3 LPA species under unadjusted assumption. This may be justified by the fact that Ang-II modulates LPA1 receptor function [22] and induces LPA production by phospholipases activation as well as the production of phosphatidic acids [23,24]. Of note, Ang-II also enhances cell proliferation and hypertrophy [25].

CHF rats exhibited an increase in LPA(16:0) and LPE(16:0-, 18:0-) associated with a decrease in LPC(18:1). All increased compounds possess an esterified saturated fatty acid (SFA) suggesting a link between CHF-mediated oxidative stress and desaturase activity [26,27]. Furthermore, it is now well known that reducing intake of saturated fatty acids may reduce LDL-c and CVD risks [28,29]. SFA possessed numerous deleterious effects such as inflammation, apoptosis, mitochondrial dysfunction and oxidative stress [30] but their impact has mainly been evaluated as non-esterified fatty acids [31]. Thus, the role of esterified SFA in LPA and LPL needs further investigations. Finally, we used an ischemic model of CHF induced by coronary ligation since myocardial infarction is the most frequent cause of CHF in humans. However, it must be noted that CHF may developed following hypertension, diabetes, kidney disease, valve diseases or congenital heart disease.

In the SuHx rat model, severe PH was induced and LPA metabolism and activity seem to be a promising pathway to target. Indeed, 80% of the monitored LPA exhibited a 2- to 3-fold increase when compared with the control group. Furthermore, exogenous addition of LPA species in a cell media of cultured human PA-SMC at a concentration of 1 µM, which is relevant when compared with circulating LPA levels [32], led to PA-SMC proliferation. In fact, accumulation of PA-SMCs within the pulmonary arterial walls is one of the most prominent features of PH, which can lead to the narrowing or occlusion of pulmonary vessels and therefore plays an important role in the occurrence and development of PH [33–36]. At this time, the role of LPA remains unclear even if their increase in
hypoxic pulmonary vascular remodeling has already been described [37]. Of note, LPA are able under physiological conditions to regulate the hypoxia-inducible factor 1α (HIF-1α) [38], that can interact with enzymes and other transcription factors in order to control vascularization and tissue growth in response to hypoxic conditions [39]. Furthermore, we observed a positive association between 16:0 LPA level, which is the most concentrated LPA specie in plasma, and the estimated PAPs in patients with cardiac dysfunction.

Although the results we obtained strengthen the hypothesis that LPA exert deleterious effects in PH, contradictory results have been reported in the literature with either protective [40,41] or harmful [42–44] effects of these lipid mediators against lung inflammation and fibrosis. Thus, blocking LPA synthesis by autotaxin or LPA actions using specific receptor antagonists in the model of SuHx represents a next step to finally demonstrate the deleterious role of LPA in PH development.

Of note, it may have been interesting to assess in parallel the evolution of LPL in a model of atherosclerosis since LPA species accumulate in high concentrations in atherosclerotic lesions contributing to thrombosis development [45,46].

### 3.3. LPA Signaling

Since LPA metabolism is complex with distinct anabolic and catabolic pathways, it may be easier to specifically act on LPA receptors rather than try to modulate enzymes involved in LPA metabolism regulation. However, it is now well admitted that LPA activate at least 6 specific G protein-coupled receptor named LPAR<sub>1-6</sub>. The downstream signals derived from activation of LPAR involve Rho, phospholipase C, phosphatidylinositol 3-kinase and adenylate cyclase intracellular pathways, thus producing diverse physio(patho)logical effects [18]. Recently, LPAR<sub>4</sub> was found to contribute to elevation of blood pressure with the more potent effect for LPA 18:1 and 20:4 species [47]. Interestingly, LPAR activation exhibited different reactivity according to the acyl chain (carbon length and number of insaturations) and the esterification position (sn-1 or sn-2) of the LPA [48]. Another key point is related to the relative expression of those receptors in lung tissue where PH takes place. All LPAR<sub>1-6</sub> were expressed in lung tissue according to https://www.genecards.org/ (accessed on 20 August 2021) and to date, it has only been demonstrated that LPAR<sub>1</sub> and LPAR<sub>2</sub> knock-down protects against pulmonary remodeling, lung inflammation and lung injury [40,49,50].

## 4. Materials and Methods

### 4.1. Experimentation

#### 4.1.1. Chemicals

Methanol (MeOH), Dichloromethan (DCM) and water of HPLC grade were purchased from Carlo Erba (Fontenayaux-Roses, France). Formic acid was purchased from VWR chemicals (Leuven, Belgium). LPA(16:0), LPA(17:1), LPC(17:1), LPE(17:1) LPA(18:1), LPC(16:0), LPC(18:1), LPE(16:0) LPS(16:0) and MAG(18:1) were purchased from Avanti Polar Lipids (Interchim, Montluçon, France). Ammonium acetate and formic acid were purchased from Carlo Erba (Fontenay-aux-Roses, France). Chromatographic Accucore XL C18 4 µm (150 × 3 mm) column was purchased from ThermoFisher Scientific (Illkirch-Graffenstaden, France)

#### 4.1.2. Animals

All the animal care and procedures were approved by French Animal Experimentation Ethics Committees and performed in accordance with the guidelines from the French National Research Council for the Care and Use of Laboratory Animals (Permit Numbers: Apafis #24107 and #11484). All experiments were performed in 10-weeks-old male wild-type Sprague-Dawley (SD/Crl) (Charles-River) rats.
## 4.2. LPL, LPA and MAG Quantitation

### 4.2.1. Blood Sampling and Processing

At the time of sacrifice, blood samples were drawn in the aorta of rats anesthetized with isoflurane using 2-mL syringes. Blood was transferred on a prechilled ethylenediaminetetraacetic acid (EDTA) tubes and immediately centrifuged 5 min at 4500× g (4 °C). Then, the plasma was frozen in liquid nitrogen and stored at -80 °C until analysis. For LPA, LPL and MAG quantitation, 300 µL of plasma were thawed at room temperature, spiked with internal standards (IS) and subjected to protein precipitation using 1 mL of MeOH. Then, liquid-liquid extraction was performed after the addition of 100 µL of HCOOH (10%), 2 mL of DCM and 350 µL of H₂O. The bottom organic phase was collected after being mixed and centrifuged (4500× g, 10 min, +4 °C), evaporated to dryness under a gentle stream of nitrogen and reconstituted in 100 µL of MeOH.

### 4.2.2. High Performance Liquid Chromatography (HPLC)

HPLC was carried out using a Prominence Shimadzu UFLC system consisting of a DGU-20A3 degasser, a LC-20AB pump, a SIL-20ACHT autosampler and a CTO20AC oven (Shimadzu, Prominence, Kyoto, Japan). Chromatographic separation was performed on an Accucore XL C18 column (4 µm particle size, 150 mm length × 3 mm inner diameter). The autosampler temperature was set at +8 °C, the column oven at 50 °C, the injected volume was 10 µL, and the flow rate was 500 µL/min. The method had the following gradient conditions using 1% CH₃COOH and 10 mM ammonium acetate in MeOH (solvent A) and 1% CH₃COOH and 10 mM ammonium acetate in water (solvent B): 0-0.5 min, 70% A, 0.5-5.5 min, 70-95% A, 5.5-10 min, 95% A, 10-10.5 min, 95-70% A, 10.5-12 min, 70% A.

### 4.2.3. Tandem Mass Spectrometry (MS²)

MS² was performed using a 4500 QTRAP operating in the positive/negative electrospray ionization (ESI) switching mode (Sciex, Toronto, Canada). Instrument control and data acquisition were performed with Analyst v1.6.3 software. The source parameters were optimized as follows; ion spray voltage: -4500 V and +4500 V for negative and positive ionization mode, respectively; nebulisation gas: 60 psi; desolvation gas: 50 psi; curtain gas: 30 psi; source temperature: 500 °C, entrance potential: -10 V and +10 V for negative and positive ionization mode, respectively, collision activation dissociation: medium. Relative quantitation was carried out using multiple reaction mode (MRM) to monitor transition from precursor to product ions. MS² parameters were optimized by direct infusion of LPA(16:0), LPA(17:1), LPA(18:1), LPC(16:0), LPC (17:1), LPC(18:1), LPE(16:0), LPE(17:1) and MAG(18:1) in pure MeOH at a concentration of 1 µg/mL. For monitored compounds without analytical standards (LPA(18:0-; 18:2-; 20:4-; 22:6-); LPC(16:1-; 18:0-; 18:2-; 20:4-; 22:6-); LPE(18:0-; 18:1-; 18:2-; 20:4-; 22:6-) and MAG(18:2-; 20:4-), transitions used were obtained from the LIPID MAPS tools "Predict MS/MS spectrum for a glycero(phospho)lipid" (https://www.lipidmaps.org/tools/) (accessed on 18 August 2021). According to typical ionization and fragmentation patterns of each LPL class, MS1 were [M-H]⁻ adduct and MS2 was set to 152.9 for LPA (Glycerol-3-phosphate ion with loss of H₂O), MS1 were [M-CH₃]⁻ adduct and MS2 was set to the carboxylate ion of the corresponding fatty acid (sn1 acyl chain [RCOO]⁻) for LPC, MS1 were [M-H]⁻ adduct and MS2 was set to 196.0 (neutral loss of fatty acid) for LPE and MS1 were [M+H]⁺ adduct and MS2 was set to the acylum form of the corresponding fatty acid (sn1 acyl chain ([RC=O]⁺)) for MAG (Table 4). Chromatograms showing the elution of the monitored compounds, calibration curves of the analytical standards used and repeatability and reproducibility analysis using pooled plasma performed in triplicate on two separate days are given in Supplementary File S1.
<table><caption>**Table 4.** MS parameters for studied analytes.</caption><thead><tr><th rowspan="2">Analyte</th><th colspan="2">Mass Transition</th><th colspan="3">MS Parameters</th><th rowspan="2">Adduct</th><th rowspan="2">Retention Time (min)</th></tr><tr><th>m/z (MS1)</th><th>m/z (MS2)</th><th>DP (V)</th><th>CE (eV)</th><th>CXP (V)</th></tr></thead><tbody><tr><td>LPA(16:0)</td><td>409.5</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.8</td></tr><tr><td>LPA(17:1)<sup>a</sup></td><td>421.1</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.5</td></tr><tr><td>LPA(18:0)<sup>b</sup></td><td>437.3</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>7.7</td></tr><tr><td>LPA(18:1)</td><td>435.2</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>7.0</td></tr><tr><td>LPA(18:2)<sup>b</sup></td><td>433.2</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.3</td></tr><tr><td>LPA(20:4)<sup>b</sup></td><td>457.2</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.3</td></tr><tr><td>LPA(22:6)<sup>b</sup></td><td>481.2</td><td>152.9</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.3</td></tr><tr><td>LPC(16:0)</td><td>480.6</td><td>255.4</td><td>-80</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>7.2</td></tr><tr><td>LPC(16:1)<sup>b</sup></td><td>478.6</td><td>253.4</td><td>-100</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>6.3</td></tr><tr><td>LPC(17:1)<sup>a</sup></td><td>492.3</td><td>267.3</td><td>-80</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>6.8</td></tr><tr><td>LPC(18:0)<sup>b</sup></td><td>508.4</td><td>283.5</td><td>-100</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>8.1</td></tr><tr><td>LPC(18:1)</td><td>506.4</td><td>281.5</td><td>-80</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>7.4</td></tr><tr><td>LPC(18:2)<sup>b</sup></td><td>504.4</td><td>279.5</td><td>-100</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>6.7</td></tr><tr><td>LPC(20:4)<sup>b</sup></td><td>528.4</td><td>303.4</td><td>-100</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>6.7</td></tr><tr><td>LPC(22:6)<sup>b</sup></td><td>552.4</td><td>327.4</td><td>-100</td><td>-40</td><td>-10</td><td>[M-CH<sub>3</sub>]<sup>-</sup></td><td>6.6</td></tr><tr><td>LPE(16:0)</td><td>452.4</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>7.2</td></tr><tr><td>LPE(17:1)</td><td>464.2</td><td>196.0</td><td>-80</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.9</td></tr><tr><td>LPE(18:0)<sup>b</sup></td><td>480.3</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>8.1</td></tr><tr><td>LPE(18:1)<sup>b</sup></td><td>478.3</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>7.4</td></tr><tr><td>LPE(18:2)<sup>b</sup></td><td>476.3</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.8</td></tr><tr><td>LPE(20:4)<sup>b</sup></td><td>500.6</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.7</td></tr><tr><td>LPE(22:6)<sup>b</sup></td><td>524.3</td><td>196.0</td><td>-50</td><td>-40</td><td>-10</td><td>[M-H]<sup>-</sup></td><td>6.7</td></tr><tr><td>MAG(18:1)</td><td>357.2</td><td>265.0</td><td>50</td><td>17</td><td>16</td><td>[M+H]<sup>+</sup></td><td>8.0</td></tr><tr><td>MAG(18:2)<sup>b</sup></td><td>355.2</td><td>263.0</td><td>50</td><td>17</td><td>16</td><td>[M+H]<sup>+</sup></td><td>7.5</td></tr><tr><td>MAG(20:4)<sup>b</sup></td><td>379.3</td><td>287.1</td><td>50</td><td>17</td><td>16</td><td>[M+H]<sup>+</sup></td><td>7.4</td></tr></tbody></table>

## 4.3. Murine Models of Cardiovascular Diseases (CVD)

### 4.3.1. Angiotensin-II Induced Hypertension (HTN)

Systemic arterial hypertension was induced using 4-week Ang-II infusion with osmotic pumps (0.25 µg/kg/day), implanted subcutaneously in isoflurane-anesthetized rats. Non-implanted rats served as controls. Tail-cuff plethysmography was used in trained conscious animals to confirm the development of systemic hypertension.

### 4.3.2. Chronic Heart Failure (CHF) Induced by Coronary Artery Ligation

For ischemic HF, myocardial infarction was induced by definitive left coronary artery ligation as previously described [51]. Briefly, rats were anesthetized using methohexthal (50 mg/kg i.p.), intubated and ventilated at 60 cycles/min (tidal volume 1 mL/100 g of body weight). A left thoracotomy was performed, and the heart exposed. A 6/0 polypropylene suture was passed around the proximal left coronary artery, which was tied in order to induce myocardial ischemia. In this case, 10 min after coronary artery ligation, the chest was closed, the pneumothorax was evacuated, and the animals were allowed to recover. Sham-operated rats, subjected to the same protocol except that the coronary artery was not occluded, served as controls. In this case, 12 weeks after surgery, the development of HF was confirmed using transthoracic echocardiography performed in rats anesthetized with isoflurane, using a Vivid 7 ultrasound echograph (GE Healthcare, Buc, France). A two-dimensional short axis view of the left ventricle was obtained at the level of the papillary muscle, in order to record M-mode tracings. Left ventricular end-diastolic (LVEDD) and systolic diameters (LVESD), allowing the determination of LV fractional shortening (FS) as $FS (\%) = ((LVEDD-LVESD)/LVEDD) \times 100$.

<footnote><sup>a</sup> Analytes used as internal standard for area under the curve normalization. <sup>b</sup> Analytes without analytical standard. CE: collision energy, CXP: collision cell exit potential, DP: declustering potential, LPA: lysophosphatidic acid, LPC: lysophosphatidylcholine, LPE: lysophosphatidylethanolamine, MAG: monoacylglycerol.</footnote>
#### 4.3.3. Pulmonary Hypertension (PH) Induced by Sugen Hypoxia (SuHx)

For PH, the SuHx model, in which occlusive neointimal lesions are observed as in humans, was used as previously described [33]. Briefly, rats were injected subcutaneously with SU5416 (a VEGF-receptor antagonist; 20 mg/kg) and exposed to hypoxia (10% $FiO_2$) for 3 weeks. Then, these rats returned to normoxia (21% $FiO_2$) for additional 5 weeks before evaluation. Control rats were not injected and remained under normoxia for 8 weeks.

To confirm the development of PH, right ventricular hemodynamic measurements were performed in anesthetized rats using a polyvinyl catheter introduced into the right external jugular vein, advanced in the right ventricle, and further, in the pulmonary artery, allowing the measurement of mean pulmonary arterial pressure.

#### 4.4. Human PA-SMC Proliferation

Lung specimens were obtained during lobectomy or pneumonectomy for localized lung cancer (control) and after lung transplant due to pulmonary hypertension (PH). For control group, lung specimens were collected at a distance from the tumor foci. This study was approved by the local ethics committee (ID CPP Est-III: 18.06.06; ID RCB: 2018-A01252-53) and all patients gave informed consent before the study. Human PA-SMC were isolated and cultured as previously described [33,34]. PA-SMC proliferation was measured by BrdU incorporation and by cell counting. BrdU staining was measured by the DELFIA® Cell proliferation kit (PerkinElmer, Courtaboeuf, France) and Time-resolved fluorometer EnVisionTM Multilabel Reader (PerkinElmer) 24 h of culture without and with LPA(18:1), LPA(18:2) and LPA(20:4) at 1 µM.

#### 4.5. Chronic Heart Failure Patients

In this case, 90 patients with cardiac dysfunction were included during a consultation visit in the Department of Cardiology of Rouen University Hospital. Blood sampling was performed allowing the quantitation of monitored compounds as described in Section 4.2. and estimated pulmonary artery systolic pressure, an indicator of PH, was obtained using transthoracic echocardiography. A value higher than 40 mm Hg is usually considered as abnormal and should prompt further investigation for PH diagnosis using right heart catheterization. The study was approved by the local ethics committee (ID RCB: 2014-A00834-43) and all patients gave written informed consent before the study.

#### 4.6. Data and Statistical Analysis

Data, statistical analysis and captions were performed using R v4.1.0 software [52] and DescTools v0.99.43 package [53]. For relative quantitation of LPL, data were expressed as median (interquartile range or IQR) and p-values were computed by *t*-test and adjusted by Benjamini and Hochberg for multiple test corrections in order to decrease the false discovery rate [19]. Captions were performed using ggplot2 v3.3.5, ggsci v2.9 and ggpubr v0.4.0 packages [54–56]. For cell culture data analysis, ANOVA was conducted first then *p*-values were computed using Dunnett's post-hoc test. For LPL, LPA and MAG quantitation, due to the lack of analytical standard for several compounds, 17:1- derivatives of each class (LPA, LPC and LPE) served as internal standard (IS) and data were expressed as compound-to-IS area under the curve (AUC) ratio that is unit-free. For MAG, 17:1-LPC was used as IS. R code for statistical analysis, rat models' data, cell culture data and human raw data are available in Supplementary Files S2–S5, respectively.

### 5. Conclusions

This study suggests that an increase in circulating LPA, observed in the relevant SuHx model and in patients with cardiac dysfunction, may contribute to the pathophysiology of PH. Additional experiments are warranted using modulators of LPA metabolism and signaling such as specific inhibitors, receptor antagonists or gene deletions in particular in the SuHx model to finally demonstrate that targeting LPA pathway may help to decrease occurrence and/or improve the outcome of PH.
**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/metabo11110784/s1, File S1: analytical method validation, File S2: R code for statistical analysis, File S3: murine models' raw data, File S4: cell culture raw data, File S5: human raw data.

**Author Contributions:** Conceptualization, C.G. and J.B.; data curation, L.T., C.G. and J.B.; formal analysis, T.D., L.T., G.F., F.B., C.G. and J.B.; funding acquisition, J.B.; investigation, T.D., L.T., M.L., H.M., D.G., G.F., S.A., R.T., F.B., M.H. and J.B.; methodology, T.D., L.T., C.G. and J.B.; project administration, V.R., C.G. and J.B.; resources, F.B., V.R., C.G. and J.B.; supervision, T.D., V.R., C.G. and J.B.; validation, T.D., C.G. and J.B.; writing—original draft, T.D., L.T., C.G. and J.B.; writing—review and editing, T.D., L.T., M.L., H.M., D.G., G.F., S.A., R.T., F.B., M.H., V.R., C.G. and J.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by grants from the French National Research Agency (ANR-16-CE17-0012).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the French Animal Experimentation Ethics Committees in accordance with the guidelines from the French National Research Council for the Care and Use of Laboratory Animals (Permit Numbers: Apafis #24107 and #11920).

**Informed Consent Statement:** For PA-SMC experiments, lung specimens' collection was approved by the local ethics committee (ID CPP Est-III: 18.06.06; ID RCB: 2018-A01252-53) and all patients gave informed consent before the study. For human blood samples from CHF patients, the collection was approved by the local ethics committee (ID CPP Nord-Ouest I: 02/014/2014; ID RCB: 2014-A00834-43) and all patients gave written informed consent before the study.

**Data Availability Statement:** The data presented in this study are available in Supplementary Materials.

**Acknowledgments:** Authors are grateful to Fabien Lamoureux, in head of the laboratory of Pharmacokinetics of the Rouen University Hospital; to M. Tony Pereira, engineer of the mass spectrometry platform of the Rouen University Hospital, for his precious help concerning the chromatographic part of the analysis; to all the staff of the laboratory of Pharmacokinetics of the Rouen University Hospital for their technical support. Authors are also grateful to Valéry Brunel and to Frederic Ziegler for the quantitation of aldosterone levels.

**Conflicts of Interest:** Over the last three years, C.G. reports grants from Acceleron and Janssen and grants and personal fees from Merck, outside the submitted work. M.H. reports grants and personal fees from Acceleron, Bayer, Merck, and Janssen and personal fees from GSK, outside the submitted work. The other authors declare no conflict of interest.

## References

1. Roth, G.A.; Mensah, G.A.; Johnson, C.O.; Addolorato, G.; Ammirati, E.; Baddour, L.M.; Barengo, N.C.; Beaton, A.Z.; Benjamin, E.J.; Benziger, C.P.; et al. Global Burden of Cardiovascular Diseases and Risk Factors, 1990–2019: Update From the GBD 2019 Study. *J. Am. Coll. Cardiol.* **2020**, *76*, 2982–3021. [CrossRef] [PubMed]

2. Soppert, J.; Lehrke, M.; Marx, N.; Jankowski, J.; Noels, H. Lipoproteins and lipids in cardiovascular disease: From mechanistic insights to therapeutic targeting. *Adv. Drug Deliv. Rev.* **2020**, *159*, 4–33. [CrossRef]

3. Eilat-Tsanani, S.; Mor, E.; Schonmann, Y. Statin Use Over 65 Years of Age and All-Cause Mortality: A 10-Year Follow-Up of 19 518 People. *J. Am. Geriatr. Soc.* **2019**, *67*, 2038–2044. [CrossRef]

4. McKenney, J.M. Optimizing LDL-C Lowering With Statins. *Am. J. Ther.* **2004**, *11*, 54–59. [CrossRef] [PubMed]

5. Chaudhary, R.; Garg, J.; Shah, N.; Sumner, A. PCSK9 inhibitors: A new era of lipid lowering therapy. *World J. Cardiol.* **2017**, *9*, 76–91. [CrossRef] [PubMed]

6. Yeang, C.; Wilkinson, M.J.; Tsimikas, S. Lipoprotein(a) and oxidized phospholipids in calcific aortic valve stenosis. *Curr. Opin. Cardiol.* **2016**, *31*, 440–450. [CrossRef]

7. Kotani, K.; Yamada, S.; Yamada, T.; Taniguchi, N.; Sakurabayashi, I. The relationship between oxidized lipoprotein(a) and carotid atherosclerosis in asymptomatic subjects: A comparison with native lipoprotein(a). *Lipids Health Dis.* **2011**, *10*, 174. [CrossRef]

8. Parthasarathy, S.; Raghavamenon, A.; Garelnabi, M.O.; Santanam, N. Oxidized Low-Density Lipoprotein. *Methods Mol Biol.* **2009**, *610*, 403–417. [CrossRef]

9. Oei, H.-H.S.; van der Meer, I.M.; Hofman, A.; Koudstaal, P.J.; Stijnen, T.; Breteler, M.M.; Witteman, J.C. Lipoprotein-Associated Phospholipase A2 Activity Is Associated with Risk of Coronary Heart Disease and Ischemic Stroke: The Rotterdam Study. *Circulation* **2005**, *111*, 570–575. [CrossRef]

10. Colley, K.J.; Wolfert, R.L.; Cobble, M.E. Lipoprotein associated phospholipase A2: Role in atherosclerosis and utility as a biomarker for cardiovascular risk. *EPMA J.* **2011**, *2*, 27–38. [CrossRef]
11. Gianazza, E.; Brioschi, M.; Fernandez, A.E.M.; Banfi, C. Lipoxidation in cardiovascular diseases. *Redox Biol.* **2019**, *23*, 101119.
[CrossRef] [PubMed]

12. Zhou, Y.; Little, P.J.; Ta, H.T.; Xu, S.; Kamato, D. Lysophosphatidic acid and its receptors: Pharmacology and therapeutic potential in atherosclerosis and vascular disease. *Pharmacol. Ther.* **2019**, *204*, 107404. [CrossRef] [PubMed]

13. Abdel-Latif, A.; Heron, P.M.; Morris, A.J.; Smyth, S.S. Lysophospholipids in coronary artery and chronic ischemic heart disease. *Curr. Opin. Lipidol.* **2015**, *26*, 432–437. [CrossRef] [PubMed]

14. Kleger, A.; Liebau, S.; Lin, Q.; von Wichert, G.; Seufferlein, T. The Impact of Bioactive Lipids on Cardiovascular Development. *Stem Cells Int.* **2011**, *2011*, 916180. [CrossRef]

15. Lee, H.; Goetzl, E.J.; An, S. Lysophosphatidic acid and sphingosine 1-phosphate stimulate endothelial cell wound healing. *Am. J. Physiol. Cell Physiol.* **2000**, *278*, C612–C618. [CrossRef] [PubMed]

16. Moolenaar, W.H. Lysophosphatidic Acid, a Multifunctional Phospholipid Messenger. *J. Biol. Chem.* **1995**, *270*, 12949–12952.
[CrossRef]

17. Smyth, S.S.; Cheng, H.-Y.; Miriyala, S.; Panchatcharam, M.; Morris, A.J. Roles of lysophosphatidic acid in cardiovascular physiology and disease. *Biochim. Biophys. Acta* **2008**, *1781*, 563–570. [CrossRef]

18. Yang, F.; Chen, G.-X. Production of extracellular lysophosphatidic acid in the regulation of adipocyte functions and liver fibrosis. *World J. Gastroenterol.* **2018**, *24*, 4132–4151. [CrossRef]

19. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. *J. R. Stat. Soc. Ser. B* **1995**, *57*, 289–300. [CrossRef]

20. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted Metabolomics. *Curr. Protoc. Mol. Biol.* **2012**, *98*, 30.2.1–30.2.24.
[CrossRef]

21. Triebl, A.; Trötzmüller, M.; Eberl, A.; Hanel, P.; Hartler, J.; Köfeler, H. Quantitation of phosphatidic acid and lysophosphatidic acid molecular species using hydrophilic interaction liquid chromatography coupled to electrospray ionization high resolution mass spectrometry. *J. Chromatogr. A* **2014**, *1347*, 104–110. [CrossRef]

22. Colín-Santana, C.C.; Avendaño-Vázquez, S.E.; Alcántara-Hernández, R.; García-Sáinz, J.A. EGF and angiotensin II modulate lysophosphatidic acid LPA1 receptor function and phosphorylation state. *Biochim. Biophys. Acta* **2011**, *1810*, 1170–1177. [CrossRef] [PubMed]

23. De Souza, A.M.; de Carvalho, T.L.G.; Lara, L.D.S.; Gomes-Quintana, E.; Gil Lopes, A.; Caruso-Neves, C. The stimulatory effect of angiotensin II on Na⁺-ATPase activity involves sequential activation of phospholipases and sustained PKC activity. *Biochim. Biophys. Acta* **2010**, *1798*, 354–359. [CrossRef] [PubMed]

24. Booz, G.W.; Taher, M.M.; Baker, K.M.; Singer, H.A. Angiotensin II induces phosphatidic acid formation in neonatal rat cardiac fibroblasts: Evaluation of the roles of phospholipases C and D. *Mol. Cell. Biochem.* **1994**, *141*, 135–143. [CrossRef]

25. Wolf, G.; Wenzel, U.O. Angiotensin II and Cell Cycle Regulation. *Hypertension* **2004**, *43*, 693–698. [CrossRef] [PubMed]

26. Tsutsui, H.; Kinugawa, S.; Matsushima, S. Oxidative stress and heart failure. *Am. J. Physiol. Heart Circ. Physiol.* **2011**, *301*, H2181–H2190. [CrossRef] [PubMed]

27. Idkowiak-Baldys, J.; Apraiz, A.; Li, L.; Rahmaniyan, M.; Clarke, C.J.; Kraveka, J.M.; Asumendi, A.; Hannun, Y.A. Dihydroceramide desaturase activity is modulated by oxidative stress. *Biochem. J.* **2010**, *427*, 265–274. [CrossRef] [PubMed]

28. Briggs, M.A.; Petersen, K.S.; Kris-Etherton, P.M. Saturated Fatty Acids and Cardiovascular Disease: Replacements for Saturated Fat to Reduce Cardiovascular Risk. *Healthcare* **2017**, *5*, 29. [CrossRef]

29. Mensink, R.P.; World Health Organization. *Effects of Saturated Fatty Acids on Serum Lipids and Lipoproteins: A Systematic Review and Regression Analysis*; World Health Organization: Geneva, Switzerland, 2016; Available online: https://apps.who.int/iris/handle/10665/246104 (accessed on 25 August 2021).

30. Beam, J.; Botta, A.; Barendregt, R.; Ghosh, S. Dietary Fatty Acids, Redox Signaling, and the Heart. *Syst. Biol. Free Radic. Antioxid.* **2014**, *1497–1522*. [CrossRef]

31. Gouaref, I.; Bouazza, A.; Abderrhmane, S.A.; Koceir, E.-A. Lipid Profile Modulates Cardiometabolic Risk Biomarkers Including Hypertension in People with Type-2 Diabetes: A Focus on Unbalanced Ratio of Plasma Polyunsaturated/Saturated Fatty Acids. *Molecules* **2020**, *25*, 4315. [CrossRef]

32. Michalczyk, A.; Budkowska, M.; Dołęgowska, B.; Chlubek, D.; Safranow, K. Lysophosphatidic acid plasma concentrations in healthy subjects: Circadian rhythm and associations with demographic, anthropometric and biochemical parameters. *Lipids Health Dis.* **2017**, *16*, 1–9. [CrossRef]

33. Eddahibi, S.; Humbert, M.; Fadel, E.; Raffestin, B.; Darmon, M.; Capron, F.; Simonneau, G.; Dartellelle, P.; Hamon, M.; Adnot, S. Serotonin transporter overexpression is responsible for pulmonary artery smooth muscle hyperplasia in primary pulmonary hypertension. *J. Clin. Investig.* **2001**, *108*, 1141–1150. [CrossRef]

34. Bordenave, J.; Thuillet, R.; Tu, L.; Phan, C.; Cumont, A.; Marsol, C.; Huertas, A.; Savale, L.; Hibert, M.; Galzi, J.-L.; et al. Neutralization of CXCL12 attenuates established pulmonary hypertension in rats. *Cardiovasc. Res.* **2019**, *116*, 686–697. [CrossRef] [PubMed]

35. Morrell, N.W.; Adnot, S.; Archer, S.L.; Dupuis, J.; Lloyd Jones, P.; MacLean, M.R.; McMurtry, I.F.; Stenmark, K.R.; Thistlethwaite, P.A.; Weissmann, N.; et al. Cellular and Molecular Basis of Pulmonary Arterial Hypertension. *J. Am. Coll. Cardiol.* **2009**, *54* (Suppl. 1), S20–S31. [CrossRef] [PubMed]
36. Xu, Y.; Bei, Y.; Shen, S.; Zhang, J.; Lu, Y.; Xiao, J.; Li, X. MicroRNA-222 Promotes the Proliferation of Pulmonary Arterial Smooth Muscle Cells by Targeting P27 and TIMP3. *Cell. Physiol. Biochem.* **2017**, *43*, 282–292. [CrossRef]

37. Shlyonsky, V.; Naeije, R.; Mies, F. Possible Role of Lysophosphatidic Acid in Rat Model of Hypoxic Pulmonary Vascular Remodeling. *Pulm. Circ.* **2014**, *4*, 471–481. [CrossRef]

38. Lee, S.-J.; No, Y.R.; Dang, D.T.; Dang, L.H.; Yang, V.W.; Shim, H.; Yun, C.C. Regulation of Hypoxia-inducible Factor 1α (HIF-1α) by Lysophosphatidic Acid Is Dependent on Interplay between p53 and Krüppel-like Factor 5. *J. Biol. Chem.* **2013**, *288*, 25244–25253. [CrossRef] [PubMed]

39. Ziello, J.E.; Jovin, I.S.; Huang, Y. Hypoxia-Inducible Factor (HIF)-1 regulatory pathway and its potential for therapeutic intervention in malignancy and ischemia. *Yale J. Biol. Med.* **2007**, *80*, 51–60.

40. Cheng, H.-Y.; Dong, A.; Panchatcharam, M.; Mueller, P.; Yang, F.; Li, Z.; Mills, G.; Chun, J.; Morris, A.J.; Smyth, S.S. Lysophosphatidic Acid Signaling Protects Pulmonary Vasculature From Hypoxia-Induced Remodeling. *Arter. Thromb. Vasc. Biol.* **2012**, *32*, 24–32. [CrossRef]

41. He, D.; Su, Y.; Usatyuk, P.V.; Spannhake, E.W.; Kogut, P.; Solway, J.; Natarajan, V.; Zhao, Y. Lysophosphatidic Acid Enhances Pulmonary Epithelial Barrier Integrity and Protects Endotoxin-induced Epithelial Barrier Disruption and Lung Injury. *J. Biol. Chem.* **2009**, *284*, 24123–24132. [CrossRef] [PubMed]

42. Li, C.; Liu, P.-P.; Tang, D.-D.; Song, R.; Zhang, Y.-Q.; Lei, S.; Wu, S.-J. Targeting the RhoA-ROCK pathway to regulate T-cell homeostasis in hypoxia-induced pulmonary arterial hypertension. *Pulm. Pharmacol. Ther.* **2018**, *50*, 111–122. [CrossRef] [PubMed]

43. Tager, A.M.; LaCamera, P.; Shea, B.S.; Campanella, G.S.; Selman, M.; Zhao, Z.; Polosukhin, V.; Wain, J.; Karimi-Shah, B.A.; Kim, N.D.; et al. The lysophosphatidic acid receptor LPA1 links pulmonary fibrosis to lung injury by mediating fibroblast recruitment and vascular leak. *Nat. Med.* **2007**, *14*, 45–54. [CrossRef] [PubMed]

44. Tang, N.; Zhao, Y.; Feng, R.; Liu, Y.; Wang, S.; Wei, W.; Ding, Q.; An, M.; Wen, J.; Li, L. Lysophosphatidic acid accelerates lung fibrosis by inducing differentiation of mesenchymal stem cells into myofibroblasts. *J. Cell. Mol. Med.* **2013**, *18*, 156–169. [CrossRef]

45. Kritikou, E.; Van Puijvelde, G.H.M.; Van Der Heijden, T.; Van Santbrink, P.J.; Swart, M.; Schaftenaar, F.H.; Kröner, M.J.; Kuiper, J.; Bot, I. Inhibition of lysophosphatidic acid receptors 1 and 3 attenuates atherosclerosis development in LDL-receptor deficient mice. *Sci. Rep.* **2016**, *6*, 37585. [CrossRef] [PubMed]

46. Cui, M.-Z. Lysophosphatidic acid effects on atherosclerosis and thrombosis. *Clin. Lipidol.* **2011**, *6*, 413–426. [CrossRef] [PubMed]

47. Kano, K.; Matsumoto, H.; Inoue, A.; Yukiura, H.; Kanai, M.; Chun, J.; Ishii, S.; Shimizu, T.; Aoki, J. Molecular mechanism of lysophosphatidic acid-induced hypertensive response. *Sci. Rep.* **2019**, *9*, 1–12. [CrossRef]

48. Bandoh, K.; Aoki, J.; Taira, A.; Tsujimoto, M.; Arai, H.; Inoue, K. Lysophosphatidic acid (LPA) receptors of the EDG family are differentially activated by LPA species. Structure-activity relationship of cloned LPA receptors. *FEBS Lett.* **2000**, *478*, 159–165. [CrossRef]

49. Zhao, J.; He, D.; Su, Y.; Berdyshev, E.; Chun, J.; Natarajan, V.; Zhao, Y. Lysophosphatidic acid receptor 1 modulates lipopolysaccharide-induced inflammation in alveolar epithelial cells and murine lungs. *Am. J. Physiol. Lung Cell Mol. Physiol.* **2011**, *301*, L547–L556. [CrossRef]

50. Huang, L.S.; Fu, P.; Patel, P.; Harijith, A.; Sun, T.; Zhao, Y.; Garcia, J.G.N.; Chun, J.; Natarajan, V. Lysophosphatidic Acid Receptor-2 Deficiency Confers Protection against Bleomycin-Induced Lung Injury and Fibrosis in Mice. *Am. J. Respir. Cell Mol. Biol.* **2013**, *49*, 912–922. [CrossRef] [PubMed]

51. Merabet, N.; Bellien, J.; Glevarec, E.; Nicol, L.; Lucas, D.; Remy-Jouet, I.; Bounoure, F.; Dreano, Y.; Wecker, D.; Thuillez, C.; et al. Soluble epoxide hydrolase inhibition improves myocardial perfusion and function in experimental heart failure. *J. Mol. Cell. Cardiol.* **2012**, *52*, 660–666. [CrossRef]

52. R Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 25 August 2021).

53. Signorell, A.; Aho, A.; Alfons, A.; Anderegg, N.; Aragon, T.; Arachchige, C.; Arppe, A.; Baddeley, A.; Barton, K.; Bolker, B.; et al. DescTools: Tools for Descriptive Statistics. R Package Version 0.99.42. 2021. Available online: https://cran.r-project.org/package=DescTools (accessed on 25 August 2021).

54. Wickham, H. *Ggplot2: Elegant Graphics for Data Analysis*; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. Available online: https://ggplot2.tidyverse.org (accessed on 25 August 2021).

55. Xiao, N. *Ggsci: Scientific Journal and Sci-Fi Themed Color Palettes for 'ggplot2'*; R Package Version 2.9. 2018. Available online: https://CRAN.R-project.org/package=ggsci (accessed on 25 August 2021).

56. Kassambara, A. *Ggpubr: 'Ggplot2' Based Publication Ready Plots*; R Package Version 0.3.0. 2020. Available online: https://CRAN.R-project.org/package=ggpubr (accessed on 25 August 2021).