ERß Regulation of Indian Hedgehog Expression in the First Wave of Ovarian Follicles.

Increased activation of ovarian primordial follicles in Erß knockout (ErßKO) rats becomes evident as early as postnatal day 8.5. To identify the ERß-regulated genes that may control ovarian primordial follicle activation, we analyzed the transcriptome profiles of ErßKO rat ovaries collected on postnatal days 4.5, 6.5, and 8.5. Compared to wildtype ovaries, ErßKO ovaries displayed dramatic downregulation of Indian hedgehog (Ihh) expression. IHH-regulated genes, including Hhip, Gli1, and Ptch1, were also downregulated in ErßKO ovaries. This was associated with a downregulation of steroidogenic enzymes Cyp11a1, Cyp19a1, and Hsd17b1. The expression of Ihh remained very low in ErßKO ovaries despite the high levels of Gdf9 and Bmp15, which are known upregulators of Ihh expression in the granulosa cells of activated ovarian follicles. Strikingly, the downregulation of the Ihh gene in ErßKO ovaries began to disappear on postnatal day 16.5 and recovered on postnatal day 21.5. In rat ovaries, the first wave of primordial follicles is rapidly activated after their formation, whereas the second wave of primordial follicles remains dormant in the ovarian cortex and slowly starts activating after postnatal day 12.5. We localized the expression of Ihh mRNA in postnatal day 8.5 wildtype rat ovaries but not in the age-matched ErßKO ovaries. In postnatal day 21.5 ErßKO rat ovaries, we detected Ihh mRNA mainly in the activated follicles in the ovaries' peripheral regions. Our findings indicate that the expression of Ihh in the granulosa cells of the activated first wave of ovarian follicles depends on ERß.

DOT1L Mediated Gene Repression in Extensively Self-Renewing Erythroblasts.

DOT1L is essential for embryonic hematopoiesis but the precise mechanisms of its action remain unclear. The only recognized function of DOT1L is histone H3 lysine 79 (H3K79) methylation, which has been implicated in both transcriptional activation and repression. We observed that deletion of the mouse Dot1L gene (Dot1L-KO) or selective mutation of its methyltransferase domain (Dot1L-MM) can differentially affect early embryonic erythropoiesis. However, both mutations result in embryonic lethality by mid-gestation and growth of hematopoietic progenitor cells (HPCs) is similarly affected in extensively self-renewing erythroblast (ESRE) cultures established from yolk sac cells. To understand DOT1L-mediated gene regulation and to clarify the role of H3K79 methylation, we analyzed whole transcriptomes of wildtype and Dot1L-mutant ESRE cells. We observed that more than 80% of the differentially expressed genes (DEGs) were upregulated in the mutant ESRE cells either lacking the DOT1L protein or the DOT1L methyltransferase activity. However, approximately 45% of the DEGs were unique to either mutant group, indicating that DOT1L possesses both methyltransferase-dependent and -independent gene regulatory functions. Analyses of Gene Ontology and signaling pathways for the DEGs were consistent, with DEGs that were found to be common or unique to either mutant group. Genes related to proliferation of HPCs were primarily impacted in Dot1L-KO cells, while genes related to HPC development were affected in the Dot1L-MM cells. A subset of genes related to differentiation of HPCs were affected in both mutant groups of ESREs. Our findings suggest that DOT1L primarily acts to repress gene expression in HPCs, and this function can be independent of its methyltransferase activity.

Granulosa cell genes that regulate ovarian follicle development beyond the antral stage: The role of estrogen receptor ß.

