# Proteomic Analysis of Maternal Urine for the Early Detection of Preeclampsia and Fetal Growth Restriction

Emmanuel Bujold, Alexandre Fillion, Florence Roux-Dalvai, Marie Pier

Scott-Boyer, Yves Giguère, Jean-Claude Forest, Clarisse Gotti, Geneviève

Laforest, Paul Guerby, Arnaud Droit

► To cite this version:

Emmanuel Bujold, Alexandre Fillion, Florence Roux-Dalvai, Marie Pier Scott-Boyer, Yves Giguère, et al. Proteomic Analysis of Maternal Urine for the Early Detection of Preeclampsia and Fetal Growth Restriction. Journal of Clinical Medicine, 2021, 10 (20), pp.4679. 10.3390/jcm10204679. hal-03771687

HAL Id: hal-03771687

https://ut3-toulouseinp.hal.science/hal-03771687v1

Submitted on 7 Sep 2022

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_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Article

# Proteomic Analysis of Maternal Urine for the Early Detection of Preeclampsia and Fetal Growth Restriction

Emmanuel Bujold<sup>1,2,*</sup>,<sup>†</sup> Alexandre Fillion<sup>1</sup>, Florence Roux-Dalvai<sup>1,‡</sup>, Marie Pier Scott-Boyer<sup>1</sup>, Yves Giguère<sup>1,3,‡</sup>, Jean-Claude Forest<sup>1,3</sup>, Clarisse Gotti<sup>1,‡</sup>, Geneviève Laforest<sup>1,‡</sup>, Paul Guerby<sup>1,4</sup> and Arnaud Droit<sup>1,‡</sup>

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

check for
updates

**Citation:** Bujold, E.; Fillion, A.; Roux-Dalvai, F.; Scott-Boyer, M.P.; Giguère, Y.; Forest, J.-C.; Gotti, C.; Laforest, G.; Guerby, P.; Droit, A. Proteomic Analysis of Maternal Urine for the Early Detection of Preeclampsia and Fetal Growth Restriction. *J. Clin. Med.* **2021**, *10*, 4679. https://doi.org/10.3390/jcm10204679

Academic Editors: Alex Heazell and
Sylvie Girard

Received: 11 September 2021

Accepted: 6 October 2021

Published: 13 October 2021

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

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

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

**Abstract:** Background: To explore the use of maternal urine proteome for the identification of preeclampsia biomarkers. Methods: Maternal urine samples from women with and without preeclampsia were used for protein discovery followed by a validation study. The targeted proteins of interest were then measured in urine samples collected at 20–24 and 30–34 weeks among nine women who developed preeclampsia, one woman with fetal growth restriction, and 20 women with uncomplicated pregnancies from a longitudinal study. Protein identification and quantification was obtained using liquid chromatography–tandem mass spectrometry (LC–MS/MS). Results: Among the 1108 urine proteins quantified in the discovery study, 21 were upregulated in preeclampsia and selected for validation. Nineteen (90%) proteins were confirmed as upregulated in preeclampsia cases. Among them, two proteins, ceruloplasmin and serpin A7, were upregulated at 20–24 weeks and 30–34 weeks of gestation (*p* < 0.05) in cases of preeclampsia, and could have served to identify 60% of women who subsequently developed preeclampsia and/or fetal growth restriction at 20–24 weeks of gestation, and 78% at 30–34 weeks, for a false-positive rate of 10%. Conclusions: Proteomic profiling of maternal urine can differentiate women with and without preeclampsia. Several proteins including ceruloplasmin and serpin A7 are upregulated in maternal urine before the diagnosis of preeclampsia and potentially fetal growth restriction.

**Keywords:** pregnancy; preeclampsia; proteomics; urine; biomarkers

## 1. Introduction

Preeclampsia (PE) affects about 2% to 5% of pregnant women in developed countries and up to 10% of pregnant women in developing countries [1–3]. PE is associated with short-term and long-term adverse outcomes in mothers and infants [4–8]. About 10% to 15% of maternal deaths and 25% of neonatal deaths are attributable to PE and its complications [8,9]. In the most severe cases, PE occurs before term and is typically associated with deep placentation disorders that can also lead to fetal growth restriction [10–12].

There is a growing body of evidence that the preterm PE, severe PE, and fetal growth restriction can be predicted in early pregnancy using a combination of biophysical, biochemical and ultrasound markers [13–15]. However, such screening tools require equipment

<footnote n="1">1. CHU de Québec—Université Laval Research Center, Université Laval, Quebec City, QC G1V 4G2, Canada; alexandre.fillion.3@ulaval.ca (A.F.); florence.roux-dalvai@crchudequebec.ulaval.ca (F.R.-D.); marie-pier.scott-boyer.1@ulaval.ca (M.P.S.-B.); Yves.Giguere@crchudequebec.ulaval.ca (Y.G.); jean-claude.forest.med@ssss.gouv.qc.ca (J.-C.F.); clarisse.gotti@crchudequebec.ulaval.ca (C.G.); genevieve.laforest@crchudequebec.ulaval.ca (G.L.); paul.guerby@gmail.com (P.G.); arnaud.droit@crchudequebec.ulaval.ca (A.D.)</footnote>

<footnote n="2">2. Department of Obstetrics and Gynecology, Faculty of Medicine, CHU de Québec—Université Laval, Quebec City, QC G1V 4G2, Canada</footnote>

<footnote n="3">3. Department of Molecular Biology, Medical Biochemistry and Pathology, Faculty of Medicine, CHU de Québec—Université Laval, Quebec City, QC G1V 4G2, Canada</footnote>

<footnote n="4">4. Department of Obstetrics and Gynecology, Paule de Viguier Hospital, CHU de Toulouse—Université de Toulouse, 31059 Toulouse, France</footnote>

