Gene expression profiling of bovine ovarian follicular and luteal cells provides insight into cellular identities and functions

Unraveling the role of lipid droplets and perilipin 2 in bovine luteal cells

Steroidogenic tissues contain cytosolic lipid droplets that are important for steroidogenesis. Perilipin 2 (PLIN2), a structural coat protein located on the surface of lipid droplets in mammalian cells, plays a crucial role in regulating lipid droplet formation and contributing to various cellular processes such as lipid storage and energy homeostasis. Herein, we examine the role that PLIN2 plays in regulating progesterone synthesis in the bovine corpus luteum. Utilizing gene array databases and Western blotting, we have delineated the expression pattern of PLIN2 throughout the follicular to luteal transition. Our findings reveal the presence of PLIN2 in both ovarian follicular and steroidogenic luteal cells, demonstrating an increase in its levels as follicular cells transition into the luteal phase. Moreover, the depletion of PLIN2 via siRNA enhanced progesterone production in small luteal cells, whereas adenovirus-mediated overexpression of both PLIN2 and Perilipin 3 (PLIN3) induced an increase in cytosolic lipid droplet accumulation and decreased hormone-induced progesterone synthesis in these cells. Lastly, in vivo administration of the luteolytic hormone prostaglandin F2α resulted in an upregulation of PLIN2 mRNA and protein expression, accompanied by a decline in serum progesterone. Our findings highlight the pivotal role of PLIN2 in regulating progesterone synthesis in the bovine corpus luteum, as supported by its dynamic expression pattern during the follicular to luteal transition and its responsiveness to luteotropic and luteolytic hormones. We suggest PLIN2 as a potential therapeutic target for modulating luteal function.

Imprinted lncRNA KCNQ1OT1 regulates CDKN1C expression through promoter binding and chromatin folding in pigs

Long noncoding RNAs (lncRNAs) are implicated in a number of regulatory functions in eukaryotic genomes. In humans, KCNQ1OT1 is a 91 kb imprinted lncRNA that inhibits multiple surrounding genes in cis. Among them, CDKN1C is closely related to KCNQ1OT1 and is involved in multiple epigenetic disorders. Here, we found that pigs also had a relatively conserved paternal allele expressing KCNQ1OT1 and had a shorter 5' end (∼27 kb) compared to human KCNQ1OT1. Knockdown of KCNQ1OT1 using antisense oligonucleotides (ASO) showed that upregulation of CDKN1C expression in pigs. However, porcine KCNQ1OT1 did not affect the DNA methylation status of the CpG islands in the promoters of KCNQ1OT1 and CDKN1C. Inhibition of DNA methyltransferase using Decitabine treatment resulted in a significant increase in both KCNQ1OT1 and CDKN1C expression, suggesting that the regulation between KCNQ1OT1 and CDKN1C may not be dependent on RNA interference. Further use of chromosome conformation capture and reverse transcription-associated trap detection in the region where CDKN1C was located revealed that KCNQ1OT1 bound to the CDKN1C promoter and affected chromosome folding. Phenotypically, inhibition of KCNQ1OT1 at the cumulus-oocyte complex promoted cumulus cell transformation, and to upregulated the expression of ALPL at the early stage of osteogenic differentiation of porcine bone marrow mesenchymal stem cells. Our results confirm that the expression of KCNQ1OT1 imprinting in pigs as well as porcine KCNQ1OT1 regulates the expression of CDKN1C through direct promoter binding and chromatin folding alteration. And this regulatory mechanism played an important role in cell differentiation.

Saturated fatty acids inhibit unsaturated fatty acid induced glucose uptake involving GLUT10 and aerobic glycolysis in bovine granulosa cells

Fatty acids have been shown to modulate glucose metabolism in vitro and in vivo. However, there is still a need for substantial evidence and mechanistic understanding in many cell types whether both saturated and unsaturated fatty acids (SFAs and UFAs) pose a similar effect and, if not, what determines the net effect of fatty acid mixes on glucose metabolism. In the present study, we asked these questions by treating granulosa cells (GCs) with the most abundant non-esterified fatty acid species in bovine follicular fluid. Results revealed that oleic and alpha-linolenic acids (UFAs) significantly increased glucose consumption compared to palmitic and stearic acids (SFAs). A significant increase in lactate production, extracellular acidification rate, and decreased mitochondrial activity indicate glucose channeling through aerobic glycolysis in UFA treated GCs. We show that insulin independent glucose transporter GLUT10 is essential for UFA driven glucose consumption, and the induction of AKT and ERK signaling pathways necessary for GLUT10 expression. To mimic the physiological conditions, we co-treated GCs with mixes of SFAs and UFAs. Interestingly, co-treatments abolished the UFA induced glucose uptake and metabolism by inhibiting AKT and ERK phosphorylation and GLUT10 expression. These data suggest that the net effect of fatty acid induced glucose uptake in GCs is determined by SFAs under physiological conditions.

Central Role for Glycolysis and Fatty Acids in LH-responsive Progesterone Synthesis

Progesterone production by the corpus luteum is fundamental for establishing and maintaining pregnancy. The pituitary gonadotropin luteinizing hormone (LH) is recognized as the primary stimulus for luteal formation and progesterone synthesis, regardless of species. Previous studies demonstrated an elevation in abundance of genes related to glucose and lipid metabolism during the follicular to luteal transition. However, the metabolic phenotype of these highly steroidogenic cells has not been studied. Herein, we determined acute metabolic changes induced by LH in primary luteal cells and defined pathways required for progesterone synthesis. Untargeted metabolomics analysis revealed that LH induces rapid changes in vital metabolic pathways, including glycolysis, tricarboxylic acid (TCA) cycle, pentose phosphate pathway, de novo lipogenesis, and hydrolysis of phospholipids. LH stimulated glucose uptake, enhanced glycolysis, and flux of [U- 13 C 6 ]-labeled glucose-derived carbons into metabolic branches associated with adenosine 5'-triphosphate (ATP) and NADH/NADPH production, synthesis of nucleotides, proteins, and lipids, glycosylation of proteins or lipids, and redox homeostasis. Selective use of small molecule inhibitors targeting the most significantly changed pathways, such as glycolysis, TCA cycle, and lipogenesis, uncovered cellular metabolic routes required for LH-stimulated steroidogenesis. Furthermore, LH via the protein kinase A (PKA) pathway triggered post- translational modification of acetyl-CoA carboxylase alpha (ACACA) and ATP citrate lyase (ACLY), enzymes involved in de novo synthesis of fatty acids. Inhibition of ACLY and fatty acid transport into mitochondria reduced LH-stimulated ATP, cAMP production, PKA activation, and progesterone synthesis. Taken together, these findings reveal novel hormone-sensitive metabolic pathways essential for maintaining LHCGR/PKA signaling and steroidogenesis in ovarian luteal cells.