Follicle development beyond the preantral stage is dependent on gonadotropins. FSH signaling is crucial for the advancement of preantral follicles to the antral stage, and LH signaling is essential for further maturation of preovulatory follicles. Estrogen is intricately tied to gonadotropin signaling during the advanced stages of folliculogenesis. We observed that Erßnull ovarian follicles fail to develop beyond the antral stage, even after exogenous gonadotropin stimulation. As ERß is primarily expressed in the granulosa cells (GCs), we explored the gonadotropin-regulated GC genes that induce maturation of antral follicles. Synchronized follicle development was induced by administration of exogenous gonadotropins to wildtype 4-wk-old female rats. The GC transcriptome was analyzed via RNA-sequencing before and after gonadotropin stimulation. An Erßnull mutant model that fails to show follicle maturation was also included in order to identify the ERß-regulated genes involved at this step. We observed that specific groups of genes were differentially expressed in response to PMSG or hCG administration in wildtype rats. While some of the PMSG or hCG-induced genes showed a similar expression pattern in Erßnull GCs, a subset of PMSG- or hCG-induced genes showed a differential expression pattern in Erßnull GCs. These latter ERß-regulated genes included previously known FSH or LH target genes including Lhcgr, Cyp11a1, Cyp19a1, Pgr, Runx2, Egfr, Kiss1, and Ptgs2, which are involved in follicle development, oocyte maturation, and ovulation. We also identified novel ERß-regulated genes including Jaml, Galnt6, Znf750, Dusp9, Wnt16, and Mageb16 that failed to respond to gonadotropin stimulation in Erßnull GCs. Our findings indicate that the gonadotropin-induced spatiotemporal pattern of gene expression is essential for ovarian follicle maturation beyond the antral stage. However, expression of a subset of those gonadotropin-induced genes is dependent on transcriptional regulation by ERß.

ERß regulated ovarian kisspeptin plays an important role in oocyte maturation.

Kisspeptin (KISS1) signaling in the hypothalamic-pituitary (H-P) axis plays an essential role in regulating gonadotropin secretion. KISS1 and KISS1 receptor (KISS1R) are also expressed in the ovary; however, the role of intraovarian KISS1 signaling remains unclear. Granulosa cell (GC)-specific expression of KISS1, and oocyte-specific expression of KISS1R indicate that GC-derived KISS1 may act on oocytes. Expression of KISS1 in GCs is induced by gonadotropins but it is absent in estrogen receptor ß knockout (Erßnull) rat ovaries. We also observed that gonadotropin stimulation failed to induce maturation of Erßnull oocytes. Interestingly, KISS1 treatment of cumulus oocyte complexes (COCs) isolated from antral follicles promotes in vitro maturation of oocytes. Treatment of oocytes with KISS1 induced intracellular Ca2+ release, and increased activation of MAP kinase ERK1/2. KISS1 treatment also induced the expression of oocyte genes that are crucial for differentiation of GCs, and maturation of oocytes. Our findings suggest that ovarian KISS1-signaling plays an important role in gonadotropin induced follicle development and oocyte maturation.

Transcriptome datasets of ESR2-regulated genes in rat granulosa cells during gonadotropin-induced follicle maturation.

Disruption of estrogen receptor beta (ESR2) dysregulates granulosa cell genes essential for follicle maturation and ovulation. The datasets presented in this article depict gonadotropin-induced genes, which are differentially expressed in Esr2-null rat granulosa cells. Synchronized follicle development was initiated in four-week-old wildtype and Esr2-null female rats by administration of PMSG. Forty-eight hours after PMSG injection, further maturation of ovarian follicles was induced by hCG treatment. Granulosa cells were collected from the ovaries before gonadotropin administration, 48 h after PMSG treatment, and 4 h after hCG injection to the PMSG-treated rats. Total RNA was purified from granulosa cells and whole transcriptome was assessed by RNA-sequencing on an Illumina HiSeq X platform. RNA-seq data of wildtype and Esr2-null granulosa cells were analyzed and differentially expressed genes were identified by CLC Genomics Workbench. Gonadotropin-induced genes were identified by comparing the transcriptome data of PMSG- or hCG-induced wildtype granulosa cells with those without gonadotropin treatment. Furthermore, differentially expressed genes in Esr2-null granulosa cells were determined by comparing the transcriptome data with that of wildtype granulosa cells. These datasets can be used to recognize the gonadotropin-induced genes in granulosa cells that are Esr2-regulated and important for ovarian follicle maturation.

A Gatekeeping Role of ESR2 to Maintain the Primordial Follicle Reserve.