<footnote>* Correspondence: emmanuel.bujold@crchudequebec.ulaval.ca</footnote>
and expertise that are not readily available throughout the world. There is evidence that
urine biomarkers could also be used for the early diagnosis of PE [16,17].

Mass-spectrometry-based proteomics analysis has gained popularity in the past decades for its ability to cover a large proportion of cellular or biological fluid proteomes [18]. Indeed, bottom-up proteomics have allowed for the identification and quantification of hundreds of proteins. This strategy relies on the identification and quantification of peptides resulting from the trypsin digestion of protein extracts by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The acquired spectra are then searched against publicly available protein databases. Protein candidates can also be accurately quantified by targeted proteomics in which up to one hundred proteins can be specifically monitored by LC-MS/MS in a single analysis [19]. The urinary proteome, which can be collected non-invasively and which contains more than 1800 protein species, is an ideal source of biomarkers for both renal/urological tract and systemic diseases [20,21].

For preeclampsia, a recent review of proteomic studies reported that at least 132 urine proteins could be used for early diagnosis [22]. More importantly, Buhimschi et al. observed that a specific urine proteomic profile combining serpin A1 and albumin could have: (1) a high accuracy in the identification of women with severe PE that required immediate delivery; and (2) could identify women at high risk to develop severe PE up to 10 to 25 weeks before its diagnosis [16]. However, these observations were limited to a small number of participants; the experiments were performed using equipment developed many years ago, and should therefore be validated in larger studies.

In the current study, we aim to assess the urine proteomic profiles in PE and to identify
biomarkers that could help in the early prediction of PE.

**2. Materials and Methods**

We performed a secondary analysis of the PEARL case-cohort study (PreEclampsia And growth Restriction: a Longitudinal study) that included both nulliparous women with PE (cases) and a cohort of nulliparous women at low risk of PE, recruited in early pregnancy and followed until delivery (controls) [23,24]. PE was defined according to the Society of Obstetricians and Gynecologists of Canada guidelines as de novo hypertension after 20 weeks of pregnancy (i.e., Systolic BP ≥ 140 mmHg and/or diastolic BP ≥ 90 mmHg) with proteinuria (≥300 mg per 24 h or ≥2 + protein on urine dipstick) or an adverse condition (thrombocytopenia, renal failure, liver injury, headache, seizures, vision problem or pulmonary edema) [25]. Fetal growth restriction was defined as a birth weight below the 10th percentile for gestational age based on a Canadian reference chart [26]. Urine samples were collected on the day of the diagnosis of PE for all cases and at several points during pregnancy for all controls. Urine samples were aliquoted and stored at −80 °C until analysis. All participants signed an informed consent form, and the study was approved by the ethics committee of the CHU de Québec—Université Laval (B14-07-2037).

For the purpose of the current study, we divided our population as follows: for the discovery study, 12 participants with PE (including 6 with fetal growth restriction) were matched to 12 controls without PE based on gestational age at urine collection, maternal age and maternal body mass index. For the validation study, another subset of 12 participants with PE (including 2 with fetal growth restriction) were matched to 12 controls also based on gestational age, maternal age and maternal body mass index. We used all samples available (n = 24) and divided by two for the discovery and validation study. For the longitudinal study, we used an independent cohort of women who were seen at each trimester of pregnancy: we selected 9 women who developed PE; 1 woman who developed fetal growth restriction without PE; and 20 women with uncomplicated pregnancies who provided urine samples at 20–24 and 30–34 weeks.

2.1. *Mass Spectrometry Analyses*

Each 500 µL sample of urine was concentrated on an Amicon Ultra-0.5 Centrifugal Filter device (Millipore, Burlington, MA, USA) by 10 min centrifugation 14,000 × g, fol-
lowed by a wash with 500 µL ammonium bicarbonate 50 mM and centrifugation in the same conditions. The protein concentration of the concentrated urine samples (volumes between 55 and 85 µL) was measured using Bradford assay. A total of 10 µg from each sample was then used for the subsequent steps. Briefly, the sample volume was adjusted to 50 µL with ammonium bicarbonate at 50 mM, and sodium deoxycholate was added to a final concentration of 1%. Protein denaturation was performed by heating at 95 °C for 5 min. The reduction and alkylation of cystein disulfide bridges was performed by the addition of 1,4 dithiothreitol (DTT) (final concentration 0.2 mM) and incubation at 37 °C for 30 min followed by the addition of iodoacetamide (final concentration 0.8 mM) and incubation at 37 °C for 30 min in the dark. An amount of 200 ng of tryspin enzyme was then added and samples were incubated at 37 °C overnight for proteolysis. Enzymatic digestion was stopped by acidification with 50% formic acid and the resulting peptides were purified on StageTip according to Rappsilber et al., using C18 Empore reverse phase (CDS) [27]. The samples were vacuum-dried and stored at −20 °C prior to mass spectrometry analysis. Each sample was resuspended at 0.2 µg/µL with 2% acetonitrile, 0.05% TFA. For validation and longitudinal studies only, iRT internal standard peptides (Biognosys, Schlieren, Switzerland) were added in each sample at 1X final concentration. An amount of 1µg from each sample was analyzed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) using a U3000 RSLCnano chromatographic system (Thermo Fisher Scientific, Waltham, MA, USA) and an Orbitrap Fusion mass spectrometer (Thermo Fisher Scientific). The chromatographic separation was performed on an Acclaim PepMap 100 C18 column (75 µm internal diameter, 3 µm particules and 500 mm length) using a 5-45% solvent B 90 min gradient (A: 5% acetonitrile, 0.1% formic acid; B: 80% acetonitrile, 0.1% formic acid). The mass spectrometer was operated in Data Dependent Acquisition mode for the discovery study, and Parallel Reaction Monitoring mode for the validation and longitudinal studies, using two peptide precursor masses for each protein selected from the discovery study.