Significance

The establishment and maintenance of pregnancy require a well-developed corpus luteum, an endocrine gland within the ovary that produces progesterone. Although there is increased awareness of intracellular signaling events initiating the massive production of progesterone during the reproductive cycle and pregnancy, there are critical gaps in our knowledge of the metabolic and lipidomic pathways required for initiating and maintaining luteal progesterone synthesis. Here, we describe rapid, hormonally triggered metabolic pathways, and define metabolic targets crucial for progesterone synthesis by ovarian steroidogenic cells. Understanding hormonal control of metabolic pathways may help elucidate approaches for improving ovarian function and successful reproduction or identifying metabolic targets for developing nonhormonal contraceptives.

Prostaglandin F2α regulates mitochondrial dynamics and mitophagy in the bovine corpus luteum

Prostaglandins are arachidonic acid-derived lipid mediators involved in numerous physiological and pathological processes. PGF2α analogues are therapeutically used for regulating mammalian reproductive cycles and blood pressure, inducing term labor, and treating ocular disorders. PGF2α exerts effects via activation of calcium and PKC signaling, however, little is known about the cellular events imposed by PGF2α signaling. Here, we explored the early effects of PGF2α on mitochondrial dynamics and mitophagy in the bovine corpus luteum employing relevant and well characterized in vivo and in vitro approaches. We identified PKC/ERK and AMPK as critical protein kinases essential for activation of mitochondrial fission proteins, DRP1 and MFF. Furthermore, we report that PGF2α elicits increased intracellular reactive oxygen species and promotes receptor-mediated activation of PINK-Parkin mitophagy. These findings place the mitochondrium as a novel target in response to luteolytic mediator, PGF2α. Understanding intracellular processes occurring during early luteolysis may serve as a target for improving fertility.

Spheroid formation and luteinization of granulosa cells of felids in a long-term 3D culture

In the present study, granulosa cells (GCs) from domestic cats and Persian leopard were cultured and characterized from selected days. The culture period was divided into two phases: maintenance, which lasted for 7 days, and luteinization, which followed for up to 11 days. Luteinization was performed on ultra-low attachment plates, supporting the formation of spheroids in a medium supplemented with insulin, forskolin, and luteinizing hormone (LH). GCs of domestic cat produced estradiol (E2) and progesterone (P4) during the maintenance phase. The gene expressions of some proteins involved in steroidogenesis were stable (STAR, HSD3B1) or decreased over time (CYP11A1, HSD17B1, CYP17A1, and CYP19A1), which was similar to the expressions of gonatropin receptors (LHCGR and FSHR). During the luteinization phase, P4 concentration significantly increased (P < 0.05), and E2, in contrast to the proliferation phase, was below detection range. The expression of genes of proteins involved in steroidogenesis (STAR, CYP11A1, HSD3B1, HSD17B1, CYP17A1, and CYP19A1) and of gonadotropin receptors (LHCGR and FSHR) significantly increased during the luteinization period, but some expressions exhibited a decrease at the end of the phase (LHCGR, FSHR, HSD17B1, CYP19A1). The morphology of the luteinized GCs of domestic cat resembled large luteal cells and had numerous vacuole-like structures. Also, the GCs of Persian leopard underwent luteinization, shown by increasing P4 production and HSD3B1 expression. This study confirms that GCs from felids can be luteinized in a 3D spheroid system which can be a basis for further studies on luteal cell function of felids. Additionally, we could show that the domestic cat can serve as a model species for establishing cell culture methods which can be transferred to other felids.

The characterization and therapeutic applications of ovarian theca cells: An update

Theca cells perform a range of roles during folliculogenesis. So far, little is known about their recruitment process and function since early research has mainly focused on the interactions between granulosa cells and the oocyte, leaving theca cells unfairly forgotten in the understanding of ovarian physiology and pathogenesis. Given that research on theca cells has greatly emerged in recent years, this review of literature aims to discuss the established theoretical concepts with the most recent findings about theca cells' characterization and origins, in vitro culture applications as models for fertility preservation and pharmacological/toxicological studies, its importance in unraveling pathogenic pathways, and stem-cell-based bioengineering for hormonal replacement therapies. Isolation and in vitro culture techniques for theca cells have led to essential advancements in their characterization as a specific cell population. Unraveling the origins of theca cells during the in vivo differentiation process in the adult ovary will assist the development of hormonal replacement therapies, reestablishment of fertility, and treatments for diseases such as premature ovarian insufficiency and polycystic ovarian syndrome, which seem to be directly influenced by theca cells.

Perfluorooctanoic acid (PFOA) promotes follicular growth and alters expression of genes that regulate the cell cycle and the Hippo pathway in cultured neonatal mouse ovaries