Over the entire reproductive lifespan in mammals, a fixed number of primordial follicles serve as the source of mature oocytes. Uncontrolled and excessive activation of primordial follicles can lead to depletion of the ovarian reserve. We observed that disruption of estrogen receptor ß (ESR2) signaling results in increased activation of primordial follicles in Esr2-null (Esr2-/-) rats. However, follicle assembly was unaffected, and the total number of follicles remained comparable between neonatal wild-type and Esr2-/- ovaries. While the activated follicle counts were increased in Esr2-/- ovary, the number of primordial follicles were markedly decreased. Excessive recruitment of primordial follicles led to premature ovarian senescence in Esr2-/- rats and was associated with reduced levels of serum AMH and estradiol. Disruption of ESR2 signaling through administration of a selective antagonist (PHTPP) increased the number of activated follicles in wildtype rats, whereas a selective agonist (DPN) decreased follicle activation. In contrast, primordial follicle activation was not increased in the absence of ESR1, indicating that the regulation of primordial follicle activation is ESR2 specific. Follicle activation was also increased in Esr2 mutants lacking the DNA binding domain, suggesting a role for the canonical transcriptional activation function. Both primordial and activated follicles express ESR2, suggesting a direct regulatory role for ESR2 within these follicles. We also detected that loss of ESR2 augmented the activation of AKT, ERK, and mTOR pathways. Our results indicate that the lack of ESR2 upregulated both granulosa and oocyte factors, which can facilitate AKT and mTOR activation in Esr2-/- ovaries leading to increased activation of primordial follicles.

SATB1-regulated transcriptome datasets of Rcho-1 rat trophoblast stem cells.

SATB homeobox 1 (SATB1) and its heterodimeric partner SATB2 play an important regulatory role in maintaining proliferation of trophoblast stem (TS) cells and in inhibiting trophoblast differentiation. To identify the SATB-regulated genes in TS cells, we studied the transcriptome changes in a 'loss of function' model of Rcho-1 rat TS cell line. Satb1 gene expression was silenced by lentiviral delivery of shRNAs targeted to exon 9 and exon 12. An Egfp shRNA was used as a non-targeted control. Total RNA was purified from shRNA-transduced Rcho-1 cells, and whole transcriptome was assessed by RNA-sequencing on an Illumina HiSeq X platform. Differentially expressed genes in Satb1 shRNA-transduced cells were identified by analyses of the RNA-sequencing data using CLC Genomics Workbench. Differentially expressed genes with each of the two different shRNAs were compared to identify SATB1-target genes and to eliminate the potential off-targets of the shRNAs. These datasets can be used to identify the SATB-regulated genes in TS cells and to understand the molecular mechanisms that regulate trophoblast proliferation and inhibit differentiation.

Transcriptome datasets of gonadotropin-induced ESR2-regulated genes in rat oocytes.

Disruption of estrogen receptor beta (ESR2) dysregulates oocyte maturation, which leads to failure of ovulation. We investigated ESR2-regulated genes during gonadotropin-induced oocyte maturation using RNA-sequencing. Through the administration of pregnant mare's serum gonadotropin (PMSG), synchronized follicle development was initiated in four-week-old wildtype and Esr2-null female rats. Forty-eight hours after the PMSG injection, human chorionic gonadotropin (hCG) was used for further maturation. Oocytes were collected from the ovaries 4 h after hCG injection. The total RNA was isolated from the oocytes and the whole oocyte transcriptome was determined by RNA-sequencing on the Illumina HiSeq4000 sequencer. RNA-sequencing data of wildtype and Esr2-null oocytes were analyzed, and differentially expressed genes were identified using the CLC Genomics Workbench. Whole oocyte transcriptome data of wildtype and Esr2-null oocytes were compared to identify the differentially expressed genes. Raw data are deposited to the NCBI Sequence Read Archive (SRA) and analyzed data are presented in this data article. These datasets can be utilized to identify the gonadotropin-induced genes in oocytes that are ESR2-regulated and important to oocyte maturation.

Disruption of ESR1 alters the expression of genes regulating hepatic lipid and carbohydrate metabolism in male rats.

