Disrupted PGR-B and ESR1 signaling underlies defective decidualization linked to severe preeclampsia

  1. Tamara Garrido-Gomez  Is a corresponding author
  2. Nerea Castillo-Marco
  3. Mónica Clemente-Ciscar
  4. Teresa Cordero
  5. Irene Muñoz-Blat
  6. Alicia Amadoz
  7. Jorge Jimenez-Almazan
  8. Rogelio Monfort-Ortiz
  9. Reyes Climent
  10. Alfredo Perales-Marin
  11. Carlos Simon  Is a corresponding author
  1. Igenomix Foundation, INCLIVA, Spain
  2. Igenomix, Spain
  3. Department of Obstetrics and Gynecology, University and Polytechnic La Fe Hospital, Spain
  4. Department of Obstetrics and Gynecology, School of Medicine, Valencia University, Spain
  5. Obstetrics & Gynecology, BIDMC Harvard University, United States
8 figures, 1 table and 4 additional files

Figures

Figure 1 with 1 supplement
Global RNA-seq transcriptomic results revealed 593 differentially expressed genes (DEGs) in severe preeclampsia (sPE) vs. control samples.

(A) Schematic drawing of the study design used to identify and validate defective decidualization (DD) fingerprinting in sPE. (B) Statistical significance (-log10 FDR) vs. gene expression log2 fold change (FC) is displayed as a volcano plot of global RNA-seq results. Label indicates: downregulated in sPE (blue dots); upregulated in sPE (red dots); not significant genes (grey dots). (C) Heatmap showing the 25 most upregulated and downregulated genes (total = 593; Figure 1—source data 1) of control vs. sPE samples. See also Figure 1—source data 1.

Figure 1—source data 1

The 593 statistically differentially expressed genes (false discovery rate [FDR] < 0.05) with at least 1.2-fold change (FC ≥ 1.2) in severe preeclampsia (sPE) vs. control cases obtained from RNA-seq analysis.

https://cdn.elifesciences.org/articles/70753/elife-70753-fig1-data1-v1.xlsx
Figure 1—figure supplement 1
Transcriptomic analysis based on gestational age at delivery of control samples.

(A) Volcano plot showing there were no significant differentially expressed genes (DEGs) between controls according to gestational age at delivery. Labels show the two criteria that we used to define the DEGs: p-value adjusted (false discovery rate [FDR] < 0.05) and fold-change (FC ≥ 1.2). Legend: not significant (FDR ≥ 0.05). (B) Principal component analysis (PCA) based on 18,476 genes after filtering out lowly expressed genes does not demonstrate clustering based on gestational age.

Defective decidualization (DD) transcriptomics in vitro vs. in vivo.

(A) Common genes between previous in vitro (left) and current in vivo approaches analyzing decidualization (right). Nine genes overlap in both approaches. (B) Box plot showing the average expression of the nine common genes between control (blue boxes) and severe preeclampsia (sPE) (orange boxes) samples. (C) From the 593 differentially expressed genes (DEGs) obtained by global RNA-seq, a subset of 263 DEGs were identified as genes with a human endometrial stromal cell (hESC) origin using the scRNA-seq data published by Wang et al., 2020. See also Figure 2—source data 1.

Figure 2—source data 1

The 593 statistically differentially expressed genes (false discovery rate [FDR] < 0.05) with at least 1.2-fold-change (FC ≥ 1.2) in severe preeclampsia (sPE) vs. control cases obtained from RNA-seq analysis.

https://cdn.elifesciences.org/articles/70753/elife-70753-fig2-data1-v1.xlsx
Severe preeclampsia defective decidualization (sPE-DD) fingerprint composed of 120 differentially expressed genes (DEGs).

(A) Volcano plot showing downregulated (blue) and upregulated (red) genes in sPE from the DD fingerprint. Each point represents one gene; gray points are the rest of the genes obtained in the global RNA-seq analysis. (B) The three most highly downregulated biological process for each major category (red, cell cycle; yellow, DNA damage response; green, cell signaling; blue, cellular response; gray, cell motility; purple, extracellular matrix; pink, immune response; brown, reproductive process). Enrichment index was calculated by -log(p-value). (C) Clustering of DD fingerprint genes shown for reproductive process, response to bacterial molecules, extracellular matrix organization, regulation of receptor signaling, and response to hormones. See also Figure 3—source data 1 and Figure 3—source data 2.

Figure 3—source data 1