Perfluorooctanoic acid (PFOA) is a synthetic chemical resistant to biodegradation and is environmentally persistent. PFOA is found in many consumer products and is a major source of water contamination. While PFOA has been identified as a contaminant of concern for reproductive health, little is known about the effects of PFOA on ovarian follicular development and growth. Recent evidence indicates that the Hippo pathway is an important regulator of ovarian physiology. Here, we investigated the effects of PFOA on ovarian folliculogenesis during the neonatal period of development and potential impacts on the Hippo signaling pathway. Post-natal day 4 (PND4) neonatal ovaries from CD-1 mice were cultured with control medium (DMSO <0.01% final concentration) or PFOA (50 μM or 100 μM). After 96 h, ovaries were collected for histological analysis of folliculogenesis, gene and protein expression, and immunostaining. Results revealed that PFOA (50 μM) increased the number of secondary follicles, which was accompanied by increases in mRNA transcripts and protein of marker of proliferation marker Ki67 with no impacts on apoptosis markers Bax, Bcl2, or cleaved caspase-3. PFOA treatment (50 μM and 100 μM) stimulated an upregulation of transcripts for cell cycle regulators Ccna2, Ccnb2, Ccne1, Ccnd1, Ccnd2, and Ccnd3. PFOA also increased abundance of transcripts of Hippo pathway components Mst1/2, Lats1, Mob1b, Yap1, and Taz, as well as downstream Hippo pathway targets Areg, Amotl2, and Cyr61, although it decreased transcripts for anti-apoptotic Birc5. Inhibition of the Hippo pathway effector YAP1 with Verteporfin resulted in the attenuation of PFOA-induced follicular growth and proliferation. Together, these findings suggest that occupationally relevant levels of PFOA (50 μM) can stimulate follicular activation in neonatal ovaries potentially through activation of the Hippo pathway.

Follicular Hyperstimulation Dysgenesis: New Explanation for Adverse Effects of Excessive FSH in Ovarian Stimulation

High follicle-stimulating hormone (FSH) doses during ovarian stimulation protocols for assisted reproductive technologies (ART) are detrimental to ovulatory follicle function and oocyte quality. However, the mechanisms are unclear. In a small ovarian reserve heifer model, excessive FSH doses lead to phenotypic heterogeneity of ovulatory size follicles, with most follicles displaying signs of premature luteinization and a range in severity of abnormalities. By performing whole transcriptome analyses of granulosa cells, cumulus cells, and oocytes from individual follicles of animals given standard or excessive FSH doses, we identified progressive changes in the transcriptomes of the 3 cell types, with increasing severity of follicular abnormality with the excessive doses. The granulosa and cumulus cells each diverged progressively from their normal phenotypes and became highly similar to each other in the more severely affected follicles. Pathway analysis indicates a possible dysregulation of the final stages of folliculogenesis, with processes characteristic of ovulation and luteinization occurring concurrently rather than sequentially in the most severely affected follicles. These changes were associated with disruptions in key pathways in granulosa and cumulus cells, which may account for previously reported reduced estradiol production, enhanced progesterone and oxytocin production and diminished ovulation rates. Predicted deficiencies in oocyte survival, stress response, and fertilization suggest likely reductions in oocyte health, which could further compromise oocyte quality and ART outcomes.

Use of Agriculturally Important Animals as Models in Biomedical Research

Livestock have contributed significantly to advances in biomedicine and offer unique advantages over rodent models. The human is the ideal biomedical model; however, ethical reasons limit the testing of hypotheses and treatments in humans. Rodent models are frequently used as alternatives to humans due to size, low cost, and ease of genetic manipulation, and have contributed tremendously to our understanding of human health and disease. However, the use of rodents in translational research pose challenges for researchers due to physiological differences to humans. The use of livestock species as biomedical models can address these challenges as livestock have several similarities to human anatomy, physiology, genetics, and metabolism and their larger size permits collection of more frequent and often larger samples. Additionally, recent advances in genetics in livestock species allow for studies in genomics, proteomics, and metabolomics, which have the added benefit of applications to both humans in biomedical research and livestock in improving production. In this review, we provide an overview of scientific findings using livestock and benefits of each model to the livestock industry and to biomedical research.

Dynamic patterns of expressed genes in granulosa cells during follicular and luteal stages in Egyptian buffaloes

A better understanding of the molecular mechanisms in granulosa cells (GC) is warranted, during different follicular and luteal developmental stages in buffalo cows. We aimed to (I) study the expression of selected genes in GC during follicular and luteal phases, (II) evaluate correlations between GC gene expression and steroid concentrations {17-beta estradiol (E2) and progesterone (P4)} in follicular fluid (FF), and (III) study effect of ovarian status on follicular population as well as follicular size frequency. Ovaries were collected in pairs from buffaloes (n = 178). Ovaries bearing corpus luteum (CL) were subdivided into hemorrhagic, developing, mature, and albicans. Follicles from luteal groups were classified only into small (< 4 mm) and large (9-20 mm), while follicles from follicular groups were classified into three subgroups: small (< 4 mm), medium (5-8 mm), and large (9-20 mm). The FF and GC were collected for steroid concentrations measurement and gene expression, respectively. In the follicular phase, luteinizing hormone/choriogonadotropin receptor (LHCGR) and cytochrome P450 aromatase (CYP19) in small follicles decreased compared to medium ones. Large follicle showed an increase in LHCGR and CYP19 compared to medium ones. Follicle-stimulating hormone receptor (FSHR) decreased in large compared to medium size follicles. Proliferating cell nuclear antigen (PCNA) increased in small and large follicles. Meanwhile, anti-Mullerian hormone (AMH) and phospholipase A2 group III (PLA2G3) decreased in small and large follicles. The different stages of luteal phase had a profound impact on GC gene expression. There were strong (positive and/or negative) correlations between gene expression and steroid hormones. The different scenarios between expressed genes in GC and steroid concentrations are required for the proper growth and development of follicles and CL.

Transcriptome analysis reveals transforming growth factor-β1 prevents extracellular matrix degradation and cell adhesion during the follicular-luteal transition in cows

