Right Ventricle Remodeling Metabolic Signature in Experimental Pulmonary Hypertension Models of Chronic Hypoxia and Monocrotaline Exposure

Introduction: Over time and despite optimal medical management of patients with pulmonary hypertension (PH), the right ventricle (RV) function deteriorates from an adaptive to maladaptive phenotype, leading to RV failure (RVF). Although RV function is well recognized as a prognostic factor of PH, no predictive factor of RVF episodes has been elucidated so far. We hypothesized that determining RV metabolic alterations could help to understand the mechanism link to the deterioration of RV function as well as help to identify new biomarkers of RV failure. Methods: In the current study, we aimed to characterize the metabolic reprogramming associated with the RV remodeling phenotype during experimental PH induced by chronic-hypoxia-(CH) exposure or monocrotaline-(MCT) exposure in rats. Three weeks after PH initiation, we hemodynamically characterized PH (echocardiography and RV catheterization), and then we used an untargeted metabolomics approach based on liquid chromatography coupled to high-resolution mass spectrometry to analyze RV and LV tissues in addition to plasma samples from MCT-PH and CH-PH rat models. Results: CH exposure induced adaptive RV phenotype as opposed to MCT exposure which induced maladaptive RV phenotype. We found that predominant alterations of arginine, pyrimidine, purine, and tryptophan metabolic pathways were detected on the heart (LV+RV) and plasma samples regardless of the PH model. Acetylspermidine, putrescine, guanidinoacetate RV biopsy levels, and cytosine, deoxycytidine, deoxyuridine, and plasmatic thymidine levels were correlated to RV function in the CH-PH model. It was less likely correlated in the MCT model. These pathways are well described to regulate cell proliferation, cell hypertrophy, and cardioprotection. These findings open novel research perspectives to find biomarkers for early detection of RV failure in PH.


Introduction
Regardless of the type of pulmonary hypertension (PH), the right ventricle (RV) initially adapts to pressure overload and is able to maintain a normal cardiac output [1][2][3]. Over time and despite optimal medical management, RV function deteriorates from adaptive to maladaptive phenotype, leading to right ventricle failure (RVF). Although RV function being a well-recognized prognostic factor of PH [4,5], no predictive factor of RVF has been elucidated so far [6]. However, several studies demonstrated that the adaptive to maladaptive phenotype shift involved defective glycolysis, RV capillary rarefaction, RV fibrosis, and RV inflammation [7][8][9]. We hypothesized that a global RV metabolic alterations study might help find new biomarkers for the early identification of PH-associated RV failure. The metabolome is the set of metabolites (small molecules, <1000 Da) found in biological media as opposed to metabolomics which deals with the large-scale detection of metabolites in biological samples. High-resolution mass spectrometry (HRMS) is preferentially used for metabolomics studies. It seems adapted to biomarkers screening as it has proved more efficient than low-resolution instruments in detecting and identifying metabolites [10]. Publications on metabolic alterations of PH patients' hearts remain scarce. In 2016, Rhodes et al. [11] investigated metabolic alterations in PH patients' plasma with mass spectrometry showing alterations of energy metabolism, tryptophan, and polyamine metabolic pathways. However, none of these alterations have been correlated to RV function or remodeling of the phenotype. Similarly, other studies limited their investigation to the plasmatic signature of PH patients considering the difficulty of collecting myocardial biopsies [12,13]. Of interest, studies on animal models focus on metabolic alterations in plasma without myocardial correlation analyses. Zheng et al. [14] analyzed plasma samples from monocrotaline (MCT)-PH rats and showed altered arginine metabolism. Zhao et al. analyzed plasma samples of both MCT-PH and Chronic-hypoxia (CH)-PH rat and showed that methionine metabolism and urea cycle metabolism were the most significant pathway involved in the pathogenesis of PH induced by MCT or CH-exposure [15].
The current study aimed at characterizing metabolic reprogramming associated with RV remodeling the phenotype during PH in order to give new insights into physiopathological processes involved in maladaptive RV remodeling and to highlight putative biomarkers for the early detection of RV failure. For that purpose, we analyzed RV and LV tissues in addition to plasma samples from MCT-PH and CH-PH rat models. We used an untargeted metabolomic approach based on liquid chromatography coupled to HRMS [16]. Plasma and myocardial metabolic signatures were compared with the goal of new RV-specific biomarker discovery in PH.

Animals
The animal facility is licensed by the French Ministry of Agriculture (agreement N • C92-019-01). This study was approved by the Committee on the Ethics of Animal Experiments (CEEA26 CAP Sud). Animal experiments were performed according to the guidelines from Directive 2010/63/EU on 22 September 2010 of the European Parliament on the protection of animals used for scientific purposes and complied with the French institution's guidelines for animal care and handling. Male Wistar rats (4 weeks old) were used in four experimental protocols: (1) PH was induced by a single MCT injection (60 mg/kg, subcutaneous) in 14 rats. MCT was dissolved in 1 N HCl and neutralized with 1 N NaOH. (2) 14 Control rats received the same volume of saline solution (0.9% NaCl) (designated here as "vehicle rats"). (3) For CH-exposure, 10 Wistar rats were placed for 3 weeks in hypoxia (10% O 2 ) and normobaric conditions. (4) Ten Normoxia Wistar rats were placed in the same room for 3 weeks as controls.

