Deciphering the preeclampsia-specific immune microenvironment and the role of pro-inflammatory macrophages at the maternal–fetal interface

  1. Haiyi Fei
  2. Xiaowen Lu
  3. Zhan Shi
  4. Xiu Liu
  5. Cuiyu Yang
  6. Xiaohong Zhu
  7. Yuhan Lin
  8. Ziqun Jiang
  9. Jianmin Wang
  10. Dong Huang
  11. Liu Liu
  12. Songying Zhang  Is a corresponding author
  13. Lingling Jiang  Is a corresponding author
  1. Assisted Reproduction Unit, Department of Obstetrics and Gynecology, Sir Run Run Shaw Hospital, School of Medicine, Zhejiang University, China
  2. Zhejiang Provincial Clinical Research Center for Reproductive Health and Disease, China
  3. Zhejiang Key Laboratory of Precise Protection and Promotion of Fertility, China
  4. Department of Obstetrics and Gynecology, Zhejiang Xiaoshan Hospital, China
7 figures, 6 tables and 1 additional file

Figures

Identification and characterization of placental immune cells using cytometry by time of flight (CyTOF).

(A) Schematic of the experimental workflow in CyTOF experiment. The placentas were obtained from individuals with normal pregnancy (NP, n = 9), preeclampsia (PE, n = 8), gestational diabetes mellitus (GDM, n = 8), or GDM&PE (n = 7). (B) t-Distributed Stochastic Neighbor Embedding (t-SNE) maps showing 8 × 104 CD45+ cells (the average cell number of the all samples) from the placenta overlaid with color-coded clusters and the distributions of B cells, CD4+ T cells, CD8+ T cells, dendritic cell (DC), γδT cells, monocytes, granulocytes, myeloid-derived suppressor cell (MDSC), and natural killer (NK) cells. (C) Percentages of each cell type of CD45+ cells in placentas. (D) t-SNE maps showing the expression of CD3, CD8, CD4, γδTCR, CD14, CD15, and CD56. (E) Heatmap showing the expression levels of markers in CD45+ cell subsets. Data were compared between NP and PE, NP and GDM, and NP and GDM&PE using the Kruskal–Wallis test and represented as mean ± SEM (*p < 0.05, NS, not significant). CyTOF soure data are stored in Dryad Digital Repository, DOI: 10.5061/dryad.4qrfj6qn0.

Figure 2 with 1 supplement
Specific altered T cell profile in the placentas of individuals with PE.

(A) Distribution of the CD4+ T cells in the group analyzed using t-Distributed Stochastic Neighbor Embedding (t-SNE). (B) Scatter dot plots showing the frequencies of cluster 8 of CD4+ T cells in the placentas of individuals with NP, PE, GDM, and GDM&PE (n = 9 in NP group; n = 8 in PE group; n = 8 in GDM group; n = 7 in GDM&PE group). (C) Expression of IL-17A and Foxp3 in CD45RO+CCR7+CD4+ T cells in placentas of individuals with NP and PE using flow cytometry (IL-17A: n = 10 in NP group, n = 15 in PE group; Foxp3: n = 7 in NP group, n = 9 in PE group). (D) Immunofluorescence co-staining of CD4 (red), CD45RO (green), and DAPI (blue) in frozen placental sections. The right panels show the fluorescence intensity of CD4 and CD45RO. Scale bar, 20 µm. (E) Scatter dot plots showing significantly altered markers of in cluster 8 of CD4+ T cells. (F) Distribution of the CD8+ T cells analyzed using t-SNE. (G) Scatter dot plots showing the frequencies of cluster 2 of CD8+ T cells in the placentas of individuals with NP, PE, GDM, and GDM&PE. (H) Scatter dot plots showing significantly altered markers in cluster 2 of CD8+ T cells. Data were compared between NP and PE, NP and GDM, and NP and GDM&PE using the Kruskal–Wallis test and represented as mean ± SEM (*p < 0.05, NS, not significant).