Ovarian angiogenesis is an extremely rapid process that occurs during the transition from follicle to corpus luteum (CL) and is crucial for reproduction. It is regulated by numerous factors including transforming growth factor-β1 (TGFB1). However, the regulatory mechanism of TGFB1 in ovarian angiogenesis is not fully understood. To address this, in this study we obtained high-throughput transcriptome analysis (RNA-seq) data from bovine luteinizing follicular cells cultured in a system mimicking angiogenesis and treated with TGFB1, and identified 455 differentially expressed genes (DEGs). Quantitative real-time PCR results confirmed the differential expression patterns of the 12 selected genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified that the MAPK and ErbB pathways, cell adhesion molecules (CAMs), and extracellular matrix (ECM)-receptor interactions may play pivotal roles in TGFB1-mediated inhibition of CL angiogenesis. TGFB1 phosphorylated ERK1/2 (MAPK1/3) and Akt, indicating that these pathways may play an important role in the regulation of angiogenesis. Several genes with specific functions in cell adhesion and ECM degradation were identified among the DEGs. In particular, TGFB1-induced upregulation of syndecan-1 (SDC1) and collagen type I alpha 1 chain (COL1A1) expression may contribute to the deposition of type I collagen in luteinizing follicular cells. These results indicate that TGFB1 inhibits cell adhesion and ECM degradation processes involving ERK1/2, ErbB, and PI3K/Akt signaling pathways, and leads to inhibition of angiogenesis during the follicular-luteal transition. Our results further reveal the molecular mechanisms underlying the actions of TGFB1 in early luteinization.

Transcriptomic data of bovine ovarian granulosa cells of control and High A4 cows