The liver helps maintain energy homeostasis by synthesizing and storing glucose and lipids. Gonadal steroids, particularly estrogens, play an important role in regulating metabolism. As estrogens are considered female hormones, metabolic disorders related to the disruption of estrogen signaling have mostly been studied in females. Estrogen receptor alpha (ESR1) is the predominant receptor in both the male and female liver, and it mediates the hepatic response to estrogens. Loss of ESR1 increases weight gain and obesity in female rats, while reducing the normal growth in males. Although Esr1-/- male rats have a reduced body weight, they exhibit increased adipose deposition and impaired glucose tolerance. We further investigated whether these metabolic disorders in Esr1-/- male rats were linked with the loss of transcriptional regulation by ESR1 in the liver. To identify the ESR-regulated genes, RNA-sequencing was performed on liver mRNAs from wildtype and Esr1-/- male rats. Based on an absolute fold change of ≥2 with a p-value ≤ 0.05, a total of 706 differentially expressed genes were identified in the Esr1-/- male liver: 478 downregulated, and 228 upregulated. Pathway analyses demonstrate that the differentially expressed genes include transcriptional regulators (Cry1, Nr1d1, Nr0b2), transporters (Slc1a2), and regulators of biosynthesis (Cyp7b1, Cyp8b1), and hormone metabolism (Hsd17b2, Sult1e1). Many of these genes are also integral parts of the lipid and carbohydrate metabolism pathways in the liver. Interestingly, certain critical regulators of the metabolic pathways displayed a sexual dimorphism in expression, which may explain the divergent weight gain in Esr1-/- male and female rats despite common metabolic dysfunctions.

Transcriptome data analyses of prostatic hyperplasia in Esr2 knockout rats.

Estrogen signaling plays an important role in the pathophysiology of prostatic hyperplasia. While signaling through estrogen receptor alpha (ESR1) increases proliferation of stromal cells, estrogen receptor beta (ESR2) plays an anti-proliferative and differentiating role in glandular epithelium. Disruption of ESR2 signaling resulted in prostatic glandular hyperplasia in the rat. To identify the ESR2-target genes, and the molecular mechanisms involved, we performed RNA-seq analyses in prostate glands of Esr2 knockout (Esr2-/-) and age matched wildtype rats. The raw data were analyzed using CLC genomics workbench. High quality RNA-seq reads were aligned to the Rattus norvegicus genome. Differentially expressed genes were identified based on an absolute fold change of 2 with pValue ≤0.05. Of the total 32,623 genes detected, 824 were differentially expressed in Esr2-/- prostate glands, 550 downregulated and 274 upregulated. Pathway analyses identified altered expression of genes involved in epithelial proliferation and benign tumor formation.

Liver transcriptome data of Esr1 knockout male rats reveals altered expression of genes involved in carbohydrate and lipid metabolism.

Estrogens are traditionally considered to be female sex steroid hormones and most of the studies examining estrogen regulation of metabolic function in the liver have been conducted in females. However, the liver expresses high levels of estrogen receptor alpha (ESR1) in both males and females, which mediates the hepatic response to estrogens. In this data article, we investigated whether metabolic disorders in Esr1 knockout (Esr1-/-) male rats were linked with loss of transcriptional regulation by ESR1 in liver. To identify the ESR1 regulated genes in the mutant liver, RNA-sequencing was performed on liver RNAs purified from young male rats. The raw data were analyzed using the CLC Genomics Workbench and high-quality RNA-sequencing reads were aligned to the Rattus norvegicus genome. Transcriptome data obtained from Esr1-/- liver RNAs were compared to that of wild type rats. Based on an absolute fold change of 2 with a p-value ≤ 0.05, a total of 618 differentially expressed genes were identified in the Esr1-/- male liver. Pathway analyses demonstrated that the majority of differentially expressed genes are regulators of carbohydrate and lipid metabolism in the liver. These differentially expressed genes and their potential roles were further examined in a companion manuscript, "Disruption of ESR1 alters the expression of genes regulating hepatic lipid and carbohydrate metabolism in male rats" (Khristi et al., 2018).

ESR2 Is Essential for Gonadotropin-Induced Kiss1 Expression in Granulosa Cells.

