Impact Factor 6.684 | CiteScore 2.7
More on impact ›

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 14 December 2021 | https://doi.org/10.3389/fcell.2021.768669

Dysregulated BMP2 in the Placenta May Contribute to Early-Onset Preeclampsia by Regulating Human Trophoblast Expression of Extracellular Matrix and Adhesion Molecules

www.frontiersin.orgYuyin Yi1, www.frontiersin.orgHua Zhu1, www.frontiersin.orgChristian Klausen1, www.frontiersin.orgHsun-Ming Chang1, www.frontiersin.orgAmy M. Inkster2, www.frontiersin.orgJefferson Terry3 and www.frontiersin.orgPeter C. K. Leung1*
  • 1Department of Obstetrics and Gynaecology, BC Children’s Hospital Research Institute, University of British Columbia, Vancouver, BC, Canada
  • 2Department of Medical Genetics, BC Children’s Hospital Research Institute, University of British Columbia, Vancouver, BC, Canada
  • 3Department of Pathology and Laboratory Medicine, BC Children’s Hospital Research Institute, University of British Columbia, Vancouver, BC, Canada

Many pregnancy disorders, including early-onset preeclampsia (EOPE), are associated with defects in placental trophoblast cell invasion and differentiation during early placental development. Bone morphogenetic protein 2 (BMP2) belongs to the TGF-β superfamily and controls various physiological and developmental processes. However, the expression of BMP2 in the placenta and underlying molecular mechanisms of how BMP2 regulates trophoblast function remain unclear. In this study, we analyzed several publicly available microarray and RNA-seq datasets and revealed differences in expression of TGF-β superfamily members between gestational age-matched non-preeclamptic control and EOPE placentas. Importantly, BMP2 levels were significantly reduced in EOPE placentas compared with controls, and RNAscope in situ hybridization further demonstrated BMP2 expression was disrupted in EOPE placental villi. To explore the molecular mechanisms of BMP2-regulated early trophoblast differentiation, we examined BMP2 expression in first-trimester human placenta and found it to be localized to all subtypes of trophoblasts and the decidua. RNA-seq analysis on control and BMP2-treated primary human trophoblast cells identified 431 differentially expressed genes, including several canonical TGF-β/BMP signaling targets (BAMBI, ID1, INHBA, IGFBP3). Gene ontology annotations revealed that differentially expressed genes were involved in cell adhesion and extracellular matrix organization. Furthermore, we identified adhesion molecule with IgG-like domain 2 (AMIGO2) as a novel target for BMP2 that contributed to BMP2-induced trophoblast invasion and endothelial-like tube formation. Overall, our findings provide insight into the molecular processes controlled by BMP2 during early placental development that may contribute to the pathogenesis of EOPE.

Introduction

Preeclampsia is a pregnancy-specific complication that affects 2–8% of pregnancies and causes a high risk of maternal and fetal morbidity and mortality (Khalil et al., 2016). Depending on the time of onset, preeclampsia can be subclassified into early-onset preeclampsia (EOPE; < 34 weeks gestation) and late-onset preeclampsia (≥34 weeks gestation). EOPE only occurs in 5–20% of preeclampsia cases but is more severe compared to late-onset disease. It is now appreciated that early- and late-onset preeclampsia have different pathophysiologies. Late-onset preeclampsia is associated to a greater degree with maternal constitutive factors or pre-existing susceptibility to vascular damage, and to a lesser extent with defective placentation (Raymond and Peterson, 2011). In contrast, EOPE has been attributed to shallow extravillous trophoblast (EVT) invasion and deficiencies in uterine spiral artery remodeling during the first-trimester of pregnancy that lead to insufficient fetal blood supply (Huppertz, 2008).

The placenta is a temporary organ that connects the fetus to the mother and provides essential supports for maintenance and protection of the developing fetus throughout pregnancy. In early pregnancy, placental cells of epithelial lineage, referred to as trophoblasts, undergo extensive proliferation and differentiation to support formation of a functional placenta. Within the placenta, the villous structure is lined with a multinucleated syncytiotrophoblast (STB) layer generated by fusion of mitotically active mononucleated villous cytotrophoblast (CTB) cells. Surrounded by maternal blood, the STB is the major site for nutrient and gas exchange and synthesis of pregnancy-related hormones (e.g. human chorionic gonadotropin, hCG) (Hamilton and Boyd, 1960; Evain-Brion and Malassine, 2003). In anchoring villi, a population of rapidly proliferating CTBs forms a cell column anchoring the placenta to the maternal uterine wall (Vicovac et al., 1995). CTBs detach from the distal ends of the cell column, gradually exit the cell cycle, and differentiate into extravillous trophoblasts (EVTs) that invade the decidua and maternal vasculature (Pijnenborg et al., 1980). Remodeling of uterine spiral arteries promotes a constant and increased maternal blood flow to the placenta, thereby sustaining the oxygen and nutrient supply for fetal growth throughout gestation (Burton et al., 2010). The process of trophoblast differentiation and invasion are precisely modulated by numerous paracrine and autocrine factors, including cytokines, growth factors, adhesion molecules, and extracellular matrix (ECM) components (Bischof et al., 2000).

In EOPE patients, many factors exhibit aberrant maternal and placental expression leading to trophoblast dysfunction and reduced placental perfusion (Nikuei et al., 2015). One group of prominent factors is the transforming growth factor-β (TGF-β) superfamily of growth factors. More than 40 TGF-β superfamily members have been identified and are largely comprised of the TGF-βs, activins and inhibins, bone morphogenic proteins (BMPs), growth differentiation factors (GDFs), and anti-Müllerian hormone (AMH). The BMP subfamily consists of over 20 members and can be classified into several subgroups based on structural homology, including the BMP2/4, BMP5/6/7/8, BMP9/10, and BMP12/13/14 (GDF5/6/7) groups (Katagiri and Watabe, 2016). Recent studies have demonstrated pro-invasive effects of BMP2 on EVTs in vitro (Zhao et al., 2018b), pointing to critical roles for BMPs in trophoblast lineage differentiation and cellular functions. However, to date, there have been no reports comparing the placental expression of BMPs in normal pregnancy or pregnancy complications. Thus, roles of the TGF-β superfamily in placental trophoblast development and in the pathophysiological events of preeclampsia await further elucidation.

Materials and Methods

Microarray Gene Expression Analysis

Microarray and RNA-seq datasets from EOPE and gestational age-matched normotensive controls were obtained through the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) under GSE75010 (Leavey et al., 2016), GSE44711 (Blair et al., 2013), and GSE114691 (Awamleh et al., 2019). Microarray data from placentas of differing gestational age was under GSE100051 (Soncin et al., 2018). Differentially expressed genes were tested with a linear model adjusted for fetal sex using limma package (Ritchie et al., 2015) in R version 3.6.1 (R Core Team, 2019). Hierarchical clustering of normalized gene expression data was performed with the hclust function (R Core Team, 2019) (method = “complete”, distance = “Euclidean””). The sigClust2 package (Huang et al., 2015) in R determined if any clusters were significantly different from one another. The pvClust package (Suzuki and Shimodaira, 2006) assessed the stability of the clusters with 1000 bootstraps.

Patient Recruitment and Tissue Collection