Figure 2—figure supplement 1
Identification of the placental T cell subsets.

(A) Heatmap showing the expression levels of markers in the CD4+ T subsets (n = 9 in NP group; n = 8 in PE group; n = 8 in GDM group; n = 7 in GDM&PE group). (B) Heatmap showing the expression levels of markers in the CD8+ T subsets. (C) Heatmap showing the expression levels of markers in the γδT subsets. (D) Distribution of the γδT cells in the group analyzed using t-Distributed Stochastic Neighbor Embedding (t-SNE). (E) Scatter dot plots showing the frequencies of cluster 15 of γδT cells in the placentas of individuals with NP, PE, GDM, and GDM&PE. (F) Scatter dot plots showing significantly altered markers in cluster 15 of γδT cells. Data were compared between NP and PE, NP and GDM, and NP and GDM&PE using the Kruskal–Wallis test and represented as mean ± SEM (*p < 0.05, **p < 0.01,NS, not significant).

Identification of the placental CD11b+ cell subsets and the interaction between placental immune cells.

(A) Distribution of the CD11b+ immune cells in group analyzed using t-Distributed Stochastic Neighbor Embedding (t-SNE). (B) Heatmap showing the expression levels of markers in the CD11b+ cells. (C) Scatter dot plots showing the frequencies of clusters 12, 23, and 25 of CD11b+ cells in the placentas. (D) Interaction between placental immune cells showed in heatmap. (E) Scatter plots of Pearson’s correlation analysis between placental immune cells. Data were compared between NP and PE, NP and GDM, and NP and GDM&PE using the Kruskal–Wallis test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001; NS, not significant).

Figure 4 with 2 supplements
The immune imbalance at the maternal–fetal interface induced by F480+CD206 pro-inflam Macs.

(A) An illustration showcasing the RNA-seq of the CD45+F4/80+CD206 pro-inflam Macs and CD45+F4/80+CD206+ anti-inflam Macs derived from the RUPP mouse model. (B) Principal component analysis (PCA) reflected the differences between the two groups of macrophages (n = 3). (C) The volcano map shows a comparison of the content and p value of gene expression between pro-inflam and anti-inflam Macs. Differential expression genes were screened out when p < 0.05. Red dots indicate genes with increased expression in pro-inflam Macs. Blue dots indicate genes with decreased expression. (D) The volcano map shows differential expression genes between pro-inflam and anti-inflam Macs. (E) Representative pathways enriched in the identified genes as determined by gene set enrichment analysis (GSEA) (p-value <0.05). (F) Embryo abortion rate of the pregnant mice, body weight and crown-rump length of pups measured on day 18.5 of gestation. Black represents mice treated with DMSO (n = 8); gray represents mice treated with PLX3397 (n = 8); blue represents mice injected with CD45+F4/80+CD206+ anti-inflammatory macrophages (n = 8); red represents mice injected with CD45+F4/80+CD206 pro-inflammatory macrophages (n = 8). (G) Systolic blood pressure (SBP) and UACR of pregnant mice in the four groups. (H) Frequencies of CD44+CD4+IL-17A+ cells, CD44+CD8+ T cells, and CD11b+Ly6G+ granulocytes analyzed by flow cytometry. Data were compared between groups using the Kruskal–Wallis test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001). RNA transcriptome soure data are stored in Dryad Digital Repository, doi:10.5061/dryad.4qrfj6qn0.

Figure 4—figure supplement 1
The construction of RUPP mouse model.