## 2.2. Bioinformatics and Statistical Treatment

For the discovery study, MaxQuant software was used to obtain protein identification by searching a Uniprot human database (Human Reference Proteome UP000005640), and quantification was obtained by the Label-Free Quantification (LFQ) method [28]. LFQ intensity values of MaxQuant were used to calculate a protein fold change (FC) between the two groups of patients, and the values were then centered by the calculation of a z-score ($z = (FC - FC average)/FC standard deviation$). Statistical analyses were then performed using R software [29]. A Limma statistical test was applied to each protein between the two groups and the corresponding *p*-values were adjusted for multiple testing using the Benjamini-Hochberg method and thus obtained *q*-values [30]. A protein was considered as significantly regulated between the two groups if it met the following criteria: $|z| > 1.96$ and $q < 0.05$. For the validation and longitudinal studies, the data analysis was performed using the Skyline software [31]. The quantification value of each targeted peptide was obtained by integration of the corresponding chromatographic peak reconstructed from the best parent/fragment transition. Two normalizations of the quantification values were then applied. The first used the median of iRT internal standard peptides intensity to correct for LC-MS/MS variability, the second used the signal of peptides from the protein biotinidase (BTD) (Uniprot accession number P43251), which displayed no variation between the two groups in our discovery analysis, to correct for Bradford protein assay. For each sample, the intensities of the two peptides of each protein were summed to obtain a protein intensity value. A modified statistical Student's test (Welch test) was performed between the two groups of patients. Statistical analyses were conducted using R software and IBM SPSS Statistics version 26.0. Interaction network analysis was performed using the STRING database [32].
## 3. Results

### 3.1. Discovery Study

The median gestational age of the 12 participants at diagnosis of PE (31.9; IQ: 28.5–34.3 weeks) was comparable to the median gestational age of the 12 participants used as controls (32.1; IQ: 29.0–35.3 weeks). Table 1 reports participant’s characteristics for the discovery, the validation and the longitudinal study.

<table><caption>**Table 1.** Characteristics of our study populations.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Discovery Study</th><th colspan="2">Validation Study</th><th colspan="2">Longitudinal Study</th></tr><tr><th>Cases (n = 12)</th><th>Controls (n = 12)</th><th>Cases (n = 12)</th><th>Controls (n = 12)</th><th>Cases (n = 10)</th><th>Controls (n = 20)</th></tr></thead><tbody><tr><td>Maternal age (years)</td><td>29 (27–31)</td><td>29 (27–33)</td><td>29 (27–31)</td><td>30 (28–33)</td><td>31 (27–35)</td><td>30 (29–33)</td></tr><tr><td>BMI (kg/m<sup>2</sup>)</td><td>33 (28–36)</td><td>28 (26–33)</td><td>30 (28–35)</td><td>30 (25–34)</td><td>31 (29–33)</td><td>24 (22–26)</td></tr><tr><td>Caucasian</td><td>12 (100%)</td><td>12 (100%)</td><td>11 (92%)</td><td>12 (100%)</td><td>9 (90%)</td><td>20 (100%)</td></tr><tr><td>Gestational age at birth</td><td>34 (29–36)</td><td>39 (38–40)</td><td>34 (30–35)</td><td>40 (40–41)</td><td>35 (33–36)</td><td>40 (39–41)</td></tr><tr><td>Preeclampsia</td><td>12 (100%)</td><td>0 (0%)</td><td>12 (100%)</td><td>0 (0%)</td><td>9 (90%)</td><td>0 (0%)</td></tr><tr><td>Fetal growth restriction</td><td>6 (50%)</td><td>1 (8%)</td><td>2 (17%)</td><td>2 (17%)</td><td>4 (40%)</td><td>0 (0%)</td></tr></tbody></table>

Median (Interquartile range) or Number (percentage).

Our proteomic analyses allowed us to quantify 1108 proteins from urine samples which were used in a Principal Component Analysis (PCA) to explore the variability of the urinary proteomic profiles of participants (Figure 1). The PCA showed two distinct groups for control and PE revealing distinct proteomic profiles.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Principal Component Analysis (PCA) of the urinary proteomic profiles in participants with and without preeclampsia. Each point represents a participant. The axes correspond to artificial axes rebuilt by the PCA to display the maximum variability between the samples. CTL: Control; PE: Preeclampsia.</figcaption></figure>

A heatmap of the intensity of the proteins quantified in this study was also generated and associated to a hierarchical clustering, which allowed us to group the most similar profiles among the urine proteomes (Figure 2).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Heatmap of the expression of the 1108 proteins extracted from the proteomic analysis. Colors correspond to the centered intensity values of each protein in each urine sample from blue (smallest intensity) to red (highest intensity). The proteins (in lines) and the samples (in columns) are presented as clustered according to a complete linkage clustering method; the corresponding dendro-gram is shown only for samples.</figcaption></figure>

Among all the quantified proteins, 62 (53 upregulated, 9 downregulated) were found statistically regulated between PE and control after filtering on z-score (obtained after centering of the protein fold change) and on the *q*-value associated to a Limma statistical test corrected for multiple testing (Figure 3). Twenty-one upregulated proteins were selected for further validation (A1BG; ALB; AFM; TTR; AZGP1; C3; CA1; CP; GC; HBA1; HBB; ITIH2; ORM1; ORM2; SERPINA1; SERPINA3; SERPINA6; SERPINA7; SERPINC1; SHBG; TF; see Table 2 for protein descriptions) based on their fold change, their total signal intensity and/or the sequence of their corresponding peptides (peptides carrying
oxidized methionine and those containing missed trypsin cleavages were excluded). Using a STRING analysis, we found that most of them have known interactions, either direct protein interaction, gene co-expression, or are cited together in the literature (Supplementary Figure S1).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Volcano plot of the 1108 proteins quantified in the proteomic analysis. The x-axis corresponds to the fold change (Log2 of PE/control protein intensity ratio), the y-axis corresponds to the statistical value (−Log10 (*p*-value)). Significantly regulated proteins are displayed in red (upregulated) or in blue (downregulated); non-significantly regulated proteins are displayed in grey.</figcaption></figure>

