# The extent of postpartum cardiac reverse remodeling is reflected in urine proteome

Ana F Ferreira, Fábio Trindade, Maria J Azevedo, Juliana Morais, Thibaut Douché, Sílvia O Diaz, Francisca A Saraiva, Carla Sousa, Ana P Machado, Mariette Matondo, et al.

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

## To cite this version:

Ana F Ferreira, Fábio Trindade, Maria J Azevedo, Juliana Morais, Thibaut Douché, et al. The extent of postpartum cardiac reverse remodeling is reflected in urine proteome. Scientific Reports, 2024, 14 (1), pp.14815. 10.1038/s41598-024-65612-1. pasteur-05076147

HAL Id: pasteur-05076147

https://pasteur.hal.science/pasteur-05076147v1

Submitted on 21 May 2025

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

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

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

OPEN

# The extent of postpartum cardiac reverse remodeling is reflected in urine proteome

Ana F. Ferreira<fnref n="1" />, Fábio Trindade<fnref n="1" />, Maria J. Azevedo²,³,⁴,⁵, Juliana Morais<fnref n="1" />, Thibaut Douché⁶, Sílvia O. Diaz<fnref n="1" />, Francisca A. Saraiva<fnref n="1" />, Carla Sousa<fnref n="1" />,⁷, Ana P. Machado⁸, Mariette Matondo⁶, Adelino Leite-Moreira<fnref n="1" />,⁹, Carla Ramalho⁹,¹⁰, Rui Vitorino<fnref n="1" />, Inês Falcão-Pires<fnref n="1" /> & António S. Barros<fnref n="1" />

The association of postpartum cardiac reverse remodeling (RR) with urinary proteome, particularly in pregnant women with cardiovascular (CV) risk factors who show long-term increased risk of cardiovascular disease and mortality is unknown. We aim to profile the urinary proteome in pregnant women with/without CV risk factors to identify proteins associated with postpartum RR. Our study included a prospective cohort of 32 healthy and 27 obese and/or hypertensive and/or diabetic pregnant women who underwent transthoracic echocardiography, pulse-wave-velocity, and urine collection at the 3rd trimester and 6 months postpartum. Shotgun HPLC-MS/MS profiled proteins. Generalized linear mixed-effects models were used to identify associations between urinary proteins and left ventricle mass (LVM), a surrogate of RR. An increase in arterial stiffness was documented from 3rd trimester to 6 months after delivery, being significantly elevated in women with CV risk factors. In addition, the presence of at least one CV risk factor was associated with worse LVM RR. We identified 6 and 11 proteins associated with high and low LVM regression, respectively. These proteins were functionally linked with insulin-like growth factor (IGF) transport and uptake regulation by IGF binding-proteins, platelet activation, signaling and aggregation and the immune system's activity. The concentration of IGF-1 in urine samples was associated with low LVM regression after delivery. Urinary proteome showed a predicting potential for identifying pregnant women with incomplete postpartum RR.

**Keywords** Pregnancy, Cardiovascular reverse remodeling, Cardiovascular risk factor, Hemodynamic overload, Urine, Proteome

**Abbreviations**

AHNAK Neuroblast differentiation-associated protein

ARF1 ADP-ribosylation factor 1

ATP1A1 Sodium/potassium-transporting ATPase subunit alpha-1

CV Cardiovascular

CDC42 Cell division control protein 42 homolog

C3 Complement C3

CTSH Pro-cathepsin H

<footnote n="1">¹Cardiovascular R&D Centre – UnIC@RISE, Department of Surgery and Physiology, Faculty of Medicine of the University of Porto, Alameda Prof. Hernâni Monteiro, 4200-319 Porto, Portugal. ²Faculdade de Medicina Dentária, Universidade do Porto, 4200-393 Porto, Portugal. ³INEB - Instituto Nacional de Engenharia Biomédica, 4200-135 Porto, Portugal. ⁴i3S - Instituto de Investigação e Inovação em Saúde, Universidade do Porto, 4200-135 Porto, Portugal. ⁵Academic Center for Dentistry Amsterdam, University of Amsterdam and Vrije Universiteit Amsterdam, 1081 LA Amsterdam, The Netherlands. ⁶Proteomic Platform, Mass Spectrometry for Biology Unit, CNRS UAR 2024, Institut Pasteur, Université Paris Cité, 75015 Paris, France. ⁷Cardiology Department, Centro Hospitalar Universitário de São João, 4200-319 Porto, Portugal. ⁸Center of Prenatal Diagnosis, Obstetrics Department, Centro Hospitalar Universitário de São João, 4200-319 Porto, Portugal. ⁹Cardiothoracic Surgery Department, Centro Hospitalar Universitário de São João, 4200-319 Porto, Portugal. ¹⁰Obstetrics, Gynaecology and Pediatrics Department, Faculty of Medicine of the University of Porto, 4200-319 Porto, Portugal. ¹¹email: asbarros@med.up.pt</footnote>
<table><tr><td>DM</td><td>Diabetes mellitus</td></tr><tr><td>DSG2</td><td>Desmoglein-2</td></tr><tr><td>FN1</td><td>Fibronectin</td></tr><tr><td>GLA</td><td>Alpha-galactosidase A</td></tr><tr><td>GLMM</td><td>Generalized linear mixed-effects model</td></tr><tr><td>GSTP1</td><td>Glutathione S-transferase P</td></tr><tr><td>HDO</td><td>Combination of hypertensive, diabetic and obese women</td></tr><tr><td>IGF</td><td>Insulin-like growth factor</td></tr><tr><td>IGF-1</td><td>Insulin-like growth factor-1</td></tr><tr><td>ITIH4</td><td>Inter-alpha-trypsin inhibitor heavy chain H4</td></tr><tr><td>LC-MS/MS</td><td>Liquid chromatography tandem mass spectrometry method</td></tr><tr><td>LA</td><td>Left atrium</td></tr><tr><td>LV</td><td>Left ventricle</td></tr><tr><td>LVM</td><td>Left ventricle mass</td></tr><tr><td>PWV</td><td>Pulse wave velocity</td></tr><tr><td>RR</td><td>Reverse remodeling</td></tr><tr><td>SIAE</td><td>Sialate O-acylesterase</td></tr><tr><td>SORT1</td><td>Sortilin</td></tr><tr><td>TF</td><td>Serotransferrin</td></tr><tr><td>TGOLN2</td><td>Trans-golgi network integral membrane protein 2</td></tr><tr><td>TKT</td><td>Transketolase</td></tr><tr><td>3T</td><td>Third trimester</td></tr><tr><td>UBA 1</td><td>Ubiquitin-like modifier-activating enzyme 1</td></tr></table>

