Introduction {#Sec1}
============

Mitogen-activated protein kinases (MAPKs) are signal transducing kinases that contribute to cell differentiation, survival and proliferation^[@CR1],[@CR2]^. MAPKs are a part of ubiquitous signaling cascades involving multiple intermediates associated with phospho-transfer including Raf-1 kinase, MEK kinase (MEKK)1, MAPK-ERK kinases (MEK) 1 and 2, and extracellular signal regulated kinases (ERK) 1 and 2 (also referred to as MAPK3 and MAPK1, respectively)^[@CR3],[@CR4]^. ERK1 and ERK2 have been long known to function integrally during embryonic development and in adult animals in a number of physiologically and pathophysiologically relevant conditions^[@CR5],[@CR6]^. During development, ERK1 and ERK2 have divergent functions since ERK1 null mice are grossly viable and fertile while ERK2 null mice are embryonic lethal^[@CR7]--[@CR9]^. Embryonic lethality of the ERK2 null mice appears to be strain-specific and/or related to the design of the knockout allele; some having reported evidence of failed early mesodermal differentiation while others reporting mis-regulation of placental development^[@CR10],[@CR11]^. Both instances of ERK2 loss lead to early embryo mortality where ERK1 does not appear to be compensatory. In adult animals, ERK1 and 2 appear largely redundant in coordinating cellular responses to growth factors, peptide hormones and other ligands via a wide array of receptor subclasses controlling proliferation, differentiation and survival depending on physiological and cellular context.

Since ERK2 null animals are embryonic lethal, a conditional approach to examining the specific role of ERK signaling in discrete cell types *in vivo* is necessary. In past studies, we reported a pituitary-specific conditional ERK double knockout (ERKdko) using the ERK1 null and ERK2 floxed alleles combined with a Cre recombinase driver regulated by the gene promoter for the α subunit to the glycoprotein hormones (αGSU)^[@CR3],[@CR4]^ and the gonadotrope-specific Cre driver (GRIC)^[@CR12],[@CR13]^. These studies were instrumental in understanding the role and requirement for GnRH-induced ERK signaling to a cohort of immediate early response genes regulating gonadotropin biosynthesis in pituitary gonadotropes. The resultant phenotype was female-specific anovulatory infertility due to a loss of luteinizing hormone (LH) β subunit expression resulting in an inability to mount a preovulatory surge of LH. Other gonadotrope-specific genes (αGSU, follicle stimulating hormone (FSH) β and the GnRHR) are also impacted when ERKs are specifically deleted in this cell type^[@CR3],[@CR12]^. Further, we examined the specific role of Raf-1 kinase in the regulation of GnRH-induced ERK activation finding that Raf-1 was dispensable for activation of this pathway. Similar studies of ERK double knockout conditionally in GnRH neurons did not impact fertility^[@CR14]^ while loss of ERKs in the granulosa cell compartment revealed a requirement for ERK signaling in the ovarian follicle and loss of LH-induced ERK activation in granulosa cells resulted in infertility^[@CR15]^. Thus, ERK signaling across the reproductive axis plays a key role in the endocrine--mediated control of fertility in female mice.

Given that in the GRIC model, Cre recombinase is expressed under the control of the endogenous *GnRHR* promoter, we also need to consider expression of GnRHR at other possible sites that could influence the overall impact of ERKdko using this Cre driver. GRIC expression has been identified in the male germ cell lineage complicating interpretation of the discrete role of ERK loss in the pituitary versus the testis. We recently described the role of ERK loss within the gonadotrope and the spermatogonial lineage, effectively parsing out the impact of ERK loss at both levels of the reproductive axis during aging^[@CR12]^. Based upon this observation of non-gonadotrope Cre expression, we were particularly vigilant in assessing GRIC-mediated gene recombination in other tissues, such as the placenta.

In the present study, we again make use of conditional deletion of ERK1 and 2 using the GRIC driver to examine the hypothesis that these hypogonadotropic-hypogonadal animals are capable of responding to exogenous gonadotropins, ovulate and establish pregnancy. A compelling part of this experimental question was whether the ERKdko females could support *corpora lutea* function following induced ovulation and during early pregnancy given the limited biosynthesis and secretion of LH in the GRIC ERKdko mouse model. We show that female mice harboring ERK gene deletions directed by the GRIC Cre driver are capable of establishing and carrying pregnancy to term. Further, Cre-mediated recombination was observed within the maternal decidua and placenta and while presently speculative, loss of ERKs within these compartments may be associated with marked abnormalities within the implantation site and a prolonged gestational phenotype and complete neonatal lethality.

Results {#Sec2}
=======

ERKdko animals can ovulate in response to exogenous gonadotropin stimulation {#Sec3}
----------------------------------------------------------------------------

We mated the ERK2^fl/fl^, ERK1^−/−^ mice with GRIC animals to improve specificity of ERK ablation to the pituitary gonadotropes compared with previous use of an αGSU-Cre driver^[@CR3],[@CR13]^. In the present studies, these animals were designated ERKdko animals while control animals were ERK2^fl/fl^, ERK1^−/−^, but Cre negative. Consistent with the αGSU model system^[@CR3]^, female ERKdko GRIC mice were infertile and anovulatory, with blunted gonadotropin secretion in response to GnRH following ovariectomy consistent with previous studies (not shown). To better understand the gonadal response to exogenous gonadotropin stimulation in this model, female control and ERKdko animals were treated with either pregnant mare serum gonadotropin (PMSG) or a standard superovulation protocol of PMSG followed by administration of human chorionic gonadotropin (hCG) to induce ovulation^[@CR16]^. Three days later, control and ERKdko females were euthanized and ovaries were collected and the presence of *corpora lutea* (CL) determined histologically. As expected, addition of hCG in the superovulation paradigm resulted in the presence of CLs in both control and ERKdko females compared with PMSG alone (Fig. [1](#Fig1){ref-type="fig"}). These data support the conclusion that the ERKdko females are fully capable of ovulation when provided exogenous gonadotropin stimulation.Figure 1ERKdko females respond to exogenous gonadotropin stimulation leading to ovulation. Photomicrographs of ovaries from Control and ERKdko female mice following administration of PMSG or PMSG and hCG. Magnification bar is shown and arrows identify some of the *corpora lutea*.

ERKdko animals can maintain pregnancy, but show prolonged gestation and parturition {#Sec4}
-----------------------------------------------------------------------------------

We next investigated whether superovulated ERKdko females could establish and maintain a pregnancy to term with viable pups. Control and ERKdko females were treated with PMSG/hCG, paired with control fertile males and were observed for copulatory plugs. In this experimental paradigm, all control and \~70% of ERKdko animals exhibited copulatory plugs. Control and ERKdko females with a copulatory plug were allowed to proceed through gestation without any human interference until near expected term. Beginning in the morning of e18.5, females were observed visually every 12 hours and assessed for behavioral signs of parturition (hiding, contractions, and visible distress) or the presence of pups. Control animals gave birth normally on gestational day 20.3; however, ERKdko females did not give birth until gestational day 23.9 (p \< 0.05; Table [1](#Tab1){ref-type="table"}). In control females, parturition was complete within a 12 hour period; while ERKdko females showed continuous behavioral signs of labor for 2.3 days (p \< 0.05; Table [1](#Tab1){ref-type="table"}) suggesting the possibility of abnormal parturition or dystocia. Finally, control animals had an average of 9 pups/litter; while litter size in the ERKdko females was 2.4 pups/litter (p \< 0.05; Table [1](#Tab1){ref-type="table"}). Pups from ERKdko females were stillborn or died soon after birth; on postnatal day 3, the percentage of live pups from the ERKdko females was zero compared with 85.7% in the control females (p \< 0.05; Table [1](#Tab1){ref-type="table"}).Table 1Reproductive characteristics of control and GRIC ERKdko female mice.Litter size\
(\# of pups)% live pups\
on PN3Gestation\
length (d)Parturition\
length (d)Control9.0 ± 1.0^a^(n = 3)85.7 ± 11.7^a^(n = 3)20.3 ± 0.25^a^(n = 5)1.0 ± 0.0^a^(n = 4)ERK DKO2.4 ± 1.3^b^(n = 7)0^b^(n = 7)23.9 ± 0.55^b^(n = 7)2.3 ± 0.42^b^(n = 4)Data are presented as mean ± standard error of the mean and sample size shown in parentheses. Differing superscripts (a, b) within columns are different (p \< 0.05).

GnRHR is preferentially expressed in the placental decidua in early gestation {#Sec5}
-----------------------------------------------------------------------------

The literature suggests healthy viable pups are routinely born to wildtype females after superovulation protocols suggesting administration of gonadotropins to induce ovulation was not likely a cause of prolonged gestation and neonatal demise in our model system^[@CR16]--[@CR18]^. Moreover, while ERKdko animals have been shown to display defects in gonadotropin production, no evidence in the literature directly links hypogonadotropism to delayed parturition and periparturient mortality^[@CR19]--[@CR22]^. We considered that the abnormal pregnancies observed may be a consequence of Cre-mediated recombination and subsequent loss of ERK signaling at the maternal/fetal interface, as GnRH and GnRHR have been detected in the placenta and decidua of several species^[@CR23]--[@CR28]^.

We used qRT-PCR to detect GnRHR mRNA in implantation sites at e10.5 and e12.5 mice and compared this directly with female pituitaries at proestrus/estrus (Fig. [2A](#Fig2){ref-type="fig"}). Levels of GnRHR mRNA in these tissues were consistent with levels of GnRHR found in adult female pituitaries at proestrus/estrus. To further validate these observations, we used the Rosa26 reporter mouse line mated to GRIC^+^ males. In this model, Cre-mediated recombination results in expression of a β-galactosidase (β-gal)-green fluorescent protein (GFP) fusion protein. *In situ* β-gal staining of whole mouse pituitaries and implantation sites confirmed recombination and β-gal activity within both tissues in ROSA^GRIC+^ mice. The highest density of β-gal activity was found at the periphery of the implantation site, indicating potential enrichment within the maternal decidua (Fig. [2B](#Fig2){ref-type="fig"}). Using placental tissues from ROSA26^GRIC−^ and ROSA26^GRIC+^ animals obtained at e12.5, we then carried out an *in vitro* β-galactosidase activity assay to quantify relative β-gal expression in a cross section of tissues (Fig. [2C](#Fig2){ref-type="fig"}). We compared ROSA26^GRIC−^ and ROSA26^GRIC+^ placental tissues along with ROSA26^GRIC+^ pituitary, hypothalamus, ovary, pancreas, liver and muscle from the same animals. Placenta lysates from ROSA26^GRIC+^ animals revealed a fourfold increase (p \< 0.05) of β-gal activity compared with the ROSA26^GRIC−^. To further parse out localization of GnRHR expression within the implantation site, e12.5 placentas were dissected to carefully separate and enrich maternal decidua from the remainder of the placental disk including the chorionic plate and labyrinth. RNA was extracted and qRT-PCR was performed, allowing us to quantify enrichment of the GnRHR mRNA in the decidua versus the remainder of the disk at e12.5. GnRHR mRNA was more abundantly expressed in the maternal decidua compared to the remainder of the placental disk where GnRHR mRNA was lower but clearly detectable (Fig. [2D](#Fig2){ref-type="fig"}).Figure 2GnRH receptor (R) mRNA is detectable in the placental disk at embryonic (e) days 10.5 and 12.5 and the GRIC Cre mediates recombination in maternal decidua and placenta. (**A**) RT-PCR was performed on placental disks from embryonic days 10.5 and 12.5. Levels of GnRHR were comparable with those seen in adult female pituitaries collected at proestrus/estrus. N = 3 placentas/gestational age. In (**B**) Adult ROSA26^GRIC+^ mice expressing a β-galactosidase (β-gal)-GFP fusion protein show β-gal activity in the female pituitary *in situ*. β-gal staining indicates the presence of GRIC Cre activity in placenta from e18.5 ROSA26^GRIC+^ embryos but not placentas from e18.5 ROSA26^GRIC−^ embryos. (**C**) Liquid β-gal activity assay confirmed GRIC Cre-mediated recombination in ROSA26^GRIC+^ pituitaries and placentas. ROSA26^GRIC−^ placental disk and ROSA26^GRIC+^ hypothalamus, ovary, pancreas, liver and muscle were used as controls. (**D**) qRT-PCR showed enrichment of GnRHR mRNA expression in placental decidua compared to whole placental disk. Differing letters (a and b) over individual bars identify differences (p \< 0.05). N = \>4 placentas/genotype.

GnRHR is expressed primarily in decidua at e12.5 and in the labyrinth at e18.5 {#Sec6}
------------------------------------------------------------------------------

Rosa26 animals were time mated with GRIC^+^ males, and euthanized at e12.5 and e18.5. Placentas were collected, snap frozen, serial sectioned and stained for GFP expression. At e12.5, GFP expression appeared to be expressed most abundantly in the decidua (Supplemental Fig. [1A,B](#MOESM1){ref-type="media"}). Punctate perinuclear GFP immunoreactivity was detected in clusters of cell within the decidua. There is modest basal expression of GFP throughout the cytoplasm of the labyrinth, but no areas of punctate expression were observed. At e18.5, punctate expression in the decidua was no longer observed; however, the labyrinth showed increased expression of GFP when compared to e12.5. Importantly, GFP was not detected within the junctional zone at e18.5 demonstrating the specificity of GFP expression within decidua and labyrinth. Negative controls (secondary antibody alone) for the immunofluorescence define the detection limit of the assay.

Cre-mediated recombination results in loss of placental ERK2 {#Sec7}
------------------------------------------------------------

Western blot analyses added additional confirmation of Cre-mediated *ERK* deletion within the decidua/placental disk at E12.5 (Supplemental Fig. [2A,B](#MOESM1){ref-type="media"}). We specifically assayed for the abundance of ERK2 since our genetic model is null for the ERK1 allele. For the remainder of the characterization, we separated the animals into three categories: control/control (control dam, control pup), ERKdko/control (ERKdko dam, control pup), and ERKdko/ERKdko (ERKdko dam, ERKdko pup). We observed \~30% reduction of ERK2 protein levels in decidua/placentas from both ERKdko groups when compared to placentas from control dams at e12.5. This suggests that the genotype of the dam may have a more significant impact on ERK expression than the genotype of the pup, again consistent with recombination ERK2 in decidua or other cells types within the uterus at this gestational age.

Loss of ERK2 in the implantation site occurred coincident with fetal growth defects and gross and histological abnormalities in the placental disk {#Sec8}
--------------------------------------------------------------------------------------------------------------------------------------------------

To examine the impact of ERK loss at the maternal fetal interface, ERKdko and control females were superovulated, mated, euthanized at e18.5 and pups and placental tissues dissected. Anecdotally, all pups from the ERKdko and control dams were noted to be alive at this gestational age and euthanized humanely. ERKdko/control (maternal/fetal) and ERKdko/ERKdko pups displayed reduced fetal weight at e18.5 compared to control/control pups (p \< 0.05; Fig. [3A](#Fig3){ref-type="fig"}). Consistent with earlier studies on ERK2 expression levels, reduced pup weight was primarily associated with loss of ERKs within the decidua (maternal ERKdko) at this time point. Intrauterine growth restriction seen in ERKdko/control or ERKdko/ERKdko pups could not be accounted for by changes in placental weight (Fig. [3B](#Fig3){ref-type="fig"}) indicating the possibility that fetal growth restriction observed was due to loss of relative placental efficiency (mass of fetus/mass of placental disk), not strictly placental size.Figure 3Pups from GRIC ERKdko females have reduced fetal weight at e18.5. Fetal (**A**) and placental (**B**) weights are shown at e18.5 from control/control, GRIC ERKdko/control and GRIC ERKdko/GRIC ERKdko pups and placentas, respectively. N = \>4 placentas/genotype. Differing letters (a and b) over individual bars identify differences (p \< 0.05).

Histological assessment of these implantation sites at e18.5 is shown in Fig. [4](#Fig4){ref-type="fig"}. For these studies, placentas were dissected free from the uterine wall, fixed in formalin, then sectioned and stained with isolectin B4 to identify the glycocalyx associated with blood vessels and lightly counter stained with hematoxylin/eosin. This approach clearly demarcates labyrinth, junctional and decidual compartments within the disk. Pups from ERKdko dams had abnormal placental architecture regardless of fetal genotype at e18.5 (Fig. [4A](#Fig4){ref-type="fig"}). Figure [4B](#Fig4){ref-type="fig"} depicts the distribution of the area of the labyrinth, junctional zone and decidua which varied with genotype. Using Aperio image analysis software, control/control (maternal/fetal) placentas were characterized by approximately 49% labyrinth, 26% junctional zone, and 17% decidua (with the remainder of the area attributed to chorionic plate). ERKdko/ERKdko placentas showed a modestly enlarged labyrinth and a smaller decidual area (Fig. [4B](#Fig4){ref-type="fig"}). ERKdko/control placentas showed an intermediate phenotype that was not significantly different from either of the other two genotypes. Pathological findings indicated that placentas from ERKdko dams and control fetuses appeared histologically disorganized compared to control placentas, displaying variability in decidual thickness. There were fingerlike projections of the junctional zone protruding into the labyrinth. GRIC-mediated recombination of ERK2 was associated with the appearance of large acellular spaces within the junctional zone. These findings were further exacerbated in the ERKdko/ERKdko placentas, with more marked histological changes, especially in the decidua and junctional zone, consistent with greater disorganization of the normal placental architecture seen in the control/control placentas.Figure 4Histological assessment of placental disks in the GRIC ERKdko model reveals abnormal architecture. (**A**) Representative sections from control/control, GRIC ERKdko/control and GRIC ERKdko/GRIC ERKdko placentas at e18.5 were stained for isolectin B4 to mark the glycocalyx associated with blood vessel. Panels on the left are low magnification while panels on the right are higher magnification. The decidua (Dec), junctional zone (JZ) and labyrinth (Lab) are indicated. The magnification bar is shown for the higher magnification sections on the right. (**B**) Quantitation of decidual, junctional zone and labyrinth are reported as % of total placental disk area in control (Con) and GRIC ERKdko (DKO) animals. N = 3 placentas/genotype. Differing letters (a and b) over individual bars identify differences (p \< 0.05).

Histological changes in the placenta occur coincident with altered placental isolectin B4 binding {#Sec9}
-------------------------------------------------------------------------------------------------

Since endothelial cell dysfunction has been associated with placental pathology and disease states^[@CR29]--[@CR32]^, we examined isolectin B4 binding to the placental glycocalyx associated with blood vessels as a proxy for endothelial cells including in decidua to assess isolectin B4 positive cell density and level of expression. Staining tissue sections to detect isolectin B4 binding highlights the extracellular matrix of endothelial cells (among other cell types), allowing visualization of highly vascularized regions within the placenta^[@CR33]^. At E18.5, both the number of isolectin B4-positive cells and the intensity of staining were blunted (p \< 0.05) in ERKdko/ERKdko placentas in both the decidua and junctional zone compared to control (Fig. [5A,B](#Fig5){ref-type="fig"}). The ERKdko/control placentas were highly variable with regard to isolectin B4 binding resulting in no statistical differences when compared to the other genotypes (not shown).Figure 5Assessment of isolectin B4 binding in GRIC ERKdko mice implantation sites. The number of isolectin B4^+^ cells (**A**) and staining intensity (**B**) are reported for control/control and GRIC ERKdko dam/GRIC ERKdko groups at e18.5 as number of cells or relative staining intensity (pixel intensity)/10 x magnification field. N = \>4 placentas/genotype. Differing letters (a and b) over individual bars identify differences (p \< 0.05).

GRIC-mediated recombination of ERK is associated with decreased placental parturition-associated apoptosis {#Sec10}
----------------------------------------------------------------------------------------------------------

We examined two possible mechanisms that may account for the prolonged gestation/delayed parturition in this model. First, the prolonged gestation interval observed in the ERKdko females was reminiscent of phenotype described for the cyclooxygenase 1 (PGHS1/COX1) knockout animal^[@CR34]--[@CR39]^. We speculated that ERK signaling within the decidua and labyrinth may be necessary for the proper expression pattern of COX1 at or near term. As an initial assessment, we used immunohistochemistry to determine if marked differences in COX1 protein expression were detectable in placental sections from control and ERKdko female mice at e18.5 (Supplemental Fig. [3](#MOESM1){ref-type="media"}). These studies demonstrated that COX1 immunoreactivity was readily detectable in control and ERKdko pregnancies suggesting that differential expression or deficiencies in COX1 was not likely causal in the context of our model of prolonged gestation.

As a second approach, we examined the progression of placental apoptosis within the maternal/fetal interface at the time of onset of parturition (Fig. [6A](#Fig6){ref-type="fig"}). Sections from control/control placentas showed focal areas of apoptosis at the interface of the decidua and uterus, putatively allowing for initial separation of the maternal/placental unit at the time of parturition^[@CR40]--[@CR43]^. The ERKdko/ERKdko placentas had markedly reduced TUNEL-positive staining around the edges of the decidua at e18.5, indicating a loss or delay of physiologic apoptosis associated with parturition (Fig. [6](#Fig6){ref-type="fig"}). Control/control placentas showed positive TUNEL staining in \~17% of the decidual area in this region of the placental disk, while ERKdko/ERKdko placentas showed only \~8% percent of the decidual area was positive for TUNEL staining at e18.5 (p \< 0.05, Fig. [6B](#Fig6){ref-type="fig"}). TUNEL staining in placentas from ERKdko/control placentas showed an intermediate phenotype and were not significantly different than either of the other two groups (data not shown). The loss of pro-apoptotic tone in ERKdko/ERKdko placentas at E18.5 occurred coincident with prolonged gestation and delayed parturition suggesting that they may be associated.Figure 6Assessment of TUNEL^+^ cells within the maternal-fetal interface at near term in GRIC ERKdko females. (**A**) Representative images of placentas from control/control and GRIC ERKdko/GRIC ERKdko placentas at e18.5 stained using TUNEL as an indicator of apoptosis (red staining). The magnification bar is shown. (**B**) TUNEL^+^ cells/μm^2^ is reported for control/control and GRIC ERKdko/GRIC ERKdko placentas at e18.5. N = 5 placentas/genotype. Differing letters (a and b) over individual bars identify differences (p \< 0.05).

Discussion {#Sec11}
==========

While loss of ERK signaling within the gonadotrope has been shown to be detrimental to the function of the hypothalamic-pituitary-gonadal axis as well as the placenta, this is the first report of a novel model of placental ERK gene deletion associated with expression of a GnRHR-Cre driver^[@CR3],[@CR12],[@CR13]^. In the present studies, the anovulatory phenotype due to loss of pituitary LH production in the GRIC-ERKdko model was overcome through exogenous gonadotropin administration, inducing a subset of ERKdko animals to ovulate and become pregnant. The resulting pregnancy in ERKdko dams was abnormal, characterized by altered placental architecture, prolonged gestation and parturition, and 100% neonatal mortality by postnatal day 3. Based upon the reporter mouse studies, GnRHR-regulated Cre expression and recombination occurred in decidua in midgestation and within the placental labyrinth during later stages of pregnancy. The altered placental architecture in ERKdko/ERKdko placentas occurred coincident with a significant reduction in physiologic apoptosis at e18.5 putatively necessary for placental separation prior to parturition, consistent with the prolonged gestation phenotype observed. Given the complexity of the expression for Cre recombinase in the GRIC model (gonadotropes and placental compartment), we cannot specifically rule out the possibility that loss of gonadotropins contributed to the placental phenotype. Moreover, an important caveat to our studies is that we have not yet demonstrated that conditional ERK deletion specifically within the implantation site (decidua and or labyrinth) is causally or directly related to the abnormal placental architecture or the delayed parturition phenotype. Our understanding of the specific impact of ERK signaling deletion within the implantation site awaits development of new mouse models. Our studies do support the hypothesis that conditional deletion of ERKs within the maternal fetal interface is likely to be associated with the phenotype observed in the present studies. Collectively, these studies shed new light on the potential role of ERK signaling within the reproductive axis including the implantation site in the mouse.

Initial analyses of human endometrium and decidua from early pregnancy revealed the presence of mRNA expression of GnRH but not the GnRHR^[@CR44]^; however, subsequent studies also examining human tissues reported that endometrial and decidual tissues were responsive the GnRH analogs as measured by induction of phosphoinositol turnover^[@CR45]^. These observations were consist with findings in the HTR-8/SVneo human extravillous trophoblast cell line where GnRH stimulated expression of several proangiogenic chemokines^[@CR46]^ and in primary human trophoblasts and decidual stromal cells from term placentas where GnRH stimulated the production of chorionic gonadotropin^[@CR23]^. Some of the most comprehensive studies of the action of GnRH on placental cell behavior come from Peter Leung's group. Using primary culture of first trimester trophoblast cells or the extravillous trophoblast model HTR-8/SVneo, this group demonstrated GnRH-induced RUNX2 expression leading to upregulation of RUNX2-dependent matrix metalloproteinases which may modulate the invasive properties of these trophoblast populations^[@CR47]^. The potential impact of GnRH/GnRHR signaling and the use of GnRHR antagonists on human endometrial stromal cells is controversial since one group suggested that these antagonists did not impact the behavior of the stromal cells^[@CR48]^; while another group suggested that antagonism of the GnRHR resulted in elevated apoptotic tone in stromal cells^[@CR49]^.

The ligands GnRH I and II localize to mononucleate villous trophoblasts, and EVTs, but GnRH I is also found in the multinucleated syncytial trophoblast layer in the chorionic villi in humans during the first trimester of pregnancy^[@CR50]^. Recently, GnRHR has been described in the canine placenta as well. Similar to the human, GnRHR was expressed at higher levels in the uteroplacental tissues compared to areas of the uterus not associated with the zonary placenta (i.e., between implantation sites in this carnivore litter-bearing species). In the canine placenta, GnRHR was found in fetal trophoblasts and maternal decidual cells, and at the surface and glandular epithelial cells in the uterus^[@CR25]^. These observations are entirely consistent with expression of Cre recombinase in the GRIC model by the endogenous *GnRHR* promoter in the mouse. In the present studies, it is important to note that we do not provide evidence for a direct and specific role for GnRH-directed GnRHR signaling to the ERK pathway within the maternal-fetal interface in mediating the potential effects of ERK deletion. Rather, we argue that this model examines the putative loss of ERK signaling within the maternal-fetal interface leading to a prolonged gestation, where regulation of the ERK signaling pathway within the placenta may be a function of numerous growth peptides and hormones (potentially including GnRH).

Other studies have shown that loss of ERK2 in null animals causes catastrophic defects early in mouse development^[@CR10]^ or at later developmental stages in the mouse placenta leading to early fetal death^[@CR8],[@CR11]^. In the former studies, specific loss of ERK2 (deletion of exon 2) resulted in embryonic demise at e6.5. Embryonic death occurred coincident with a lack of mesoderm formation and a marked increase in apoptotic tone within the developing embryo. These observations are difficult to reconcile with delayed or absence apoptotic tone in the GRIC ERKdko model at term, but may reflect the timing of ERK loss in the individual models. In the later studies, ERK2 deletion was mediated by disruption of exon 3^[@CR8]^ or introduction of a neomycin cassette between exons 2 and 3 disrupting splicing of the mature ERK2 mRNA^[@CR11]^. In both cases, placental competence was compromised resulting in embryonic death later in development (e8.5 to 12.5). All of these studies point to the specific importance of ERK2 during development since in each experimental paradigm, ERK1 did not compensate for the loss of ERK2. Thus it may not be too surprising that loss of ERK signaling within the implantation site in the GRIC model was detrimental to normal pregnancy. However, our studies did not phenocopy the ERK2 null mouse (early fetal demise) suggesting that the timing of ERK2 loss within the placenta and/or the cell types impacted by such loss help to explain the apparent differences in phenotypes comparing the ERK 2 null mouse with the conditional GRIC ERKdko in the present studies.

Loss of placental ERK2 resulted in reduced intensity of isolectin B4 binding within the blood vessel associated glycocalyx, and placental abnormalities that may have contributed to reduced fetal weight in the current studies. We speculate that expression of GnRHR (and Cre) and subsequent loss of ERKs during the formation/expansion of the decidua prior to e12.5 is a likely cause of the changes in placental histology seen at e18.5^[@CR51]^. Either the change in architecture or loss of ERK2 proteins could impact trophoblast invasion and subsequent maternal spiral artery remodeling, reducing the number of vessels observed in the decidua^[@CR52]^. Additionally, matrix metalloproteases, which are regulated by GnRH/GnRHR in the placenta, are involved in maternal spiral artery remodeling, leading to another hypothesis for reduced decidual vascularization; however, the specific role of ERK signaling in this situation remains to be elucidated^[@CR52],[@CR53]^. Although speculative, the reduction in the number of isolectin B4-positive cells and the magnitude of isolectin B4 binding may reflect a reduction in endothelial cell populations or functional competence of the microenvironment within the decidua, which may lead to reduced placental transfer between mother and fetus and impact fetal size. Loss of normal placental architecture, general disorganization of the decidua and junctional zone may also contribute to loss of efficient nutrient exchange between the fetus and dam. The direct cause(s) of the alteration in placenta architecture is currently unclear. Interestingly, the most obvious changes are attributable to the finger-like projections of the junctional zone into the labyrinth. Recall that our Rosa26 reporter mice and subsequent immunofluorescence studies revealed that cells within the junctional zone did not undergo apparent genetic recombination. This may support the conclusion that a secreted product of either the maternal decidua or the labyrinth (or both) may be impacting the organization of the cells within the junctional zone. Interestingly, similar junctional zone histological abnormalities have been demonstrated in studies of implantations sites from various forms of assisted reproductive technologies from the Bartolomei group^[@CR54]^. In these studies, this abnormal placental architecture was associated with disruptions in genomic imprinting. In other studies, genetic disruption of Plac1 (a maternally derived X chromosome-linked gene expressed in trophoblasts) caused similar junctional zone architectural abnormalities^[@CR55]^. At present, the potential relationship between ERK signaling within the implantation site and aspects of assisted reproduction, genomic imprinting or Plac1 expression is unknown but may provide some important insights into possible mechanisms.

Prolonged gestation and parturition in this model may reflect the local impact of loss of ERKs within the maternal-fetal interface. Multiple factors initiate parturition, including fetal cortisol, maternal drop in progesterone, oxytocin, prostaglandin F2α and other endocrine factors^[@CR56]--[@CR59]^. Our initial examination that ERK signaling may have been required for normal expression of COX1 was a logical hypothesis given the prolonged gestation phenotype of the COX1 deficient mouse^[@CR34],[@CR35],[@CR37]--[@CR39]^. However, we did not detect appreciable differences in COX1 protein expression using immunohistochemistry. Alternatively, distinct changes in apoptotic tone at the time of normal parturition were observed suggesting that local loss of ERK signaling may be critical for establishing this aspect of initiation of parturition. Other alternatives include any one of a number of paracrine and/or endocrine signaling events described above. Progesterone withdrawal is a contributing factor to initiation of parturition and the possibility exists that the baseline reduction in progesterone in our model system does not occur in a timely manner or at a level to initiate parturition^[@CR60]^. The loss of ERK signaling in GnRHR-expressing cells within the maternal-fetal interface may blunt or abrogate the release of PGF2α thus contributing to the delay in onset of parturition^[@CR61]^. Alteration in oxytocin release is likely not related to the delayed parturition since this peptide is not implicated as an important factor in the onset of parturition in mice^[@CR62]^. Future studies will examine these possibilities.

This model holds interesting potential as a model of human disease. GnRHR has been shown to be important in placental development and function; it is even expressed at the implantation sites of ectopic pregnancy^[@CR63]^. Although multiple mouse models have been used to recapitulate idiopathic hypogonadal hypogonadism (IHH), none have focused on loss of function in cells expressing GnRHR and how ERK signaling may be functioning in these tissue compartments^[@CR3],[@CR64],[@CR65]^. Up to 40 percent of IHH patients suffer from a GnRHR mutation, which could alter placental function as well^[@CR66]^. Work by Janet Hall's lab provides important evidence that although some women with IHH respond appropriately to gonadotropin supplementation and are able to carry a normal pregnancy, others are not^[@CR67],[@CR68]^. Some of these women suffer from multiple miscarriages, most often early in gestation. Not all patients with these mutations showed early pregnancy losses, indicating a non-genetic component, such as alterations in placental function as a potentially contributing factor^[@CR67],[@CR69],[@CR70]^. Further investigation of this model system may elucidate missing links and help improve assisted reproduction in women with IHH.

Materials and Methods {#Sec12}
=====================

Animals {#Sec13}
-------

Animals were handled in compliance with a protocol approved by the Cornell University Institutional Animal Care and Use Committee. ERK1 null (*ERK1*^−/−^), ERK2 floxed (*ERK2*^*fl/fl*^) and GnRH receptor (GnRHR) IRES Cre (GRIC) mice have been described previously^[@CR3],[@CR12],[@CR13],[@CR64]^. The presence of the internal ribosome entry site (IRES) allows for co-expression of the GnRHR and Cre recombinase regulated by the endogenous *GnRHR* promoter. For experiments involving analyses of Cre activity, ROSA26-GNZ KI mice were purchased from Jackson Laboratory. This reporter strain was crossed with GRIC animals, and designated ROSA26^GRIC**+**^ or ROSA26^GRIC−^. For timed matings, mice were paired and checked for copulatory plugs daily. Pairs were separated after a copulatory plug was observed (indicating embryonic day (e)0.5). On e10.5, e12.5 and e18.5, the dams were euthanized, and maternal-fetal implantation sites were carefully dissected. For superovulation studies, females were injected with 100 IU of pregnant mare serum gonadotropin (PMSG) intraperitoneally. Approximately 46--48 hours later, they were injected with 100 IU human chorionic gonadotropin (hCG) intraperitoneally. For assessing *corpora lutea* (CL) and follicle formation, superovulated animals were euthanized, and ovaries were collected for histological examination. For induction of pregnancy and assessment of parturition, superovulated control and ERKdko females were placed with fertile males and checked for copulatory plugs the following morning day. Following evidence of copulation, females were minimally handled until e18.5, at which point they were checked every 12 hours for behavioral signs of initiation of parturition, contractions, distress, hiding behaviors, blood on the cage bedding, and pups.

Genotyping {#Sec14}
----------

Genomic DNA was isolated from tail snips, using an E-Z Tissue DNA Kit (Omega Biotek, Norcross, GA) per the manufacturer's instructions. Routine PCR genotyping was performed on animals as previously described^[@CR71]^. PCR confirmation of *ERK1* knockout, *ERK2* floxed, *GRIC* and *ROSA26* alleles were performed with primers reported in Table [2](#Tab2){ref-type="table"}.Table 2PCR primer sets used in genotyping mice.Primer setSequence (5′ to 3′)Rosa26 Reporter*Rosa26* ForwardTAA GCC TGC CCA GAA GAC TC*Rosa26* ReverseAAA GTC GCT CTG AGT TGT TAT*Rosa26* CommonTCC AGT TCA ACA TCA GCC GCT ACAERK1*Erk1* ForwardAAG GTT AAC ATC CGG TCC AGC A*Erk1*ReverseAAG CAA GGC TAA GCC GTA CCERK2*Erk2* ForwardAGC CAA CAA TCC CAA CCC TG*Erk2* ReverseGGC TGC AAC CAT CTC ACA ATGnRH Receptor*Gnrhr* ForwardGAA CTA CAG CTG AAT CAG TC*Gnrhr* ReverseCTC TAA CAA ACT CTG TAC A*Gnrhr* HomozygousCGG AAT TCA TCG ATC ATA TCA GAT CC

Immunoblotting {#Sec15}
--------------

Placentas were dissected such that a portion was placed in formalin for histological studies, snap frozen for RNA analysis, or homogenized in lysis buffer containing 20 mM Tris-HCl (pH 8.0), 130 mM NaCl, 10% glycerol, 1% Nonidet P-40, 0.1% sodium dodecyl sulfate, 0.5% deoxycholate, 2 mM EDTA, 5 mM sodium vanadate, 0.2 mM phenylmethysulfonylfluoride, and 5 mM benzamidine. Protein concentrations of lysates were determined by Bradford assay. Samples were boiled for 5 minutes in sodium dodecyl sulfate loading buffer (2x contains 100 mM Tris pH 6.8, 4% sodium dodecyl sulfate, 0.2% bromophenol blue and 200 mM dithiothreitol), resolved by SDS-PAGE, and transferred to polyvinylidine difluoride membranes by electroblotting. Membranes were blocked with 5% nonfat dry milk in TBST (10 mM Tris-HCl, pH 7.5; 150 mM NaCl; 0.05% Tween 20) and then incubated with specified antisera (anti-ERK2, anti-β actin and horseradish peroxidase-conjugated secondary antibodies from Santa Cruz Biotech, Dallas, TX)^[@CR3]^. Protein bands were visualized using enhanced chemiluminescence according to manufacturer's instructions (BioRad, Berkeley, CA) and imaged on a ChemiDoc XRS (BioRad, Berkeley, CA) Band densities were analyzed using Image Lab software (BioRad, Berkeley, CA).

RNA isolation and quantitative PCR {#Sec16}
----------------------------------

Placentas were subjected to Trizol (ThermoFisher, Waltham, MA) extraction per the manufacturer's instructions. Reverse transcription of RNA was performed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA) according to the manufacturer's directions. The sequence of the primers used in these amplification studies are reported in Table [3](#Tab3){ref-type="table"}. Amplifications were carried out using a BioRad CFX96 Touch Real-Time OCR Detection System (BioRad, Berkeley, CA). RNA levels were standardized using the internal control Gapdh and calibrated to corresponding transcript levels of a control group by the ΔCt method^[@CR3]^.Table 3Quantitative PCR primer sets for assessing expression of GnRH receptor and GAPDH.Primer setSequence (5′ to 3′)Glyceraldehydes-3-phosphate dehydrogenase*Gapdh* forwardATGTTTGTGATGGGTGTGAA*Gapdh* reverseATGCCAAAGTTGTCATGGATGnRH Receptor*Gnrhr* forwardTGCTCGGCCATCAACAACA*Gnrhr* reverseGGCAGTAGAGAGTAGGAAAAGGA

Histology, immunofluorescence, Immunohistochemistry and β-galactosidase assays {#Sec17}
------------------------------------------------------------------------------

Mice were paired and checked for copulatory plugs daily. Pairs were separated after a copulatory plug was observed (e0.5). On e12.5 and e18.5, the dams were euthanized, and placentas were carefully dissected. For histological examination, tissues were fixed in 10% neutral-buffered formalin, paraffin embedded, serially sectioned at 4 µm, and stained with hematoxylin and eosin, terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL), smooth muscle actin or isolectin B4 using standard techniques as described previously^[@CR72]^. Sections were scanned and digitized using an Aperio Scanscope (Vista CA). Placental area, areas of decidua, junctional zone and labyrinth layers were quantitated using Aperio software. TUNEL was quantified using the Positive Pixel Algorithm on ImageScope (Leica Biosystems, Buffalo Grove, IL).

For immunofluorescence labeling of GFP, placentas were embedded in Tissue-Tek OCT media (Sakura Finetek, Torrence, CA) and maintained at −80 °C in 2-methylbutane for 24 hours. Sections (10 µM) were cut using a cryotome and stored at −80 °C. Frozen sections were fixed with 4% paraformaldehyde, blocked for one hour, and stained with anti-GFP antiserum (Abcam, Cambridge, United Kingdom) at 1:100 dilution overnight at 4 C followed with a FITC goat-anti-rabbit secondary antibody (Vector Laboratories, Burlingame, CA). Sections were cover slipped with SlowFade Gold with DAPI (Life Technologies, Carlsbad, CA) and imaged on an AxioVision fluorescent microscope with Zen software (Zeiss, Oberkochen, Germany). Sections were also scanned and digitized using an Aperio Scanscope (Vista CA). For immunohistochemistry studies, placental sections were blocked using an antibody-specific HRP/DAB detection kit (per the manufacturer's instructions; Abcam, Cambridge, United Kingdom) and stained with the PGHS1/COX1 antibody (Abcam, Cambridge, United Kingdom) at 1:50 dilution overnight. Negative controls included sections treated identically except for addition of the primary antibody. Detection of epitopes was carried out using the kit described above. Sections were cover slipped as described above and imaged and digitized using an Aperio Scanscope (Vista, CA).

For β-galactosidase *in vitro* assays, tissues were fixed in 4% paraformaldehyde/PBS for 1 hour at 4 °C, then rinsed 3 times for 30 minutes each in a rinse buffer (100 mM sodium phosphate (pH 7.3), 2 mM MgCl~2~, 0.01% sodium deoxycholate, 0.02% NP-40 (by volume)). Whole tissues or sections were stained overnight in staining buffer (rinse buffer containing 5 mM potassium ferricyanide, 5 mM potassium ferrocyanide, 1 mg/ml X-gal). Tissues or sections were cleared using methyl salicylate, fixed overnight in 10% formalin, then washed with distilled water twice for 30 minutes. They were then dehydrated by sequential ethanol washes (70%, 95%, twice in 100%), then washed in methyl salicylate until the tissue cleared. For liquid β-Galactosidase Enzyme Assay System with Reporter Lysis Buffer (PROMEGA, Madison, WI), tissues were collected and assayed according to the manufacturer's instructions.

Statistics {#Sec18}
----------

Pairwise comparisons were made by Student's *t*-test. For comparisons of more than two means, a one way ANOVA was performed and differences were determined by a Duncan's test. All data are expressed as means ± standard error of the mean. A *p* value of \<0.05 was considered statistically significant.

Supplementary information
=========================

 {#Sec19}

Dataset 1

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information
=========================

**Supplementary information** accompanies this paper at 10.1038/s41598-019-45997-0.

We thank Dr. Lynn Dong for technical assistance with immunohistochemistry. These studies were supported by NIH/NICHD R01 HD034722 and R21 HD082780. We thank members of the Roberson lab for helpful discussions in the development of this manuscript.

J. Brown, J. Sones and M. Roberson designed the research; J. Brown, J. Sones, C. Angulo, K. Abbott and M. Roberson performed the research; U. Boehm contributed new reagents; A. Miller analyzed pathology data; and J. Brown and M. Roberson wrote the manuscript.

The authors declare no competing interests.