Microarray analysis using Affymetrix Bovine GeneChip 1.0 ST Array to determine RNA expression analysis was performed on somatic granulosa cells from two different groups of cows classified based on androstenedione concentration within the follicular fluid (Control vs High A4) of estrogen-active dominant follicles. The normalized linear microarray data was deposited to the NCBI GEO repository (GSE97017 - RNA Expression Data from Bovine Ovarian Granulosa Cells from High or Low Androgen-Content Follicles). Subsequent ANOVA determined genes that were enriched (≥ 1.5 fold more) or decreased (≤ 1.5 fold less) in the High A4 granulosa cells compared to Control granulosa cells and analyzed filtered datasets of these differentially expressed genes are presented as tables. MicroRNAs that are differentially expressed in Control and High A4 granulosa cells are also reported in tables. The standard deviation of the analyzed array data in relation to the log of the expression values are shown as a figure. Ingenuity Pathway Analysis determined upstream regulators of differently expressed genes as presented in a table. These data have been further analyzed and interpreted in the companion article "A High-Androgen Microenvironment Inhibits Granulosa Cell Proliferation and Alters Cell Identity" (McFee et. al., 2021 [1].

A high-androgen microenvironment inhibits granulosa cell proliferation and alters cell identity

A naturally occurring bovine model with excess follicular fluid androstenedione (High A4), reduced fertility, and polycystic ovary syndrome (PCOS)-like characteristics has been identified. We hypothesized High A4 granulosa cells (GCs) would exhibit altered cell proliferation and/or steroidogenesis. Microarrays of Control and High A4 GCs combined with Ingenuity Pathway Analysis indicated that High A4 GCs had cell cycle inhibition and increased expression of microRNAs that inhibit cell cycle genes. Granulosa cell culture confirmed that A4 treatment decreased GC proliferation, increased anti-Müllerian hormone, and increased mRNA for CTNNBIP1. Increased CTNNBIP1 prevents CTNNB1 from interacting with members of the WNT signaling pathway thereby inhibiting the cell cycle. Expression of CYP17A1 was upregulated in High A4 GCs presumably due to reduced FOS mRNA expression compared to Control granulosa cells. Furthermore, comparisons of High A4 GC with thecal and luteal cell transcriptomes indicated an altered cellular identity and function contributing to a PCOS-like phenotype.

Attainment and maintenance of pubertal cyclicity may predict reproductive longevity in beef heifers†

We hypothesized the manner that heifers achieve puberty may indicate their future reproductive longevity. Heifers with discontinued or delayed cyclicity during puberty attainment may have irregular reproductive cycles, anovulation, and infertility in their first breeding season contributing to a shorter reproductive lifespan. Therefore, plasma progesterone (P4) was measured from weaning to breeding on 611 heifers born 2012-2017 and four pubertal classifications were identified: (1) Early; P4 ≥ 1 ng/ml < March 12 with continued cyclicity, (2) Typical; P4 ≥ 1 ng/ml ≥ March 12 with continued cyclicity, (3) Start-Stop; P4 ≥ 1 ng/ml but discontinued cyclicity, and (4) Non-Cycling; no P4 ≥ 1 ng/ml. Historical herd records indicated that 25% of heifers achieved puberty prior to March 12th in the 10 years prior to the study. Start-Stop and Non-Cycling yearling heifers were lighter indicating reduced growth and reproductive maturity traits compared with Early/Typical heifers. In addition, Non-Cycling/Start-Stop heifers were less responsive to prostaglandin F2 alpha (PGF2α) to initiate estrous behavior and ovulation to be artificially inseminated. Non-Cycling heifers had fewer reproductive tract score-5 and reduced numbers of calves born in the first 21-days-of-calving during their first breeding season. Within the Start-Stop classification, 50% of heifers reinitiated cyclicity with growth traits and reproductive parameters that were similar to heifers in the Early/Typical classification while those that remained non-cyclic were more similar to heifers in the Non-Cycling group. Thus, heifers with discontinued cyclicity or no cyclicity during puberty attainment had delayed reproductive maturity resulting in subfertility and potentially a shorter reproductive lifespan.

Prostaglandin F2α regulation and function during ovulation and luteinization in cows

Although prostaglandins are important in the ovulation process, a precise role for prostaglandin F2α (PGF) has not been elucidated. This study aimed to evaluate the regulation of PGF receptor mRNA (PTGFR) in granulosa cells and the local effect of PGF on ovulation and luteinization. In Experiment 1, using samples collected in vivo before (Day 2), during (Day 3) and after (Day 4) follicular deviation, expression of PTGFR in bovine granulosa cells was more abundant in the dominant follicle after deviation than in subordinates (P < 0.05). However, the expression of PTGFR was not regulated (P = 0.1) in preovulatory follicles at different time-points (0, 3, 6, 12 and 24 h) after ovulation induction with GnRH. In Experiment 2, to assess the role of systemic PGF treatment on luteinization and vascularization of preovulatory follicles, flunixin meglumine (FM), a nonsteroidal anti-inflammatory drug, was used to inhibit endogenous prostaglandin synthesis. Cows with preovulatory follicles were induced to ovulate with GnRH (0 h) and allocated to three groups: Control, with no further treatment; FM, treated with 2.2 mg/kg FM im 17 h after GnRH treatment; and FM + PGF, treated with FM 17 h after GnRH, followed by 25 mg dinoprost tromethamine (PGF) 23 h after GnRH treatment. FM injection was able to reduce the concentration of PGF in the follicular fluid (FF) (P < 0.001). However, contrary to our hypothesis, color Doppler ultrasound evaluations revealed decreased vascular flow in FM + PGF group (P < 0.05), and no effect of the treatments on intrafollicular P4 and E2 concentrations 24 h after GnRH. The prostaglandin metabolite (PGFM) concentrations in the FF were greater in cows receiving systemic PGF (P < 0.001), which prompted us to further check its role on ovulation. Therefore, in Experiment 3, in a final attempt to demonstrate the local effect of PGF on ovulation, cows with preovulatory follicles received an intrafollicular injection (IFI) of PBS (Control) or 100 ng/mL purified PGF (PGF group). PGF treatment did not affect the time of ovulation after IFI (66 ± 6.4 and 63 ± 8.5 h for control and PGF, respectively; P > 0.05), further suggesting that it has no direct effect in the ovulatory process. Based on our findings, we concluded that FM decreased PGF synthesis within the follicle, whereas PGF treatment decreased follicular vascularization. In addition, the in vivo model of intrafollicular injection evidenced that PGF alone is not able to locally induce ovulation.

Prostaglandin F2 Alpha Triggers the Disruption of Cell Adhesion with Cytokeratin and Vimentin in Bovine Luteal Theca Cells

Simple Summary

Luteolysis is an important event in the control of the corpus luteum function in bovines. However, some aspects of the luteolytic mechanism remain unclear. We evaluated changes in cell adhesion in luteal cells during regression of corpus luteum. Bovine luteal theca cells (LTCs) were treated in vitro with Prostaglandin F2 alpha (PGF2α). Cytokeratin, vimentin and desmoplakin proteins in LTCs were disrupted by PGF2α, affecting cell adhesion. These results suggest that PGF2α plays an important function in cell adhesion during the regression of corpus luteum.

Abstract

Intermediate filaments (IFs) maintain cell–cell adhesions and are involved in diverse cellular processes such as cytokinesis, cell migration and the maintenance of cell structure. In this study, we investigated the influence of prostaglandin F2 alpha (PGF2α) on cytokeratin and vimentin IFs, Rho-associated protein kinase (ROCK), and cell-cell adhesion in bovine luteal theca cells (LTCs). The luteal cells were isolated from bovine corpus luteum (CL), and the LTCs were treated with 0, 0.01, 0.1 and 1.0 mM PGF2α. Cytokeratin, vimentin and desmoplakin proteins were disrupted and the ROCK protein was significantly increased in PGF2α-treated LTCs. In addition, cell–cell adhesion was significantly (p < 0.05) decreased in the PGF2α-induced LTCs compared to control group (0 mM PGF2α). In conclusion, PGF2α affected the adhesion of cell to cell via disruption of desmoplakin, cytokeratin and vimentin, additionally increasing ROCK in bovine LTCs. These results may provide a better understanding of the mechanism of bovine CL regression.

Quercetin alleviated H2 O2 -induced apoptosis and steroidogenic impairment in goat luteinized granulosa cells

Quercetin (Que) is a natural flavonoid in most plants. Luteinized granulosa cell (LGC) culture is necessary for the study of follicle growth/differentiation. In the present study, we analyzed the role of Que in steroid production and apoptosis in hydrogen peroxide (H2 O2 )-treated goat LGCs. The results showed that treatment with H2 O2 induced apoptosis in goat LGCs, and treatment with Que decreased LGC apoptosis induced by H2 O2 (P < .05), accompanied with the different expressions of BAX, BCL-2, Caspase 3, and Cleaved caspase 3. Meanwhile, the messenger RNA expressions of nuclear factor erythroid 2 like 2 (Nrf2) and its downstream genes were upregulated with H2 O2 +Que treatment, accompanied by the increased cellular viability (P < .05). Furthermore, Que alleviated H2 O2 -induced reduction in the secretion of progesterone (P4 ) (P < .05); however, it had no effect on the secretion of estrogen (E2 ). Simultaneously, the expressions of StAR and P450scc were increased when treated with Que +H2 O2 , compared with the group treated with only H2 O2 (P < .05). In conclusion, it is observed that Que could alleviate the H2 O2 -induced apoptosis and steroidogenic impairment in goat LGCs, which might be mediated by the Nrf2 pathway.

miR-4463 regulates aromatase expression and activity for 17β-estradiol synthesis in response to follicle-stimulating hormone

Objective

The aim of this study was to investigate microRNAs (miRNAs) related to follicle-stimulating hormone (FSH) responsiveness using miRNA microarrays and to identify their target genes to determine the molecular regulatory pathways involved in FSH signaling in KGN cells.

Methods

To change the cellular responsiveness to FSH, KGN cells were treated with FSH receptor (FSHR)-specific small interfering RNA (siRNA) followed by FSH. miRNA expression profiles were determined through miRNA microarray analysis. Potential target genes of selected miRNAs were predicted using bioinformatics tools, and their regulatory function was confirmed in KGN cells.

Results

We found that six miRNAs (miR-1261, miR-130a-3p, miR-329-3p, miR-185-5p, miR-144-5p and miR-4463) were differentially expressed after FSHR siRNA treatment in KGN cells. Through a bioinformatics analysis, we showed that these miRNAs were predicted to regulate a large number of genes, which we narrowed down to cytochrome P450 family 19 subfamily A member 1 (CYP19A1) and estrogen receptor alpha (ESR1) as the main targets for miR-4463. Functional analysis revealed that miR-4463 is a regulatory factor for aromatase expression and function in KGN cells.

Conclusion

In this study, we identified differentially expressed miRNAs related to FSH responsiveness. In particular, upregulation of miR-4463 expression by FSHR deficiency in human granulosa cells impaired 17β-estradiol synthesis by targeting CYP19A1 and ESR1. Therefore, our data might provide novel candidates for molecular biomarkers for use in research into poor responders.

Molecular profiling demonstrates modulation of immune cell function and matrix remodeling during luteal rescue†

The corpus luteum (CL) is essential for maintenance of pregnancy in all mammals and luteal rescue, which occurs around day 16-19 in the cow, is necessary to maintain luteal progesterone production. Transcriptomic and proteomic profiling were performed to compare the day 17 bovine CL of the estrous cycle and pregnancy. Among mRNA and proteins measured, 140 differentially abundant mRNA and 24 differentially abundant proteins were identified. Pathway analysis was performed using four programs. Modulated pathways included T cell receptor signaling, vascular stability, cytokine signaling, and extracellular matrix remodeling. Two mRNA that were less in pregnancy were regulated by prostaglandin F2A in culture, while two mRNA that were greater in pregnancy were regulated by interferon tau. To identify mRNA that could be critical regulators of luteal fate, the mRNA that were differentially abundant during early pregnancy were compared to mRNA that were differentially abundant during luteal regression. Eight mRNA were common to both datasets, including mRNA related to regulation of steroidogenesis and gene transcription. A subset of differentially abundant mRNA and proteins, including those associated with extracellular matrix functions, were predicted targets of differentially abundant microRNA (miRNA). Integration of miRNA and protein data, using miRPath, revealed pathways such as extracellular matrix-receptor interactions, abundance of glutathione, and cellular metabolism and energy balance. Overall, this study has provided a comprehensive profile of molecular changes in the corpus luteum during maternal recognition of pregnancy and has indicated that some of these functions may be miRNA-regulated.

Trafficking of cholesterol from lipid droplets to mitochondria in bovine luteal cells: Acute control of progesterone synthesis

The corpus luteum is a transient endocrine gland that synthesizes and secretes the steroid hormone, progesterone, which is vital for establishment and maintenance of pregnancy. Luteinizing hormone (LH) via activation of protein kinase A (PKA) acutely stimulates luteal progesterone synthesis via a complex process, converting cholesterol via a series of enzymatic reactions, into progesterone. Lipid droplets in steroidogenic luteal cells store cholesterol in the form of cholesterol esters, which are postulated to provide substrate for steroidogenesis. Early enzymatic studies showed that hormone sensitive lipase (HSL) hydrolyzes luteal cholesterol esters. In this study, we tested whether HSL is a critical mediator of the acute actions of LH on luteal progesterone production. Using LH-responsive bovine small luteal cells our results reveal that LH, forskolin, and 8-Br cAMP-induced PKA-dependent phosphorylation of HSL at Ser563 and Ser660, events known to promote HSL activity. Small molecule inhibition of HSL activity and siRNA-mediated knock down of HSL abrogated LH-induced progesterone production. Moreover, western blotting and confocal microscopy revealed that LH stimulates phosphorylation and translocation of HSL to lipid droplets. Furthermore, LH increased trafficking of cholesterol from the lipid droplets to the mitochondria, which was dependent on both PKA and HSL activation. Taken together, these findings identify a PKA/HSL signaling pathway in luteal cells in response to LH and demonstrate the dynamic relationship between PKA, HSL, and lipid droplets in luteal progesterone synthesis.

Formation and characterization of lipid droplets of the bovine corpus luteum

Establishment and maintenance of pregnancy depends on progesterone synthesized by luteal tissue in the ovary. Our objective was to identify the characteristics of lipid droplets (LDs) in ovarian steroidogenic cells. We hypothesized that LDs are a major feature of steroidogenic luteal cells and store cholesteryl esters. Whole bovine tissues, isolated ovarian steroidogenic cells (granulosa, theca, small luteal, and large luteal), and isolated luteal LDs were assessed for LD content, LD-associated proteins and lipid analyses. Bovine luteal tissue contained abundant lipid droplets, LD-associated perilipins 2/3/5, hormone-sensitive lipase, and 1-acylglycerol-3-phosphate O-acyltransferase ABHD5. Luteal tissue was enriched in triglycerides (TGs) compared to other tissues, except for adipose tissue. Luteal cells were distinguishable from follicular cells by the presence of LDs, LD-associated proteins, and increased TGs. Furthermore, LDs from large luteal cells were numerous and small; whereas, LDs from small luteal cells were large and less numerous. Isolated LDs contained nearly all of the TGs and cholesteryl esters present in luteal tissue. Isolated luteal LDs were composed primarily of TG, with lesser amounts of cholesteryl esters, diglyceride and other phospholipids. Bovine luteal LDs are distinct from LDs in other bovine tissues, including follicular steroidogenic cells.

Yes-associated protein 1 is required for proliferation and function of bovine granulosa cells in vitro†

Yes-associated protein 1 (YAP1) is a major component of the Hippo signaling pathway. Although the exact extracellular signals that control the Hippo pathway are currently unknown, increasing evidence supports a critical role for the Hippo pathway in embryonic development, regulation of organ size, and carcinogenesis. Granulosa cells (GCs) within the ovarian follicle proliferate and produce steroids and growth factors, which facilitate the growth of follicle and maturation of the oocyte. We hypothesize that YAP1 plays a role in proliferation and estrogen secretion of GCs. In the current study, we examined the expression of the Hippo signaling pathway in bovine ovaries and determined whether it was important for GC proliferation and estrogen production. Mammalian STE20-like protein kinase 1 (MST1) and large tumor suppressor kinase 2 (LATS2) were identified as prominent upstream components of the Hippo pathway expressed in granulosa and theca cells of the follicle and large and small cells of the corpus luteum. Immunohistochemistry revealed that YAP1 was localized to the nucleus of growing follicles. In vitro, nuclear localization of the downstream Hippo signaling effector proteins YAP1 and transcriptional co-activator with PDZ-binding motif (TAZ) was inversely correlated with GC density, with greater nuclear localization under conditions of low cell density. Treatment with verteporfin and siRNA targeting YAP1 or TAZ revealed a critical role for these transcriptional co-activators in GC proliferation. Furthermore, knockdown of YAP1 in GCs inhibited follicle-stimulating hormone (FSH)-induced estradiol biosynthesis. The data indicate that Hippo pathway transcription co-activators YAP1/TAZ play an important role in GC proliferation and estradiol synthesis, two processes necessary for maintaining normal follicle development.

Identification of a New Theca/Interstitial Cell-Specific Gene and Its Biological Role in Growth of Mouse Ovarian Follicles at the Gonadotropin-Independent Stage

Theca/interstitial cells are responsible for the growth and maturation of ovarian follicles. However, little is known about the theca/interstitial cell-specific genes and their functions. In this study, we explored transcriptomes of theca/interstitial cells by RNA-seq, and the novel biological roles of a theca cell marker, asporin (Aspn)/periodontal ligament-associated protein 1 (PLAP-1). RNA-seq detected 432 and 62 genes expressed specifically in theca/interstitial cells and granulosa cells isolated from 3-weeks old mouse ovaries. Gene ontology analysis demonstrated that these genes were largely categorized into four major groups: extracellular matrix organization-related terms, chemotaxis-related terms, the angiogenesis-related terms, and morphogenesis-related terms. In situ hybridization demonstrated that the newly detected representative gene, Aspn/PLAP-1, was detected specifically in the outer layer of theca cells in contrast with the expression of the basal lamina-specific gene, Nidgen-1. Intriguingly, an Aspn/PLAP-1 antibody completely arrested the growth of secondary follicles that is the gonadotropin-independent follicle developmental stage. Furthermore, transforming growth factor-β (TGF-β)-triggered signaling was induced by the Aspn/PLAP-1 antibody treatment, which is consistent with the inhibitory effect of Aspn/PLAP-1 on TGF-β. Altogether, these results suggest that theca cells are classified into subpopulations on the basis of new marker genes and their biological functions, and provide evidence that Aspn/PLAP-1 is expressed exclusively in the outer layer of theca cells and plays a pivotal role in the growth of secondary follicles via downregulation of the canonical TGF-β signaling cascade.

Melatonin Stimulates STAR Expression and Progesterone Production via Activation of the PI3K/AKT Pathway in Bovine Theca Cells

Melatonin is present in mammalian follicular fluid and plays an important role in regulating steroidogenesis in follicular cells. In this study, we report the effect of melatonin on steroidogenesis in the theca interna (TI) in small bovine follicles and theca cells (TCs) cultured in vitro. Treatment with melatonin significantly increased the expression of steroidogenic acute regulatory protein (STAR) and the production of progesterone in both TI and in TCs. Melatonin stimulated the phosphorylation of AKT but not ERK1/2, and the addition of luzindole (a nonspecific MT1 and MT2 inhibitor) or 4P-PDOT (specific MT2 inhibitor) reduced melatonin-induced STAR expression, progesterone secretion, and PI3K/AKT pathway activation. The effect of melatonin on the TI in follicles was more obvious than on the TCs in vitro. Results indicate that melatonin stimulates the steroidogenesis of TCs mainly via the activation of the PI3K/AKT pathway by MT1 and MT2.

Perturbations in Lineage Specification of Granulosa and Theca Cells May Alter Corpus Luteum Formation and Function

Anovulation is a major cause of infertility, and it is the major leading reproductive disorder in mammalian females. Without ovulation, an oocyte is not released from the ovarian follicle to be fertilized and a corpus luteum is not formed. The corpus luteum formed from the luteinized somatic follicular cells following ovulation, vasculature cells, and immune cells is critical for progesterone production and maintenance of pregnancy. Follicular theca cells differentiate into small luteal cells (SLCs) that produce progesterone in response to luteinizing hormone (LH), and granulosa cells luteinize to become large luteal cells (LLCs) that have a high rate of basal production of progesterone. The formation and function of the corpus luteum rely on the appropriate proliferation and differentiation of both granulosa and theca cells. If any aspect of granulosa or theca cell luteinization is perturbed, then the resulting luteal cell populations (SLC, LLC, vascular, and immune cells) may be reduced and compromise progesterone production. Thus, many factors that affect the differentiation/lineage of the somatic cells and their gene expression profiles can alter the ability of a corpus luteum to produce the progesterone critical for pregnancy. Our laboratory has identified genes that are enriched in somatic follicular cells and luteal cells through gene expression microarray. This work was the first to compare the gene expression profiles of the four somatic cell types involved in the follicle-to-luteal transition and to support previous immunofluorescence data indicating theca cells differentiate into SLCs while granulosa cells become LLCs. Using these data and incorporating knowledge about the ways in which luteinization can go awry, we can extrapolate the impact that alterations in the theca and granulosa cell gene expression profiles and lineages could have on the formation and function of the corpus luteum. While interactions with other cell types such as vascular and immune cells are critical for appropriate corpus luteum function, we are restricting this review to focus on granulosa, theca, and luteal cells and how perturbations such as androgen excess and inflammation may affect their function and fertility.

Genome-wide differences in DNA methylation changes in caprine ovaries between oestrous and dioestrous phases

Background

DNA methylation plays a vital role in reproduction. Entire genome DNA methylation changes during the oestrous phase (ES) and dioestrous phase (DS) in the ovaries of Guanzhong dairy goats were investigated using bisulphite sequencing to understand the molecular biological mechanisms of these goats’ oestrous cycle.

Results

We discovered distinct genome-wide DNA methylation patterns in ES and DS ovaries. A total of 26,910 differentially methylated regions were upregulated and 21,453 differentially methylated regions were downregulated in the ES samples compared with the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Differentially methylated region analysis showed hypomethylation in the gene body regions and hypermethylation in the joining region between upstream regions and gene bodies. The methylation ratios of the STAR, FGF2, FGF12, BMP5 and SMAD6 genes in the ES samples were lower than those of the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Conversely, the methylation ratios of the EGFR, TGFBR2, IGF2BP1 and MMD2 genes increased in the ES samples compared with the DS samples. In addition, 223 differentially methylated genes were found in the GnRH signalling pathway (KO04912), ovarian steroidogenesis pathway (KO04913), oestrogen signalling pathway (KO04915), oxytocin signalling pathway (KO04921), insulin secretion pathway (KO04911) and MAPK signalling pathway (KO04010).

Conclusions

This study is the first large-scale comparison of the high-resolution DNA methylation landscapes of oestrous and dioestrous ovaries from dairy goats. Previous studies and our investigations have shown that the NR5A2, STAR, FGF2 and BMP5 genes might have potential application value in regulating caprine oestrus.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0301-x) contains supplementary material, which is available to authorized users.