(A) Experimental design of mice model with PE by reducing uterine perfusion pressure (RUPP) on day 12.5 of gestation. Mice with sham operation were considered as controls. Systolic blood pressure (SBP) and urine albumin creatine ratio (UACR) were measured on days 12.5 and 16.5 of gestation, respectively. Mice were sacrificed on day 18.5 of gestation. (B) The embryo abortion rate of pregnant mice, body weight, and crown-rump length of pups measured on day 18.5 of gestation in the Sham and RUPP group. Embryo abortion rate = number of absorbed embryos/total number of embryos. Blue represents mice in the Sham group (n = 9); red represents mice in the RUPP group (n = 9). (C) SBP and UACR of pregnant mice in Sham and RUPP group. (D) The embryo abortion rate of pregnant mice, body weight, and crown-rump length of pups measured on day 18.5 of gestation in the mice injected with PBS or RUPP-derived macrophages. Embryo abortion rate = number of absorbed embryos/total number of embryos. Blue represents mice in the mice injected with PBS group (n = 6); red represents mice in the mice injected with RUPP-derived macrophages group (n = 6). (E) SBP and UACR of pregnant mice in PBS and RUPP-derived macrophages group. (F) Experimental design of GFP pregnant mice model. (G) Frequencies of GFP+CD11b+ fetal-derived macrophages and GFPCD11b+ maternal-derived macrophages were analyzed by flow cytometry. Data were compared between the two groups using the Student’s t-test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 4—figure supplement 2
Clodronate liposomes were used to deplete the macrophages of pregnant mice to demonstrate that pro-inflam Macs lead to immune imbalance.

(A) Frequencies of F4/80+7AAD macrophages in PLX3397-treated mouse were analyzed by flow cytometry. (B) Frequencies of F4/80+7AAD macrophages in clodronate liposomes-treated mouse were analyzed by flow cytometry. (C) Embryo abortion rate of the pregnant mice, body weight and crown-rump length of pups measured on day 18.5 of gestation. Black represents mice treated with control liposomes (n = 6); gray represents mice treated with clodronate liposomes (n = 6); blue represents mice injected with CD45+F4/80+CD206+ anti-inflam Macs (n = 6); red represents mice injected with CD45+F4/80+CD206 pro-inflam Macs (n = 6). (D) Systolic blood pressure (SBP) and UACR of pregnant mice in the four groups. (E) Frequencies of CD44+CD4+IL-17A+ cells, CD44+CD8+ T cells, and CD11b+Gr1+ granulocytes analyzed by flow cytometry. Data were compared between groups using one-way ANOVA and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 5 with 1 supplement
Identification of the macrophage subsets in mice injected with pro-inflam and anti-inflam Macs by using single-cell RNA sequencing (scRNA-seq).

(A) Uniform Manifold Approximation and Projection (UMAP) maps showing the eight types of mouse immune cells at the maternal–fetal interface. (B) Heatmap showing clustering analysis for markers distinguished different type of immune cells. (C) UMAP maps showing the 15 clusters of mouse macrophages was listed. (D) Bar graph showing the frequencies of clusters of macrophages in the two groups of mice was listed in the right panel. (E) UMAP maps showing the distribution of specific markers of clusters 0 and 1. (F) Dot plot depicting GO enrichment terms that were significantly enriched in the differentially expressed genes in clusters 0 and 1 from the pro-inflam Macs group and the control group. (G) Violin plot of specific differential gene expression in clusters 0 and 1 between the pro-inflam Macs and the control groups. Single-cell RNA transcriptome soure data are stored in Dryad Digital Repository, 10.5061/dryad.4qrfj6qn0.

Figure 5—figure supplement 1
Heatmap of different clusters of macrophages.

(A) Heatmap showing clustering analysis for markers distinguished different clusters of macrophages.

Figure 6 with 2 supplements
Identification of the T/NK cells subsets in mice injected pro-inflam and anti-inflam Macs using single-cell RNA sequencing (scRNA-seq).