This study was approved by the Research Ethics Board of the University of British Columbia (certificate #H07-01149). First-trimester placental villi and corresponding decidua samples (6–10 weeks) were collected from women undergoing elective surgical termination of pregnancy in the CARE Program at the BC Women’s Hospital & Health Centre. Informed written consent was obtained from participating women. Deidentified term formalin fixed paraffin embedded (FFPE) samples were obtained from placentas submitted for examination between 2014 and 2019 in the BC Children’s Hospital Anatomical Pathology archives. The pre-eclampsia group was defined according to the following clinical features: systolic blood pressure (BP) ≥ 140 mmHg and/or diastolic BP ≥ 90 mmHg after 20 weeks gestation plus proteinuria >300 mg/day or ≥2 + by dipstick. Control samples were derived from otherwise normal preterm multiple gestation deliveries.

Human Primary Trophoblast and HTR-8/SVneo Immortalized Human Extravillous Trophoblasts Cells

Primary human trophoblast cells were isolated from chorionic villous explants and cultured at 37°C in a 5% CO2 air atmosphere as previously described (Irving et al., 1995). Briefly, chorionic villi tips were rinsed in cold Dulbecco’s Modified PBS (Gibco, Thermo Fisher Scientific) and mechanically minced into 1–2 mm fragments. Villous tissue fragments were cultured in Dulbecco’s Minimum Essential Medium (DMEM; Gibco) supplemented with 10% fetal bovine serum (FBS; Gibco), 100  U/ml penicillin, and 100 μg/ml streptomycin for 3–4 days and non-attached fragments were removed. After an additional 10–14 days of culture to allow trophoblast outgrowth, cells were collected from attached fragments by trypsinization.

The HTR-8/SVneo human immortalized extravillous trophoblast cell line was a generous gift from Dr. P. K. Lala (Western University, London, Ontario, Canada) (Graham et al., 1993). Cells were cultured in a 1:1 (v/v) mixture of DMEM/F12 medium (Gibco) supplemented with 10% FBS. Cultures were maintained at 37°C in a humidified atmosphere of 5% CO2 in air.

RNAScope

RNA in situ hybridization was performed using RNAscope 2.5 HD Assay-BROWN (Advanced Cell Diagnostics) following the manufacturer’s instructions (Wang et al., 2012). Briefly, formalin-fixed, paraffin-embedded sections were deparaffinized and pretreated with RNAscope Hydrogen Peroxide (room temperature for 10 min), RNAscope Target Retrieval Reagent (95°C for 15 min), and RNAscope Protease Plus Reagent (40°C for 30 min). RNAScope probes targeting BMP2 (Probe-Hs-BMP2; 430641) and negative control (Probe-dapB; 310043) were incubated on sections for 2 h at 40°C and hybridized with RNAscope 2.5 AMP 1-6. RNAscope signal was detected for 10 min at room temperature with a 1:1 mixture of RNAscope 2.5 HD DAB-A and RNAscope 2.5 HD DAB-B. Sections were counterstained with Harris hematoxylin (Sigma-Aldrich) and mounted in a xylene-based mounting medium. Images were captured with an Olympus BX61 light microscope (Olympus) and quantitative assessment of BMP2 mRNA expression was performed utilizing CellProfiler (www.cellprofiler.org). Nuclei were grouped into 5 bins based on the number of dots per nucleus (Bin 0, 0 dots/nucleus; Bin 1, 1–3 dots/nucleus; Bin 2, 4–9 dots/nucleus; Bin 3, 10–15 dots/nucleus; Bin 4, >15 dots/nucleus). H-scores were calculated by totaling the percentage of nuclei in each bin, according to the weighted formula:

Hscore=0×(% of nuclei in bin0)+1×(% of nuclei in bin1)+2×(% of nuclei in bin2)+3×(% of nuclei in bin3)+4×(% of nuclei in bin4)

RNA-Sequencing and Bioinformatic Analysis

Total RNA from human primary trophoblasts was extracted using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s instructions. Sample quality control was performed with the 2100 Bioanalyzer (Agilent Technologies) and qualifying samples (RNA integrity number >7.5) were prepared following the standard protocol for Next Ultra II Stranded mRNA (New England Biolabs). Sequencing was performed on the Illumina NextSeq 500 with paired-end 42 bp × 42 bp reads at the Biomedical Research Centre (University of British Columbia). Reads were aligned with STAR to the UCSC Homo sapiens hg19 (GRCh37) reference sequence (Dobin et al., 2013). Fragments per kilobase of exon per million reads (FPKMs) and differential expression of novel and reference transcripts were calculated using Cufflinks2 (Trapnell et al., 2012). For downstream analyses, genes had to be ≥10 FPKM and ≥1.5-fold change. Gene Ontology (GO) enrichment analysis and KEGG pathway analysis were performed using the enrichGO package (Yu et al., 2012) in R3.6.1. Data have been deposited in NCBI’s Gene Expression Omnibus repository (Edgar et al., 2002) and are accessible through GEO Series accession number GSE172040 with a reviewer token wdojmqssxbqnvcf. More detailed information of the materials and methods are provided in the Supplementary Material and Methods.

Reverse Transcription Quantitative Real-Time PCR

Total RNA from cells was isolated using TRIzol (Invitrogen, Thermo Fisher Scientific) according to the manufacturer’s instructions. cDNA was synthesized using Moloney murine leukemia virus (M-MLV) reverse transcriptase (Promega, Madison, WI). Quantitative real-time PCR was performed with SYBR Green PCR Master Mix (Applied Biosystems, Thermo Fisher Scientific) as per the manufacturer’s instructions on a QuantStudio 3 instrument and data were analyzed with QuantStudio Design and Analysis software (Applied Biosystems). Relative quantification of mRNA levels was performed by the comparative Cq method with GAPDH as the reference gene and using the formula 2−ΔΔCq. The primers used in this study are listed in Supplementary Table S1.

Western Blots

Cells were lysed in lysis buffer (Cell Signaling Technology) supplemented with protease inhibitor cocktail (Sigma). Equal amounts of protein were separated by 10% SDS-PAGE and transferred to PVDF membranes (Bio-Rad). After blocking with 5% non-fat dry milk in Tris-buffered saline for 1 h at room temperature, membranes were incubated at 4°C overnight with anti-AMIGO2 (1:1000, HPA054004, Sigma-Aldrich) or anti-GAPDH (1:2000, sc-365062, Santa Cruz Biotechnology) primary antibodies diluted in TBS with 5% milk. Membranes were then washed and probed for 1 h at room temperature with horseradish peroxidase (HRP)-conjugated secondary. Immunoreactive proteins were detected using the enhanced chemiluminescence kit (Pierce, Thermo Fisher Scientific) and X-ray film.

Immunohistochemistry and Immunofluorescence

First-trimester human placenta samples (8–10 weeks, n = 5) were fixed in 4% formaldehyde at 4°C overnight. Tissues were dehydrated with 70% ethanol, embedded in paraffin, and serially sectioned at 5 µm onto glass slides. Sections were deparaffinized in xylene, rehydrated through graded ethanol, and processed for wet heat-induced antigen retrieval in a steamer for 20 min with a modified citrate buffer (pH6.0; Dako). Endogenous peroxidase activity was quenched with 3% H2O2 for 10 min at room temperature. After blocking with serum-free protein block (Dako) for 1 h at room temperature, sections were incubated overnight at 4°C with the following primary antibodies: mouse monoclonal anti-human leukocyte antigen-G (HLA-G, 1:100, clone 4H84; ExBio), mouse monoclonal anti-cytokeratin 7 (1:100, RCK105, Santa Cruz Biotechnology), and rabbit polyclonal anti-AMIGO2 (1:100, HPA054004, Sigma-Aldrich). Immunoreactivity was detected with the LSAB 2 System-HRP kit (K0690, Dako) and 3, 3′-diaminobenzidine (DAB) chromogen solution (Dako). Slides were counterstained with Harris hematoxylin (Sigma-Aldrich), dehydrated through graded ethanol to xylene, and mounted in a xylene-based mounting medium.

For immunofluorescence staining, cells were cultured on coverslips, fixed with 4% paraformaldehyde for 10 min at room temperature, and permeabilized with 0.1% Triton X-100 (Sigma) for 5 min. Coverslips were blocked with 5% BSA in PBS for 1 h at room temperature and then incubated with primary antibodies overnight at 4°C. Coverslips were incubated for 1 h at room temperature with Alexa Fluor goat anti-mouse 594 conjugated secondary antibodies (Life Technologies) and mounted using ProLong Gold mounting media containing DAPI (Life Technologies). Images were captured with an Olympus BX61 light microscope (Olympus Canada, Richmond Hill, ON) and analyzed using ImageJ software.

Small Interfering RNA Transfection

To knockdown endogenous AMIGO2, HTR-8/SVneo cells (1.5 × 105/well) or primary trophoblasts (3 × 105/well) were seeded in a 6-well plate 24 h before transfection. Cells were transfected with 50 nM ONTARGETplus SMARTpool siRNA (Dharmacon, Lafayette, CO) targeting AMIGO2 using 5 µl Lipofectamine RNAiMAX (Invitrogen, Thermo Fisher Scientific). The siCONTROL non-targeting pool (Dharmacon) was used as the transfection control. Cells were then cultured for 48 h and harvested for assessment of knockdown efficiency (Western blot) and functional analyses (invasion and endothelial-like tube formation).

Invasion and Endothelial-like Tube Formation Assays

For invasion assays, transwell cell culture inserts (8-μm pore size, 24 wells, BD Biosciences) were coated with 1 mg/ml of growth factor-reduced Matrigel (BD Biosciences). Treated cells (1 × 105 cells/insert) in medium supplemented with 0.1% FBS were incubated for 48 h against a gradient of 10% FBS. Non-invading cells were removed with a cotton swab from the upper side of the membrane. Cells that penetrated the membrane were fixed with cold methanol and stained with crystal violet (0.5%, Sigma) for 30 min and subsequently washed thoroughly with tap water.

For endothelial-like tube formation assays, individual wells of a 96-well plate were each coated with 30 μl growth factor-reduced Matrigel (BD Biosciences, 5 mg/ml). Coated wells were seeded with 4 × 104 HTR-8/SVneo cells suspended in 100 μl DMEM/F12 medium supplemented with 0.1% FBS and incubated at 37°C for 24 h. Digital images were taken with a light microscope and data was analyzed with Image-J software.

Statistical Analysis

Results are presented as the mean ± SEM from a minimum of three independent experiments unless otherwise indicated. Student’s t test was used for comparisons between two groups. Multiple group comparisons were analyzed by one-way analysis of variance (ANOVA) followed by Tukey’s test using GraphPad Prism version 7.00 for Windows (GraphPad Software, La Jolla California United States, www.graphpad.com). Means were considered significantly different from each other if p < 0.05.

Results

TGF-β Pathway-Associated Genes Are Differentially Expressed in Placentas From Healthy and Early-Onset Preeclampsia Pregnancies

Microarray data from 157 placenta samples was obtained from GSE75010(Leavey et al., 2016). EOPE (diagnosis <34 weeks) tends to be more severe and associated with the defects in early placental development and trophoblast function, whereas late-onset preeclampsia (≥34 weeks) may arise from a maternal genetic predisposition to cardiovascular and metabolic disease (Burton et al., 2019). Thus, we filtered to cases delivered on or before 34 weeks of gestation, as this is implicitly EOPE (n = 49). Patients with preterm labour and no other complications before 34 weeks of gestation were included as controls (n = 35). We first investigated the association between EOPE and 84 TGF-β pathway-associated genes encoding TGF-β superfamily ligands, receptors, SMADs, and other relevant proteins. Differential expression analysis was conducted focused on these 84 candidates and identified 17 TGF-β pathway genes that were significantly altered in the EOPE placentas (adjusted p < 0.05, 1.2-fold change cut-off). Among them, 8 genes were upregulated and 9 genes were downregulated in EOPE compared with the matched normal placental tissues (Figure 1A and Supplementary Table S2). We applied hierarchical clustering to all samples on the expression values of the 17 differentially expressed TGF-β pathway genes. Unsurprisingly, samples were clustered into two stable and statistically different clusters. Cluster A included almost all preterm control samples whereas Cluster B contained most EOPE samples, suggesting differentially expressed TGF-β pathway genes successfully distinguished EOPE from gestational age-matched non-preeclamptic controls (Figure 1B). Additionally, we confirmed that maternal ethnicity and fetal sex had no significant impact on the bioinformatic clustering of the study (Table 1). Previous meta-analyses of placental gene expression using the same cohort have identified 5 distinct molecular subclasses (Leavey et al., 2016), including control pregnancies or maternal preeclampsia patients with similar placental gene expression profiles (Cluster 1); “canonical” preeclampsia with placental dysfunction markers (Cluster 2); “immunological” preeclampsia patients associated with maternal-placental incompatibility (Cluster 3); the preterm control samples (Cluster 4); and a Cluster 5 with chromosomal abnormalities. We sought to further compare the clusters identified by differentially expressed TGF-β pathway genes with the clusters in the meta-analyses. As shown in Figure 1C, 24 out of 26 Cluster A samples were in Cluster 1 or 4 reported in the previous study (Leavey et al., 2016), implying Cluster A represents the preterm control with placental gene profile similar to healthy pregnancies. In contrast, 43 out of 58 Cluster B samples corresponded with Cluster 2 samples, the “canonical” preeclampsia subclass that likely arises from a placental origin. Leavey et al. has reported that the preeclampsia disease appears to be more severe in Cluster 2 and more frequently associated with delivered preterm and coexisting diagnoses of IUGR. Interestingly, our analysis showed that all non-EOPE IUGR cases cluster with the Cluster B containing most EOPE cases (Table 1 and Supplementary Figure S1), although we observed no statistically significant distinctions within Cluster B in the diagnosis with coexisting IUGR.

FIGURE 1
www.frontiersin.org

FIGURE 1. Microarray analysis of TGF-β pathway-associated gene expression in healthy and EOPE placentas. (A), Heatmap of 17 differentially expressed TGF-β pathway genes from GSE75010 whole-genome mRNA microarray analysis between EOPE (n = 49) and matched normal control placentas (n = 35). (B), Hierarchical clustering of patients in the GSE75010 dataset using 17 differentially expressed genes associated with the TGF-β pathway and the results were plotted as phylogenetic trees. Samples plotting closer together demonstrated more similar gene expression features, with tips colored based on clinical outcome. The sigClust2 and pvCluster packages in R were used to determine stable and significantly different clusters. p values signify clusters are significantly different from one another. (C), Principal component analysis (PCA) plots for the visualization. Upper panel was modified from original PCA plot in prior gene expression study (Leavey et al., 2016) and showed the 5 Clusters of samples based on expression of 3663 genes. Bottom panel: Two major clusters (Cluster A and Cluster B) identified in current study corresponded with original Cluster 1/4 and Cluster 2

TABLE 1
www.frontiersin.org

TABLE 1. Clinical information Across Clusters (generated from GSE57010).

Among the TGF-β superfamily ligands and receptors, TGFB1, INHBA, INHA, and BMPR2 were significantly upregulated in EOPE, whereas BMP4, BMP5, TGFB3, and TGFBR2 were significantly downregulated, although the differences between groups were small (Supplementary Table S2 and Supplementary Figure S2). Although some studies have shown that plasma levels of TGF-β1 and TGF-β2 are elevated in patients with preeclampsia (Djurovic et al., 1997; Shaarawy et al., 2001; Peracoli et al., 2008), others have not found any differences in serum or placental levels of TGF-βs (Lyall et al., 2001; Perucci et al., 2014). Interestingly, our analysis showed that BMP family ligands – BMP4, BMP5 – were decreased in EOPE. However, the importance of the BMP family in trophoblast differentiation and function during early pregnancy and their relevance to the pathogenesis of placenta-related disorders remains unclear.

Bone Morphogenetic Protein Family Ligands Are Dysregulated in Early-Onset Preeclampsia Placentas

To further investigate the clinical significance of BMP family ligands in EOPE, we focused our subsequent analyses on five ligands with relatively high expression in the human placenta; BMP2, BMP4, BMP5, BMP6 and BMP7 (Figure 2A) (Zhu et al., 2016). In addition, we introduced another two datasets which were analyzed in parallel with GSE75010 (Control, n = 35; EOPE, n = 49) – a microarray dataset GSE44711(Blair et al., 2013) with placenta samples from EOPE pregnancies (n = 8) and gestational age-matched control pregnancies (n = 8), and an RNA-seq dataset (GSE114691(Awamleh et al., 2019)) generated from human EOPE placentas (n = 20) and matched control samples (n = 21). As shown in Figure 2B, BMP4 and BMP5 showed a significant downregulation in the placentas from patients with EOPE compared with controls in all three datasets. In addition, analyses from GSE44711 and GSE114691 datasets showed that BMP2 was significantly downregulated in EOPE placentas compared with controls. Of note, BMP2 showed a decreased expression in EOPE placentas in GSE75010, but this difference was not significant (Figure 2B). In contrast, we only observed a significant but relatively small decrease of BMP6 and BMP7 in EOPE placentas from the GSE114691 RNA-seq dataset. Together these analyses suggest the downregulation of BMP2, BMP4 and BMP5 might be associated with EOPE.

FIGURE 2
www.frontiersin.org

FIGURE 2. BMP family ligands expression in non-preeclamptic and EOPE placentas. (A), RNA levels of BMP family ligands in the human placenta. Data was generated from GSE83115. (B), Expression levels of BMP family ligands in placental tissues from EOPE vs. matched controls. Data were generated from microarray datasets GSE75010 (Control, n = 35; EOPE, n = 49) and GSE44711 (Control, n = 8; EOPE, n = 8) as well as RNA-seq dataset GSE114691 (Control, n = 21; EOPE, n = 20). Expression values were represented as log2 expression microarray data or log2 (counts+1) for RNA-seq data.

Our previous work has suggested that BMP2 contributes to the regulation of invasive trophoblast differentiation and spiral artery remodelling (Zhao et al., 2018a; Zhao et al., 2018b; Zhao et al., 2020a; Zhao J. et al., 2020b). Therefore, we next sought to confirm the dysregulation of placental BMP2 expression in preeclampsia samples. Lacking a BMP2 antibody with appropriate specificity, we used RNAscope in situ hybridization to examine the localization of BMP2 mRNA. A retrospective cohort study was conducted to compare the expression and localization of BMP2 in placentas from EOPE (n = 20) and non-preeclampsia preterm controls (n = 20). The expression of BMP2 in term placental villi was detected in the STB layer and some villous stromal cells. Importantly, the BMP2 mRNA level (H-score) significantly decreased in the EOPE group compared with the gestational age-matched control group (Figure 3 and Supplementary Figure S3). The clinical characteristics of the preterm control pregnancies and EOPE patients are shown in Table 2. The decreased BMP2 levels in the EOPE placenta link the critical role of BMP2 in early trophoblast development to the pathogenesis of EOPE.

FIGURE 3
www.frontiersin.org

FIGURE 3. BMP2 is aberrantly expressed in EOPE placenta. Representative image showing in situ localization of BMP2 mRNA transcript (brown) in control and EOPE placentas (200×). Bottom left; Negative control probe (DapB) as a negative control for the RNAscope in situ hybridization assay (200×). Bottom right; Box and whisker plot of BMP2 H-score (the weighted number of dots/cell, see Materials and Methods) in control and EOPE placentas (p = 0.0051 by Students’t test). Scale bars = 100 µm.

TABLE 2
www.frontiersin.org

TABLE 2. Clinical information of EOPE RNAscope study.

Localization of Bone Morphogenetic Protein 2 in Human First-Trimester Placental Villi

To examine the expression of BMP2 in the human placenta throughout pregnancy, we re-analyzed another microarray-based transcriptome dataset comprised of normal human placenta samples from different gestational age time-points (4–39 weeks gestation, n = 54, GSE100051(Soncin et al., 2018)). As shown in Figure 4A, BMP2 transcript was expressed in placentas across all three trimesters, with significantly higher levels in the first- and second-trimesters than at term. Further analysis of the first-trimester showed an increasing trend in BMP2 levels, peaking at 10–12 weeks gestation (Figure 4B). Despite the increasing evidence indicating the importance of BMP2 in placental trophoblast functions, the exact expression pattern and localization of BMP2 in first-trimester human placenta has yet to be described. Next, we used in situ hybridization to examine the localization of BMP2 mRNA in first-trimester human placental tissues (n = 5; 8–10 weeks gestation). To identify different subtypes of trophoblast cells, we also performed immunohistochemistry with antibodies against the pan-trophoblast marker cytokeratin 7 (KRT7) and the EVT specific marker HLA-G. RNAscope in situ hybridization showed specific and intense BMP2 mRNA localization to cells within the stroma of placental villi as well as to CTBs, STBs and EVTs (Figure 4C and Supplementary Figure S4). Notably, BMP2 mRNA was localized to EVTs within the anchoring columns and decidua as indicated by positive staining for the HLA-G antibody (Figure 4C). BMP2 mRNA was also detected in maternal decidua cells whereas no specific positive signal was detected in serial sections incubated with the negative control probe (Figure 4C).

FIGURE 4
www.frontiersin.org

FIGURE 4. BMP2 expression and localization at the early human maternal-fetal interface. (A,B) Microarray analysis of BMP2 levels in placentas from GSE100051. (A), BMP2 levels in placentas across gestation (1st trimester, n = 42; 2nd trimester, n = 7; term, n = 5). Results are expressed as the mean ± SD and statistical analysis was performed using the Kruskal-Wallis test with Dunn’s correction. (B), Microarray analysis shows an increase in BMP2 levels in human placentas during the first trimester. n = 5 for each gestational age indicated, except 4 weeks gestation (n = 3) and 8 weeks gestation (n = 4). (C), Top left: Representative image showing in situ localization of BMP2 mRNA transcript (brown) within a first-trimester placental villus (9 weeks +5 days gestation; 100×). Top centre: Magnification of the boxed area in the top left panel (200×). Top right: Negative control probe (DapB) incubation of a serial section as a negative control for the RNAscope in situ hybridization assay (100×). Bottom panels: Hematoxylin and eosin (H&E) staining (left) as well as immunohistochemical staining for cytokeratin 7 (KRT7; centre) and human leukocyte antigen-G (HLA-G; right) was performed on serial sections of the same placental villi (100×). CTB, cytotrophoblast; STB, syncytiotrophoblast; EVT: extravillous trophoblast. Scale bars = 100 µm.

Bone Morphogenetic Protein 2 Orchestrates a Gene Network of Cell Adhesion and Extracellular Matrix Molecules in Primary Human Trophoblasts

Primary human trophoblasts were isolated from first-trimester placental villi as previously described (Irving et al., 1995). To confirm the characteristics of these primary human trophoblasts, cells were immunofluorescently stained with KRT7 and HLA-G and assessed microscopically to confirm purity of invasive trophoblasts >90% (Figure 5A). To gain mechanistic insight into how BMP2 modulates trophoblast development and function, RNA-seq was used to analyze global gene expression. In this study, primary trophoblast cells from five first-trimester placentas (n = 5) were treated with vehicle control or 25 ng/ml BMP2 for 24 h as previously described (Zhao et al., 2018b). RNA-seq analysis identified 431 differentially expressed genes between control and BMP2 treated cells (253 genes upregulated, 178 genes downregulated; false discovery rate <0.05; Supplementary Table S3). A heat-map highlighting patterns of differentially expressed genes across control and BMP2 treated trophoblasts is shown in Figure 5B. RNA-seq analysis identified several canonical TGF-β/BMP signaling targets that were differentially expressed between control and BMP2-treated trophoblasts, including several that have previously been identified as critical for the regulation of human trophoblast invasion and differentiation by BMP2 (BAMBI, CDH2, ID1, INHBA, IGFBP3; Figure 5C) (Zhao et al., 2018a; Zhao et al., 2018b; Zhao et al., 2020a; Zhao et al., 2020b). RT-qPCR was used to validate the differential expression of a subset of TGF-β/BMP2 target genes (ID1, ID2, ID3, ID4, ENG and IL6; Figure 5D). Gene Ontology (GO) analysis revealed that BMP2-regulated genes are involved in cell adhesion, cell migration, cell proliferation and ECM organization (Figure 5E and Supplementary Table S4). Unsurprisingly, KEGG pathway analysis demonstrated that total differentially expressed genes were enriched in TGF-β signaling and pluripotency stem cell pathways (Figure 5F). Interestingly, BMP2 also regulates a set of genes enriched in FoxO signaling, insulin resistance, and fluid shear stress and atherosclerosis, suggesting potential links to trophoblast metabolic regulation and trophoblast-endothelial cell remodelling (Figure 5F).

FIGURE 5
www.frontiersin.org

FIGURE 5. Comparison of differentially expressed genes in human primary trophoblasts treated with or without BMP2. (A). Purified primary human trophoblast cells were immunofluorescently stained with antibodies directed against cytokeratin 7 (KRT7) and human leukocyte antigen-G (HLA-G), and nuclei were stained with DAPI (200 × magnification, scale bars = 100 µm). (B), Hierarchical clustering of control and BMP2 treated trophoblasts using a total of 431 differentially expressed genes from mRNA-seq analysis (row Z-score). (C), RNA-seq data of differentially expressed BMP signaling targets in primary trophoblasts treated with or without 25 ng/ml BMP2 (FPKM, fragments per kilobase of exon per million reads). False discovery rate values for all the genes presented were <0.05 from CuffDiff2. (D), RT-qPCR validation of the differential expression of BMP signaling target genes in primary trophoblasts treated with or without 25 ng/ml BMP2, n = 5. Data were normalized to GAPDH and values represent the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001 from Student’s t test. (E), Gene Ontology (GO) analysis of differentially expressed genes indicating the top 20 enriched GO terms in BMP2 treated primary trophoblasts. (F), Top 5 KEGG pathways enriched in BMP2 treated primary trophoblasts.

To verify the functional connectivity of differentially expressed genes, we performed functional annotation clustering using Cytoscape software with the ClueGo plugin. Network mapping identified regulation of cell adhesion and ECM-associated processes as the most enriched clusters, but also showed enrichment of pathways related to circulatory system development, muscle structure development, and cellular response to TGF-β stimulation (Figure 6A). Genes associated with cellular adhesion included CDH2, SERPINE2, SPON2, CD9, ID1, AMIGO2, WNT5A, CCN2, CCL2, CAV1 (Figures 6B,C). ECM-associated genes regulated by BMP2 included ECM components (COL6A1, COL7A1, COL11A1, COL13A1, FBLN1, FBLN2, FBLN5), ECM receptors (ITGA2, ITGA6, DDR1), and proteinases (ADAM12, ADAMTS2, MMP11; Figures 6B,C). RT-qPCR was used to validate the regulation of a subset of these cell adhesion and ECM-associated genes by BMP2 (Figure 6D). Our RNA-seq analysis provides strong evidence that BMP2 regulates gene programmes and processes that are critical for cell adhesion and cell-ECM interaction, thereby promoting trophoblast differentiation along the extravillous invasive pathway.

FIGURE 6
www.frontiersin.org

FIGURE 6. BMP2 regulates the expression of cell adhesion and ECM-associated genes in primary human placental trophoblasts. (A), Functional clustering map for genes differentially expressed in control and BMP2 treated trophoblasts. The map was generated using Cytoscape with ClueGo plug-in and each node represents a pathway or GO term. Links between nodes represent the genes shared by both pathways (p cutoff 0.01, Kappa score threshold 0.5). Clusters were manually drawn and annotated. Groupings with less than 4 connections were excluded from the final list of networks. (B), Gene interaction networks of differentially expressed genes within the top 4 enriched GO terms (biological adhesion, cell adhesion, extracellular matrix organization, and extracellular structure organization). Each small dot represents a gene and large dots represent GO terms. Links indicate the association of each gene with the GO terms. Log-transformed fold change for each gene is shown by a colored scale. (C), Relative expression levels (row Z-score) of significantly altered genes belonging to the cell-adhesion and ECM-associated categories from RNA-seq data. (D), RT-qPCR validation of differential expression of cell adhesion and ECM-associated genes in primary trophoblasts treated with or without 25 ng/ml BMP2, n = 5. Data were normalized to GAPDH and values represent the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001 from Student’s t test.

Identification of Adhesion Molecule With IgG-Like Domain 2 as a Novel Mediator Required for the Pro-Invasive Effects of Bone Morphogenetic Protein 2 on Human Trophoblast Cell Invasion

Among the transcriptional targets of BMP2 in primary human trophoblasts, we identified AMIGO2 as a novel adhesion molecule associated with cell adhesion and ECM formation (Rabenau et al., 2004) based on its high expression in placenta among human organs and tissues (Figure 7A). In addition, analysis of RNA-seq data of first-trimester human CTB, EVT, STB, and placenta-derived stromal cells (Okae et al., 2018) found that AMIGO2 expression was highest in EVTs (Figure 7B). Further immunohistochemical analysis of first-trimester placentas showed strong expression of AMIGO2 in the cytoplasm and nucleus of villous CTBs and EVTs within anchoring columns as well as decidua cells, whereas expression was weak in STBs (Figure 7C and Supplementary Figure S5). In both RNA-seq and RT-qPCR analyses, AMIGO2 expression was elevated by BMP2 treatment in primary human trophoblasts (Figures 6C,D). Western blot results confirmed the upregulation of AMIGO2 protein levels following treatment of primary and immortalized human trophoblasts HTR-8/SVneo with BMP2 (Figures 8A,B). Together with our demonstration of AMIGO2 expression in trophoblasts, we hypothesized that AMIGO2 may play a role in regulating trophoblast invasiveness. To study the function of AMIGO2 in human trophoblasts, we used Matrigel-coated transwell assay to model the process of trophoblast invasion through ECM during early placental development. As shown in Figure 8C, siRNA-mediated knockdown of AMIGO2 significantly reduced primary human trophoblast cell invasion, indicating endogenous AMIGO2 is required for trophoblast invasion. We next examined the involvement of AMIGO2 in BMP2-induced trophoblast invasion. Consistent with our previous studies (Zhao et al., 2018b), treatment with BMP2 significantly increased primary human trophoblasts cell invasion and, importantly, these effects were attenuated by AMIGO2 knockdown (Figure 8C). BMP2 has also been shown to promote endothelial-like tube formation of HTR-8/SVneo cells (Zhao et al., 2020a). We showed that knockdown of AMIGO2 dramatically impaired basal and BMP2-mediated endothelial-like tube formation in trophoblastic HTR-8/SVneo cells (Figure 8D and Supplementary Figure S6). Taken together, these results suggest that AMIGO2 contributes to BMP2-induced invasive and endovascular differentiation of human trophoblasts during early placental development. Analysis of microarray datasets showed that AMIGO2 levels were significantly reduced in placenta samples from EOPE patients compared with gestational age-matched non-preeclamptic controls (GSE75010 and GSE114691), suggesting that AMIGO2 could be relevant to the pathogenesis of EOPE (Figure 8E).

FIGURE 7
www.frontiersin.org

FIGURE 7. AMIGO2 expression and localization in human placenta. (A), Highest RNA levels of AMIGO2 detected in placenta among the human organs and tissues. Data was generated from GSE83115. (B), The RNA-seq data of AMIGO2 levels in primary CTB (n = 3), EVT (n = 3), STB (n = 3), and placenta-derived stromal cells (n = 2). Data were generate from Okae H et al.,and presented as log2(FPKM+1) (Okae et al., 2018). (C), Representative immunohistochemical image showing AMIGO2 localization (brown) within a first-trimester placental villus (9 weeks + 4 days gestation; 100 × magnification, scale bar = 100 µm). Cytokeratin 7 (KRT7; centre) and human leukocyte antigen-G (HLA-G; right) was performed on serial sections of the same placental villi. CTB, cytotrophoblast; SCT, syncytiotrophoblast; EVT, extravillous trophoblast; CC, cell column.

FIGURE 8
www.frontiersin.org

FIGURE 8. AMIGO2 mediates BMP2-induced human trophoblast invasion and endothelial-like tube formation. (A,B), Representative Western blot images showing the effect of exogenous BMP2 treatment on AMIGO2 protein levels in primary human trophoblasts (n = 5) and HTR-8/SVneo cells (n = 3). (C,D), Cells were transfected for 48 h with non-targeting control siRNA (siCtrl) or siRNA targeting AMIGO2 (siAMIGO2) prior to treatment with vehicle (Ctrl) or 25 ng/ml BMP2. AMIGO2 knockdown was confirmed by Western blot (left panels). (C), Representative images from Matrigel-coated transwell invasion assays examining the effect of AMIGO2 knockdown on BMP2-induced invasion of primary human trophoblasts (n = 5; centre panel). Summarized quantitative results are expressed as the mean ± SD of five independent experiments and values without common letters are significantly different (p < 0.05; right panel). (D), Representative images of the effect of AMIGO2 knockdown on BMP2-induced endothelial-like tube formation in HTR-8/SVneo cells (n = 3; centre panel). Summarized quantitative results are expressed as the mean ± SD of three independent experiments and values without common letters are significantly different (p < 0.05; right panel). (E), AMIGO2 expression is downregulated in EOPE placentas compared with non-preeclamptic placenta tissues. Data were generated from microarray datasets GSE75010 (Control, n = 35; EOPE, n = 49) and GSE44711 (Control, n = 8; EOPE, n = 8) as well as RNA-seq dataset GSE114691 (Control, n = 21; EOPE, n = 20).

Discussion

Our data highlight the disrupted expression of TGF-β superfamily members as a potential etiologic factor in the development of EOPE. We further emphasize the differential expression and putative role of one TGF-β superfamily member, BMP2, in the pathologic condition of EOPE. For the first time, we profile the expression of BMP2 at the fetal-maternal interface in early pregnancy. In addition, our data provide evidence that BMP2 transcriptionally regulates cell adhesion and ECM molecules to promote trophoblast differentiation along the extravillous pathway. Our study also uncovers AMIGO2 as a novel target for BMP2 and sheds new light on the expression and the pro-invasive effects of AMIGO2 in human trophoblasts.

EOPE is associated with shallow EVT invasion and deficiencies in uterine spiral artery remodeling during the first-trimester of pregnancy, and dysregulation of many factors resulting in the trophoblast dysfunction has been implicated in EOPE pathogenesis (Nikuei et al., 2015). Here, we have generated a comprehensive expression profile of TGF-β superfamily-associated genes in the placentas of EOPE versus normotensive women. We found that TGFB1 was upregulated in EOPE, which supports the elevated expression and serum levels of TGF-β1 in patients with preeclampsia (Djurovic et al., 1997; Peracoli et al., 2008). Blocking TGF-β3 has been shown to restore the invasive capability of EVTs in preeclamptic pregnancies, revealing its suppressive role for EVT invasion (Caniggia et al., 1999). In this context, the downregulation of TGFB3 in our analysis could suggest some type of compensatory effect with respect to EVT invasion. Supported by our data, activin A (INHBA) and inhibin A (INHA), but not inhibin B (INHB), are increased in patients with preeclampsia (Laivuori et al., 1999; Yair et al., 2001; Silver et al., 2002; Bersinger et al., 2003). Given that TGF-β ligands have been implicated in the control of trophoblast invasion and migration (Graham et al., 1992; Jones et al., 2006; Zhao et al., 2018b), our data support placental contributions to the pathology of EOPE and highlight the TGF-β superfamily as a potential etiologic factor in EOPE. Interestingly, our results have identified that three BMPs subfamily ligand genes (BMP2, BMP4, BMP5) are downregulated in placentas from women with EOPE compared to controls, suggesting that BMP pathways are relevant to the pathogenesis of EOPE. Increased expression of BMPR2, which encodes type II receptors for BMPs, could constitute a compensatory mechanism in response to the reduced levels of BMP ligands in EOPE. The dysregulation of BMP pathway components in placental tissues from women with EOPE may provide novel diagnostic and/or therapeutic strategies to manage placenta-related pregnancy disorders.

Previous studies of several mouse models have demonstrated that BMP signaling is involved in the regulation of blastocyst implantation, endometrial stromal cell decidualization, and placental development. BMP2 has been detected in the uterus of mice and its expression is spatiotemporally correlated with embryo implantation during pregnancy (Paria et al., 2001). Female mice with conditional deletion of Bmp2 in the uterus are infertile due to impaired decidualization of the endometrial stroma (Lee et al., 2007). Moreover, conditional ablation of Bmpr2 in the uterus causes defective spiral artery remodeling at mid-gestation that results in IUGR and fetal death (Nagashima et al., 2013). In humans, relatively little is known about the expression and roles of BMPs in the uterus and placenta during pregnancy. BMP2 is induced during in vitro decidualization of human endometrial stromal cells in response to steroids and cAMP, and addition of BMP2 stimulates decidualization (Li et al., 2007). Additionally, BMP2 has been implicated in regulating human trophoblast invasion and spiral artery remodeling using in vitro cell models (Zhao et al., 2018b; Zhao et al., 2020b). Our study provides the first evidence of BMP2 localization in all trophoblast cell populations of human first-trimester placenta and decidua tissue. Microarray analysis shows an increasing level of BMP2 in the human placenta during the first trimester of pregnancy. These data suggest that autocrine/paracrine BMP2 produced by both placenta trophoblasts and decidua promotes human trophoblast differentiation along the EVT pathway. Although we did not observe differential expression of BMP2 in GSE75010, analyses from two other datasets (GSE44711, and GSE114691) showed reduced levels of BMP2 in EOPE. The difference across datasets may be due to differences in patients’ ethnicity, platforms used, and study design between different cohorts. Importantly, the RNAscope data support our bioinformatic analysis and reveal that the mRNA expression of BMP2 is abnormal in the EOPE placentas. Our finding of altered expression of BMP2 and receptors in EOPE placentas, combined with the aforementioned mouse work, indicates that BMP2 signaling might be important in human placental development and trophoblast biology. A recent study has shown that the serum concentration of BMP2 is significantly higher in first-trimester pregnancies compared with non-pregnant individuals, suggesting that the developing placenta may be the source of the increased circulating BMP2(You et al., 2021). Future research designed to examine the maternal circulating levels of BMPs in routine blood samples and to compare the BMPs serum levels between EOPE and non-preeclamptic women will provide the theoretical basis for clinical decision-making regarding the use of BMP2 as a potential biomarker for EOPE diagnosis.

Transcriptome analysis has shown BMP2 targets are mainly involved in regulating cell adhesion and ECM organization. The process of trophoblasts invading into the maternal decidua and eventually to the endometrial spiral arteries is accompanied by altered expression of cell adhesion molecules and remodelling of the ECM by various proteinases, including serine proteases, metalloproteinases and collagenases (McEwan et al., 2009). Indeed, we show that BMP2 induces expression of cellular adhesion genes such as CDH2 and AMIGO2 as well as ECM-associated genes including ECM components (COL6A1, COL7A1), ECM receptors (ITGA2, ITGA6), and proteinases (ADAM12, MMP11). Combined with our previous studies of BMP2 expression and function in human trophoblast cells, our findings suggest that BMP2 promotes differentiation along the EVT pathway and enhances molecular pathways related to cell adhesion and ECM remodeling.

AMIGO2 is a member of the adhesion molecule with IgG-like domain family of proteins (AMIGO1, 2, 3), which were identified as type I transmembrane proteins induced by the neurite-promoting protein amphoterin in neurons (Kuja-Panula et al., 2003). AMIGO2 is preferentially expressed in the central nervous system and plays an important role in neuronal development and survival (Kuja-Panula et al., 2003; Ono et al., 2003). AMIGO2 has also been linked to the processes of tumorigenicity and metastasis due to its ability to regulate tumour cell adhesion, migration and survival (Rabenau et al., 2004; Fontanals-Cirera et al., 2017; Kanda et al., 2017). In addition, recent studies suggest that AMIGO2 contributes to vascular development and angiogenesis by regulating endothelial cell adhesion, migration, survival, and angiogenesis (Park et al., 2015). Our discovery that AMIGO2 is preferentially expressed in progenitor CTBs and invasive EVTs combined with our finding that disruption of AMIGO2 impairs BMP2-regulated trophoblast cell invasion and endothelial-like tube formation suggest it is required for BMP2-mediated trophoblast invasive and endovascular differentiation. Moreover, the downregulation of placental AMIGO2 in EOPE suggests its involvement in the pathophysiology of EOPE, which is associated with aberrant trophoblast functions. Future studies are required to determine the precise mechanism of how AMIGO2 modulates trophoblast functions. AMIGO2 contains six leucine-rich repeats (LRRs), a single IgG-like domain, a transmembrane domain, and a cytosolic domain (Kuja-Panula et al., 2003). It is possible that AMIGO2 may regulate trophoblast adhesion to ECM and decidua cells via its LRR domains.

Our study combined transcriptomic datasets analyses, in situ gene expression profiling, and in vitro functional evaluation to elucidate the novel mechanism of BMP2 action in the first-trimester maternal-fetal interface and links the critical role of BMP2 in early trophoblast development to the potential pathological mechanisms for poor placentation in EOPE.

Data Availability Statement

Data have been deposited in NCBI’s Gene Expression Omnibus repository and are accessible through GEO Series accession number GSE172040.

Ethics Statement

The studies involving human participants were reviewed and approved by the Research Ethics Board of the University of British Columbia (certificate #H07-01149). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YY, HZ, CK, H-MC, and PL conceived the research. AI contributed to the clinical data and microarray analyses. YY and HZ conducted RNA-seq and analyses of sequence data. JT contributed to the clinical sample collection and critical discussion. YY designed and performed research, analyzed and interpreted the data, and prepared the manuscript in partnership with CK and PL.

Funding

This study was supported by a Canadian Institutes of Health Research Foundation Grant to PL (#143317).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank the patients and staff of the CARE program at the BC Women’s Hospital & Health Centre for providing the placental samples used in this study. We thank Dr. Wendy P. Robinson for discussions on microarray analysis. We thank Ryan Vander Werff at the Biomedical Research Centre (University of British Columbia) for his assistance with RNA-seq analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.768669/full#supplementary-material

References

Awamleh, Z., Gloor, G. B., and Han, V. K. M. (2019). Placental microRNAs in Pregnancies With Early Onset Intrauterine Growth Restriction and Preeclampsia: Potential Impact on Gene Expression and Pathophysiology. BMC Med. Genomics. 12 (1), 91. doi:10.1186/s12920-019-0548-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bersinger, N. A., Smárason, A. K., Muttukrishna, S., Groome, N. P., and Redman, C. W. (2003). Women With Preeclampsia Have Increased Serum Levels of Pregnancy‐Associated Plasma Protein A (PAPP‐A), Inhibin A, Activin A and Soluble E‐selectin. Hypertens. Pregnancy. 22 (1), 45–55. doi:10.1081/PRG-120016794

PubMed Abstract | CrossRef Full Text | Google Scholar

Bischoff, P., Meisser, A., and Campana, A. (2000). Paracrine and Autocrine Regulators of Trophoblast Invasion- A Review. Placenta. 21 (Suppl. A), S55–S60. doi:10.1053/plac.2000.0521

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, J. D., Yuen, R. K. C., Lim, B. K., McFadden, D. E., von Dadelszen, P., and Robinson, W. P. (2013). Widespread DNA Hypomethylation at Gene Enhancer Regions in Placentas Associated With Early-Onset Pre-Eclampsia. Mol. Hum. Reprod. 19 (10), 697–708. doi:10.1093/molehr/gat044

PubMed Abstract | CrossRef Full Text | Google Scholar

Burton, G. J., Jauniaux, E., and Charnock-Jones, D. S. (2010). The Influence of the Intrauterine Environment on Human Placental Development. Int. J. Dev. Biol. 54 (2-3), 303–312. doi:10.1387/ijdb.082764gb

CrossRef Full Text | Google Scholar

Burton, G. J., Redman, C. W., Roberts, J. M., and Moffett, A. (2019). Pre-Eclampsia: Pathophysiology and Clinical Implications. BMJ. 366, l2381. doi:10.1136/bmj.l2381

PubMed Abstract | CrossRef Full Text | Google Scholar

Caniggia, I., Grisaru-Gravnosky, S., Kuliszewsky, M., Post, M., and Lye, S. J. (1999). Inhibition of TGF-Β3 Restores the Invasive Capability of Extravillous Trophoblasts in Preeclamptic Pregnancies. J. Clin. Invest. 103 (12), 1641–1650. doi:10.1172/JCI6380

CrossRef Full Text | Google Scholar

Djurovic, S., Schjetlein, R., Wisløff, F., Haugen, G., Husby, H., and Berg, K. (1997). Plasma Concentrations of Lp(a) Lipoprotein and TGF-Β1 Are Altered in Preeclampsia. Clin. Genet. 52 (5), 371–376. doi:10.1111/j.1399-0004.1997.tb04356.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics. 29 (1), 15–21. doi:10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res. 30 (1), 207–210. doi:10.1093/nar/30.1.207

PubMed Abstract | CrossRef Full Text | Google Scholar

Evain-Brion, D., and Malassine, A. (2003). Human Placenta as an Endocrine Organ. Growth Horm. IGF Res. 13, S34–S37. doi:10.1016/S1096-6374(03)00053-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Fontanals-Cirera, B., Hasson, D., Vardabasso, C., Di Micco, R., Agrawal, P., Chowdhury, A., et al. (2017). Harnessing BET Inhibitor Sensitivity Reveals AMIGO2 as a Melanoma Survival Gene. Mol. Cell. 68 (4), 731–744. doi:10.1016/j.molcel.2017.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, C. H., Hawley, T. S., Hawley, R. C., MacDougall, J. R., Kerbel, R. S., Khoo, N., et al. (1993). Establishment and Characterization of First Trimester Human Trophoblast Cells With Extended Lifespan. Exp. Cell Res. 206 (2), 204–211. doi:10.1006/excr.1993.1139

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, C. H., Lysiak, J. J., McCrae, K. R., and Lala, P. K. (1992). Localization of Transforming Growth Factor-β at the Human Fetal-Maternal Interface: Role in Trophoblast Growth and Differentiation1. Biol. Reprod. 46 (4), 561–572. doi:10.1095/biolreprod46.4.561

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamilton, W. J., and Boyd, J. D. (1960). Development of the Human Placenta in the First Three Months of Gestation. J. Anat. 94, 297–328.

Google Scholar

Huang, H., Liu, Y., Yuan, M., and Marron, J. S. (2015). Statistical Significance of Clustering Using Soft Thresholding. J. Comput. Graphical Stat. 24 (4), 975–993. doi:10.1080/10618600.2014.948179

PubMed Abstract | CrossRef Full Text | Google Scholar

Huppertz, B. (2008). Placental Origins of Preeclampsia. Hypertension. 51 (4), 970–975. doi:10.1161/HYPERTENSIONAHA.107.107607

PubMed Abstract | CrossRef Full Text | Google Scholar

Irving, J. A., Lysiak, J. J., Graham, C. H., Hearn, S., Han, V. K. M., and Lala, P. K. (1995). Characteristic's of Trophoblast Cells Migrating from First Trimester Chorionic Villus Explants and Propagated in Culture. Placenta. 16 (5), 413–433. doi:10.1016/0143-4004(95)90100-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, R. L., Findlay, J. K., Farnworth, P. G., Robertson, D. M., Wallace, E., and Salamonsen, L. A. (2006). Activin A and Inhibin A Differentially Regulate Human Uterine Matrix Metalloproteinases: Potential Interactions During Decidualization and Trophoblast Invasion. Endocrinology. 147 (2), 724–732. doi:10.1210/en.2005-1183

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanda, Y., Osaki, M., Onuma, K., Sonoda, A., Kobayashi, M., Hamada, J., et al. (2017). Amigo2-Upregulation in Tumour Cells Facilitates Their Attachment to Liver Endothelial Cells Resulting in Liver Metastases. Sci. Rep. 7, 43567. doi:10.1038/srep43567

PubMed Abstract | CrossRef Full Text | Google Scholar

Katagiri, T., and Watabe, T. (2016). Bone Morphogenetic Proteins. Cold Spring Harb Perspect. Biol. 8 (6), a021899. doi:10.1101/cshperspect.a021899

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalil, A., O'Brien, P., and Townsend, R. (2016). Current Best Practice in the Management of Hypertensive Disorders in Pregnancy. Integr. Blood Press. Control. 9, 79–94. doi:10.2147/IBPC.S77344

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuja-Panula, J., Kiiltomäki, M., Yamashiro, T., Rouhiainen, A., and Rauvala, H. (2003). AMIGO, a Transmembrane Protein Implicated in Axon Tract Development, Defines a Novel Protein Family With Leucine-Rich Repeats. J. Cell Biol. 160 (6), 963–973. doi:10.1083/jcb.200209074

CrossRef Full Text | Google Scholar

Laivuori, H., Kaaja, R., Turpeinen, U., Stenman, U.-H., and Ylikorkala, O. (1999). Serum Activin A and Inhibin A Elevated in Pre-Eclampsia: No Relation to Insulin Sensitivity. BJOG:An Int. 106 (12), 1298–1303. doi:10.1111/j.1471-0528.1999.tb08185.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Leavey, K., Benton, S. J., Grynspan, D., Kingdom, J. C., Bainbridge, S. A., and Cox, B. J. (2016). Unsupervised Placental Gene Expression Profiling Identifies Clinically Relevant Subclasses of Human Preeclampsia. Hypertension. 68 (1), 137–147. doi:10.1161/HYPERTENSIONAHA.116.07293

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, K. Y., Jeong, J.-W., Wang, J., Ma, L., Martin, J. F., Tsai, S. Y., et al. (2007). Bmp2 Is Critical for the Murine Uterine Decidual Response. Mol. Cell Biol. 27 (15), 5468–5478. doi:10.1128/MCB.00342-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Kannan, A., Wang, W., Demayo, F. J., Taylor, R. N., Bagchi, M. K., et al. (2007). Bone Morphogenetic Protein 2 Functions via a Conserved Signaling Pathway Involving Wnt4 to Regulate Uterine Decidualization in the Mouse and the Human. J. Biol. Chem. 282 (43), 31725–31732. doi:10.1074/jbc.M704723200

CrossRef Full Text | Google Scholar

Lyall, F., Simpson, H., Nicola Bulmer, J., Barber, A., and Courtenay Robson, S. (2001). Transforming Growth Factor-β Expression in Human Placenta and Placental Bed in Third Trimester Normal Pregnancy, Preeclampsia, and Fetal Growth Restriction. Am. J. Pathol. 159 (5), 1827–1838. doi:10.1016/s0002-9440(10)63029-5

CrossRef Full Text | Google Scholar

McEwan, M., Lins, R. J., Munro, S. K., Vincent, Z. L., Ponnampalam, A. P., and Mitchell, M. D. (2009). Cytokine Regulation During the Formation of the Fetal-Maternal Interface: Focus on Cell-Cell Adhesion and Remodelling of the Extra-Cellular Matrix. Cytokine Growth Factor. Rev. 20 (3), 241–249. doi:10.1016/j.cytogfr.2009.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagashima, T., Li, Q., Clementi, C., Lydon, J. P., DeMayo, F. J., and Matzuk, M. M. (2013). BMPR2 Is Required for Postimplantation Uterine Function and Pregnancy Maintenance. J. Clin. Invest. 123 (6), 2539–2550. doi:10.1172/JCI65710

CrossRef Full Text | Google Scholar

Nikuei, P., Malekzadeh, K., Rajaei, M., Nejatizadeh, A., and Ghasemi, N. (2015). The Imbalance in Expression of Angiogenic and Anti-Angiogenic Factors as Candidate Predictive Biomarker in Preeclampsia. Iran J. Reprod. Med. 13 (5), 251–262.

Google Scholar

Okae, H., Toh, H., Sato, T., Hiura, H., Takahashi, S., Shirane, K., et al. (2018). Derivation of Human Trophoblast Stem Cells. Cell Stem Cell. 22 (1), 50–63. e56. doi:10.1016/j.stem.2017.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ono, T., Sekino-Suzuki, N., Kikkawa, Y., Yonekawa, H., and Kawashima, S. (2003). Alivin 1, a Novel Neuronal Activity-dependent Gene, Inhibits Apoptosis and Promotes Survival of Cerebellar Granule Neurons. J. Neurosci. 23 (13), 5887–5896. doi:10.1523/jneurosci.23-13-05887.2003

CrossRef Full Text | Google Scholar

Paria, B. C., Ma, W.-g., Tan, J., Raja, S., Das, S. K., Dey, S. K., et al. (2001). Cellular and Molecular Responses of the Uterus to Embryo Implantation Can Be Elicited by Locally Applied Growth Factors. Proc. Natl. Acad. Sci. 98 (3), 1047–1052. doi:10.1073/pnas.98.3.1047

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, H., Lee, S., Shrestha, P., Kim, J., Park, J. A., Ko, Y., et al. (2015). AMIGO2, a Novel Membrane Anchor of PDK1, Controls Cell Survival and Angiogenesis via Akt Activation. J. Cell Biol. 211 (3), 619–637. doi:10.1083/jcb.201503113

PubMed Abstract | CrossRef Full Text | Google Scholar

Peraçoli, M. T. S., Menegon, F. T. F., Borges, V. T. M., de Araújo Costa, R. A., Thomazini-Santos, I. A., and Peraçoli, J. C. (2008). Platelet Aggregation and TGF-Beta1 Plasma Levels in Pregnant Women With Preeclampsia. J. Reprod. Immunol. 79 (1), 79–84. doi:10.1016/j.jri.2008.08.001

CrossRef Full Text | Google Scholar

Perucci, L. O., Gomes, K. B., Freitas, L. G., Godoi, L. C., Alpoim, P. N., Pinheiro, M. B., et al. (2014). Soluble Endoglin, Transforming Growth Factor-Beta 1 and Soluble Tumor Necrosis Factor Alpha Receptors in Different Clinical Manifestations of Preeclampsia. PLoS One. 9 (5), e97632. doi:10.1371/journal.pone.0097632

PubMed Abstract | CrossRef Full Text | Google Scholar

Pijnenborg, R., Dixon, G., Robertson, W. B., and Brosens, I. (1980). Trophoblastic Invasion of Human Decidua From 8 to 18 Weeks of Pregnancy. Placenta. 1 (1), 3–19. doi:10.1016/s0143-4004(80)80012-9

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.

Google Scholar

Rabenau, K. E., O'Toole, J. M., Bassi, R., Kotanides, H., Witte, L., Ludwig, D. L., et al. (2004). DEGA/AMIGO-2, a Leucine-Rich Repeat Family Member, Differentially Expressed in Human Gastric Adenocarcinoma: Effects on Ploidy, Chromosomal Stability, Cell Adhesion/Migration and Tumorigenicity. Oncogene. 23 (29), 5056–5067. doi:10.1038/sj.onc.1207681

PubMed Abstract | CrossRef Full Text | Google Scholar

Raymond, D., and Peterson, E. (2011). A Critical Review of Early-Onset and Late-Onset Preeclampsia. Obstet. Gynecol. Surv. 66 (8), 497–506. doi:10.1097/OGX.0b013e3182331028

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaarawy, M., Meleigy, M. E., and Rasheed, K. (2001). Maternal Serum Transforming Growth Factor Beta-2 in Preeclampsia and Eclampsia, a Potential Biomarker for the Assessment of Disease Severity and Fetal Outcome. J. Soc. Gynecol. Investig. 8 (1), 27–31. doi:10.1177/107155760100800105

CrossRef Full Text | Google Scholar

Silver, H. M., Lambert-Messerlian, G. M., Reis, F. M., Diblasio, A. M., Petraglia, F., and Canick, J. A. (2002). Mechanism of Increased Maternal Serum Total Activin a and Inhibin a in Preeclampsia. J. Soc. Gynecol. Investig. 9 (5), 308–312. doi:10.1016/s1071-5576(02)00165-x10.1177/107155760200900509

CrossRef Full Text | Google Scholar

Soncin, F., Khater, M., To, C., Pizzo, D., Farah, O., Wakeland, A., et al. (2018). Comparative Analysis of Mouse and Human Placentae Across Gestation Reveals Species-Specific Regulators of Placental Development. Development. 145 (2), dev156273. doi:10.1242/dev.156273

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, R., and Shimodaira, H. (2006). Pvclust: an R Package for Assessing the Uncertainty in Hierarchical Clustering. Bioinformatics. 22 (12), 1540–1542. doi:10.1093/bioinformatics/btl117

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential Gene and Transcript Expression Analysis of RNA-Seq Experiments With TopHat and Cufflinks. Nat. Protoc. 7 (3), 562–578. doi:10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Vento-Tormo, R., Efremova, M., Botting, R. A., Turco, M. Y., Vento-Tormo, M., Meyer, K. B., et al. (2018). Single-cell Reconstruction of the Early Maternal-Fetal Interface in Humans. Nature. 563 (7731), 347–353. doi:10.1038/s41586-018-0698-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vicovac, L., Jones, C. J., and Aplin, J. D. (1995). Trophoblast Differentiation During Formation of Anchoring Villi in a Model of the Early Human Placenta In Vitro. Placenta. 16 (1), 41–56. doi:10.1016/0143-4004(95)90080-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F., Flanagan, J., Su, N., Wang, L.-C., Bui, S., Nielson, A., et al. (2012). RNAscope. J. Mol. Diagn. 14 (1), 22–29. doi:10.1016/j.jmoldx.2011.08.002

CrossRef Full Text | Google Scholar

Yair, D., Eshed-Englender, T., Kupferminc, M. J., Geva, E., Frenkel, J., and Sherman, D. (2001). Serum Levels of Inhibin B, unlike Inhibin A and Activin A, Are Not Altered in Women With Preeclampsia. Am. J. Reprod. Immunol. 45 (3), 180–187. doi:10.1111/j.8755-8920.2001.450310.x

CrossRef Full Text | Google Scholar

You, J., Wang, W., Chang, H.-M., Yi, Y., Zhao, H., Zhu, H., et al. (2021). The BMP2 Signaling Axis Promotes Invasive Differentiation of Human Trophoblasts. Front. Cell Dev. Biol. 9, 607332. doi:10.3389/fcell.2021.607332

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118

CrossRef Full Text | Google Scholar

Zhao, H.-J., Chang, H.-M., Zhu, H., Klausen, C., Li, Y., and Leung, P. C. K. (2018a). Bone Morphogenetic Protein 2 Promotes Human Trophoblast Cell Invasion by Inducing Activin A Production. Endocrinology. 159 (7), 2815–2825. doi:10.1210/en.2018-00301

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H.-J., Klausen, C., Li, Y., Zhu, H., Wang, Y.-L., and Leung, P. C. K. (2018b). Bone Morphogenetic Protein 2 Promotes Human Trophoblast Cell Invasion by Upregulating N-Cadherin via Non-Canonical SMAD2/3 Signaling. Cell Death Dis. 9 (2), 174. doi:10.1038/s41419-017-0230-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H.-J., Chang, H.-M., Klausen, C., Zhu, H., Li, Y., and Leung, P. C. K. (2020a). Bone Morphogenetic Protein 2 Induces the Activation of WNT/β-Catenin Signaling and Human Trophoblast Invasion Through Up-Regulating BAMBI. Cell Signal. 67, 109489. doi:10.1016/j.cellsig.2019.109489

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H. J., Klausen, C., Zhu, H., Chang, H. M., Li, Y., and Leung, P. C. K. (2020b). Bone Morphogenetic Protein 2 Promotes Human Trophoblast Cell Invasion and Endothelial‐Like Tube Formation Through ID1‐mediated Upregulation of IGF Binding Protein‐3. FASEB j. 34, 3151–3164. doi:10.1096/fj.201902168RR

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J., Chen, G., Zhu, S., Li, S., Wen, Z., Bin Li, L., et al. (2016). Identification of Tissue-Specific Protein-Coding and Noncoding Transcripts Across 14 Human Tissues Using RNA-Seq. Sci. Rep. 6, 28400. doi:10.1038/srep28400

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: TGF-β superfamily, early-onset preeclampsia, trophoblast, mRNA-seq, BMP2 (bone morphogenetic protein 2), AMIGO

Citation: Yi Y, Zhu H, Klausen C, Chang H-M, Inkster AM, Terry J and Leung PCK (2021) Dysregulated BMP2 in the Placenta May Contribute to Early-Onset Preeclampsia by Regulating Human Trophoblast Expression of Extracellular Matrix and Adhesion Molecules. Front. Cell Dev. Biol. 9:768669. doi: 10.3389/fcell.2021.768669

Received: 01 September 2021; Accepted: 24 November 2021;
Published: 14 December 2021.

Edited by:

Denny Sakkas, Boston IVF, United States

Reviewed by:

Jaime Agustín Gutiérrez, San Sebastián University, Chile
Shawn L. Chavez, Oregon Health and Science University, United States

Copyright © 2021 Yi, Zhu, Klausen, Chang, Inkster, Terry and Leung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Peter C. K. Leung, peter.leung@ubc.ca