Echocardiographic Measurement
Evaluation by trans-thoracic echocardiography (TTE) was performed with a digital ultrasound system (Vivid E9, GE Healthcare, Buc, France) by using a high-frequency phased array transducer (12 S-D 4-12MHz, GE Healthcare, Buc, France). The echocardiographic evaluation procedure was performed under general anesthesia and spontaneous breathing with an Isoflurane Rodent Anesthesia System (Minerve, Esternay, France) (maintenance isoflurane 2% at room air). Rats were shaved, and body temperature was controlled during experiments. All analyses were performed in a blinded manner: rats' experimental conditions were unknown by the operator during TTE examination and data interpretation. In the parasternal short-axis view, we measured pulmonary artery acceleration time (PAAT; ms), pulmonary artery velocity time integral (VTI-PA), RV ejection time (RVET; ms), cycle length (ms) as previously described [17]. In the 4-cavity view, performed subcostally, we measured ascending aorta velocity time integral (VTI-Aorta (VTI-Ao)) reflecting LV cardiac output [18,19], VTI-Ao value reproducible in heart failure [20,21], Tricuspid annular plane systolic excursion (TAPSE; mm), RV and LV thickness (mm). RV and LV fractional shortening (RV FS and LV FS%) corresponded to the percent change in LV and RV cavity diameters, respectively. LV FS = ((LV EDd-LV EDs/LV EDd)*100 or RV FS = (RV EDd-RV EDs/RV EDd)*100, as previously described [17]. Echocardiographic measurements were illustrated in Supplementary Figure S1.

Right Heart Catheterization and Tissues Harvest
Rats were placed under general anesthesia and spontaneous breathing with an isoflurane Rodent Anesthesia System (Minerve Esternay, France) (maintenance: isoflurane 2% at room air). Hemodynamic measurements were performed, such as cardiac output (CO; mL/min), heart rate (HR; beat/min), and Min and Max dPdt (MindPdt and MaxdPdt, mmHg/s). Recording of the right ventricular pressure (RVP) allowed the measurement of the right ventricular systolic pressure (RVS.P; mmHg). Pulmonary vascular resistance (PVR) was assessed by the ratio of RVS.P and cardiac output. Hemodynamic parameters were blindly analyzed in unventilated anesthetized rats using LabChart software. LabChart allowed us to calculate several ventricular wave form parameters, including End diastolic pressure (EDP), Contractility index which corresponds to Max dP/dt divided by the RV pressure (RVP), at the time of Max dP/dt. Four rats from the MCT sample sets died before catheterization. Therefore, their RVS.P, HR, Max and MindPdt, CO, and PVR could not be assessed. After the catheterization, plasma, RV and LV were collected and stocked at −80 • C before metabolomic analyses. Fulton's index was calculated before sample congelation: the ratio of RV weight to LV plus septal (S) weight (RV/LV+S).

Pulmonary Vessel Remodeling Analysis
The lungs were fixed in 4% paraformaldehyde, embedded in paraffin, and serially sectioned (5 µm). The histopathological evaluation of the lungs was performed by a pulmonary pathologist who was blinded to the genotypes and treatment group assignments of the rats. The morphometric analyses were performed on sections stained with hematein-eosin-safran (HES). The lung samples from all the conditions were analyzed by conventional light microscopy using quantitative semi-automated software (NIS-BR; Nikon, France). Vascular remodeling was assessed in all the pulmonary vessels larger than 50 µm and less than 100 µm that were identified in 20 randomly selected microscopic fields. The wall thickness was calculated according to the following equation: ((External diameter − Internal diameter)/(External diameter)) × 100.