## 3.2. *Validation Study*

The median gestational age at urine collection of our PE cases and term uncomplicated delivery controls was 31.5 (IQ: 29.0–34.3) weeks and 32.5 weeks (IQ: 29.6–33.9) weeks, respectively (Table 1). The 21 proteins previously selected were monitored by targeted proteomics to obtain an accurate protein quantification. Table 2 and Figure 4 report the average and distribution of intensity values for each group, PE or control, of each of the selected 21 proteins monitored by targeted proteomics. We observed that 19 (90%) of the 21 proteins were significantly upregulated in the urine of PE patients when compared to controls with a *p*-value < 0.05. Among them, 13 have a *p*-value < 0.001.
<table><caption>**Table 2.** Validation study by targeted proteomics of 21 urinary proteins in women with and without preeclampsia.</caption><thead><tr><th>Protein Accession</th><th>Gene Name</th><th>Protein Description</th><th>Ratio</th><th>p-Value</th><th>Significance</th></tr></thead><tbody><tr><td>P04217</td><td>A1BG</td><td>Alpha-1B-glycoprotein</td><td>4.85</td><td>0.00000</td><td>***</td></tr><tr><td>P43652</td><td>AFM</td><td>Afamin</td><td>6.89</td><td>0.00001</td><td>***</td></tr><tr><td>P25311</td><td>AZGP1</td><td>Zinc-alpha-2-glycoprotein</td><td>14.39</td><td>0.00005</td><td>***</td></tr><tr><td>P01024</td><td>C3</td><td>Complement C3</td><td>1.92</td><td>0.01029</td><td>*</td></tr><tr><td>P00915</td><td>CA1</td><td>Carbonic anhydrase 1</td><td>6.20</td><td>0.03347</td><td>*</td></tr><tr><td>P00450</td><td>CP</td><td>Ceruloplasmin</td><td>6.55</td><td>0.00005</td><td>***</td></tr><tr><td>P02774</td><td>GC</td><td>Vitamin D-binding protein</td><td>11.31</td><td>0.00001</td><td>***</td></tr><tr><td>P19823</td><td>ITIH2</td><td>Inter-alpha-trypsin inhibitor heavy chain 2</td><td>6.35</td><td>0.01660</td><td>*</td></tr><tr><td>P02763</td><td>ORM1</td><td>Alpha-1-acid-glycoprotein 1</td><td>2.96</td><td>0.06649</td><td></td></tr><tr><td>P19652</td><td>ORM2</td><td>Alpha-1-acid-glycoprotein 2</td><td>1.66</td><td>0.28634</td><td></td></tr><tr><td>P01009</td><td>SERPINA1</td><td>Alpha-1-antitrypsin (serpin A1)</td><td>9.50</td><td>0.01026</td><td>*</td></tr><tr><td>P01011</td><td>SERPINA3</td><td>Alpha-1-antichymotrypsin (serpin A3)</td><td>3.83</td><td>0.00011</td><td>**</td></tr><tr><td>P08185</td><td>SERPINA6</td><td>Corticosteroid-binding globulin (serpin A6)</td><td>2.27</td><td>0.00039</td><td>**</td></tr><tr><td>P05543</td><td>SERPINA7</td><td>Thyroxine-binding globulin (serpin A7)</td><td>16.42</td><td>0.00013</td><td>**</td></tr><tr><td>P01008</td><td>SERPINC1</td><td>Antithrombin-III (serpin C1)</td><td>4.38</td><td>0.00000</td><td>***</td></tr><tr><td>I3L145</td><td>SHBG</td><td>Sex-hormone-binding globulin</td><td>4.35</td><td>0.00000</td><td>***</td></tr><tr><td>P02766</td><td>TTR</td><td>Transthyretin</td><td>4.32</td><td>0.00000</td><td>***</td></tr><tr><td>P02768</td><td>ALB</td><td>Albumin</td><td>2.57</td><td>0.00001</td><td>***</td></tr><tr><td>P69905</td><td>HBA1</td><td>Hemoglobin subunit alpha 1</td><td>3.74</td><td>0.02510</td><td>*</td></tr><tr><td>P68871</td><td>HBB</td><td>Hemoglobin subunit beta</td><td>24.62</td><td>0.00004</td><td>***</td></tr><tr><td>P02787</td><td>TF</td><td>Serotransferrin</td><td>2.82</td><td>0.02372</td><td>*</td></tr></tbody></table>

For each protein (Uniprot accession number given in the first column), average intensity ratio of PE over control groups (PE/CTL) is shown associated to its Student's test statistical *p*-value and the corresponding significance. * *p* < 0.05; ** *p* < 0.001; *** *p* < 0.0001.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Validation study: Distribution of intensities measured by targeted proteomics of 21 urinary proteins in women with and without preeclampsia. The boxplots represent the distribution of intensities in each group (grey: control, yellow: PE) as well as the median (line), the interquartile (box) and the maximum and minimum values within 1.5 times the interquartile range (whisker) of 21 proteins monitored by targeted proteomics in urine samples of women with preeclampsia or women with uncomplicated pregnancies. Statistical significance is displayed above each graph: * *p* < 0.05; ** *p* < 0.01; *** *p* < 0.001.</figcaption></figure>
### 3.3. Longitudinal Study