(A) Uniform Manifold Approximation and Projection (UMAP) maps showing the 12 clusters of mouse T/NK cells was listed. Heatmap showing clustering analysis for markers distinguished 12 different clusters of T/NK cells. (B) Bar graph showing the frequencies of clusters of T/NK cells in the two groups of mice was listed in the down panel. (C) Dot plot depicting GO enrichment terms that were significantly enriched in the differentially expressed genes in cluster 0 from the pro-inflam Macs and the control groups. (D) Violin plot of specific differential gene expression in cluster 0 between the pro-inflam Macs and the control groups. (F) Frequencies of CD4+CD44+ T cells and the percentages of Foxp3+ or IL-17A+ cells in CD4+CD44+ T cells at the maternal–fetal interface in Sham and RUPP group analyzed by flow cytometry. (G) The embryo abortion rate of pregnant mice, body weight, and crown-rump length of pups measured on day 18.5 of gestation in mice injected PBS, Sham mouse-derived or RUPP mouse-derived CD4+CD44+ T cells. Black represents mice injected with PBS (n = 6); blue represents mice injected with Sham mouse-derived CD4+CD44+ T cells (n = 6); red represents mice injected with RUPP mouse-derived CD4+CD44+ T cells (n = 6). (H) Systolic blood pressure (SBP) and UACR of pregnant mice injected with PBS, Sham mouse-derived or RUPP mouse-derived CD4+CD44+ T cells. Data were compared between groups using one-way ANOVA and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001). Data were compared between the two groups using the Student’s t-test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 6—figure supplement 1
Memory-like Th17 cells may be associated with the recurrence of PE.

(A) Heatmap showing clustering analysis for markers distinguished 12 different clusters of T/NK cells. (B) Dot plot depicting gene that were significantly enriched in the differentially expressed genes in cluster 1 and cluster 2. (C) Embryo abortion rate of pregnant mice, body weight and crown-rump length of pups measured on day 18.5 of gestation. Blue represents mice with previous normal pregnancy (n = 9); red represents mice with previous pregnancy with PE (n = 9). (D) Systolic blood pressure (SBP) and UACR of second pregnant mice with a history of PE or NP in the first pregnancy measured on day 16.5 of gestation. (E) Frequencies of CD4+ CD44+ T cells and the levels of IL-17A in CD4+ CD44+ T cells in mice with a pregnancy history with NP and PE analyzed by flow cytometry. Data were compared using the Student’s t-test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001; NS, not significant). Data were compared between the two groups using the Student’s t-test and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 6—figure supplement 2
Identification of the granulocytes subsets in mice injected pro-inflam and anti-inflam Macs use single-cell RNA sequencing (scRNA-seq).

(A) Uniform Manifold Approximation and Projection (UMAP) maps showing the nine clusters of mouse granulocytes. Bar graph showing the frequencies of clusters of granulocytes in the two groups of mice. (B) Heatmap showing clustering analysis for markers distinguished nine different clusters of granulocytes cells. (C) UMAP maps showing the distribution of specific markers of cluster 3. (D) Dot plot depicting GO enrichment terms that were significantly enriched in the differentially expressed genes in cluster 3 from the pro-inflam Macs and the control groups. (E) Violin plot of specific differential gene expression in cluster 3 between the pro-inflam Macs and the control groups. (F) Immunofluorescence co-staining of CD15 (red), CD66b (green), and DAPI (blue) in frozen placental sections. Scale bar, 20 µm. Data were compared between the two groups using the Student’s t-test and represented as mean ± SEM (*p < 0.05).

Figure 7 with 1 supplement
Pro-inflam Macs induce the generation of memory-like Th17 cells via IGF1–IGF1R.