Hypothalamic expression of Kiss1 plays an essential role in the onset of puberty, gonadal development, and ovulation. Estrogens regulate the expression of Kiss1 in the hypothalamus through estrogen receptor-α. Kiss1 is also expressed in the ovary, where its expression correlates with the onset of puberty and progression of the estrous cycle. To date, estrogen regulation of Kiss1 expression in the ovary has not been investigated. We recently observed that gonadotropin-induced Kiss1 expression was absent in Esr2-null rat ovaries even though Esr1 was present. Wild-type granulosa cells abundantly expressed Kiss1 and oocytes expressed the Kiss1 receptor. We characterized estrogen receptor-ß (ESR2) regulation of Kiss1 expression in granulosa cells by identifying granulosa cell-specific transcript variants and potential regulatory regions. The Kiss1 promoter, an upstream enhancer, and a downstream enhancer all possessed conserved estrogen response elements (EREs) and showed active histone marks in gonadotropin-stimulated granulosa cells. The transcriptionally active Kiss1 promoter, as well as the enhancers, also revealed enrichment for ESR2 binding. Furthermore, activity of a Kiss1 promoter construct was induced after overexpression of ESR2 and was blocked upon mutation of an ERE within the promoter. Finally, pregnant mare serum gonadotropin and human chorionic gonadotropin administration induced phosphorylation of ESR2 and upregulated the AP-1 proteins FOSL2 and JUNB in granulosa cells. Activated MAPK ERK2 was associated with the ESR2 phosphorylation in granulosa cells, and AP-1 factors could synergistically activate the Kiss1 promoter activity. These gonadotropin-induced changes paralleled Kiss1 expression in granulosa cells. We conclude that gonadotropin-stimulated Kiss1 expression in granulosa cells is dependent on both the activation of ESR2 and the upregulation of AP-1.

Differentially regulated genes in Esr2-mutant rat granulosa cells.

RNA seq analyses were performed in granulosa cells (GCs) collected from gonadotropin treated ESR2 mutant rats. Data obtained from a null mutant with Esr2 exon 3 deletion (∆3) and another DNA binding domain (DBD) mutant with exon 4 deletion (∆4) were compared to that of wildtype (WT) rats. The raw data were analyzed using CLC genomics workbench. High quality RNA-sequencing reads were aligned to the Rattus norvegicus genome. Differentially expressed genes in ∆3 or ∆4 Esr2-mutant GCs were identified based on the following criteria: FDR p-Value ≤0.05 and an absolute fold change of 2. Fewer differentially expressed genes were identified in ∆3 compared to the ∆4 mutant group. As both mutant groups demonstrated a common phenotype of ovulation failure, differentially expressed genes common to both in ∆3 and ∆4 mutant rats were emphasized and further analyzed in the companion article "ESR2 regulates granulosa cell genes essential for follicle maturation and ovulation" [1].

ESR2 regulates granulosa cell genes essential for follicle maturation and ovulation.

Estrogen receptor 2 (ESR2) plays a critical role in folliculogenesis and ovulation. Disruption of ESR2-function in the rats results in female infertility due to failure of ovulation. Ovulation failure occurred in two distinct rat models, a null mutant and a DNA binding domain (DBD) mutant of ESR2, indicating that transcriptional regulation by ESR2 is indispensable for ovulation. To define the regulatory role of ESR2 in preovulatory follicular maturation and ovulation, we investigated ovarian responsiveness to exogenous gonadotropins in prepubertal females. Granulosa cells (GCs) play a vital role in follicle maturation and ovulation, and ESR2-dependent estrogen signaling is predominant in GCs, therefore, we examined the differential expression of gonadotropin-induced genes in GCs. Of 32,623 genes detected by RNA-sequencing, 1696 were differentially expressed in Esr2-mutant rats (789 downregulated, and 907 upregulated, absolute fold change 2, FDR p < 0.05). Molecular pathway analyses indicated that these differentially expressed genes are involved in steroidogenesis, follicle maturation, and ovulation. Many of these genes are known regulators of ovarian function and a subset were also disrupted in Esr2-mutant mice. Interestingly, Kiss1 was identified as one of the differentially expressed genes implicating a potential role within the follicle and its regulation by ESR2. Our findings indicate that ESR2 regulates key genes in GCs that are essential for follicle maturation and ovulation in the rat.