List of genes selected as defective decidualization signature in severe preeclampsia (sPE) (120 differentially expressed genes [DEGs] with false discovery rate [FDR] < 0.05 and fold-change [FC] ≥ 1.4).

https://cdn.elifesciences.org/articles/70753/elife-70753-fig3-data1-v1.xlsx
Figure 3—source data 2

Biological process Gene Ontology (GO) terms computed by the 120 genes included in the defective decidualization (DD) fingerprinting in severe preeclampsia (sPE) (N, number of genes associated to GO term; DE, number of genes differentially expressed in this GO term; P.DE, unadjusted p-value; FDR, adjusted p-value).

https://cdn.elifesciences.org/articles/70753/elife-70753-fig3-data2-v1.xlsx
Validation of the defective decidualization (DD) fingerprint in severe preeclampsia (sPE).

(A) Principal component analysis (PCA) based on 120 genes included in the fingerprinting in the training set. Each sample is represented as a colored point (blue, control; orange, sPE). (B) Heatmap dendrogram of expression of the 120 genes included in the final fingerprinting for each sample of the training set (control, n = 12; sPE, n = 17). (C) PCA based on the fingerprinting in the test set. Each sample is represented as a colored point (blue, control; orange, sPE). (D) Heatmap dendrogram of expression of the 120 genes included in the final fingerprinting for each sample of the test set (control, n = 4; sPE, n = 7). See also Figure 4—source data 1.

Figure 4—source data 1

List of genes selected as defective decidualization signature in severe preeclampsia (sPE) (120 differentially expressed genes [DEGs] with false discovery rate [FDR] < 0.05 and fold-change [FC] ≥ 1.4).

https://cdn.elifesciences.org/articles/70753/elife-70753-fig4-data1-v1.xlsx
Estrogen receptor 1 (ER1) and progesterone receptor-B (PR-B) are linked to defective decidualization (DD) fingerprinting in severe preeclampsia (sPE).

(A) Venn diagram displaying genes included in the fingerprinting (120) predominantly expressed in the endometrium based on Human Protein Atlas data that overlap with genes modulated by ESR1 described by Hewitt et al., 2010 and genes associated with PGR silencing described by Mazur et al., 2015. (B) Network showing the connections between proteins codified by DD fingerprinting and the hormonal receptors, ER1 and PR. Shapes indicate different clusters established by String k-means method. Squares, cluster involved in gland morphogenesis and cell migration; circles, cluster involved in extracellular matrix organization and stromal cell differentiation; hexagons; cluster involved in cellular response to DNA damage and regulation of cell cycle. Color gradient indicate gene expression in terms of log2FC. Hub genes are shown with an asterisk. (C-H) Gene expression levels of IHH, MSX2, ESR1, PGR, PGR-A, and PGR-B assessed for sPE (n=13) vs. controls (n=9) by RT-qPCR (gray bars, control; green bars, sPE). RT-qPCR values are expressed as mean± SE. *** p<0.001, ** p<0.01, *p<0.05. (I-J) Tissue sections of control (n=4) and sPE (n=4) endometrium during late secretory phase were immunostained with antibody against ER1 or PR. Nuclei were visualized with DAPI. Scale bar: 50 µM.

Modeling of the molecular mechanism for defective decidualization (DD) in severe preeclampsia (sPE).

(A) Decidualization induced by P4 and E2 in control pregnancy including the interaction of immune response and endothelium. (B) Hypothetical network that could link DD and dysregulated hormone signaling in sPE. All genes were downregulated. Biological processes specified are candidates to be impaired based on functions associated with the observed dysregulation. Red arrows show the downregulation of decidualization modulators.

Author response image 1
Author response image 2
Validation of the DD fingerprint in sPE.