(A) Signaling modules indicated by ligand–receptor pairing between macrophages and other types of immune cells at the maternal–fetal interface using CellPhoneDB. (B) Frequencies of IGF1+CD14+ and IGFIR+CD4+ cells in placentas of individuals with NP and PE (n = 12 in NP group, n = 4 in PE group). (C) Schematic of the experimental workflow to induce memory-like T cells in vitro. Macrophages, after incubating with PBS, NP-EVs or PE-EVs, were co-cultured with CD4+ naive T cells treated with DMSO or BMS-754807. Cells were isolated from human peripheral blood. (D) Frequencies of CD45RO+CCR7+Th17 cells. Black represents CD4+ naive T cells treated with DMSO; gray represents CD4+ naive T cells co-cultured with NP-EV-treated macrophages; blue represents CD4+ naive T cells co-cultured with PE-EV-treated macrophages; red represents CD4+ naive T cells treated with BMS-754807 before co-cultured with PE-EV-treated macrophages (n = 10 in each group). (E) Schematic of mice transferred CD4+ T cells treated with BMS-754807 or PBS. Anti-CD4 antibody was used to deplete CD4+ T cells in mice on day 10.5 of gestation. CD4+ T cells were transferred into mice on day 11.5 of gestation. 20 µg/kg lipopolysaccharide (LPS) was intraperitoneally injected on days 12.5 and 15.5 of gestation to induce a PE-like pregnant mice model. Mice were sacrificed on day 18.5 of gestation. (F) Embryo abortion rate of pregnant mice, body weight and crown-rump length of pups were measured on day 18.5 of gestation. Black represents the control group mice (n = 6); gray represents mice treated with LPS (20 μg/kg) to construct an animal model of PE (n = 6); blue represents anti-CD4 antibody treated PE mice (n = 6); red represents anti-CD4 antibody treated PE mice injected with CD4+ T cells with DMSO treatment (n = 7); orange represents anti-CD4 antibody treated PE mice injected with CD4+ T cells with BMS754807 treatment (n = 7). (G) Systolic blood pressure (SBP) and UACR of pregnant mice in the five groups. (H) The frequencies of CD4+ CD44+ IL-17A+ cells analyzed by flow cytometry. Data were compared between groups using one-way ANOVA and represented as mean ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 7—figure supplement 1
Cell–cell communications in immune cells when pro-inflam Macs accumulated at the maternal–fetal interface.

(A) Abundance of connections between different cell types at the maternal–fetal interface analyzed by CellPhoneDB. (B) Capacities for interactions between immune cells are showed. Each line indicates the ligands expressed by the cell population represented by the same color (labeled). The lines connect to cell types that express cognate receptors. Line thickness is proportional to the number of ligands when cognate receptors are present in the recipient cell type.

Tables