RV Fibrosis Assessment
RV were fixed in 4% paraformaldehyde, paraffin-embedded, and serially sectioned (5 µm). The sections were stained with trichrome Masson (#MST-100T Biognost) to evaluate fibrosis as a percentage of the total tissue area % of fibrosis was measured using Image J software (Bethesda, Rockville, MD, USA).

Sample Preparation for LC-HRMS Analysis
Please find more details about methods in Supplementary Material-From plasma: Metabolites were extracted from 50 µL of plasma as previously described [16]. A quality control (QC) sample was obtained by pooling 20 µL of each sample. Briefly, 50 µL of plasma were treated with 200 µL of cold methanol containing internal standards and left on ice for 90 min to allow protein precipitation. Supernatants were then collected and dried. Prior to LC-HRMS analysis, dried extracts were dissolved in 150 µL of H2O/Acetonitrile (95:5, v/v), containing 0.1% formic acid + external standards or 10 mM ammonium carbonate pH 10.5 + external standards/Acetonitrile (40:60, v/v) for metabolite analysis using C18 and ZIC-pHILIC columns, respectively. QC samples were injected every 10 samples in order to evaluate the signal variations of any metabolite. From heart tissues: Metabolites were extracted from 30 mg of ground ventricle resuspended in 170 µL of ultrapure water by sonication 5 times for 10 s. At this step, the total protein concentration was determined. A volume of 350 µL of methanol containing internal standards was added to the tissue lysate, which was then left on ice for 90 min for protein precipitation. Cell debris were then removed by centrifugation, and supernatants were split into two equal aliquots for C18 and HILIC analyses. The resulting aliquots were then dried and stored at −80 • C until analysis. Prior to LC-MS analysis, dried extracts were resuspended to reach a fixed protein concentration (equivalent to 20 mg/mL) using variable volumes of H2O/Acetonitrile (95:5, v/v), containing 0.1% formic acid + external standards or 10 mM ammonium carbonate pH 10.5 + external standards/Acetonitrile (40:60, v/v) for metabolite analysis using C18 and ZIC-pHILIC columns, respectively. A quality control (QC) sample was obtained by pooling 20 µL of each sample preparation.

LC-HRMS Conditions
Ten microliters of extracts were analyzed by LC-MS using an Ultimate 3000 chromatographic system (ThermoFisher Scientific, Courtaboeuf, 91140 France) coupled to an Exactive mass spectrometer (ThermoFisher Scientific Courtaboeuf, 91140 France) fitted with an electrospray source and operating in the positive and negative ion modes for metabolite separations on C18 (C18(+) conditions) and ZIC-pHILIC (HILIC(-) conditions) columns, respectively. Ions were detected from m/z 85 to 1000 [16]. Ultra-high-performance LC (UHPLC) separation was performed on a Hypersil GOLD C18 (1.9 µm, 2.1 mm × 150 mm) column at 30 • C (ThermoFisher Scientific, les Ulis, France), using a gradient program of 0.1% formic acid in water (phase A) and 0.1% formic acid in acetonitrile (phase B). High-performance LC (HPLC) separation was performed on a Sequant ZICpHILIC column (5 µm, 2.1 × 150 mm) at 15 • C (Merck, Darmstadt, Germany), using a gradient program of an aqueous buffer of 10 mM ammonium carbonate pH 10.5 (phase A) and acetonitrile (phase B). After LC-HRMS analysis, QC samples were re-injected for performing additional higher-energy collisional dissociation (HCD) tandem mass spectrometry experiments (MS/MS) in both negative and positive ionization modes, with the instrument set in the targeted acquisition mode, using inclusion lists. HCD mass spectra were inspected manually to confirm annotations of detected metabolites.

Data Processing and Statistical Analysis
Data processing workflow and statistical analyses were performed on the open webbased platform workfow4metabolomics, a collaborative research infrastructure for computational metabolomics (W4M: https://workflow4metabolomics.usegalaxy.fr [22], version 4.0.0, accessed on 8 July 2020). Automatic peak detection and integration were performed using the matched filter algorithm in the W4M pre-processing package (including XCMS software). XCMS features were thereafter filtered according to the following criteria: (i) a correlation between dilution factors of QC samples and areas of chromatographic peaks superior to 70%, (ii) repeatability (the coefficient of variations obtained for chromatographic peak areas of QC samples should be below 30%), and (iii) ratio of chromatographic peak area of biological to blank samples above a value of 3. Features were also annotated by matching their accurate measured masses and chromatographic retention times with our spectral database obtained from~1000 pure, authentic standards [16]. In addition, using the present untargeted metabolomics approach, some isomeric metabolites cannot be resolved. Thus, some chromatographic peaks could correspond to more than 1 metabolite. This is typically the case for glucose, galactose, and other hexose isomers or also for pentose phosphates. MS/MS patterns of metabolites that are particularly mentioned in our work were compared with MS/MS patterns of authentic standards to validate their identification, as proposed by the Metabolomics Standards Initiative [23]. Statistical analyses were performed with W4M. Normal distribution of was evaluated with a Shapiro-Wilk normality test. Echocardiography and catheterization parameters were considered as discriminant when p-value < 0.05 with a t-test or Mann-Whitney test, applied depending on normality of measures distribution. Metabolites were considered as discriminant when showing both p-value < 0.05 with Mann-Whitney statistical test and VIP > 1 in multivariate analysis (OPLS-DA).

RV Remodeling after CH and MCT-Exposure in Rats
CH and MCT-induced PH models were compared with their respective controls (normoxia and vehicle) ( Figure 1A). Principal component analysis (PCA) based on echocardiographic and right heart catheterization parameters showed similarities between normoxia and vehicle rats compared with both CH-and MCT-PH rats. Of interest, PCA showed strong differences between CH and MCT-PH rats ( Figure 1B). However, both CH-and MCT-PH rats developed severe PH, while normoxia and vehicle rats showed normal pulmonary pressure. This was highlighted by the strong increase in RV systolic pressure (RVS.P) ( Figure 1C,E), pulmonary vascular resistance ( Figure 1D), end-diastolic pressure (EDP) ( Table 1) as well as by the decrease in the pulmonary artery acceleration time (PAAT) in CH-PH and MCT-PH rats compared to their respective controls ( Table 1). As a consequence, CH-and MCT-exposed rats developed RV hypertrophy as indicated by the increased Fulton's index and RV free wall thickness ( Figure 1F,G). Interestingly, CH exposure induced an adaptive RV phenotype as opposed to MCT exposure which induced a maladaptive RV phenotype. Indeed, in the CH-PH model, the RV contractility was preserved as indicated by RV contractility index values (121.7 ± 13.63 vs. 75.25 ± 14.7 for CH and MCT-PH rats, respectively) ( Table 1 and Figure 1H,I) and RV hypertrophy (calculated by Fulton's index as the weight ratio of RV and (LV + septum)) (0.63 ± 0.12 vs. 0.49 ± 0.11) was more pronounced than in MCT-PH rats (Table 1 and Figure 1F). On the opposite, RV was dysfunctional in MCT-PH rats. As a matter of fact, the elevation of RVS.P and PVR was more pronounced in MCT-PH (89.46 ± 26.02) rat group than in CH-PH rat group (68.26 ± 10.83), suggesting that PH was more severe in MCT-PH rats compared to CH-PH rats (Table 1). Moreover, cardiac output (CO) was more markedly decreased in MCT-PH-rats ( Figure 1H), which could be a consequence of compromised LV filling or associated with RV dysfunction in MCT-PH rats. To support this hypothesis, we found that the RV contractility index ( Figure 1I) was reduced in MCT-PH compared to CH-PH rats and control animals. We also found that the tricuspid annular plane systolic excursion (TAPSE) was decreased in MCT-PH rats and not decreased in CH-vs. normoxia rats (Table 1). No significant changes were observed in RVET between CH and MCT rats. LV echocardiography analysis revealed that LV fractional shortening (LV FS) was similar in CH-PH and MCT-PH compared to their respective controls, while the Aorta velocity integral time (VTI-Ao) was reduced in MCT-PH rats compared to vehicle-exposed rats ( Table 1). significant changes were observed in RVET between CH and MCT rats. LV echocardiography analysis revealed that LV fractional shortening (LV FS) was similar in CH-PH and MCT-PH compared to their respective controls, while the Aorta velocity integral time (VTI-Ao) was reduced in MCT-PH rats compared to vehicle-exposed rats (Table 1).  Supplementary Table S1. Significance: * p ≤ 0.05; ** p ≤ 0.01; *** p ≤ 0.001, **** p ≤ 0.0001.  In addition to hemodynamic characterization of PH, we also analyzed pulmonary vascular remodeling in MCT-PH and CH-PH rats using hematoxylin and eosin staining of paraffin-embedded lung section. We found a significant increase in the vascular wall thickness in MCT-PH rats and in CH-PH rats compared with their respective controls ( Figure 2A,B). Importantly, we measured that the vascular wall thickness was more important in lung from MCT-PH rats compared with lung from CH-PH rats (Figure 2A,B).
Three weeks after PH induction by MCT or CH exposure, we also measured the presence of RV fibrosis using trichrome Masson staining ( Figure 2). We found a strong increase in RV fibrosis in CH-exposed rats and MCT-PH rats ( Figure 2C,D), and the RV fibrosis was more pronounced in RV from MCT-PH compared with RV from CH-PH rats ( Figure 2C,D).
In addition to hemodynamic characterization of PH, we also analyzed pulmonary vascular remodeling in MCT-PH and CH-PH rats using hematoxylin and eosin staining of paraffin-embedded lung section. We found a significant increase in the vascular wall thickness in MCT-PH rats and in CH-PH rats compared with their respective controls (Figure 2A,B). Importantly, we measured that the vascular wall thickness was more important in lung from MCT-PH rats compared with lung from CH-PH rats (Figure 2A,B).  . Pulponary vessel occlusion (%) in the control, the control, MCT, Vehicle, CH and Nx groups (n = 4 different rats per condition). (C). Interstitial fibrosis was identified with trichrome Masson staining in the RV compartments of control, the control, MCT, Vehicle, CH, and Nx groups. Scale bar = 100 µm (D). Quantification of the percentage of fibrosis in RV tissues from control and MCT-exposed rats (n = 20 images per rat from 4 rats). t-tests were used after verification of normal distribution of values (Shapiro-Wilk normality test). Significance: ** p ≤ 0.01; *** p ≤ 0.001.

PH-Induced Cardiac Remodeling Associated with Alteration of Four Metabolic Pathways
KEGG database [24] was consulted to determine involved metabolic pathways (Supplementary Table S2). Arginine, purine, pyrimidine, and tryptophan metabolic pathways were altered in both models ( Table 2). The arginine pathway was the most altered pathway with 18 altered metabolites in CH-PH rats (12 in plasma, 14 in RV and LV tissues) and 12 altered metabolites in MCT rats (5 in plasma, 11 in heart tissues). This pathway showed the strongest alterations with seven metabolites with fold changes greater than 1.8 in heart tissues and significant statistics with p-values < 0.001 for half of the significant results. Pyrimidine pathway was the second most altered pathway with 11 altered metabolites in CH-PH rats (10 in plasma and 8 in heart tissues) and 12 altered metabolites in MCT rats (7 in plasma and 10 in heart tissues). An overall increase in metabolites from the pyrim-idine pathway was observed in both MCT and CH-PH rats and systematically detected in RV of both CH-PH rats and MCT-PH rats (unless two metabolites were specifically altered in plasma). As opposed to the arginine and pyrimidine pathways that were globally altered in both heart and plasma, alterations of metabolites from purine metabolism were mostly found in heart tissues. Over the 10 metabolites altered in the CH-PH model, 9 were altered in heart tissues versus 3 in plasma. Similarly, from the eight metabolites altered in MCT rats, all were altered in the heart but only two in plasma. Alteration of tryptophan metabolism followed another different pattern since half of tryptophan metabolism showed a specific plasmatic alteration. The tryptophan pathway was represented by six altered metabolites in CH-PH rats (five in plasma and four in heart tissues) and eight alterations in MCT rats (six in plasma and five in heart tissues). The results from the tryptophan pathway were less significant and difficult to interpret due to a mismatch between MCT-PH and CH-PH rats.

Five Metabolites as PH-Induced Cardiac Remodeling Metabolic Signature
Alteration of polyamines (i.e., spermidine, acetyl-spermidine, and putrescine) was only observed in LV and RV tissues with p-values < 0.01, and fold changes were among the most important FC observed in the current study (Table 2). Particularly, acetyl-spermidine was found to be elevated more than twice in the heart of PH-models (Figure 3(B1,C1)). In addition, arginine and guanidinoacetate were altered in plasma and heart (LV and RV from both CH and MCT rats). Guanidinoacetate was the most altered metabolite of the arginine metabolic pathway with FC > 3 and p < 0.0001 in CH-PH RV compared to normoxia rats (Table 3 and Figure 3(A2,B2,C2)). Arginine decreased in both plasma and LV and RV tissues ( Figure 3A3,B3,C3), in accordance with the broad literature about arginine in PH and RV failure [14,[25][26][27][28][29]. Orotate was the most impacted metabolite from the pyrimidine metabolic pathway, with fold changes >2 in RV of both MCT and CH models and with highly significant p-values in CH-PH heart and plasma (p < 0.0001) ( Table 2 and Figure 3(A4,B4,C4)). Concerning the purine metabolic pathway, allantoin was the only metabolite altered in plasma and LV and RV tissues of both CH-PH and MCT-PH rats (Figure 3(A5,B5,C5)). Furthermore, xanthosine, xanthine, guanosine, and inosine were specifically altered in RV of both CH-PH and MCT-PH models.

Focus on Arginine Bioavailability and Enzymatic Activities
Arginine bioavailability was estimated by the arginine-to-ornithine ratio (which reflects arginase activity) and arginine-over-guanidinoacetate ratio (which reflects arginine:glycine amidinotransferase activity). A significant decrease in the arginine-overornithine ratio illustrated a stimulation of arginase in plasma and RV of both CH-PH and MCT-PH rat models (Supplementary Figure S2). Similarly, arginine-over-guanidinoacetate ratios were reduced in plasma, RV and LV of CH-PH rat, reflecting an increase in arginine:glycine amidinotransferase activity. These ratios evaluation underlined a double source of arginine consumption, reducing arginine bioavailability.

Correlation between Metabolic Signature and RV Function
(arg pathway) Levels of acetylspermidine and putrescine in heart tissues were correlated to RVS.P, Fulton's index, PVR, and PAAT in both models but particularly in CH-PH rats with overall absolute values of correlation coefficients superior to 0.7 (Table 3). Guanidinoacetate correlated with 11 echocardiography and hemodynamic parameters. (pyruvate pathway) Moreover, orotate RV levels were highly correlated with clinical parameters, since they showed R > 0.8 and p < 0.0001 in each organ for PVR correlation and R > 0.8 and p < 0.00001 in RV for RVS.P, RV weight, and Fulton's index correlations in CH-PH rats (Table 3). Cytosine, deoxycytidine, deoxyuridine, and plasmatic thymidine levels were also highly correlated with RVS.P, RV weight, and Fulton's index. (purine pathway) Allantoin levels were the most correlated to RVS.P, RV weight, Fulton's index, and PVR in CH-PH rats (Table 3). Xanthosine, xanthine, guanosine, and inosine levels were not correlated with echocardiography and catheterization parameters.

Discussion
To the best of our knowledge, this is the first time that metabolic disorders have been studied from heart tissues in PH animal models. We analyzed plasma and heart tissues with an untargeted metabolomic approach to detect and compare variations of more than 280 metabolites in two different rat models (MCT-PH and CH-PH). Metabolomic analysis showed that four main metabolic pathways (arginine, pyrimidine, purine, and tryptophan) were altered in both models. Among these metabolic pathways, many metabolite levels from arginine, pyrimidine, tryptophan, and purine pathways were highly correlated with cardiac morphometric, RV, and closed-chest heart catheterization parameters in the adaptive phenotype of CH-PH models, suggesting their involvement in PH-induced cardiac remodeling and their potential role as specific biomarkers for RV in PH.
RV remodeling and dysfunction in MCT-PH rats have been often associated with maladaptive RV (reduction in exercise capacity, reduction in RV function, capillary rarefaction, RV inflammation). In CH-PH rats, the pulmonary arterial medial hypertrophy leads to elevated RVS.P and RV hypertrophy without RV failure [30]. In the present study, we confirmed that MCT exposure induces severe PH in rats with maladaptive RV while rats exposed to CH develop severe PH with adaptive RV. Therefore, we confirmed from plasma and heart biopsy samples that MCT and CH-rats are complementary models as suggested by Zhao et al. [15].

Arginine and Tryptophan Metabolic Pathways, a Duel for the Control of Pulmonary Arterial Tone in PH Plasma
Arginine metabolism was the most impacted pathway in accordance with previous studies that described such alteration in plasma [12,15] and lungs of MCT-PH rats [25]. In the urea cycle, arginine can be metabolized competitively by nitrite oxide synthase (NOS) to produce nitric oxide (NO) or by arginase to synthesize ornithine and then polyamines (Figure 4(A1,B1)). This competitive metabolism of arginine regulates the NO/polyamines balance [27]. The role of polyamines in the physiopathology of PH is unclear. Indeed, while some studies linked spermidine and putrescine with an increase in RVS.P in CH-induced pulmonary vascular (PV) remodeling [31], others demonstrated their anti-inflammatory effects [32] and, therefore, their putative benefit against PH [33]. On the opposite, the critical role of NO for vasodilatation is well-referenced [34]. Numerous studies linked the progression of PH with NO alteration [35][36][37] and demonstrated that the urea cycle played a major role in PH severity [14,38]. Alteration of arginine metabolism leads to regulation of tryptophan metabolism as already described in the literature about PH [11,15]: Indeed, Lewis et al. reported that indoleamine 2,3-dioxygenase (IDO), the heme-containing enzyme metabolizing tryptophan into kynurenine, was over-activated by a decrease in NO as a compensatory vasodilator pathway [12,39] (Figure 4A3,B3). IDO-deficient mice developed exaggerated hypoxic pulmonary vasoconstriction [40]. Moreover, lung mRNA expression of IDO was increased in PH models, and IDO-tryptophan metabolites were correlated with RV dysfunction [12]. However, alteration of tryptophan metabolism has specifically been found altered in pulmonary circulation [12], explaining why kynurenate and xanthurenate were not altered in heart tissues in our study. Tryptophan can also be metabolized by tryptophan hydroxylase (TPH), which yields serotonin. Our results demonstrated alterations of serotonin in the plasma of CH-PH rats. Serotonin has been involved in the pathophysiology of PH, contributing to vasoconstriction and vascular remodeling in PH lungs [41,42]. In our study, a decrease in serotonin was observed in the plasma of the carotid artery, suggesting that tryptophan is mainly metabolized by IDO in our model. Indeed, Lewis et al. reported that indoleamine 2,3-dioxygenase (IDO), the heme-containing enzyme metabolizing tryptophan into kynurenine, was over-activated by a decrease in NO as a compensatory vasodilator pathway [12,39] (Figure 4A3,B3). IDO-deficient mice developed exaggerated hypoxic pulmonary vasoconstriction [40]. Moreover, lung mRNA expression of IDO was increased in PH models, and IDO-tryptophan metabolites were correlated with RV dysfunction [12]. However, alteration of tryptophan metabolism has specifically been found altered in pulmonary circulation [12], explaining why kynurenate and xanthurenate were not altered in heart tissues in our study. Tryptophan can also be metabolized by tryptophan hydroxylase (TPH), which yields serotonin. Our results demonstrated alterations of serotonin in the plasma of CH-PH rats. Serotonin has been involved in the pathophysiology of PH, contributing to vasoconstriction and vascular remodeling in PH lungs [41,42]. In our study, a decrease in serotonin was observed in the plasma of the carotid artery, suggesting that tryptophan is mainly metabolized by IDO in our model.

A Balance between Vasoconstriction and Vasodilatation in the Heart
For the first time, we showed stimulation of arginase activity in heart biopsies from PH models, which suggests that vasoconstriction induced by NO depletion could also affect cardiac arterial tone. Moreover, an alteration of purine metabolism was mainly spe-

A Balance between Vasoconstriction and Vasodilatation in the Heart
For the first time, we showed stimulation of arginase activity in heart biopsies from PH models, which suggests that vasoconstriction induced by NO depletion could also affect cardiac arterial tone. Moreover, an alteration of purine metabolism was mainly specific to heart tissues ( Figure 4(A4,B4)). Previous studies analyzing the role of purine metabolic pathways in PH mostly focused on serum urate, and cardiac alteration of purines has been weakly investigated in this disease [43,44]. Apart from the PH context, cardiac efflux of purine metabolites triggered by ischemia or CH are well referenced [45]. Several studies correlated alteration of purine metabolism with a decrease in functional capacities in patients with cardiac failure [46,47] and deficient oxidative metabolism during hypoxic states [12,13,46]. As an end-product of purine metabolism, allantoin was the most impacted metabolite and the most correlated with RV-PV function parameters, as previously reported [12]. Allantoin has been associated with central antihypertensive effects in rats [48]. Our results suggest that allantoin would be a putative biomarker of cardiac hypoxia during PH. In addition, kynurenine, from the tryptophan pathway, was increased in the heart of CH-PH rats and correlated with PVR, RVS.P, RV weight, and Fulton's index (Figure 4(A3,B3)). This metabolite has been described as a putative endothelium-derived vasodilator in several models [49]. Thus, our results might suggest a protective role of purine and tryptophan metabolic pathways in response to vasoconstriction in heart of PH rats, suggesting that these pathways could be important for RV remodeling during PH due to hypoxemia.

Polyamine Mediated Cardiac Hypertrophy
We demonstrated that putrescine and acetyl-spermidine were increased in both functional and dysfunctional hypertrophied RV in PH. Indeed, super activated arginase led to a strong increased in polyamines (putrescine, spermidine, and acetyl-spermidine), specifically in the heart of CH and MCT-PH rats and particularly in RV (Figure 4(A1,B1)). This alteration was highly correlated to RVS.P, RV weight, Fulton's index, and PVR. Yet, polyamines have been known to contribute to pathological processes in the heart, including cardiac hypertrophy [50][51][52][53]. Particularly, the increase in polyamines is mediated by pulmonary artery constriction, enhanced cardiac work consecutive to hypoxic PH, and reduced availability of oxygen after cardiac overload [53]. Moreover, several studies reported that polyamines, and particularly spermidine, are thought to be responsible for the induction of autophagy [54], a critical process involved in cardiac hypertrophy [55]. Furthermore, polyamines can affect the function of many ion channels, including the inward rectifying K + current, NMDA receptor, cyclic nucleotide-gated channels, voltage-gated sodium channel in different cell types [56][57][58]. By interaction with 14-3-3 proteins, increased polyamines may also downregulate TASK-1 channels [59][60][61], which contributes to the reduction in RV cardiomyocytes' excitability [30]. Moreover, we previously demonstrated that loss of TASK-1 function was a hallmark of RV hypertrophy/dysfunction associated with PH, including in RV from MCT-PH and CH-PH models [62]. Thus, inhibition of these ion channel families by increased polyamines could have profound consequences on RV function, including resting membrane potential, action potential duration, and excitation, which are observed in isolated cardiomyocytes from different PH animal models [30]. Therefore, significant evidences emerge for the role of polyamines in cardiac pathophysiology and point out these metabolites as putative targets for the treatment of cardiovascular diseases [52].

Metabolic Alteration Related to Cardioprotection
Pyrimidine metabolism was the second most impaired pathway in our study (Figure 4(A2,B2)), with a major increase in uridine and orotate, as already reported in PH and heart ischemia [12,63]. Due to the high myocardial turnover of pyrimidine metabolites, any alteration in their utilization or supply could have severe metabolic consequences [64]. In case of alteration, cardiac de novo synthesis of pyrimidines is inefficient, and precursors are mainly released by the liver and transferred through the bloodstream [65]. This fact might explain why we detected almost all pyrimidine alterations in both plasma and heart of CH-PH rats. Cardioprotection of pyrimidine metabolites, such as uridine and orotate, has been demonstrated for more than twenty years [63]. An international symposium about the cardioprotective effect of orotate highlighted that this metabolite (i) improves the tolerance of infarcted heart to ischemia, (ii) reduces the severity of chronic myocardial dysfunction, and (iii) improves contractile function in rats with hypertrophying hearts [64]. These effects could be explained by the fact that orotate increases pyrimidines' base de novo-production in the liver, which consequently activates a metabolic cascade preventing ATP depletion in the myocardium [64,66,67]. We suggest an orotate increase in MCT-and CH-PH rats as a cardioprotective response to the cardiac vasoconstriction and hypertrophy.
Our results contribute to the understanding of PH pathogenesis as well as RV remodeling. Thanks to metabolomics, some hypotheses were formulated about the communication between four main metabolic pathways ( Figure 5). Arginine metabolism is supposed to be central in the disease outbreak through the depletion in NO and, therefore, the induction of vasoconstriction at both peripheral and cardiac levels. Some compensatory reactions are probably promoted to slow down the rise in arterial tone thanks to the stimulation of IDO activity (in plasma and heart) and synthesis of allantoin (specifically in heart). We assume that increased polyamines are involved in the development of cardiac hypertrophy and that cardiac injury is partly offset by several cardioprotective pyrimidines. These different factors stand as putative new therapeutic targets or candidate biomarkers for the early diagnostic of cardiac impairment secondary to PH. Such biomarkers could be detected in blood samples or used as PET tracers when limited to the heart compartment. orotate, has been demonstrated for more than twenty years [63]. An international symposium about the cardioprotective effect of orotate highlighted that this metabolite (i) improves the tolerance of infarcted heart to ischemia, (ii) reduces the severity of chronic myocardial dysfunction, and (iii) improves contractile function in rats with hypertrophying hearts [64]. These effects could be explained by the fact that orotate increases pyrimidines' base de novo-production in the liver, which consequently activates a metabolic cascade preventing ATP depletion in the myocardium [64,66,67]. We suggest an orotate increase in MCT-and CH-PH rats as a cardioprotective response to the cardiac vasoconstriction and hypertrophy.
Our results contribute to the understanding of PH pathogenesis as well as RV remodeling. Thanks to metabolomics, some hypotheses were formulated about the communication between four main metabolic pathways ( Figure 5). Arginine metabolism is supposed to be central in the disease outbreak through the depletion in NO and, therefore, the induction of vasoconstriction at both peripheral and cardiac levels. Some compensatory reactions are probably promoted to slow down the rise in arterial tone thanks to the stimulation of IDO activity (in plasma and heart) and synthesis of allantoin (specifically in heart). We assume that increased polyamines are involved in the development of cardiac hypertrophy and that cardiac injury is partly offset by several cardioprotective pyrimidines. These different factors stand as putative new therapeutic targets or candidate biomarkers for the early diagnostic of cardiac impairment secondary to PH. Such biomarkers could be detected in blood samples or used as PET tracers when limited to the heart compartment.

Limitations
Using experimental PH animal models allowed myocardium metabolic disturbances to be studied in comparison to plasmatic alterations, but there is no animal model fully mimicking pathological features of humans, and clinical trials are necessary to confirm

Limitations
Using experimental PH animal models allowed myocardium metabolic disturbances to be studied in comparison to plasmatic alterations, but there is no animal model fully mimicking pathological features of humans, and clinical trials are necessary to confirm our observations. Moreover, the current study did not allow gender influence in RV dysfunction to be considered since only male rats were analyzed for their higher sensitivity to experimental PH. In our study, MCT rats gave less homogenous and interpretable results than CH rats. Moreover, the CH model is considered as inducing a weaker intensity of obstructive remodeling and modest hemodynamic abnormalities, which can make the alterations less representative of severe human cardiac remodeling. As we know that a greater proportion of female PAH patients have been consistently observed, the metabolic signature of RV failure could be investigated in females in future studies. As RV and LV tissues are composed of several cell types, including endothelial cells, fibroblasts, cardiomyocytes, inflammatory cells, we plan to analyze the metabolic signature of adult rat ventricular myocytes in different experimental conditions in future studies. At least, the study explored the surface of metabolite change in RV; neither a theory nor validation of these findings was conducted. One should emphasize that metabolomic differences between groups may be due to rat model differences rather than RV phenotype.
In conclusion, this study investigated the metabolomic signature of cardiac remodeling on two PH rat models using a robust untargeted metabolomic platform coupled to an inhouse database for the identification of hundreds of metabolites. Predominant alterations of arginine, pyrimidine, purine, and tryptophan metabolic pathways were detected on both plasma and heart tissue samples. These pathways are known to be involved in vasodilatation, cell proliferation, hypertrophy, and cardioprotection in either PH or cardiac ischemia. Here, for the first time, we gave new insights into the organ origin of their alteration and brought new insights about their implication, specifically in cardiac failure during PH. Metabolites from arginine, pyrimidine, tryptophan, and purine pathways that are highly correlated with cardiac clinical parameters were highlighted as candidate biomarkers of RV remodeling or as putative therapeutic targets. These finding open novel research perspectives to find biomarkers for early detection of RV failure in PH, and further clinical trials are needed to confirm our hypotheses in patients.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/cells10061559/s1, Figure S1: Representative echocardiographic illustrations allowing us to measure all parameters presented in the study, Figure S2: Multivariate analyses on metabolomic analyses of CH-PH and MCT-PH sample sets, Figure S3: Arginine over ornithine ratios and arginine over guanidinoacetate ratios in CH-PH and MCT-PH rats and their respective controls. Table S1: Summary of metabolomic analyses of plasma and heart of CH-PH and MCT-PH rats.

Informed Consent Statement: Not Applicable.
Data Availability Statement: Row data were submitted for paper review process and are available on demand via the corresponding author.

Conflicts of Interest:
The sponsors had no role in the design, collection, analysis, or interpretation of data, the writing of this article, or the decision to submit it for publication. M.H. has relationships with drug companies, including Actelion, Bayer, GSK, Novartis, and Pfizer. In addition to being an