(A) PCA based on 166 genes included in the fingerprinting in the training set. Each sample is represented as a colored point (blue, control; orange, sPE). (B) Heatmap dendogram of expression of the 166 genes included in the final fingerprinting for each sample of the training set (control, n=12; sPE, n=17). Sex of the fetus is represented by color (yellow, female; green, male). (C) PCA based on the fingerprinting in the test set. Each sample is represented as a colored point (blue, control; orange, sPE). (D) Heatmap dendogram of expression of the 166 genes included in the final fingerprinting for each sample of the test set (control, n=4; sPE, n=7). Sex of the fetus is represented by color (yellow, female; green, male).

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Biological sample (Homo sapiens)Endometrial biopsiesUniversity and Polytechnic La Fe Hospital (Valencia, Spain)Freshly isolated from human donors
AntibodyAnti-progesterone receptor antibody [YR85] (rabbit monoclonal anti-human)AbcamCat: AB32085RRID:AB_777452Dilution: (1:50)
AntibodyAnti-estrogen receptor alpha antibody (mouse monoclonal antibody)Santa CruzCat: sc-8002RRID:AB_627558Dilution: (1:50)
AntibodyGoat anti-rabbit IgG H&L (Alexa Fluor 488) (goat polyclonal)AbcamCat: ab150077RRID:AB_2630356Dilution: (1:1000)
AntibodyGoat anti-mouse IgG (H + L) Cro Alexa Fluor 488 (goat polyclonal)InvitrogenCat: A-11001RRID:AB_2534069Dilution: (1:1000)
Sequence-based reagentRT-qPCR primersThis paperSupplementary file 3
Commercial assay or kitQIAsymphony RNA KitQiagen931636Global RNA-seq library preparation
Commercial assay or kitIllumina TruSeq Stranded mRNA sample prep kitIllumina20020595Global RNA-seq library preparation
Commercial assay or kitKapa SYBR fast qPCR kitKapa Biosystems IncKK4602Global RNA-seq library preparation
Commercial assay or kitTruSeq RNA CD Index Plate (96 indexes, 96 samples)Illumina20019792RNA sequencing
Commercial assay or kitNextSeq 500/550 cartridge of 150 cyclesIlluminaFC-404-2002RNA sequencing
Commercial assay or kitSuperScript VILO cDNA Synthesis KitThermo Fisher Scientific11754250RT-qPCR. cDNA preparation
Software, algorithmSTARDobin et al., 2013URL: http://code.google.com/p/rna-star/RRID:SCR_004463RNA-seq analysisRead alignerVersion 2.4.2a
Software, algorithmFastQCURL: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/RRID:SCR_014583RNA-seq analysisQuality of FASTQ file determinationVersion 0.11.2
Software, algorithmSAMtoolsLi et al., 2009URL: http://htslib.org/RRID:SCR_002105RNA-seq analysisSAM and BAM manipulation files Version 1.1
Software, algorithmHTSeqAnders et al., 2015URL: http://htseq.readthedocs.io/en/release_0.9.1/RRID:SCR_005514RNA-seq analysisTo count the number of reads per geneVersion 0.6.1p1
Software, algorithmBEDtoolsQuinlan and Hall, 2010URL: https://github.com/arq5x/bedtools2RRID:SCR_006646RNA-seq analysisTo obtain gene coverageVersion 2.17.0
Software, algorithmedgeRRobinson et al., 2010URL: http://bioconductor.org/packages/edgeR/RRID:SCR_012802RNA-seq analysisTo analyze differentially expressed genesVersion 3.24.3
Software, algorithmStringJensen et al., 2009URL: http://string.embl.de/RRID:SCR_005223Interaction Network.
Software, algorithmCytoscapeShannon et al., 2003URL: http://cytoscape.orgSCR_003032Interaction Network
Software, algorithmCytoHubbaChin et al., 2014URL: http://apps.cytoscape.org/apps/cytohubbaRRID:SCR_017677Interaction Network
OtherCustom scriptsURL: https://github.com/mclemente-igenomix/garrido_et_al_2021The specific script to run RNA-seq analysis

Additional files

Supplementary file 1

Maternal and neonatal characteristics of endometrial donors.

https://cdn.elifesciences.org/articles/70753/elife-70753-supp1-v1.docx
Supplementary file 2

Biological and technical variables of interest for controlling confounding effects in the RNA-seq analysis.

https://cdn.elifesciences.org/articles/70753/elife-70753-supp2-v1.xlsx
Supplementary file 3

RT-qPCR primers list.

https://cdn.elifesciences.org/articles/70753/elife-70753-supp3-v1.doc
Transparent reporting form
https://cdn.elifesciences.org/articles/70753/elife-70753-transrepform1-v1.docx

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. Tamara Garrido-Gomez
  2. Nerea Castillo-Marco
  3. Mónica Clemente-Ciscar
  4. Teresa Cordero
  5. Irene Muñoz-Blat
  6. Alicia Amadoz
  7. Jorge Jimenez-Almazan
  8. Rogelio Monfort-Ortiz
  9. Reyes Climent
  10. Alfredo Perales-Marin
  11. Carlos Simon
(2021)
Disrupted PGR-B and ESR1 signaling underlies defective decidualization linked to severe preeclampsia
eLife 10:e70753.
https://doi.org/10.7554/eLife.70753