In the longitudinal study, we observed that 6 of the 21 proteins monitored by targeted proteomics were significantly upregulated at 30–34 weeks in women who subsequently developed PE in comparison to the control group (serpin A7; ceruloplasmin; afamin; inter-alpha-trypsin inhibitor heavy chain; transferrin; alpha-1B-glycoprotein). Two, serpin A7 and ceruloplasmin, were also found significantly upregulated at 20–24 weeks in women who subsequently developed PE in comparison to the control group (Table 3 and Figure 5).

<table><caption>**Table 3.** Longitudinal study by targeted proteomics of 21 urinary proteins in women with and without preeclampsia.</caption><thead><tr><th rowspan="2">Gene Name</th><th colspan="3">20-24 Weeks</th><th colspan="3">30-34 Weeks</th></tr><tr><th>Ratio PE/CTL</th><th>p-Value</th><th>Significance</th><th>Ratio PE/CTL</th><th>p-Value</th><th>Significance</th></tr></thead><tbody><tr><td>SERPINA7</td><td>1.42</td><td>0.04247</td><td>*</td><td>1.62</td><td>0.00966</td><td>**</td></tr><tr><td>CP</td><td>1.81</td><td>0.04515</td><td>*</td><td>4.87</td><td>0.01784</td><td>*</td></tr><tr><td>AFM</td><td>1.29</td><td>0.33745</td><td></td><td>2.48</td><td>0.01821</td><td>*</td></tr><tr><td>ITIH2</td><td>1.14</td><td>0.66046</td><td></td><td>3.67</td><td>0.02808</td><td>*</td></tr><tr><td>TF</td><td>1.86</td><td>0.27263</td><td></td><td>8.42</td><td>0.03446</td><td>*</td></tr><tr><td>A1BG</td><td>1.19</td><td>0.50528</td><td></td><td>2.22</td><td>0.04651</td><td>*</td></tr><tr><td>SERPINA3</td><td>1.39</td><td>0.22978</td><td></td><td>1.95</td><td>0.05920</td><td></td></tr><tr><td>GC</td><td>1.20</td><td>0.28638</td><td></td><td>1.97</td><td>0.06296</td><td></td></tr><tr><td>ALB</td><td>1.49</td><td>0.20602</td><td></td><td>5.54</td><td>0.06954</td><td></td></tr><tr><td>SERPINA1</td><td>1.36</td><td>0.24731</td><td></td><td>4.11</td><td>0.07165</td><td></td></tr><tr><td>C3</td><td>1.07</td><td>0.85042</td><td></td><td>1.63</td><td>0.07776</td><td></td></tr><tr><td>ORM1</td><td>1.69</td><td>0.42028</td><td></td><td>2.52</td><td>0.07781</td><td></td></tr><tr><td>SERPINA6</td><td>1.05</td><td>0.84649</td><td></td><td>1.97</td><td>0.08995</td><td></td></tr><tr><td>ORM2</td><td>1.35</td><td>0.49657</td><td></td><td>1.52</td><td>0.19926</td><td></td></tr><tr><td>TTR</td><td>1.23</td><td>0.29396</td><td></td><td>1.48</td><td>0.21370</td><td></td></tr><tr><td>HBA1</td><td>5.27</td><td>0.01524</td><td>*</td><td>55.40</td><td>0.30338</td><td></td></tr><tr><td>HBB</td><td>5.43</td><td>0.02956</td><td>*</td><td>30.91</td><td>0.33029</td><td></td></tr><tr><td>SERPINC1</td><td>0.99</td><td>0.97216</td><td></td><td>1.28</td><td>0.34586</td><td></td></tr><tr><td>AZGP1</td><td>1.44</td><td>0.59655</td><td></td><td>1.19</td><td>0.61454</td><td></td></tr><tr><td>CA1</td><td>0.62</td><td>0.55159</td><td></td><td>0.70</td><td>0.63019</td><td></td></tr><tr><td>SHBG</td><td>0.92</td><td>0.58530</td><td></td><td>1.07</td><td>0.72452</td><td></td></tr></tbody></table>

For each protein (Uniprot accession number given in the first column) at two-time points (20–24 or 30–34 weeks of gestation), average intensity of PE over control groups (PE/CTL) is shown associated to its Student's test statistical *p*-value and the corresponding significance.
* *p* < 0.05; ** *p* < 0.01.

Figure 6A,B reports the receiver operating characteristic (ROC) curves for the prediction of PE or fetal growth restriction using ceruloplasmin and serpin A7 protein concentrations in urine at 20–24 weeks and 30–34 weeks of gestation. At a false positive rate of 10%, ceruloplasmin could have predicted between 60% and 78% of the PE and/or fetal growth restriction cases at 20–24 weeks and 30–34 weeks of gestation, respectively. Interestingly, the case of fetal growth restriction without preeclampsia would have been detected with ceruloplasmin but not with serpin A7.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Longitudinal study: Distribution of intensities measured by targeted proteomics of 21 urinary proteins in women with and without preeclampsia at 20–24 or 30–34 weeks of gestation. The boxplots represent the distribution of intensities in each group as well as the median (line), the interquartile (box) and the maximum and minimum values within 1.5 times the interquartile range (whisker) of 21 proteins monitored by targeted proteomics in urine samples collected at 20–22 weeks (blue) and 30–32 weeks (green) of gestation in women with preeclampsia (n = 9) (dark color) or women with uncomplicated pregnancies (n = 20) (light color). Statistical significance is displayed above each graph: \* *p* < 0.05; \*\* *p* < 0.01.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Receiver operating characteristic (ROC) curves for the prediction of preeclampsia and/or fetal growth restriction based on ceruloplasmin and serpin A7 proteins measured by targeted proteomics in urine of women at 20–24 and 30–34 weeks of gestation. The ROC curves present the predictive values of ceruloplasmin and serpin A7 measurements using targeted proteomic on maternal urine samples collected at 20–24 weeks of gestation (A) and at 30–34 weeks of gestation (B). The area under the ROC curves (AUC) were significant for ceruloplasmin at 20–24 weeks (AUC: 0.78; 95%CI: 0.58–0.97, *p* = 0.016); serpin A7 at 20–24 weeks (AUC: 0.75; 95%CI: 0.57–0.92, *p* = 0.028); ceruloplasmin at 30–34 weeks (AUC: 0.92; 95%CI: 0.82–1.00, *p* < 0.001); serpin A7 at 30–34 weeks (AUC: 0.77; 95%CI: 0.56–0.98, *p* = 0.024).</figcaption></figure>
## 4. Discussion