Global gene expression analysis indicates that small luteal cells are involved in extracellular matrix modulation and immune cell recruitment in the bovine corpus luteum

Genome wide mRNA expression analysis of small and large luteal cells, isolated from the mature staged corpora lutea (CL), was not performed in any species. In the current study, we have isolated bovine small and large luteal cells from mid-cycle (day 10-11) animals and characterized their transcriptomes using "GeneChip™ Bovine Gene 1.0 ST Arrays". A total of 1276 genes were identified to be differentially expressed between small and large luteal cells. Data evaluation revealed that novel functions, extracellular matrix synthesis and immune cell recruitment, were enriched in small luteal cells. On contrary, functions regarding the regulation of folliculogenesis, luteal regression, fatty acid and branched chain amino acid metabolism were differentially enriched in large luteal cells. Overall, the current data offer a first and detailed insight into the functional roles of small and large luteal cells in the mature bovine corpus luteum.

Transcriptome analysis of ovine granulosa cells reveals differences between small antral follicles collected during the follicular and luteal phases

Ovarian follicular growth occurs in both the follicular and luteal phases of the estrous cycle but in very different endocrine contexts. In both phases, many small antral follicles with similar morphologic and histologic characteristics are present within the ovaries as a reserve for the terminal folliculogenesis. However, there are several gaps in our molecular knowledge of the gene expression profiles of small antral follicles in the follicular and luteal phases. The aim of the present study was to use RNA sequencing to compare and analyze the global transcriptional profile of ovine granulosa cells collected from small antral follicles (1-3 mm) either during the follicular or the luteal phase of the estrous cycle, with the hypothesis that they should be differential. We identified 663 genes whose mRNA was differentially expressed or accumulated in the granulosa cell layer of small antral follicles in the two phases. A comprehensive interpretation of these data was performed through integrative analyses (Gene Ontology, Ingenuity Pathway Analysis) and the exploitation of already available transcriptomic data on follicular growth and atresia. In particular, we observed that the contrasted endocrine context between follicular and luteal phases may have an impact on estradiol, follicle-stimulating hormone (FSH), and on the activin/inhibin signaling pathways. Furthermore, we reveal the possible initiation of early follicular atresia in small antral follicles during the follicular phase in interaction with the presence of immune cells. This study provides new insights into the gene expression profile in ovine granulosa cells, and we suggest that these molecular changes may have an implication at the time of follicle selection.