Table 1
Details of the individuals included in the cytometry by time of flight (CyTOF).
ParametersNP (n = 9)PE (n = 8)GDM (n = 8)GDM&PE (n = 7)p value
Age (years)31.33 ± 2.10831.75 ± 4.20635.25 ± 3.41932.43 ± 2.3210.1008
BMI (kg/m2)27.35 ± 1.74030.84 ± 2.88829.56 ± 2.68930.16 ± 2.2800.0518
Gestational age (weeks)38.89 ± 0.55837.39 ± 1.25338.73 ± 0.41936.22 ± 1.372<0.001
Number of living children0.4444 ± 0.4970.500 ± 0.5001.125 ± 0.3310.429 ± 0.4950.0236
Previous abortions0.7778 ± 1.0301.000 ± 0.8661.000 ± 1.0001.714 ± 2.0500.5848
Mean systolic blood pressure109.7 ± 5.400151 ± 7.225106.9 ± 2.848150.7 ± 9.161<0.0001
Table 2
Cytometry by time of flight (CyTOF) antibody panel used for analyzing placentas from individuals with normal pregnancy (NP), preeclampsia (PE), gestational diabetes mellitus (GDM), and GDM complicated with PE (GDM&PE).
TargetMetal tag
CD45HI30
CD3UCHT1
CD68Y1/82A
CD56NCAM16.2
gd TCR5A6.E9
CD19HIB19
CCR6G034E3
CD38HIT2
CD103B-Ly7
CD39A1
CXCR3G025H7
PD-L129E.2A3
PD-1EH12.2H7
CD11CBU15
CD107aH4A3
pAKTD9E
CCR4L291H4
TIGITA15153G
CD20615-2
Helios22F6
CD28CD28.2
pSTAT34/p-Stat3
GITR110416
CD33WM53
CTLA-4BN13
FOXP3PCH101
CD163GHI/61
CD45RAHI100
ICOSC398.4A
CD69FN50
CCR7G043H7
NKG2A131411
CD15W6D3
TIM-3F38-2E2
CD86Fun-1
HLA-DRL243
Granzyme BQA16A02
CD14M5E2
pS6A17020B
CD127A019D5
CD4RPA-T4
CD8RPA-T8
CD11bM1/70
Table 3
Expression levels of markers identified in each immune subset in the placentas.
B cellsCD4+ TCD8+ TgdTCRNK cellsMacrophageGranulocyteMDSCDendritic cells
CD45++++++++++++++++++
CD3++++
CD68++++++++++++
CD56+
gd TCR++++
CD19++++
CCR6+
CD38+++++++++++
CD103
CD39++++++++
CXCR3+++++++++
PD-L1+
PD-1
CD11C+++++++++
CD107a++++
pAKT++++++++++
CCR4
TIGIT
CD206++
Helios+++++++++
CD28+++
pSTAT3+++++
GITR
CD33++
CTLA-4
FOXP3
CD163
CD45RA+++++++
ICOS++
CD69+++++++++
CCR7++++++
NKG2A
CD15++
TIM-3
CD86+++
HLA-DR+++++++++
Granzyme B++++++
CD14++
pS6
CD127+++
CD4++
CD8++
CD11b+++++++
Table 4
The marker profile of PE-specific immune subsets.
Memory-like CD4+ T cellsMemory-like CD8+ T cellsgMDSCsPro-inflammatory macrophagesAnti-inflammatory macrophages
Cluster 08 of CD4+ T cellsCluster 02 of CD8+ T cellsCluster 12 ofCD3CD11b+ cellsCluster 23 ofCD3CD11b+ cellsCluster 25 ofCD3CD11b+ cells
CD45++++++++++
CD3++++
CD68-+++++
CD56
gd TCR-
CD19-
CCR6-
CD38++++++
CD103++
CD39++
CXCR3++
PD-L1
PD-1
CD11C-++++++
CD107a-++
pAKT++
CCR4
TIGIT
CD206++
Helios++
CD28++
pSTAT3-
GITR
CD33+
CTLA-4
FOXP3
CD163
CD45RA-
ICOS++
CD69+
CCR7++++
NKG2A
CD15++
TIM-3
CD86++
HLA-DR-++++
Granzyme B-+
CD14++
pS6
CD127+
CD4++
CD8++
CD11b++++++
Table 5
Details of the individual with NP or PE included in the study.
ParametersNP (n = 30)PE (n = 30)p value
Age (years)31.13 ± 2.974731.50 ± 4.16930.0741
BMI (kg/m2)27.24 ± 2.804829.48 ± 3.92950.0155
Gestational age (weeks)38.68 ± 0.603835.53 ± 3.29<0.0001
Number of living children0.6000 ± 0.49000.4667 ± 0.61820.3665
Table 6
Details of the individual with early- or late-onset preeclampsia included in the study.
ParametersEarly-onset preeclampsia (n = 11)Late-onset preeclampsia (n = 19)p value
The gestational week when blood pressure begins to rise28.82 ± 2.822036.05 ± 1.129<0.0001
Age (years)32.45 ± 4.676832.36 ± 3.84720.9568
BMI (kg/m2)29.64 ± 2.684028.71 ± 4.45060.5388
Gestational age (weeks)33.45 ± 3.788736.77 ± 1.21710.0014
Number of living children0.36 ± 0.50450.36 ± 0.49560.98

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Haiyi Fei
  2. Xiaowen Lu
  3. Zhan Shi
  4. Xiu Liu
  5. Cuiyu Yang
  6. Xiaohong Zhu
  7. Yuhan Lin
  8. Ziqun Jiang
  9. Jianmin Wang
  10. Dong Huang
  11. Liu Liu
  12. Songying Zhang
  13. Lingling Jiang
(2025)
Deciphering the preeclampsia-specific immune microenvironment and the role of pro-inflammatory macrophages at the maternal–fetal interface
eLife 13:RP100002.
https://doi.org/10.7554/eLife.100002.3