We observed that women with preeclampsia have a different urine proteomic profile than women without preeclampsia and such differences can be present up to 12 weeks before the first signs and symptoms of preeclampsia. This observation, consistent with previous studies, suggests that women at high-risk of preeclampsia could be detected from urine biomarkers. Since fetal growth restriction shares a common mechanism of disease, it is possible that it could be detected in early pregnancy as well.

More specifically, the PCA performed on our dataset discriminates the two groups of patients (PE vs. controls). Moreover, the hierarchical clustering associated with the heatmap of protein intensities generates two main clusters containing either control or PE samples. To reveal the proteins dysregulated between PE and controls, we centered the control/PE ratios by calculating a z-score. Indeed, the proteinuria associated to PE strongly affects the protein content of the analyzed samples. Therefore, lower mass spectrometry signals obtained in PE cases, as can be observed on the heatmap of protein expressions, might result in a misinterpretation of the dysregulated proteins between the two groups. By applying a z-score correction and using a statistical test corrected for multiple testing (q-values), we could confidently identify 62 proteins with statistical differences in intensity between control and PE groups. Using targeted proteomics, known for their high accuracy in protein quantification, and another subset of patients from the same cohort, we confirmed our findings from the discovery study. The monitoring of these 21 proteins in our longitudinal study revealed that 6 of these proteins were significantly upregulated at 30–34 weeks of gestation and 2 of them (ceruloplasmin and serpin A7) were significantly upregulated at 20–24 weeks of gestation, up to 12 weeks before the clinical onset of PE. Ceruloplasmin is a copper-binding protein involved in iron transport across cellular membranes, and has antioxidant ferroxidase properties. It could be upregulated in PE as a response to placental ischemia [33,34]. As for serpin A7, also known as thyroxine-binding globulin (TBG), the major thyroid hormone transport protein in serum, studies have observed lower serum TBG in women with PE [35,36]. However, little is known about how and why serpin A7 is upregulated in the urine of PE cases.

Overall, our observations suggest that the urine proteomic profile could be used to predict preeclampsia several weeks before the first signs and symptoms manifest themselves. More specifically, we observed that at least two proteins, namely, ceruloplasmin and serpin A7, are increased in maternal urine at 20–24 and 30–34 weeks of gestation in most women who subsequently develop PE.

Our study is in agreement with the study of Buhimschi et al., who observed that serpin A1, another member of the serpin family, was significantly increased in PE cases [16]. Their study also reported misfolded proteins in the urine of women with preeclampsia bound to Congo Red dye (urine congophilia or affinity for the amyloidophilic dye Congo Red) [37]. Ultimately, their study reported promising results from a prototype point-of-care test for the detection of urine congophilia [38]. A review by Navajas et al. found nine publications from 2008 to 2020 that observed 132 proteins that were differently expressed in urine in PE cases compared to controls [22]. Nineteen of these showed high potential for PE prediction as they were consistently higher or lower in PE, including ceruloplasmin, serpin A1, serpin A5, C3, ALB, TF and HBB. Starodubtseva et al. reported that the estimation of serpin A1 peptides in urine was also related to the severity of PE [39]. Placental growth factor (PlGF), a proangiogenic protein commonly used for the prediction of PE and other placental-mediated outcomes of pregnancy is also predictive of PE and PE-related adverse outcomes when measured in urine [17,23,40]. Altogether, these studies, including ours, provide hope that rapid identification of PE and potentially early prediction of PE is possible using urine studies.

The small number of cases available at each step is a limit of our study. However, we used standardized methods for proteomic analysis and, using high-resolution LC-MS/MS, we obtained the largest urinary proteome coverage ever published on preeclampsia samples (1108 proteins quantified). One of the major strengths of the current study is the validation
of our findings in an independent subset of patients and subsequent confirmation in a
prospective study, which includes the collection of urine samples up to three months before
the first signs and symptoms of PE.

**5. Conclusions**

PE is a major cause of perinatal morbidity and mortality around the world, primarily
in developing countries. PE is commonly associated with fetal growth restriction, as they
share a common mechanism of disease. Our study and current literature strongly suggest
that PE, and potentially fetal growth restriction, are syndromes that are highly detectable
in their early phases using the proteomic analysis of maternal urine. More efforts should
be devoted to the development of rapid point-of-care tests on maternal urine that could
help in the prevention of PE-related adverse outcomes of pregnancy.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/jcm10204679/s1, Figure S1: STRING interaction network of the 21 proteins selected for validation study.

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

**Funding:** This research was funded by Jeanne et Jean-Louis Lévesque, Perinatal Research Chair—Université Laval.

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Centre de recherché du CHU de Québec (protocol code B14-07-2037 on 7 August 2014).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** We acknowledge the work of M. Michael Spanier, for the editing and review of the document.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**References**

1. Mol, B.W.J.; Roberts, C.T.; Thangaratinam, S.; Magee, L.A.; de Groot, C.J.M.; Hofmeyr, G.J. Pre-eclampsia. *Lancet* **2016**, *387*, 999–1011. [CrossRef]