Biological functions and role of CCN1/Cyr61 in embryogenesis and tumorigenesis in the female reproductive system (Review)

Cysteine-rich angiogenic inducer 61 (CCN1/Cyr61) is a prompt response transcription product activated by growth factors. As a member of the CCN family, it mediates cell survival, proliferation, differentiation, migration, adhesion and synthesis of the extracellular matrix by binding directly to the integrins and heparin sulfate proteoglycans or activating multiple signaling transduction pathways. It has previously been demonstrated that CCN1/Cyr61 exhibits an important role in the female reproductive system during embryogenesis and tumorigenesis. However, the functions of CCN1/Cyr61 in the female reproductive system have not been systematically investigated, therefore, the primary aim of the present review is to introduce the role and function of CCN1/Cyr61 in the female reproductive system. The current review presents the molecular structure and biological function of CCN1/Cyr61 and provides detailed data on its expression pattern and contribution to the female reproductive system, including the role in embryogenesis and tumorigenesis.

Early transcriptome responses of the bovine midcycle corpus luteum to prostaglandin F2α includes cytokine signaling

In ruminants, prostaglandin F2alpha (PGF2α)-mediated luteolysis is essential prior to estrous cycle resumption, and is a target for improving fertility. To deduce early PGF2α-provoked changes in the corpus luteum a short time-course (0.5-4 h) was performed on cows at midcycle. A microarray-determined transcriptome was established and examined by bioinformatic pathway analysis. Classic PGF2α effects were evident by changes in early response genes (FOS, JUN, ATF3) and prediction of active pathways (PKC, MAPK). Several cytokine transcripts were elevated and NF-κB and STAT activation were predicted by pathway analysis. Self-organizing map analysis grouped differentially expressed transcripts into ten mRNA expression patterns indicative of temporal signaling cascades. Comparison with two analogous datasets revealed a conserved group of 124 transcripts similarly altered by PGF2α treatment, which both, directly and indirectly, indicated cytokine activation. Elevated levels of cytokine transcripts after PGF2α and predicted activation of cytokine pathways implicate inflammatory reactions early in PGF2α-mediated luteolysis.

Transcriptomes of bovine ovarian follicular and luteal cells

Affymetrix Bovine GeneChip® Gene 1.0 ST Array RNA expression analysis was performed on four somatic ovarian cell types: the granulosa cells (GCs) and theca cells (TCs) of the dominant follicle and the large luteal cells (LLCs) and small luteal cells (SLCs) of the corpus luteum. The normalized linear microarray data was deposited to the NCBI GEO repository (GSE83524). Subsequent ANOVA determined genes that were enriched (≥2 fold more) or decreased (≤−2 fold less) in one cell type compared to all three other cell types, and these analyzed and filtered datasets are presented as tables. Genes that were shared in enriched expression in both follicular cell types (GCs and TCs) or in both luteal cells types (LLCs and SLCs) are also reported in tables. The standard deviation of the analyzed array data in relation to the log of the expression values is shown as a figure. These data have been further analyzed and interpreted in the companion article “Gene expression profiling of ovarian follicular and luteal cells provides insight into cellular identities and functions” (Romereim et al., 2017) [1].