Pregnancy induces significant multisystemic changes to ensure a successful adaptation to fetal growth and development. Pregnancy's cardiovascular response to blood volume expansion triggers an increase in preload accompanied by a reduction of peripheral vascular stiffness and a subsequent compensatory decline in afterload¹,². These circulatory adaptations induce cardiac remodeling, involving non-pathological left ventricle (LV) eccentric hypertrophy and slight diastolic dysfunction with left atrium (LA) enlargement¹,³ but with preserved ejection fraction¹,³. This physiological cardiovascular (CV) remodeling may be influenced by CV risk factors, such as arterial hypertension, gestational diabetes and obesity. Indeed, the CV remodeling in pregnant women with any of these CV risk factors is characterized by diastolic dysfunction combined with an increased left ventricular mass (LVM) and increased systemic vascular resistances⁴,⁵. Some authors suggest that women's heart fully recovers without cardiovascular complications during gestation and that global and segmental myocardial performances return to their pre-gravid state/condition in a process known as reverse remodeling (RR)²,³. Although it has been reported that pregnant women with CV risk factors have an increased risk for potential cardiovascular disease and mortality in the long term, few studies have explored their impact on cardiac RR during postpartum⁶,⁷. We hypothesized that cardiac RR after delivery is associated with distinct proteomic traits in pregnant women with and without CV risk factors and that these are mirrored in urine. As a plasma filtrate, urine can be viewed as a vast collection of products from cardiovascular activity. Therefore, it should harbor many proteins potentially differentiating physiological from pathological cardiac RR. Besides, urine analysis brings additional advantages when compared to plasma/serum, including sample stability (the bulk of proteolytic activity is completed in the bladder), lower matrix complexity (protein dynamic range in urine is about half of plasma's), and the non-invasive collection⁸. Urine proteomics allows for uncovering potential markers and for developing predictive tools, such as the recently reported Heart Failure Predictor, a classifier based on 96 urinary peptides⁹. Identifying pregnant women with a higher risk of a pathological cardiac RR is important to establish a tight monitoring program and eventually extend the follow-up until full cardiac RR is attained. Thus, this study aimed to profile the urinary proteome in pregnant women with or without CV risk factors and identify putative proteins associated with cardiac postpartum RR.

# Methods

## Study design and participants recruitment

We designed a prospective cohort study¹⁰, whose participants were recruited at the Department of Obstetrics of Centro Hospitalar Universitário São João (CHUS) and Unidade Local de Saúde de Matosinhos—Hospital Pedro Hispano (ULSM-HPH). The Ethics Committees of CHUS and ULSM-HPH approved the study (ID 201/18 and ID 154/20/RS, respectively). All included participants provided written informed consent. The confidentiality and data anonymity complied with the Declaration of Helsinki.

Participants were recruited from March 2019 to March 2021 at their first medical appointment (in the 1st or 3rd trimester of pregnancy) or voluntarily through online forms available at https://perimyrobb.wordpress.com/.

The inclusion criteria were adult pregnant women (> 18 years old) with or without CV risk factors, namely chronic and/or gestational hypertension and/or gestational diabetes mellitus (DM) or type-2 DM, and/or obesity. Arterial hypertension was defined as systolic blood pressure (SBP) ≥ 140 mmHg and/or diastolic blood pressure (DBP) ≥ 90 mmHg measured in office or in-hospital before 20 weeks of gestation. Gestational hypertension was defined as arterial hypertension diagnosed after 20 weeks of gestation and resolving within 42 days postpartum. Gestational DM was diagnosed if 92 ≤ fasting-glucose ≥ 126 mg/dL at 1st trimester or ≥ 180 mg/dL or ≥ 153 mg/dL, respectively, 1 or 2 h after an oral glucose tolerance test (75 g oral glucose load) performed at 24–28 pregnancy weeks. Type-2 DM was diagnosed through the following criteria: glycated haemoglobin ≥ 6.5% or fasting plasma
glucose ≥ 126 mg/dL or two-hour post-load venous plasma glucose ≥ 200 mg/dL after oral glucose tolerance test. Obesity was considered if the body mass index ≥ 30 kg/m² before pregnancy.

Women with twin pregnancies, pre-existing cardiomyopathy, renal disease, chronic obstructive airway disease, active systemic infection, genetic syndromes or type-1 diabetes mellitus were excluded.

Two study groups were set according to the presence or absence of CV risk factors defined above, resulting in the definition of a healthy group and an HDO group (which included any combination of hypertensive, diabetic and obese women).

## Measurements

In the 3rd trimester (30–35 weeks), at the peak of cardiac remodeling (when cardiovascular adaptations are noticeable) and 6–7 months after delivery (during RR), participants underwent clinical characterization (maternal cardiovascular health, health-related habits, smoking habits, parity, medical history, demographics, obstetric and perinatal outcomes were obtained from questionnaires or electronic medical records), transthoracic echocardiography, pulse wave velocity (PWV) measurement and urine collection. The LVM was defined as the surrogate for assessing cardiac RR.

## Echocardiographic assessment

Conventional transthoracic echocardiography evaluation was performed with a 3 MHz phased-array probe (ACUSON SC2000 PRIME™) by a single operator and measurements obtained from standard views according to European Society of Cardiology recommendations for chamber quantification and diastolic function evaluation (Supplementary Table 1)¹¹,¹². At least three cardiac cycle images were acquired for data analysis. Two certified cardiologists independently analyzed, interpreted and harmonized the results. Myocardial deformation was assessed in LV through strain and strain rate analysis by Syngo Velocity Vector Imaging software, version 3.5 (Siemens Healthcare, Erlangen, Germany).

## Vascular stiffness evaluation

After 10 min of rest, brachial blood pressure was measured. Complior® device (Alam Medical, France) quantified arterial stiffness by carotid-femoral PWV, calculated from carotid-femoral distance/transit time. The mean value of two PWV measurements (whose difference was ≤ 0.5 m/s) were performed in each session for each participant.

## Urinary proteomic profiling

Fifty milliliters of first-morning urine samples were collected, centrifuged (2370 × g, 15 min, 4 °C), and stored at –80 °C in each visit (3 T and 6 months after delivery). Protein was concentrated using 10 kDa-cutoff concentrators (Vivaspin 500, Sartorius Biotech, Göttingen, Germany). The protein retentate was solubilized in 0.3 M Tris pH 6.8 and 4% sodium dodecyl sulphate and quantified by a standard bicinchoninic acid assay (Pierce, ThermoFisher™, Rockford, IL, USA).

The urinary proteome was profiled by a label-free shotgun approach. Briefly, 100 µg of protein was incubated 1 h at RT in the dark for denaturation, reduction, and cysteine alkylation using a homogenization buffer composed of SDS 1%, TCEP 5 mM and chloroacetamide 20 mM in ammonium bicarbonate 50 mM, pH 8.

Protein digestion was carried out on an auto Single-Pot Solid-Phase-enhanced Sample Preparation (auto-SP3) on the Bravo (Agilent technologies) platform, with trypsin, for 12 h at 37 °C. Peptides were cleaned up using C18 cartridge (Agilent Technologies,), vacuum-dried and resuspended in acetonitrile (ACN) 2%, formic acid (FA) 0.1% buffer.

The peptides were separated and identified through a reverse-phase nanochromatographic system (Proxeon EASY-nLC 1200-Thermo Fisher Scientific) coupled to a Q Exactive Plus mass spectrometer (Thermo Fisher Scientific). Peptides were eluted with a multi-step gradient over 118 min. Mass spectra were acquired using Xcalibur software using a data-dependent Top 10 method with a survey scan (300–1700 m/z) at a resolution of 70,000 and MS/MS scans at a resolution of 17,500.

Raw data were analyzed using MaxQuant software version 2.0.3.0 using the Andromeda search engine. The MS/MS spectra were searched against a UniProt *Homo sapiens* database (download in 13/12/2021). One unique peptide to the protein group was required for the protein identification. For quantification, we only considered protein groups with two unique peptides. A false discovery rate (FDR) cutoff of 1% was applied at the peptide and protein levels.

Proteins showing to be independently associated with cardiac RR were subjected to protein-protein interaction and functional enrichment analyses with STRING (version 11.5) webtool (accession date: 28th March 2023). A score of 0.4 was set as the minimum threshold for consideration of the validated and putative interactions.

## Quantification of insulin-like growth factor-1 through ELISA assay kit

Insulin-like growth factor-1 (IGF-1, DY291, R&D Systems) was quantified in pairwise urine samples collected 3 T and 6 months after delivery. Three operators, blinded to clinical and echocardiographic information, assessed this ELISA assay kit on each urine sample using a commercially available kit following the manufacturer's instructions.
## Statistical analysis

Clinical parameters, physical assessment, echocardiographic and PWV values, and urinary proteins detected and quantified were gathered on a database and analyzed using *lmer*, *ggplot2* and *gtsummary* R packages, in the R statistical software environment version 4.2.1.

Continuous variables are expressed as median, 1st and 3rd quartiles (Q1, Q3), as appropriate. Data normality was checked after examining the histograms and Q-Q plots. Absolute values and relative frequencies are presented for categorical variables. Protein intensity and IGF-1 quantification were log2-transformed and scaled.

A multivariable generalized linear mixed-effects model (GLMM) was used to explore pregnancy-associated cardiac RR by including the 3rd trimester of pregnancy and 6 months after delivery. The LVM (including all LVM measured in 3rd trimester and 6 months after delivery) was defined as the only outcome in all GLMM performed. The follow-up time (3rd trimester and 6 months after delivery) was taken as a covariate, allowing the longitudinal analysis of cardiac (reverse) remodeling (using the left ventricular mass as surrogate). The analysis of LVM progression during pregnancy and postpartum was adjusted for CV risk factors, age at recruitment moment, parity, weight gain until 3rd trimester and evaluation moments (3rd trimester and 6 months after delivery) due to their clinical relevance. Only adjusted statistically significant proteins with an absolute effect size (beta coefficient) higher than 6 standard deviations and IGF-1 quantification were considered for exploratory analysis. A threshold of 6 standard deviations emerged from data-driven exploration of the volcano plot. This methodological approach prioritised minimising false positives given the sample size limitations and maximising the signal/noise ratio often present in proteomics analysis. *P* < 0.05 was deemed significant.

## Results

### Baseline demographic and clinical characterization

We included 59 participants that had paired urine samples at 32 [31; 34] weeks of gestation and 6 [6;7] months after delivery, with a median age of 34 [32; 37] years-old (50.8% nulliparous). The demographic and clinical characteristics of participants are shown in Table 1. Twenty-seven (45.8%) women were diagnosed with at least one CV risk factor prior to or during gestation (Table 1). The study groups were matched for smoking habits (*p* = 0.477), parity (*p* = 0.604) and obstetric outcomes (type of delivery, *p* = 0.279; gestational age at delivery, *p* = 0.291) (Table 1). In addition, women with CV risk factors were older (36 [33; 39] years versus 34 [30; 36]

<table><caption>**Table 1.** Demographic and clinical baseline characteristics. Values expressed by median [interquartile range]. Significant values are in bold.</caption><thead><tr><th></th><th>Healthy group n = 32</th><th>HDO group n = 27</th><th>p value</th></tr></thead><tbody><tr><td>Age, years</td><td>34 [30;36]</td><td>36 [33;39]</td><td><strong>0.049</strong></td></tr><tr><td>Height, meters</td><td>1.64 [1.60;1.69]</td><td>1.66 [1.60;1.70]</td><td>0.492</td></tr><tr><td>Pre-pregnancy body mass index, kg/m²</td><td>22.3 [20.8;24.3]</td><td>27.2 [22.2;31.7]</td><td><strong>0.001</strong></td></tr><tr><td>Cardiovascular risk factors</td><td></td><td></td><td></td></tr><tr><td>Arterial hypertension n (%)</td><td>n/a</td><td>16 (59.3)</td><td>n/a</td></tr><tr><td>Obesity n (%)</td><td>n/a</td><td>11 (40.7)</td><td>n/a</td></tr><tr><td>Gestational diabetes n (%)</td><td>n/a</td><td>12 (44.4)</td><td>n/a</td></tr><tr><td>Type 2 Diabetes n (%)</td><td>n/a</td><td>1 (3.7)</td><td>n/a</td></tr><tr><td>Pharmacological therapy during pregnancy</td><td></td><td></td><td></td></tr><tr><td>Aspirin</td><td>4 (12.5)</td><td>14 (51.9)</td><td><strong>0.002</strong></td></tr><tr><td>Metformin</td><td>0 (0.0)</td><td>11 (20.0)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Insulin</td><td>0 (0.0)</td><td>10 (18.2)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Beta-blockers</td><td>0 (0.0)</td><td>2 (7.4)</td><td>0.205</td></tr><tr><td>Nifedipine</td><td>0 (0.0)</td><td>1 (3.7)</td><td>0.458</td></tr><tr><td>Smoking habits</td><td></td><td></td><td></td></tr><tr><td>Non-smoker n (%)</td><td>18 (56.3)</td><td>20 (74.1)</td><td rowspan="5">0.477</td></tr><tr><td>Smoker n (%)</td><td>3 (9.4)</td><td>2 (7.4)</td></tr><tr><td>Stopped only during pregnancy n (%)</td><td>1 (3.1)</td><td>1 (3.7)</td></tr><tr><td>Stopped at pregnancy beginning n (%)</td><td>2 (6.3)</td><td>2 (7.4)</td></tr><tr><td>Ex-smoker n (%)</td><td>8 (25.0)</td><td>2 (7.4)</td></tr><tr><td>Primiparous women n (%)</td><td>15 (46.9)</td><td>15 (55.6)</td><td>0.604</td></tr><tr><td>Pregnancy weeks at delivery</td><td>39 [31;33]</td><td>39 [38;40]</td><td>0.291</td></tr><tr><td>Type of delivery</td><td></td><td></td><td></td></tr><tr><td>Vaginal delivery n (%)</td><td>19 (59.4)</td><td>20 (74.1)</td><td rowspan="2">0.279</td></tr><tr><td>Caesarean delivery n (%)</td><td>13 (40.6)</td><td>7 (25.9)</td></tr><tr><td>Evaluation moment</td><td></td><td></td><td></td></tr><tr><td>Pregnancy, weeks</td><td>31 [32;39]</td><td>32 [32;34]</td><td>0.146</td></tr><tr><td>Postpartum, months</td><td>7 [6;8]</td><td>6 [6;7]</td><td><strong>0.033</strong></td></tr></tbody></table>
years, *p* = 0.049) than healthy participants. Aspirin was the only drug prescribed to 4 healthy pregnant women. Other pharmacological classes were prescribed exclusively to the participants in the HDO group. Regarding vascular reverse remodeling, a significant increase in PWV and, consequently, mean arterial pressure was noticed from 3rd trimester to 6 months after delivery in both study groups (Table 2). Women in the HDO group had higher blood pressure and arterial stiffness in both evaluation moments (Table 2). Three pregnancies developed pre-eclampsia (two hypertensive and one without any CV risk factor), and one from healthy group developed cholestasis of pregnancy. Two experienced other peripartum complications (arrhythmia was documented by one participant of the healthy group, and uncontrolled arterial hypertension was reported in hypertensive woman).

## Characterization of cardiac reverse remodeling

Both groups similarly decreased heart rate and cardiac output from the 3rd trimester to 6 months after delivery (Table 2). Pregnancy impaired diastolic function in both groups, as assessed by a significant decrease of E/A and an increase in E/e' and left atrial volume (LAV) compared to the postpartum period (Table 2). At postpartum, participants showed a significant reduction in LVM, volume and relative wall thickness, irrespective of the presence of CV risk factors (Table 2).

Interestingly, global longitudinal strain was the only parameter significantly reduced in the 3rd trimester of HDO pregnant women (−23.3[−26.2; −21.3]% versus −21.0[−22.1; −20.1]%,*p* = 0.019). After delivery, the cardiac anatomic and functional normalization was much less noticeable in the HDO group, which maintained values significantly higher than those in the healthy group (Table 2). The diastolic changes observed in the HDO group become significantly different from those in the healthy group at 6 months postpartum (E/e', 6.1[5.4;7.4] versus 5.6[4.9;6.0], *p* = 0.008; LAV, 43[37;51]mL versus 38[31;43]mL, *p* = 0.045).

## Identification of proteins in urine sample

Over 2700 proteins were identified in urine samples through LC-MS/MS, after excluding reverse sequences, potential contaminants and proteins only identified by site. To explore the predictive power of the urinary proteins, we only considered proteins with two or more peptides identified and proteins with less than 20% of missing values, resulting in a selection of 324 proteins (Supplemental Table 2; Supplemental file 1).

## Urinary proteins showing an association with cardiac reverse remodeling

The LVM was defined as the outcome in a multivariable generalized linear mixed-effects model (GLMM) to identify urinary proteins associated with postpartum cardiac RR. In the pre-specified GLMM containing only clinical variables, the presence of at least one CV risk factor (β: 1.4, 95% CI [3.30; 24.0], *p* = 0.011) was associated with low LVM regression after delivery (Table 3). Upon addition of each one of 324 proteins on TOP to the pre-specified model, seventeen urinary proteins were associated with LVM (Fig. 1A and B). According to STRING network configuration, many of these proteins interact amongst them, participating in common pathways. This includes the regulation of insulin-like growth factor (IGF) transport and uptake by IGF binding proteins (involving Complement C3 [C3], Fibronectin [FN1], Serotransferrin [TF] and Trans-Golgi network integral membrane protein 2 [TGOLN2]); platelet activation, signaling and aggregation (including Inter-alpha-trypsin inhibitor heavy chain H4 [ITIH4], Serotransferrin [TF], Fibronectin [FN1] and Cell division control protein 42 homolog [CDC42]) and immune system (namely, Complement C3 [C3], Fibronectin [FN1], Cell division control protein 42 homolog [CDC42], ADP-ribosylation factor 1 [ARF1], Alpha-galactosidase A [GLA], Ubiquitin-like modifier-activating enzyme 1 [UBA 1], Pro-cathepsin H [CTSH] and Glutathione S-transferase P [GSTP1]) (Fig. 1C).

We then evaluated the added value of the proteins to the clinical model to predict LVM regression. For this, we only considered relevant urinary proteins with an absolute magnitude (expressed by the beta coefficient) higher than 6 (Inter-alpha-trypsin inhibitor heavy chain H4 [ITIH4], Glutathione S-transferase P [GSTP1], Fibronectin [FN1], ADP-ribosylation factor 1 [ARF1], Alpha-galactosidase A [GLA], Sialate O-acylesterase [SIAE], Sortilin [SORT1], Desmoglein-2 [DSG2]). Nevertheless, the model performance improved significantly (*p* < 0.001, Table 4). Of note, the combination of the 8 proteins, or the single addition of the proteins GSTP1, SORT1, or DSG2 increased model performance regarding LVM regression (lowest AIC [Akaike Information Criterion] and higher R²).

## Insulin-like growth factor-1 showing an association with cardiac reverse remodeling

Insulin-like growth factor-1 (IGF-1) is involved in the pathway of the regulation of IGF transport and uptake by IGF-binding proteins. In our cohort, increased urinary levels of IGF-1, quantified through ELISA Assay Kit were associated with low LVM regression after delivery (Table 3).

## Discussion

This work is the first exploratory prospective cohort study that included healthy pregnant women and with CV risk factors to identify urinary proteins associated with cardiac RR, as assessed through LVM variation, from the 3rd trimester to 6 months after delivery.

Our cohort showed a significant regression of LVM and volume, and LAV, associated with an improvement of diastolic function (E/e') from the 3rd trimester to 6 months after delivery, confirming previous reports²,³. Although pregnant women with and without CV risk factors have similar LVM in the 3rd trimester, pregnant women with CV risk factors showed lower LVM regression, that is, incomplete cardiac RR at 6 months postpartum compared to healthier ones. The same holds true for diastolic function in the CVR factors group. Indeed, pregnancy seems to mask the diastolic changes observed in the HDO group, as these differences only become evident and significantly different from those in the healthy group at 6 months postpartum. The presence of at least one CV risk factor (arterial hypertension, obesity, gestational diabetes or type-II diabetes) was associated
<table><caption>**Table 2.** Physical, vascular and echocardiographic assessment. Values expressed by median [interquartile range]. Significant values are in bold.</caption><thead><tr><th rowspan="2">Group</th><th colspan="3">Healthy group</th><th colspan="3">HDO group</th><th colspan="3">Healthy group versus HDO group</th><th colspan="3">Healthy group</th></tr><tr><th>Parameter</th><th>3rd trimester</th><th>6 months postpartum</th><th>Variation</th><th>3rd trimester</th><th>6 months Postpartum</th><th>Variation</th><th>3<sup>rd</sup> trimester<br>p value</th><th>6 months postpartum<br>p value</th><th>Variation<br>p value</th><th>3<sup>rd</sup> trimester<br>versus 6 months<br>postpartum<br>p value</th><th>3<sup>rd</sup> trimester<br>versus 6 months<br>postpartum<br>p value</th></tr></thead><tbody><tr><td colspan="13"><strong>Physical assessment</strong></td></tr><tr><td>Maternal body weight, kg</td><td>70 [67-81]</td><td>61 [57-69]</td><td>-10 [-11; -8]</td><td>82 [73-95]</td><td>76 [64-89]</td><td>-8 [-10; -3.8]</td><td>0.004</td><td>0.003</td><td>0.043</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Weight variation, kg</td><td>11 [8;14]</td><td>-10 [-11; -8]</td><td>-20 [-25; -17]</td><td>9 [7;12.5]</td><td>-8 [-9.5; -3.8]</td><td>-15 [-21; -13]</td><td>0.089</td><td>0.043</td><td>0.036</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td colspan="13"><strong>Vascular assessment</strong></td></tr><tr><td>Systolic blood pressure, mmHg</td><td>110 [104;120]</td><td>110 [102;118]</td><td>0 [-10; 9]</td><td>110 [110;128]</td><td>120 [114;140]</td><td>10 [0; 15]</td><td>0.160</td><td>&lt;0.001</td><td>0.022</td><td>0.753</td><td>0.009</td></tr><tr><td>Diastolic blood pressure, mmHg</td><td>65 [60;70]</td><td>76 [70;80]</td><td>10 [0; 17]</td><td>80 [70;80]</td><td>82 [80-93]</td><td>10 [0; 20]</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.565</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Mean arterial pressure, mmHg</td><td>80 [76;86]</td><td>88 [84;91]</td><td>7 [-2; 13]</td><td>90 [82;97]</td><td>95 [90;108]</td><td>7 [3; 18]</td><td>0.003</td><td>&lt;0.001</td><td>0.188</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Pulse wave velocity, cm/s</td><td>60 [5;36;3]</td><td>6.7 [6.2;7.3]</td><td>0.7 [0.5; 1.3]</td><td>6.7 [5.9;7.3]</td><td>7.2 [6.7;8.0]</td><td>0.7 [0.3; 1.2]</td><td>&lt;0.001</td><td>0.021</td><td>0.451</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td colspan="13"><strong>Echocardiographic assessment</strong></td></tr><tr><td>Heart rate, beats/min</td><td>78 [70;85]</td><td>62 [58;67]</td><td>-14 [-21; -9]</td><td>82 [72;86]</td><td>66 [64;72]</td><td>-12 [-22; -6]</td><td>0.364</td><td>0.027</td><td>0.690</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>LVES, %</td><td>60 [57;62]</td><td>62 [60;65]</td><td>3.0 [0.6; 7.8]</td><td>61 [58;62]</td><td>61 [58;63]</td><td>0.4 [-20; 4.0]</td><td>0.622</td><td>0.163</td><td>0.032</td><td>0.008</td><td>0.493</td></tr><tr><td>Global longitudinal strain, %</td><td>-23.3 [-26.2; -21.3]</td><td>-23.7 [-25.4; -21.5]</td><td>-1.2 [-2.7; 3.7]</td><td>-21.0 [-22.1; -20.1]</td><td>-22.3 [-23.5; -20.4]</td><td>-1.6 [-4.4; 1.3]</td><td>0.019</td><td>0.073</td><td>0.245</td><td>0.753</td><td>0.140</td></tr><tr><td colspan="13"><strong>Mitral inflow parameters</strong></td></tr><tr><td>E, cm/s</td><td>88 [77;96]</td><td>82 [72;90]</td><td>-7 [-17; 2]</td><td>92 [78;101]</td><td>80 [72;88]</td><td>-2 [-19; 3]</td><td>0.867</td><td>0.690</td><td>0.606</td><td>0.009</td><td>0.061</td></tr><tr><td>A, cm/s</td><td>59 [52;74]</td><td>48 [42;54]</td><td>-14 [-21; -6]</td><td>66 [59;4]</td><td>52 [44;71]</td><td>-13 [-18; -7]</td><td>0.171</td><td>0.097</td><td>0.600</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>E/A</td><td>1.43 [1.26;1.70]</td><td>1.69 [1.43;1.94]</td><td>0.19 [0.01; 0.51]</td><td>1.33 [1.09;1.53]</td><td>1.47 [1.31;1.66]</td><td>0.22 [-0.05; 0.35]</td><td>0.161</td><td>0.062</td><td>0.542</td><td>&lt;0.001</td><td>0.007</td></tr><tr><td colspan="13"><strong>Pulmonary vein inflow parameters</strong></td></tr><tr><td>S, cm/s</td><td>55 [48;59]</td><td>48 [44;56]</td><td>-4 [-12; 0]</td><td>56 [50;67]</td><td>51 [45;60]</td><td>-8 [-12; 2]</td><td>0.522</td><td>0.089</td><td>0.777</td><td>0.028</td><td>0.076</td></tr><tr><td>D, cm/s</td><td>50 [47;57]</td><td>52 [41;60]</td><td>-1 [-10; 5]</td><td>52 [43;58]</td><td>49 [41;53]</td><td>-7 [-11; 3]</td><td>0.089</td><td>0.402</td><td>0.349</td><td>0.637</td><td>0.088</td></tr><tr><td>S/D</td><td>0.93 [0.83;1.21]</td><td>1.08 [0.78;1.27]</td><td>-0.11 [-0.34; 0.11]</td><td>1.17 [0.99;1.25]</td><td>1.15 [1.06;1.28]</td><td>0.04 [-0.11; 0.14]</td><td>0.332</td><td>0.152</td><td>0.325</td><td>0.122</td><td>0.557</td></tr><tr><td colspan="13"><strong>TDI mitral annulus</strong></td></tr><tr><td colspan="13"><strong>Septal</strong></td></tr><tr><td>e' sept, cm/s</td><td>10.8 [9.5;12.6]</td><td>12.2 [10.8;13.3]</td><td>1.1 [-0.1; 2.0]</td><td>10.1 [8.6;11.2]</td><td>10.8 [9.6;11.9]</td><td>1.1 [-0.2; 1.8]</td><td>0.103</td><td>0.006</td><td>0.508</td><td>0.001</td><td>0.088</td></tr><tr><td colspan="13"><strong>Lateral</strong></td></tr><tr><td>e' lat, cm/s</td><td>17.2 [15.2;19.3]</td><td>17.9 [16.0;19.8]</td><td>1.3 [-1.1; 2.6]</td><td>15.5 [12.6;18.3]</td><td>15.3 [14.4;17.1]</td><td>-0.6 [-1.1; 1.6]</td><td>0.054</td><td>0.005</td><td>0.219</td><td>0.079</td><td>0.732</td></tr><tr><td>E/e' average (E/e')</td><td>6.3 [5.6; 7.4]</td><td>5.6 [4.9;6.0]</td><td>-1.1 [-1.7; -0.1]</td><td>7.0 [5.7;8.2]</td><td>6.1 [5.4;7.4]</td><td>-0.6 [-1.4; -0.1]</td><td>0.122</td><td>0.008</td><td>0.552</td><td>&lt;0.001</td><td>0.008</td></tr><tr><td>Stroke volume, mL</td><td>62 [53;69]</td><td>60 [51;65]</td><td>-3 [-11; 3]</td><td>60 [53;65]</td><td>60 [46;68]</td><td>-2 [-12; 2]</td><td>0.486</td><td>0.847</td><td>0.895</td><td>0.131</td><td>0.166</td></tr><tr><td>Cardiac output, L/min</td><td>5.0 [4.2;5.5]</td><td>3.7 [3.2;4.1]</td><td>-1 [-2; -1]</td><td>4.9 [4.1;5.5]</td><td>3.8 [2.8;4.5]</td><td>-1 [-2; -1]</td><td>0.865</td><td>0.905</td><td>0.841</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Relative wall thickness</td><td>0.34 [0.31;0.38]</td><td>0.31 [0.29;0.33]</td><td>-0.03 [-0.05; 0.01]</td><td>0.37 [0.33;0.40]</td><td>0.34 [0.30;0.37]</td><td>-0.05 [-0.10; 0.01]</td><td>0.061</td><td>0.093</td><td>0.553</td><td>0.007</td><td>0.017</td></tr><tr><td>LV mass, g</td><td>129 [114;135]</td><td>99 [94;109]</td><td>-22 [-35; -15]</td><td>128 [109;158]</td><td>112 [101;133]</td><td>-12 [-33; 3]</td><td>0.475</td><td>0.014</td><td>0.207</td><td>&lt;0.001</td><td>0.003</td></tr><tr><td>Maximum LA volume, mL</td><td>49 [45;53]</td><td>38 [31;43]</td><td>-11 [-18; -6]</td><td>51 [47;56]</td><td>43 [37;51]</td><td>-10 [-14; -5]</td><td>0.317</td><td>0.045</td><td>0.258</td><td>&lt;0.001</td><td>0.006</td></tr><tr><td>LV end-diastolic volume, mL</td><td>94 [82;100]</td><td>78 [71;86]</td><td>-14 [-18; -7]</td><td>101 [85;108]</td><td>88 [75;100]</td><td>-7 [-20; 1]</td><td>0.209</td><td>0.029</td><td>0.128</td><td>&lt;0.001</td><td>0.006</td></tr></tbody></table>
<table><caption>**Table 3.** Multivariable generalized linear mixed models for LVM. Significant values are in bold.</caption><thead><tr><th>Parameter</th><th>β, 95% CI</th><th>p value</th></tr></thead><tbody><tr><td colspan="3"><strong>Pre-specified clinical model</strong></td></tr><tr><td>Follow-up time</td><td>-20.00 [-26.00; -15.00]</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Age</td><td>-0.79 [-2.10;0.50]</td><td>0.200</td></tr><tr><td>Cardiovascular risk group</td><td>14.00 [3.30;24.00]</td><td><strong>0.011</strong></td></tr><tr><td>Nulliparous</td><td>-5.70 [-17.00;5.10]</td><td>0.300</td></tr><tr><td>Weight gain during pregnancy</td><td>0.28 [-0.75;1.30]</td><td>0.600</td></tr><tr><td colspan="3"><strong>Pre-specified clinical model with inclusion of IGF-1 quantification</strong></td></tr><tr><td>Follow-up time</td><td>-13.00 [-22.00; -4.60]</td><td><strong>0.003</strong></td></tr><tr><td>Age</td><td>-0.82 [-2.10;0.50]</td><td>0.200</td></tr><tr><td>Cardiovascular risk group</td><td>10.00 [-0.47;21.00]</td><td>0.060</td></tr><tr><td>Nulliparous</td><td>-6.90 [-18.00, 4.10]</td><td>0.200</td></tr><tr><td>Weight gain during pregnancy</td><td>0.09 [-0.93;1.10]</td><td>0.900</td></tr><tr><td>Log-transformed insulin-like growth factor-1</td><td>6.10 [0.18;12.00]</td><td><strong>0.044</strong></td></tr></tbody></table>

with worse LVM regression. Additionally, women with these cardiometabolic comorbidities had higher postpartum cardiac volumes and LV filling pressure than those in the healthy group, as confirmed by previous studies focusing on hypertensive diseases of pregnancy¹³⁻¹⁵. In addition, both study groups showed an increase in arterial stiffness between evaluation moments, being significantly elevated in women with CV risk factors, which is in line with previous reports¹⁶⁻¹⁹.

To identify proteins associated with higher cardiac RR, we profiled the urine proteome before and after delivery in pregnant women with and without CV risk factors. We then assessed whether these proteins could independently predict postpartum LVM. Alpha-galactosidase A (GLA), sialate O-acetylesterase (SIAE), trans-golgi network integral membrane protein 2 (TGOLN2), sortilin (SORT1), desmoglein-2 (DSG2) and neuroblast differentiation-associated protein (AHNAK) were found to be potential indicators of a higher LVM regression (β < 0). In contrast, serotransferrin (TF), inter-α-trypsin inhibitor heavy chain 4 (ITIH4), complement C3 (C3), transketolase (TKT), fibronectin (FN1), ubiquitin-like modifier-activating enzyme 1 (UBA1), cell division control protein 42 homolog (CDC42), glutathione s-transferase P1 (GSTP1), sodium/potassium-transporting ATPase subunit alpha-1 (ATP1A1), ADP-ribosylation factor 1 (ARF1) and pro-cathepsin H (CTSH) were associated with a lower LVM regression after delivery (β > 0).

Many putative proteins associated with LVM were found to be functionally integrated with biological pathways, such as the regulation of insulin-like growth factor (IGF) transport and uptake by IGF-binding proteins; platelet activation, signaling and aggregation and immune system activity. These signaling pathways evidence a complex systemic interplay between organs during pregnancy, including women's response to partum-induced cardiac RR.

For instance, insulin-like growth factor-I (IGF-1) can bind to IGF-1 receptors and activate the Akt pathway²⁰, thereby promoting growth. IGF-1/PI3K/Akt pathway is involved in LV hypertrophy²¹, including pregnancy-induced cardiac hypertrophy²². Thus, IGF bioavailability regulation, through binding to IGF-binding proteins, may impact IGF-induced pro-hypertrophic cardiac activity. Interestingly, three out of four proteins associated with IGF transport and uptake (TF, C3 and FN1) were associated with worse LVM regression, which suggests that IGF signaling blockage can be a protective mechanism against cardiac hypertrophy in the pregnancy setting. Corroborating it, the IGF-1 levels quantified in urine samples of our cohort were associated with low LVM regression after delivery. Among the proteins mapped to this pathway, we reported for the first time that FN1 is associated with low LVM regression. Indeed, FN1 has been correlated with cardiac remodeling induced by volume overload and arterial hypertension²³. Fibronectin accelerates pathological hypertrophy, leading to the expression of fetal genes, such as ANP, BNP, SkA, and MHC²³,²⁴. Thereby, FN1 is a potential target for modification of hypertrophic phenotype²⁴. Intimately correlated with an increase in FN1 is the ubiquitin-proteasome pathway, responsible for the degradation of most intracellular proteins in the heart²⁵,²⁶. For instance, UBA1, a protein involved in the initial step of ubiquitylation²⁵ whose role in cardiac remodeling and RR regulation remains unclear, was associated with lower LVM regression in our cohort. In our previous study, ubiquitin correlated positively with interventricular septum and posterior wall thickness (surrogates of cardiac hypertrophy) in patients with aortic stenosis undergoing valve replacement²⁷. Moreover, the inhibition of UBA1 by PYR-41 was linked to a significant decrease in blood pressure, cardiac hypertrophy, fibrosis, oxidative stress and inflammation and improved angiotensin II-induced cardiac contractile dysfunction²⁵.

According to STRING's functional analysis, some proteins associated with postpartum LVM, namely FN1, TF, ITIH4, and CDC42, are involved in platelet activation, signaling and aggregation. Platelets are actively involved in regulating vascular tone and inflammation and participate in the immune system, showing an important role in hemostasis, thrombosis and the pathophysiology of uteroplacental diseases²⁸,²⁹. Increased platelet aggregation accompanied by a reduction of circulating platelets has been reported throughout pregnancy²⁹⁻³¹. The decrease in platelet count may result from physiologic hemodilution or from platelet aggregation. In this case, it may lead to uteroplacental diseases such as recurrent miscarriage and pre-eclampsia²⁹. The low number of preeclampsia cases present in our cohort may be explained since approximately half of the participants with CV risk factors
<figure><img src="image_4.png" /></figure>

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

<figure><img src="image_6.png" /><figcaption>**Figure 1.** (A) Volcano Plot from multivariable generalized linear mixed-effects models (adjusted to cardiovascular risk factors, age at recruitment moment, parity, weight gain until 3rd trimester, and evaluation moments) showing the distribution of the proteins according to their effect (beta coefficient) and significance level (−log₁₀(*p*-value)) regarding LVM regression. Proteins are identified with the respective gene name. Red labels are non-significant associations whereas blue labels signify significant associations. (B) Urine proteins showing a significant association with LVM, according to the GLMM. (C) STRING Protein–Protein interaction network of the 17 most relevant proteins associated with postpartum LVM, according to a GLMM. Each node represents a protein, identified with the respective gene name. The most relevant biological processes found through gene ontology enrichment analysis are depicted through node colors. Each type of protein interaction is represented by a different edge color (cyan: known interaction from curated databases, pink: experimentally determined; yellow: text mining; black: co-expression; green: predicted interaction from gene neighborhood).</figcaption></figure>

were taking acetylsalicylic acid, a platelet aggregation inhibitor, according to their risk for preeclampsia estimated by *The Fetal Medicine Foundation* algorithm. Indeed, the American College of Obstetricians and Gynecologists (ACOG) recommended low-dose aspirin prophylaxis for women with a high risk of pre-eclampsia³². Sero-transferrin (TF) is involved in iron transfer from the placenta to the fetus through the transferrin receptor in the syncytiotrophoblast membrane. Therefore, it rises towards the end of pregnancy to meet the elevated iron requirements for fetus development³³. Our findings suggest that CV risk factors mediate the positive association between TF and LVM (i.e., higher gestational TF is associated with worse cardiac mass regression), as observed in our study and supported by studies showing TF overexpression in women with CV risk factors³⁴. Regarding *platelet activation, signaling and aggregation*, ITIH4 belongs to the ITIH family and stabilizes the extracellular matrix³⁵. Therefore, high levels of ITIH proteins may have an inhibitory effect on the degradation of extracellular matrix, thereby maintaining anti-angiogenic³⁶ and hypertrophy effects. This might explain its inverse association with LVM regression. CDC42 holds a critical role in cardiomyocyte proliferation, sarcomere organization and cell–cell adhesion³⁷. Furthermore, CDC42 is a member of the Rho GTPase family and may regulate LV hypertrophy, directly controlling non-contractile actin and microtubules cytoskeleton formation through JNK
<table><caption>**Table 4.** Performance of multivariable generalized linear mixed-effects model. AIC, akaike information criterion; ARF1, ADP-ribosylation factor 1; cond, conditional; DSG2, desmoglein-2; FN1, fibronectin; GLA, alpha-galactosidase A; GSTP1, glutathione S-transferase P1; ITIH4, inter-α-trypsin inhibitor heavy chain 4; marg., marginal; R², coefficient of determination; RMSE, root mean square error; SORT1, sortilin; SIAE, sialate O-acylesterase. The pre-specified clinical variables included in multivariable generalized linear mixed-effects model (GLMM) were: CV risk factors, age at recruitment moment, parity, weight gain until 3rd trimester and evaluation moments (3rd trimester and 6 months after delivery) due to their clinical relevance. The LVM was defined as the outcome in all GLMMs.</caption><thead><tr><th></th><th>AIC</th><th>R<sup>2</sup> (cond)</th><th>R<sup>2</sup> (marg.)</th><th>RMSE</th></tr></thead><tbody><tr><td>Pre-specified clinical variables multivariable GLMM</td><td>1032.8</td><td>0.63</td><td>0.24</td><td>12.087</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + GSTP1</td><td>861.9</td><td>0.70</td><td>0.25</td><td>10.015</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + ARF1</td><td>864.6</td><td>0.68</td><td>0.25</td><td>10.622</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + FN1</td><td>1024.3</td><td>0.67</td><td>0.28</td><td>11.368</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + SORT1</td><td>870.1</td><td><strong>0.72</strong></td><td>0.26</td><td>9.829</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + SIAE</td><td>889.8</td><td>0.61</td><td>0.29</td><td>12.535</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + GLA</td><td>959.2</td><td>0.62</td><td>0.28</td><td>12.226</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + ITIH4</td><td>1027.8</td><td>0.61</td><td>0.27</td><td>12.309</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + DSG2</td><td>837.2</td><td><strong>0.70</strong></td><td>0.25</td><td>10.482</td></tr><tr><td>Pre-specified clinical variables multivariable GLMM + panel of last 8 urinary proteins</td><td>621.0</td><td><strong>0.70</strong></td><td>0.28</td><td>9.638</td></tr></tbody></table>

anti-hypertrophic pathway³⁸. However, CDC42 shows controversial results, having been described as either pro- or anti-cardiac hypertrophic signaling actions in response to physiological and pathological conditions, depending on the experimental setting³⁹,⁴⁰. In our cohort, which included diabetic and hypertensive women, CDC42 demonstrated a positive correlation with LVM. This observation is consistent with the hormone-induced proliferative state of pregnancy. Of note, if and how the regulation of platelet-mediated hemostasis through FN1/TF/ITIH4/CDC42 impacts cardiovascular remodeling during pregnancy remains elusive and deserves further scrutiny.

According to STRING, the immune system is also an important player among the proteins deemed to be more strongly associated with cardiac RR. All phases of pregnancy, from implantation to fetal growth and parturition, require a fine-tuning of immunological processes⁴¹, consistent with pregnancy's cardiovascular adaptations. We found the complement system protein C3 to be positively correlated with LVM. Indeed, in coronary heart disease, C3 inhibition seems to restrain hypertrophy and improve cardiac function and survival⁴². In contrast, Shahini et al. showed that C3 inhibition did not affect the cardiac maladaptive remodeling in response to pressure overload⁴³. Among the 8 proteins mapped to the immune system, GSTP1, an anti-inflammatory protein, was upregulated in the hypertrophied LV to balance elevated levels of reactive oxygen species in heart failure patients, being an independent predictor of adverse cardiac remodeling⁴⁴,⁴⁵.

The pre-specified clinical model to assess the probability of postpartum RR found that 8 proteins (ITIH4, GSTP1, FN1, ARF1, GLA, SIAE, SORT1, DSG2) could improve the model's performance individually (GSTP1, DSG2 and SORT-1) or combined. Therefore, we developed the first model, including clinical variables and urinary proteins, to assess cardiac RR, through LVM quantification during and after pregnancy. Regarding the proteins relevant to improving the robustness of the GLM model, we found an inverse association between GLA, DSG2, SIAE and SORT1 expression and LVM. Curiously, despite most of these not being described previously in the context of cardiac remodeling or RR⁴⁶⁻⁴⁹, they do predict postpartum cardiac RR.

This highlights the relevance of the present study and the need to validate and unravel the involved molecular mechanisms and pathways. In this context, we desire to validate the eight proteins (ITIH4, GSTP1, FN1, ARF1, GLA, SIAE, SORT1, DSG2) in an independent cohort of pregnant women. Beyond the functional validation through IGF-1 in the present work, we want to explore the mechanisms of the most relevant biological pathways associated with incomplete cardiac RR using an already established model of pregnancy in Wistar rats⁵⁰.

Pregnancy becomes an important window to predict and prevent future cardiovascular events and heart failure progression in women with cardiometabolic comorbidities⁵¹. The incomplete cardiac RR pattern that we showed in the present study has been reported in women with a history of pregnancy hypertensive diseases, who show a higher probability of developing concentric left ventricular (LV) remodeling accompanied by LV diastolic dysfunction, substrates for Heart Failure with Preserved Ejection Fraction (HFpEF) development later in life⁵¹. For this reason, we believe identifying potential urinary proteins (ITIH4, GSTP1, FN1, ARF1, GLA, SIAE, SORT1, DSG2) as biomarkers of incomplete postpartum RR will be promising for signaling women with a higher risk of developing HFpEF at the long-term. Identifying these women during pregnancy may promote a close medical follow-up after delivery and better pharmacological therapy optimisation, reducing HFpEF incidence. Regarding study limitations, COVID-19 pandemic strongly impacted the progression of this study by limiting participants' recruitment and cardiovascular evaluation. This resulted in follow-up losses and, consequently, a reduced sample size, compromising the proteomics analysis according to the study group. This also precluded the comparison of cardiac RR and its association with a urinary proteomic profile according to each cardiovascular risk (chronic arterial hypertension versus obesity versus type 2 diabetes versus gestational hypertension versus gestational diabetes). The lack of adjustment for multiple testing performed in this study may also be considered a limitation. However, we used stringent effects size criteria to select the putative proteins to inform the extent
of cardiac RR in the model. Furthermore, the participants were stratified at the recruitment moment according to their baseline conditions (presence or absence of cardiovascular risk factors, namely chronic hypertension, gestational hypertension, gestational diabetes and obesity), not taking into account potential subsequent crossovers (such as the diagnosis of pre-eclampsia at the end of pregnancy).

In conclusion, we found that cardiac RR was characterized by significant LV chamber regression and diastolic function improvement. However, pregnant women with CV risk factors revealed compromised LVM regression and diastolic function improvement at postpartum. The proteins GSTP1, ARF1, FN1, SORT1, SIAE, GLA, ITIH4 and DGS2 were independently associated with LVM regression, and all improved the performance of a model for predicting cardiac RR. Indeed, the proteomic profile of urine reflected the cardiac reverse remodeling process, endorsing urine protein analysis as an alternative or complementary approach to assess the risk of postpartum incomplete cardiac RR. In addition, the association between IGF and postpartum LVM regression was validated.

## Data availability

The mass spectrometry proteomics data have been deposited to the PRIDE Archive (http://www.ebi.ac.uk/pride/archive/) via the PRIDE partner repository with the data set identifier PXD042655.

Received: 26 January 2024; Accepted: 21 June 2024

Published online: 27 June 2024

## References

1. Sengupta, S. P. et al. Gestational changes in left ventricular myocardial contractile function: new insights from two-dimensional speckle tracking echocardiography. *Int. J. Cardiovasc. Imaging* **33**(1), 69–82 (2017).

2. Meah, V. L. et al. Cardiac output and related haemodynamics during pregnancy: A series of meta-analyses. *Heart* **102**(7), 518–526 (2016).

3. Cong, J. et al. Structural and functional changes in maternal left ventricle during pregnancy: A three-dimensional speckle-tracking echocardiography study. *Cardiovasc. Ultrasound* **13**, 6 (2015).

4. Castleman, J. S. et al. Echocardiographic structure and function in hypertensive disorders of pregnancy: A systematic review. *Circ. Cardiovasc. Imaging* https://doi.org/10.1161/CIRCIMAGING.116.004888 (2016).

5. Buddeberg, B. S. et al. Impact of gestational diabetes mellitus on maternal cardiac adaptation to pregnancy. *Ultrasound Obstet. Gynecol.* **56**(2), 240–246 (2020).

6. Oliver-Williams, C. et al. Association between hypertensive disorders of pregnancy and later risk of cardiovascular outcomes. *BMC Med.* **20**(1), 19 (2022).

7. Ying, W., Catov, J. M. & Ouyang, P. Hypertensive disorders of pregnancy and future maternal cardiovascular risk. *J. Am. Heart Assoc.* **7**(17), e009382 (2018).

8. Decramer, S. et al. Urine in clinical proteomics. *Mol. Cell Proteom.* **7**(10), 1850–1862 (2008).

9. Zhang, Z. Y. et al. Novel urinary peptidomic classifier predicts incident heart failure. *J. Am. Heart Assoc.* https://doi.org/10.1161/JAHA.116.005432 (2017).

10. Ferreira, A. F. et al. The PERInatal MYocardial Remodeling (PERIMYR) cohort study protocol: A prospective study of cardiac remodeling and “recovery” in pregnancy as a model to understand the impact of comorbidities in cardiac remodeling and reverse remodeling. *Rev. Port Cardiol.* https://doi.org/10.1016/j.repc.2022.08.015 (2023).

11. Lang, R. M. et al. Recommendations for chamber quantification: a report from the American society of Echocardiography's guidelines and standards committee and the chamber quantification writing group, developed in conjunction with the European association of echocardiography, a branch of the European society of cardiology. *J. Am. Soc. Echocardiogr.* **18**(12), 1440–1463 (2005).

12. Nagueh, S. F. et al. Recommendations for the evaluation of left ventricular diastolic function by echocardiography: An update from the American society of echocardiography and the European association of cardiovascular imaging. *J. Am. Soc. Echocardiogr.* **29**(4), 277–314 (2016).

13. Simmons, L. A., Gillin, A. G. & Jeremy, R. W. Structural and functional changes in left ventricle during normotensive and preec- lamptic pregnancy. *Am. J. Physiol. Heart Circ. Physiol.* **283**(4), H1627–H1633 (2002).

14. Scantlebury, D. C. et al. Left ventricular hypertrophy after hypertensive pregnancy disorders. *Heart* **101**(19), 1584–1590 (2015).

15. Yu, L. et al. Left ventricular function of patients with pregnancy-induced hypertension evaluated using velocity vector imaging echocardiography and N-terminal pro-brain natriuretic peptide. *Echocardiography* **35**(4), 459–466 (2018).

16. Lekva, T. et al. Aortic stiffness and cardiovascular risk in women with previous gestational diabetes mellitus. *PLoS One* **10**(8), e0136892 (2015).

17. Castleman, J. S. et al. Altered cardiac and vascular stiffness in pregnancy after a hypertensive pregnancy. *J. Hum. Hypertens.* https://doi.org/10.1038/s41371-022-00662-4 (2022).

18. Orabona, R. et al. Endothelial dysfunction and vascular stiffness in women with previous pregnancy complicated by early or late pre-eclampsia. *Ultrasound Obstet. Gynecol.* **49**(1), 116–123 (2017).

19. Tihtonen, K. et al. Hemodynamic adaptation during pregnancy in chronic hypertension. *Hypertens. Pregnancy* **26**(3), 315–328 (2007).

20. Aguirre, G. A. et al. Insulin-like growth factor-1 deficiency and metabolic syndrome. *J. Transl. Med.* **14**(1), 3 (2016).

21. Stitt, T. N. et al. The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. *Mol. Cell* **14**(3), 395–403 (2004).

22. Aryan, L. et al. Pregnancy-associated cardiac dysfunction and the regulatory role of microRNAs. *Biol. Sex Differ.* **11**(1), 14 (2020).

23. Konstandin, M. H. et al. Fibronectin contributes to pathological cardiac hypertrophy but not physiological growth. *Basic Res. Cardiol.* **108**(5), 375 (2013).

24. Chen, H. et al. Gene expression changes associated with fibronectin-induced cardiac myocyte hypertrophy. *Physiol. Genom.* **18**(3), 273–283 (2004).

25. Shu, Q. et al. Administration of ubiquitin-activating enzyme UBA1 inhibitor PYR-41 attenuates angiotensin II-induced cardiac remodeling in mice. *Biochem. Biophys. Res. Commun.* **505**(1), 317–324 (2018).

26. Chen, C. et al. Ubiquitin-activating enzyme E1 inhibitor PYR41 attenuates angiotensin II-induced activation of dendritic cells via the IκBa/NF-κB and MKP1/ERK/STAT1 pathways. *Immunology* **142**(2), 307–319 (2014).

27. Trindade, F. et al. Preoperative myocardial expression of E3 ubiquitin ligases in aortic stenosis patients undergoing valve replacement and their association to postoperative hypertrophy. *PLoS One* **15**(9), e0237000 (2020).

28. Valero, P. et al. Platelets and endothelial dysfunction in gestational diabetes mellitus. *Acta Physiol.* **237**(4), e13940 (2023).

29. Su, X. & Zhao, W. Platelet aggregation in normal pregnancy. *Clin. Chim. Acta* **536**, 94–97 (2022).

30. Juan, P. et al. Platelets in pregnancy. *J. Prenat. Med.* **5**(4), 90–92 (2011).
31. Park, Y. H. Diagnosis and management of thrombocytopenia in pregnancy. *Blood Res.* **57**(S1), 79–85 (2022).

32. ACOG Committee Opinion No. 743: Low-dose aspirin use during pregnancy. *Obstet. Gynecol.* **132**(1), e44–e52 (2018).

33. Léger, D. et al. Physiological significance of the marked increased branching of the glycans of human serotransferrin during pregnancy. *Biochem. J.* **257**(1), 231–238 (1989).

34. Kralova, A. et al. Differential transferrin expression in placentae from normal and abnormal pregnancies: A pilot study. *Reprod. Biol. Endocrinol.* **6**(1), 27 (2008).

35. Zhao, S. et al. Serum peptidomic screening identified circulating peptide biomarkers predictive for preeclampsia. *Front. Cardiovasc. Med.* **9**, 946433 (2022).

36. Chen, H. et al. Maternal plasma proteome profiling of biomarkers and pathogenic mechanisms of early-onset and late-onset preeclampsia. *Sci. Rep.* **12**(1), 19099 (2022).

37. Li, J. et al. Essential role of Cdc42 in cardiomyocyte proliferation and cell-cell adhesion during heart development. *Dev.Biol.* **421**(2), 271–283 (2017).

38. Kehat, I. & Molkentin, J. D. Molecular pathways underlying cardiac remodeling during pathophysiological stimulation. *Circulation* **122**(25), 2727–2735 (2010).

39. Ye, H. et al. Nebivolol induces distinct changes in profibrosis MicroRNA expression compared with atenolol, in salt-sensitive hypertensive rats. *Hypertension* **61**(5), 1008–1013 (2013).

40. Raut, S. K. et al. miR-30c mediates upregulation of Cdc42 and Pak1 in diabetic cardiomyopathy. *Cardiovasc. Ther.* **33**(3), 89–97 (2015).

41. Weng, J., Couture, C. & Girard, S. Innate and adaptive immune systems in physiological and pathological pregnancy. *Biology* https://doi.org/10.3390/biology12030402 (2023).

42. Lappegård, K. T. et al. A vital role for complement in heart disease. *Mol. Immunol.* **61**(2), 126–134 (2014).

43. Shahini, N. et al. Complement component C3 and the TLR co-receptor CD14 are not involved in angiotensin II induced cardiac remodelling. *Biochem. Biophys. Res. Commun.* **523**(4), 867–873 (2020).

44. Ponikowska, B. et al. Biomarkers of myocardial injury and remodeling in heart failure. *J. Pers. Med.* https://doi.org/10.3390/jpm12050799 (2022).

45. Berezin, A. E. & Berezin, A. A. Biomarkers in heart failure: From research to clinical practice. *Ann. Lab. Med.* **43**(3), 225–236 (2023).

46. Akhtar, M. M. & Elliott, P. M. Anderson-Fabry disease in heart failure. *Biophys. Rev.* **10**(4), 1107–1119 (2018).

47. Chatterjee, D. et al. An autoantibody identifies arrhythmogenic right ventricular cardiomyopathy and participates in its pathogenesis. *Eur. Heart J.* **39**(44), 3932–3944 (2018).

48. Han, W. et al. The association between sortilin and inflammation in patients with coronary heart disease. *J. Inflamm. Res.* **13**, 71–79 (2020).

49. Tsai, S. et al. Transcriptional profiling of human placentas from pregnancies complicated by preeclampsia reveals disregulation of sialic acid acetylesterase and immune signalling pathways. *Placenta* **32**(2), 175–182 (2011).

50. Bracchi, I. et al. The cardiometabolic impact of rebaudioside a exposure during the reproductive stage. *Biology (Basel)* https://doi.org/10.3390/biology13030163 (2024).

51. Kittleson Michelle, M. et al. 2023 ACC expert consensus decision pathway on management of heart failure with preserved ejection fraction. *J. Am. Coll. Cardiol.* **81**(18), 1835–1878 (2023).

## Acknowledgements

All authors takes responsibility for all aspects of the reliability and freedom from bias of the data presented and their discussed interpretation.

## Author contributions

I.F.P., R.V., and A.S.B. conceived and designed research; A.F.F., F.T., M.J.A., J.M., A.P.M. and C.R. performed experiments; A.F.F., F.T., T.D., M.M., R.V., I.F.P. and A.S.B. analyzed data; A.F.F., F.T., F.A.S., T.D., S.O.D., C.S., M.M., I.F.P. and A.S.B. interpreted results of experiments; A.F.F., F.T. and A.S.B. prepared figures; A.F.F. drafted manuscript; F.T., M.J.A., J.M., F.S., T.D., S.O.D., C.S., A.P.M., A.L.M., B.S.M., C.R. and I.F.P. edited and revised manuscript; A.F.F., M.J.A., J.M., F.T., S.O.D., C.S., A.L.M., B.S.M., C.R., A.S.B. and I.F.P. approved final version of manuscript.

## Funding

This work has been supported by EPIC-XS, project number 318, funded by the Horizon 2020 program of the European Union; Bolsa de Estudo João Porto and Prémio Saúde Cardiovascular na Mulher da Sociedade Portuguesa de Cardiologia; by RTP Maratona da Saúde 2017; and by national funds through FCT—Portuguese Foundation for Science and Technology, under the scope of the Cardiovascular R&D Center – UnIC (UIDB/00051/2020 and UIDP/00051/2020). European Regional Development Fund (ERDF), through the North Regional Operational Program in the framework of the project HEALTH-UNORTE: Setting-up biobanks and regenerative medicine strategies to boost research in cardiovascular, musculoskeletal, neurological, oncological, immunological and infectious diseases (reference NORTE-01-0145-FEDER-000039). Ana F. Ferreira, Maria J. Azevedo and Juliana Morais are supported by Foundation for Science and Technology (SFRH/BD/138925/2018, SFRH/BD/144982/2019 and UIBD1523062021, respectively). Fábio Trindade is supported by a post-doctoral research grant by FCT through UnIC (UIDP/00051/2020).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-65612-1.

**Correspondence** and requests for materials should be addressed to A.S.B.

**Reprints and permissions information** is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

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

© The Author(s) 2024