2. Rana, S.; Lemoine, E.; Granger, J.P.; Karumanchi, S.A. Preeclampsia: Pathophysiology, Challenges, and Perspectives. *Circ. Res.* **2019**, *124*, 1094–1112. [CrossRef] [PubMed]

3. Nakimuli, A.; Chazara, O.; Byamugisha, J.; Elliott, A.M.; Kaleebu, P.; Mirembe, F.; Moffett, A. Pregnancy, parturition and preeclampsia in women of African ancestry. *Am. J. Obstet. Gynecol.* **2014**, *210*, 510–520. [CrossRef] [PubMed]

4. Murphy, M.S.; Smith, G.N. Pre-eclampsia and Cardiovascular Disease Risk Assessment in Women. *Am. J. Perinatol.* **2016**, *33*, 723–731. [CrossRef] [PubMed]

5. Wu, P.; Haththotuwa, R.; Kwok, C.S.; Babu, A.; Kotronias, R.A.; Rushton, C.; Zaman, A.; Fryer, A.A.; Kadam, U.; Chew-Graham, C.A.; et al. Preeclampsia and Future Cardiovascular Health: A Systematic Review and Meta-Analysis. *Circ. Cardiovasc. Qual. Outcomes* **2017**, *10*, e003497. [CrossRef]

6. Pinheiro, T.V.; Brunetto, S.; Ramos, J.G.; Bernardi, J.R.; Goldani, M.Z. Hypertensive disorders during pregnancy and health outcomes in the offspring: A systematic review. *J. Dev. Orig. Health Dis.* **2016**, *7*, 391–407. [CrossRef]

7. Davis, E.F.; Lazdam, M.; Lewandowski, A.J.; Worton, S.A.; Kelly, B.; Kenworthy, Y.; Adwani, S.; Wilkinson, A.R.; McCormick, K.; Sargent, I.; et al. Cardiovascular risk factors in children and young adults born to preeclamptic pregnancies: A systematic review. *Pediatrics* **2012**, *129*, e1552–e1561. [CrossRef]
8. Duley, L. The global impact of pre-eclampsia and eclampsia. *Semin. Perinatol.* **2009**, *33*, 130–137. [CrossRef]

9. Steegers, E.A.; von Dadelszen, P.; Duvekot, J.J.; Pijnenborg, R. Pre-eclampsia. *Lancet* **2010**, *376*, 631–644. [CrossRef]

10. Burton, G.J.; Redman, C.W.; Roberts, J.M.; Moffett, A. Pre-eclampsia: Pathophysiology and clinical implications. *BMJ* **2019**, *366*, 12381. [CrossRef]

11. Phipps, E.A.; Thadhani, R.; Benzing, T.; Karumanchi, S.A. Pre-eclampsia: Pathogenesis, novel diagnostics and therapies. *Nat. Rev. Nephrol.* **2019**, *15*, 275–289. [CrossRef]

12. Hod, T.; Cerdeira, A.S.; Karumanchi, S.A. Molecular Mechanisms of Preeclampsia. *Cold Spring Harb. Perspect. Med.* **2015**, *5*, a023473. [CrossRef]

13. Boutin, A.; Gasse, C.; Guerby, P.; Giguere, Y.; Tetu, A.; Bujold, E. First-Trimester Preterm Preeclampsia Screening in Nulliparous Women: The Great Obstetrical Syndrome (GOS) Study. *J. Obstet. Gynaecol. Can.* **2021**, *43*, 43–49. [CrossRef]

14. Rolnik, D.L.; Wright, D.; Poon, L.C.Y.; Syngelaki, A.; O’Gorman, N.; de Paco Matallana, C.; Akolekar, R.; Cicero, S.; Janga, D.; Singh, M.; et al. ASPRE trial: Performance of screening for preterm pre-eclampsia. *Ultrasound Obstet. Gynecol.* **2017**, *50*, 492–495. [CrossRef]

15. Mosimann, B.; Amylidi-Mohr, S.K.; Surbek, D.; Raio, L. First Trimester Screening for Preeclampsia—A Systematic Review. *Hypertens. Pregnancy* **2020**, *39*, 1–11. [CrossRef] [PubMed]

16. Buhimschi, I.A.; Zhao, G.; Funai, E.F.; Harris, N.; Sasson, I.E.; Bernstein, I.M.; Saade, G.R.; Buhimschi, C.S. Proteomic profiling of urine identifies specific fragments of SERPINA1 and albumin as biomarkers of preeclampsia. *Am. J. Obstet. Gynecol.* **2008**, *199*, 551.e1–551.e16. [CrossRef] [PubMed]

17. Lecarpentier, E.; Gris, J.C.; Cochery-Nouvellon, E.; Mercier, E.; Abbas, H.; Thadhani, R.; Karumanchi, S.A.; Haddad, B. Urinary Placental Growth Factor for Prediction of Placental Adverse Outcomes in High-Risk Pregnancies. *Obstet. Gynecol.* **2019**, *134*, 1326–1332. [CrossRef] [PubMed]

18. Parker, C.E.; Pearson, T.W.; Anderson, N.L.; Borchers, C.H. Mass-spectrometry-based clinical proteomics—A review and prospective. *Analyst* **2010**, *135*, 1830–1838. [CrossRef]

19. Uzozie, A.C.; Aebersold, R. Advancing translational research and precision medicine with targeted proteomics. *J. Proteom.* **2018**, *189*, 1–10. [CrossRef]

20. Beasley-Green, A. Urine Proteomics in the Era of Mass Spectrometry. *Int. Neurourol. J.* **2016**, *20*, S70–S75. [CrossRef]

21. Marimuthu, A.; O’Meally, R.N.; Chaerkady, R.; Subbannayya, Y.; Nanjappa, V.; Kumar, P.; Kelkar, D.S.; Pinto, S.M.; Sharma, R.; Renuse, S.; et al. A comprehensive map of the human urinary proteome. *J. Proteome Res.* **2011**, *10*, 2734–2743. [CrossRef] [PubMed]

22. Navajas, R.; Corrales, F.; Paradela, A. Quantitative proteomics-based analyses performed on pre-eclampsia samples in the 2004-2020 period: A systematic review. *Clin. Proteom.* **2021**, *18*, 6. [CrossRef] [PubMed]

23. Fillion, A.; Guerby, P.; Lachance, C.; Comeau, M.-P.; Bussières, M.-C.; Doucet-Gingras, F.-A.; Zérounian, S.; Demers, S.; Laforest, G.; Menzies, D.; et al. Placental Growth Factor and soluble Fms-like tyrosine kinase-1 in Preeclampsia: A case-cohort (PEARL) study. *J. Obstet. Gynaecol. Can.* **2020**, *42*, 1235–1242. [CrossRef] [PubMed]

24. Fillion, A.; Guerby, P.; Menzies, D.; Lachance, C.; Comeau, M.P.; Bussieres, M.C.; Doucet-Gingras, F.A.; Zerounian, S.; Bujold, E. Pathological investigation of placentas in preeclampsia (the PEARL study). *Hypertens. Pregnancy* **2021**, *40*, 56–62. [CrossRef] [PubMed]

25. Magee, L.A.; Pels, A.; Helewa, M.; Rey, E.; von Dadelszen, P. Diagnostic, evaluation et prise en charge des troubles hypertensifs de la grossesse: Resume directif. *J. Obstet. Gynaecol. Can.* **2016**, *38*, S426–S452. [CrossRef]

26. Kramer, M.S.; Platt, R.W.; Wen, S.W.; Joseph, K.S.; Allen, A.; Abrahamowicz, M.; Blondel, B.; Breart, G.; Fetal/Infant Health Study Group of the Canadian Perinatal Surveillance, S. A new and improved population-based Canadian reference for birth weight for gestational age. *Pediatrics* **2001**, *108*, E35. [CrossRef]

27. Rappsilber, J.; Mann, M.; Ishihama, Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. *Nat. Protoc.* **2007**, *2*, 1896–1906. [CrossRef]

28. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. *Nat. Protoc.* **2016**, *11*, 2301–2319. [CrossRef]

29. RCoreTeam. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2020.

30. Smyth, G. *Limma: Linear Models for Microarray Data*; Springer: New York, NY, USA, 2005; pp. 397–420.

31. Pino, L.K.; Searle, B.C.; Bollinger, J.G.; Nunn, B.; MacLean, B.; MacCoss, M.J. The Skyline ecosystem: Informatics for quantitative mass spectrometry proteomics. *Mass Spectrom. Rev.* **2020**, *39*, 229–244. [CrossRef]

32. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. *Nucleic Acids Res.* **2021**, *49*, D605–D612. [CrossRef]

33. Serdar, Z.; Gur, E.; Develioglu, O. Serum iron and copper status and oxidative stress in severe and mild preeclampsia. *Cell Biochem. Funct.* **2006**, *24*, 209–215. [CrossRef] [PubMed]

34. Guller, S.; Buhimschi, C.S.; Ma, Y.Y.; Huang, S.T.; Yang, L.; Kuczynski, E.; Zambrano, E.; Lockwood, C.J.; Buhimschi, I.A. Placental expression of ceruloplasmin in pregnancies complicated by severe preeclampsia. *Lab. Investig.* **2008**, *88*, 1057–1067. [CrossRef] [PubMed]
35. Alemu, A.; Terefe, B.; Abebe, M.; Biadgo, B. Thyroid hormone dysfunction during pregnancy: A review. *Int. J. Reprod. Biomed.* **2016**, *14*, 677-686. [CrossRef] [PubMed]

36. Kaya, E.; Sahin, Y.; Ozkececi, Z.; Pasaoglu, H. Relation between birth weight and thyroid function in preeclampsia-eclampsia. *Gynecol. Obstet. Investig.* **1994**, *37*, 30-33. [CrossRef]

37. Buhimschi, I.A.; Nayeri, U.A.; Zhao, G.; Shook, L.L.; Pensalfini, A.; Funai, E.F.; Bernstein, I.M.; Glabe, C.G.; Buhimschi, C.S. Protein misfolding, congophilia, oligomerization, and defective amyloid processing in preeclampsia. *Sci. Transl. Med.* **2014**, *6*, 245ra292. [CrossRef]

38. Bracken, H.; Buhimschi, I.A.; Rahman, A.; Smith, P.R.S.; Pervin, J.; Rouf, S.; Bousieguez, M.; Lopez, L.G.; Buhimschi, C.S.; Easterling, T.; et al. Congo red test for identification of preeclampsia: Results of a prospective diagnostic case-control study in Bangladesh and Mexico. *EClinicalMedicine* **2021**, *31*, 100678. [CrossRef]

39. Starodubtseva, N.; Nizyaeva, N.; Baev, O.; Bugrova, A.; Gapaeva, M.; Muminova, K.; Kononikhin, A.; Frankevich, V.; Nikolaev, E.; Sukhikh, G. SERPINA1 Peptides in Urine as A Potential Marker of Preeclampsia Severity. *Int. J. Mol. Sci.* **2020**, *21*, 914. [CrossRef]

40. Levine, R.J.; Thadhani, R.; Qian, C.; Lam, C.; Lim, K.H.; Yu, K.F.; Blink, A.L.; Sachs, B.P.; Epstein, F.H.; Sibai, B.M.; et al. Urinary placental growth factor and risk of preeclampsia. *JAMA* **2005**, *293*, 77-85. [CrossRef]