Changes in granulosa cells gene expression associated with growth, plateau and atretic phases in medium bovine follicles

WGBS of embryonic gonads revealed that long non-coding RNAs in the MHM region might be involved in cell autonomous sex identity and female gonadal development in chickens

DNA methylation plays a key role in sex determination and differentiation in vertebrates. However, there are few studies on DNA methylation involved in chicken gonad development, and most focused on male hypermethylated regions (MHM). It is unclear whether there are specific differentially methylated regions (DMRs) in chicken embryonic gonads regulating sex determination and differentiation. Here, the DNA methylation maps showed that the difference of DNA methylation level between sexes was much higher at embryonic day 10 (E10) than that at embryonic day 6 (E6), and the significant differentially methylated regions at both stages were mainly distributed on the Z chromosome, including MHM1 and MHM2. The results of bisulphite sequencing PCR (BSP) and qRT-PCR showed hypomethylation of female MHM and upregulation of long non-coding RNAs (lncRNAs) whose promoter in the MHM region was consistent with the sequencing results, and similar results were in brain and muscle. In female sex-reversed gonads, the methylation pattern of MHM remained unchanged, and the expression levels of the three candidate lncRNAs were significantly decreased compared with those in females, but were significantly increased compared to males. The fluorescence in situ hybridization (FISH) results also showed that these lncRNAs were highly expressed in female embryonic gonads. The results of methyltransferase inhibitor and dual-luciferase reporter assay suggest that lncRNA expression may be regulated by DNA methylation within their promoters. Therefore, we speculated that MHM may be involved in cell-autonomous sex identity in chickens, and that lncRNAs regulated by MHM may be involved in female sexual differentiation.

Decoding molecular features of bovine oocyte fate during antral follicle growth via single-cell multi-omics analysis†

Antral follicle size is a useful predictive marker of the competency of enclosed oocytes for yielding an embryo following in vitro maturation and fertilization. However, the molecular mechanisms underpinning oocyte developmental potential during bovine antral follicle growth are still unclear. Here, we used a modified single-cell multi-omics approach to analyze the transcriptome, DNA methylome, and chromatin accessibility in parallel for oocytes and cumulus cells collected from bovine antral follicles of different sizes. Transcriptome profiling identified three types of oocytes (small, medium, and large) that underwent different developmental trajectories, with large oocytes exhibiting the largest average follicle size and characteristics resembling metaphase-II oocytes. Differential expression analysis and real-time polymerase chain reaction assay showed that most replication-dependent histone genes were highly expressed in large oocytes. The joint analysis of multi-omics data revealed that the transcription of 20 differentially expressed genes in large oocytes was associated with both DNA methylation and chromatin accessibility. In addition, oocyte-cumulus interaction analysis showed that inflammation, DNA damage, and p53 signaling pathways were active in small oocytes, which had the smallest average follicle sizes. We further confirmed that p53 pathway inhibition in the in vitro maturation experiments using oocytes obtained from small antral follicles could improve the quality of oocytes and increased the blastocyte rate after in vitro fertilization and culture. Our work provides new insights into the intricate orchestration of bovine oocyte fate determination during antral folliculogenesis, which is instrumental for optimizing in vitro maturation techniques to optimize oocyte quality.

Granulosa cells provide transcriptomic information on ovarian follicle dynamics in southern white rhinoceros

Much remains unknown about the reproductive physiology of southern white rhinoceros (SWR) and the effect of ovarian stimulation prior to ovum pickup (OPU) have not been fully elucidated. Granulosa cells (GC) provide valuable insight into follicle growth and oocyte maturation status. The goals of this study were to evaluate transcriptomic changes in GC from three stages of follicle development and to identify biomarkers possibly associated with follicular growth and maturation as a result of ovarian stimulation. GC collected from SWRs following OPU were assigned stages based upon follicle size. Total RNA was isolated, and cDNA libraries were prepared and sequenced on a NovaSeq 6000. All bioinformatics analyses were performed utilizing the Galaxy web platform. Reads were aligned to CerSimCot1.0, and the manual curation was performed with EquCab3.0. Overall, 39,455 transcripts (21,612 genes) were identified across follicle stages, and manual curation yielded a 61% increase in gene identification from the original annotation. Granulosa cells from preovulatory follicles expressed the highest number of unique transcripts. The following seven biomarkers were determined based upon cluster analysis and patterns of expression: COL1A1, JMY, FBXW11, NRG1, TMPO, MACIR and COL4A1. These data can be used to potentially evaluate the effects of different ovarian stimulation protocols on follicle dynamics, improve OPU results, and support conservation efforts in this species.

Developmental, cytogenetic and epigenetic consequences of removing complex proteins and adding melatonin during in vitro maturation of bovine oocytes

Background

In vitro maturation (IVM) of germinal vesicle intact oocytes prior to in vitro fertilization (IVF) is practiced widely in animals. In human assisted reproduction it is generally reserved for fertility preservation or where ovarian stimulation is contraindicated. Standard practice incorporates complex proteins (CP), in the form of serum and/or albumin, into IVM media to mimic the ovarian follicle environment. However, the undefined nature of CP, together with batch variation and ethical concerns regarding their origin, necessitate the development of more defined formulations. A known component of follicular fluid, melatonin, has multifaceted roles including that of a metabolic regulator and antioxidant. In certain circumstances it can enhance oocyte maturation. At this stage in development, the germinal-vesicle intact oocyte is prone to aneuploidy and epigenetic dysregulation.

Objectives

To determine the developmental, cytogenetic and epigenetic consequences of removing CP and including melatonin during bovine IVM.

Materials and methods

The study comprised a 2 x 2 factorial arrangement comparing (i) the inclusion or exclusion of CP, and (ii) the addition (100 nM) or omission of melatonin, during IVM. Cumulus-oocyte complexes (COCs) were retrieved from stimulated cycles. Following IVM and IVF, putative zygotes were cultured to Day 8 in standard media. RNAseq was performed on isolated cumulus cells, cytogenetic analyses (SNP-based algorithms) on isolated trophectoderm cells, and DNA methylation analysis (reduced representation bisulfite sequencing) on isolated cells of the inner-cell mass.

Results

Removal of CP during IVM led to modest reductions in blastocyst development, whilst added melatonin was beneficial in the presence but detrimental in the absence of CP. The composition of IVM media did not affect the nature or incidence of chromosomal abnormalities but cumulus-cell transcript expression indicated altered metabolism (primarily lipid) in COCs. These effects preceded the establishment of distinct metabolic and epigenetic signatures several days later in expanded and hatching blastocysts.

Conclusions

These findings highlight the importance of lipid, particularly sterol, metabolism by the COC during IVM. They lay the foundation for future studies that seek to develop chemically defined systems of IVM for the generation of transferrable embryos that are both cytogenetically and epigenetically normal.

Lin28B overexpression decreases let-7b and let-7g levels and increases proliferation and estrogen secretion in Dolang sheep ovarian granulosa cells

Although ovine puberty initiation has been previously studied, the mechanism by which the RNA-binding protein Lin28B affects this process has not been investigated. The present study aimed to investigate the effects of Lin28B overexpression on let-7b, let-7g, cell proliferation, and estrogen secretion in Dolang sheep ovine ovarian granulosa cells. In this study, a Lin28B vector was constructed and transfected into ovarian granulosa cells using liposomes. After 24, 48, and 72 h of overexpression, quantitative real-time PCR (qRT-PCR) was used for measuring let-7b and let-7g microRNA (miRNA) levels, and estrogen secretion was measured using the enzyme-linked immunosorbent assay (ELISA). A CCK-8 (Cell Counting Kit-8) kit was used for evaluating cell viability and proliferation in response to Lin28B overexpression at 24 h. The results showed that the expression of let-7b and let-7g decreased significantly after Lin28B overexpression, and the difference was consistent over different periods. The result of ELISA showed that estradiol (E2) levels significantly increased following Lin28B overexpression. Additionally, Lin28B overexpression significantly increased the cell viability and proliferation. Therefore, the Lin28B-let-7 family axis may play a key role in the initiation of female ovine puberty.

Proteomic approach-based comparison of metabolic pathways and functional activities of whey proteins derived from Guishan and Saanen goat milk

Guishan goats, a unique goat breed in Yunnan Province, have a long history and representation, but their whey protein and function remain unclear. In this study, we carried out a quantitative analysis of the Guishan and Saanen goat whey proteome using a label-free proteomic approach. A total of 500 proteins were quantified from the 2 kinds of goat whey proteins, including 463 common proteins, 37 uniquely expressed whey proteins (UEWP), and 12 differentially expressed whey proteins (DEWP). Bioinformatics analysis indicated that UEWP and DEWP were mainly involved in cellular and immune system processes, membrane, and binding. In addition, UEWP and DEWP in Guishan goats participated primarily in metabolism and immune-related pathways, whereas Saanen goat whey proteins were associated mostly with environmental information processing-related pathways. Guishan goat whey promoted the growth of RAW264.7 macrophages more than Saanen goat whey, and significantly reduced the production of nitric oxide in lipopolysaccharide-stimulated RAW264.7 cells. This study provides a reference for further understanding these 2 goat whey proteins and finding functional active substances from them.

Cumulus Cell Transcriptome after Cumulus-Oocyte Complex Exposure to Nanomolar Cadmium in an In Vitro Animal Model of Prepubertal and Adult Age

Cadmium (Cd), a highly toxic pollutant, impairs oocyte fertilization, through oxidative damage on cumulus cells (CCs). This study analysed the transcriptomic profile of CCs of cumulus-oocyte complexes (COCs) from adult and prepubertal sheep, exposed to Cd nanomolar concentration during in vitro maturation. In both age-groups, CCs of matured oocytes underwent RNA-seq, data analysis and validation. Differentially expressed genes (DEGs) were identified in adult (n = 99 DEGs) and prepubertal (n = 18 DEGs) CCs upon Cd exposure. Transcriptomes of adult CCs clustered separately between Cd-exposed and control samples, whereas prepubertal ones did not as observed by Principal Component Analysis. The transcriptomic signature of Cd-induced CC toxicity was identified by gene annotation and literature search. Genes associated with previous studies on ovarian functions and/or Cd effects were confirmed and new genes were identified, thus implementing the knowledge on their involvement in such processes. Enrichment and validation analysis showed that, in adult CCs, Cd acted as endocrine disruptor on DEGs involved in hormone biosynthesis, cumulus expansion, regulation of cell signalling, growth and differentiation and oocyte maturation, whereas in prepubertal CCs, Cd affected DEGs involved in CC development and viability and CC-oocyte communications. In conclusion, these DEGs could be used as valuable non-invasive biomarkers for oocyte competence.

The genomic response of human granulosa cells (KGN) to melatonin and specific agonists/antagonists to the melatonin receptors

Melatonin is a known modulator of follicle development; it acts through several molecular cascades via binding to its two specific receptors MT1 and MT2. Even though it is believed that melatonin can modulate granulosa cell (GC) functions, there is still limited knowledge of how it can act in human GC through MT1 and MT2 and which one is more implicated in the effects of melatonin on the metabolic processes in the dominant follicle. To better characterize the roles of these receptors on the effects of melatonin on follicular development, human granulosa-like tumor cells (KGN) were treated with specific melatonin receptor agonists and antagonists, and gene expression was analyzed with RNA-seq technology. Following appropriate normalization and the application of a fold change cut-off of 1.5 (FC 1.5, p ≤ 0.05) for each treatment, lists of the principal differentially expressed genes (DEGs) are generated. Analysis of major upstream regulators suggested that the MT1 receptor may be involved in the melatonin antiproliferative effect by reprogramming the metabolism of human GC by activating the PKB signaling pathway. Our data suggest that melatonin may act complementary through both MT1 and MT2 receptors to modulate human GC steroidogenesis, proliferation, and differentiation. However, MT2 receptors may be the ones implicated in transducing the effects of melatonin on the prevention of GC luteinization and follicle atresia at the antral follicular stage through stimulating the PKA pathway.

Glucose consumption and gene expression in granulosa cells collected before and after in vitro oocyte maturation in the southern white rhinoceros (Ceratotherium simum simum)

CONTEXT

With two northern white rhinos (NWR) remaining, the continued existence of this species relies on studying their relative, the southern white rhino (SWR).

AIMS

(1) Characterise gene expression in granulosa cells (GC) from SWR cumulus oocyte complexes (COCs) prior to (Pre-) and after (Post-) in vitro maturation (IVM), comparing culture media and oocytes from donors treated with or without gonadotropin stimulation prior to ovum recovery; and (2) evaluate COC glucose consumption in spent media.

METHODS

COCs were retrieved from four SWRs. Granulosa cells were collected before and after IVM in SDZ or IZW medium. Total RNA was evaluated by qPCR.

KEY RESULTS

Oocyte maturation was greater in SDZ than IZW media. Expression of genes associated with follicle development increased in Pre-IVM GC. Six genes were differentially expressed in Post-IVM GC from stimulated compared to unstimulated donors. COCs from stimulated animals consumed more glucose. Fifty seven percent of oocytes in SDZ medium consumed all available glucose.

CONCLUSIONS

Gene expression changed upon in vitro maturation and gonadotropin stimulation. Higher glucose availability might be needed during IVM.

IMPLICATIONS

This is the first study examining GC gene expression and COC metabolic requirements in rhinoceros, which are critical aspects to optimise IVM of rhinoceros oocytes.

Melatonin Signaling Pathways Implicated in Metabolic Processes in Human Granulosa Cells (KGN)

Female reproduction depends on the metabolic status, especially during the period of folliculogenesis. Even though it is believed that melatonin can improve oocyte competence, there is still limited knowledge of how it can modulate metabolic processes during folliculogenesis and which signaling pathways are involved in regulating gene expression. To investigate the effects of melatonin on metabolic signals during the antral stage of follicular development, human granulosa-like tumor cells (KGN) were treated with melatonin or forskolin, and gene expression was analyzed with RNA-seq technology. Following appropriate normalization and the application of a fold change cut-off of 1.5 (FC 1.5, p ≤ 0.05), 1009 and 922 genes were identified as differentially expressed in response to melatonin and forskolin, respectively. Analysis of major upstream regulators suggested that melatonin may activate PKB/mTOR signaling pathways to program the metabolism of KGN cells to support slower growth and differentiation and to prevent follicular atresia. Similarly, PKA activation through stimulation of cAMP synthesis with FSK seemed to exert the same effects as melatonin in reducing follicular growth and regulating differentiation. This study suggests that melatonin may act through PKA and PKB simultaneously in human granulosa cells to prevent follicular atresia and early luteinization at the antral stage.

The two-step process of ovarian follicular growth and maturation in mammals can be compared to a fruit ripening where quality depends on the second step

In human IVF, the main uncertainty factor impacting on success is oocyte quality, which largely depends on the follicular status at the time of collection. Decades of debate ensued to find the perfect stimulation protocol demonstrated the complexity of the ovarian response to exogenous gonadotropins and the dynamic nature of late folliculogenesis. Although several follicular markers, proteins, RNA from granulosa cells or microRNA and follicular fluid metabolites have been associated with outcome, the possibility to influence them during stimulation remains elusive. The heterogeneity of the follicle's maturity following control ovarian stimulation is also an important factor to explain average poor oocyte quality still observed today. In this review, the analogy between the apple ripening on the tree and follicular development is presented to focus the attention on a biphasic process: growth and differentiation. The molecular analysis of the progressive follicular differentiation indicates 2 competing phenomena: growth and differentiation where a delicate balance must operate from one to the other to ensure proper maturity at ovulation. As long as FSH stimulates growth, follicles remain green, and it is only when FSH is replaced by LH that the ripening process begins, and "apples" become red. Both fruits, follicles and apples, depend on a perfect timing of events to generate offspring.

Gene cascade analysis in human granulosa tumor cells (KGN) following exposure to high levels of free fatty acids and insulin

Maternal metabolic disorders such as obesity and diabetes are detrimental factors that compromise fertility and the success rates of medically assisted procreation procedures. During metabolic stress, adipose tissue is more likely to release free fatty acids (FFA) in the serum resulting in an increase of FFA levels not only in blood, but also in follicular fluid (FF). In humans, high concentrations of palmitic acid and stearic acid reduced granulosa cell survival and were associated with poor cumulus-oocyte complex (COC) morphology. Obesity and high levels of circulating FFA were also causatively linked to hampered insulin sensitivity in cells and compensatory hyperinsulinemia. To provide a global picture of the principal upstream signaling pathways and genomic mechanisms involved in this metabolic context, human granulosa-like tumor cells (KGN) were treated with a combination of palmitic acid, oleic acid, and stearic acid at the higher physiological concentrations found in the follicular fluid of women with a higher body mass index (BMI) (≥ 30.0 kg/m2). We also tested a high concentration of insulin alone and in combination with high concentrations of fatty acids. Transcription analysis by RNA-seq with a cut off for fold change of 1.5 and p-value 0.05 resulted in thousands of differentially expressed genes for each treatment. Using analysis software such as Ingenuity Pathway Analysis (IPA), we were able to establish that high concentrations of FFA affected the expression of genes mainly related to glucose and insulin homoeostasis, fatty acid metabolism, as well as steroidogenesis and granulosa cell differentiation processes. The combination of insulin and high concentrations of FFA affected signaling pathways related to apoptosis, inflammation, and oxidative stress. Taken together, our results provided new information on the mechanisms that might be involved in human granulosa cells exposed to high concentrations of FFA and insulin in the contexts of metabolism disorders.

Metabolic exchanges between the oocyte and its environment: focus on lipids

Finely regulated fatty acid (FA) metabolism within ovarian follicles is crucial to follicular development and influences the quality of the enclosed oocyte, which relies on the surrounding intra-follicular environment for its growth and maturation. A growing number of studies have examined the association between the lipid composition of follicular compartments and oocyte quality. In this review, we focus on lipids, their possible exchanges between compartments within the ovarian follicle and their involvement in different pathways during oocyte final growth and maturation. Lipidomics provides a detailed snapshot of the global lipid profiles and identified lipids, clearly discriminating the cells or fluid from follicles at distinct physiological stages. Follicular fluid appears as a main mediator of lipid exchanges between follicular somatic cells and the oocyte, through vesicle-mediated and non-vesicular transport of esterified and free FA. A variety of expression data allowed the identification of common and cell-type-specific actors of lipid metabolism in theca cells, granulosa cells, cumulus cells and oocytes, including key regulators of FA uptake, FA transport, lipid transformation, lipoprotein synthesis and protein palmitoylation. They act in harmony to accompany follicular development, and maintain intra-follicular homeostasis to allow the oocyte to accumulate energy and membrane lipids for subsequent meiotic divisions and first embryo cleavages.

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.

Flow cytometric evaluation of canine follicular cell apoptosis during the oestrous cycle

Apoptosis is the cellular mechanism of ovarian follicular atresia in mammals; the aim of this study was to examine the apoptosis-related cyclic changes in follicular cells of different-sized antral follicles throughout the oestrous cycle in canines. Ovaries were collected from 26 adult female dogs (1-4 years) following routine ovariohysterectomy. Antral follicles were classified as small, medium or large antral or preovulatory. The percentage of apoptotic cells was determined flow cytometrically using the terminal deoxynucleotidyl transferase-mediated deoxyuridine triphosphate (dUTP) DNA nick end labelling (TUNEL) assay. Apoptosis rate was quantified as the percentage of TUNEL-positive cells on a logarithmic scale. Percentages of TUNEL-positive cells obtained in the flow cytometric assay were compared among oestrous phases and follicular sizes using analysis of variance. Apoptotic follicles were observed in all types of canine follicles in different cycle phases and stages of development, possibly corresponding to the physiological process of the oestrous cycle. Both the oestrous phase and follicular size significantly influenced the apoptosis rate (p < .05). Apoptosis rate increased significantly (p < .05) as follicular development progressed. Apoptosis rate was the highest in large follicles during the oestrous phase (9.2%; p < .05) and the lowest in small follicles during the anestrus period (1.8%; p < .05). In conclusion, our results demonstrate significant differences in the apoptosis rate during the oestrous cycle related to follicle development in the canine ovary. Furthermore, flow cytometry using the TUNEL assay was found to be an effective method for detecting apoptosis in canine follicles.

Effects of bone morphogenetic protein 4, gremlin, and connective tissue growth factor on estradiol and progesterone production by bovine granulosa cells

Bone morphogenetic proteins (BMPs) are members of the transforming growth factor-β family of proteins that have been implicated in the paracrine regulation of granulosa cell (GC) function, but whether responses to BMPs change with follicular size or interact with connective tissue growth factor (CTGF) or BMP antagonists (e.g., gremlin [GREM]) to directly affect GC function of cattle is unknown. Therefore, to determine the effects of BMP4 on proliferation and steroidogenesis of GCs and its interaction with GREM or CTGF, experiments were conducted using bovine GC cultures. In vitro, BMP4 (30 ng/mL) inhibited (P < 0.05) follicle-stimulating hormone (FSH) plus insulin-like growth factor 1 (IGF1)-induced progesterone and estradiol production by large- and small-follicle GCs, but the inhibitory effect of BMP4 on estradiol production was much more pronounced in large-follicle GCs. In small-follicle GCs, BMP4 had no effect (P > 0.10) on IGF1-induced proliferation, but GREM inhibited (P < 0.05) cell proliferation and estradiol and progesterone production in IGF1 plus FSH-treated GCs. In large-follicle GCs, BMP4 (10 to 30 ng/mL) increased (P < 0.05) GC numbers and GREM (100 ng/mL) blocked this effect. In large-follicle GCs, CTGF inhibited (P < 0.05) FSH plus IGF1-induced progesterone and estradiol production, and CTGF blocked the stimulatory effect of BMP4 on GC proliferation. These results indicate that BMP4, GREM, and CTGF inhibit GC aromatase activity and progesterone production. Also, the stimulatory effect of BMP4 on GC proliferation and the inhibitory effects of BMP4 on GC steroidogenesis are more pronounced in large vs. small follicles.

FMR1 expression in human granulosa cells and variable ovarian response: control by epigenetic mechanisms

In humans, FMR1 (fragile X mental retardation 1) is strongly expressed in granulosa cells (GCs) of the female germline and apparently controls efficiency of folliculogenesis. Major control mechanism(s) of the gene transcription rate seem to be based on the rate of CpG-methylation along the CpG island promoter. Conducting CpG-methylation-specific bisulfite-treated PCR assays and subsequent sequence analyses of both gene alleles, revealed three variably methylated CpG domains (FMR1-VMR (variably methylated region) 1, -2, -3) and one completely unmethylated CpG-region (FMR1-UMR) in this extended FMR1-promoter-region. FMR1-UMR in the core promoter was exclusively present only in female GCs, suggesting expression from both gene alleles, i.e., escaping the female-specific X-inactivation mechanism for the second gene allele. Screening for putative target sites of transcription factors binding with CpG methylation dependence, we identified a target site for the transcriptional activator E2F1 in FMR1-VMR3. Using specific electrophoretic mobility shift assays, we found E2F1 binding efficiency to be dependent on CpG-site methylation in its target sequence. Comparative analysis of these CpGs revealed that CpG 94-methylation in primary GCs of women with normal and reduced efficiency of folliculogenesis statistically significant differences. We therefore conclude that E2F1 binding to FMR1-VMR3 in human GCs is part of an epigenetic mechanism regulating the efficiency of human folliculogenesis. Our data indicate that epigenetic mechanisms may control GC FMR1-expression rates.

STAT4 targets KISS1 to promote the apoptosis of ovarian granulosa cells

Background

In mammals, it is known that the estradiol-17β (E2) is mainly synthetized in ovarian granulosa cells (GCs), and the excessive apoptosis of GCs induces the follicular atresia. Many studies have implicated the essential role of KISS1, with the pro-synthetic effect of E2 and the anti-apoptotic effect on GCs, in the mammalian folliculogenesis, and several STAT4 potential binding sites were previously predicted on the promoter of KISS1 in pigs. However, the biological effects of STAT4 on GCs and the molecular regulation between STAT4 and KISS1 remained largely unknown.

Methods

Using the porcine GCs as the cellular model, the overexpression plasmid, small interfering RNA, 5′-deletion and luciferase assay were applied to investigate the molecular mechanisms for STAT4 regulating the expression of KISS1.

Results

In this study, the STAT4 negatively regulated the mRNA and protein levels of KISS1 in porcine GCs, and the mRNA level of STAT4 was observed to significantly decrease from immature to mature follicles, which was inversed with that of KISS1. The relative luciferase activity of KISS1 promoter was significantly increased with deletion of the fourth potential binding site (− 305/− 295), and ChIP further confirmed that the STAT4 bound at − 305/− 295 region of KISS1. Besides, the STAT4 significantly regulated the mRNA levels of PDK1, FOXO3 and TSC2 of PI3K signaling pathway to promote the cell apoptosis and the percentage of cells at G0/G1 phase of cell cycle in GCs. Alternatively, the STAT4 significantly decreased the mRNA levels of CYP17, 3B-HSD, 17B-33 HSD, ESR1, and ESR2, as well as the concentration of E2 in GCs. Furthermore, interfering with the expression of STAT4 was observed to significantly stimulate the pro-synthetic effect of E2 and anti-apoptotic effect of KISS1 in GCs.

Conclusions

Collectively, the STAT4 might directly target at − 305/− 295 region of KISS1 to negatively regulate the transcription of KISS1, promote the cell apoptosis via PI3K signaling pathway, suppress the synthesis of E2 through the estrogen signaling pathway in porcine GCs. These proposed works could provide useful insight in further investigations on the molecular functionalities of STAT4 and KISS1 in the folliculogenesis of mammals.

Gene analysis of major signaling pathways regulated by gonadotropins in human ovarian granulosa tumor cells (KGN)†

The female reproductive function largely depends on timing and coordination between follicle-stimulating hormone (FSH) and luteinizing hormone. Even though it was suggested that these hormones act on granulosa cells via shared signaling pathways, mainly protein kinases A, B, and C (PKA, PKB, and PKC), there is still very little information available on how these signaling pathways are regulated by each hormone to provide such differences in gene expression throughout folliculogenesis. To obtain a global picture of the principal upstream factors involved in PKA, PKB, and PKC signaling in granulosa cells, human granulosa-like tumor cells (KGN) were treated with FSH or specific activators (forskolin, SC79, and phorbol 12-myristate 13-acetate) for each pathway to analyze gene expression with RNA-seq technology. Normalization and cutoffs (FC 1.5, P ≤ 0.05) revealed 3864 differentially expressed genes between treatments. Analysis of major upstream regulators showed that PKA is a master kinase of early cell differentiation as its activation resulted in the gene expression profile that accompanies granulosa cell differentiation. Our data also revealed that the activation of PKC in granulosa cells is also a strong differentiation signal that could control "advanced" differentiation in granulosa cells and the inflammatory cascade that occurs in the dominant follicle. According to our results, PKB activation provides support for PKA-stimulated gene expression and is also involved in granulosa cell survival throughout follicular development. Taken together, our results provide new information on PKA, PKB, and PKC signaling pathways and their roles in stimulating a follicle at the crossroad between maturation/ovulation and atresia.

Regulation of the transcription factor E2F1 mRNA in ovarian granulosa cells of cattle

The E2F family of transcription factors plays an important role in the control of the cell cycle, cell proliferation, and differentiation, and their role in ovarian function is just emerging. Although some evidence suggests a possible role of E2F1 in ovarian follicular development, what regulates its production in ovarian cells is unknown. Objectives of this study were to determine whether: (i) E2F1 gene expression in granulosa cells (GCs) and theca cells (TCs) change with follicular development and (ii) E2F1 mRNA abundance in TC and GC is hormonally regulated. Using real-time PCR, E2F1 mRNA abundance in GC was 5.5-fold greater (P < 0.05) in small (SM; 1 to 5 mm) than large (LG; >8 mm) follicles, but in TC, E2F1 expression did not differ among follicle sizes. SM-follicle GC had 2.1-fold greater (P < 0.05) E2F1 mRNA than TC. In SM-follicle GC, FGF9 induced a 7.6-fold increase in E2F1 mRNA abundance; however, FGF9 did not affect (P > 0.10) abundance of E2F1 mRNA in LG-follicle TC or GC. Follicle-stimulating hormone (FSH) had no effect (P > 0.10) on E2F1 gene expression in SM- or LG-follicle GC. SM-follicle GC were concomitantly treated with insulin-like growth factor 1 (30 ng/mL), FSH (30 ng/mL), and either 0 or 30 ng/mL of FGF9 with or without 50 µM of an E2F inhibitor (E2Fi; HLM0064741); FGF9 alone increased (P < 0.05) GC numbers, whereas E2Fi alone decreased (P < 0.05) GC numbers, and concomitant treatment of E2Fi with FGF9 blocked (P < 0.05) this stimulatory effect of FGF9. Estradiol production was inhibited (P < 0.05) by FGF9 alone and concomitant treatment of E2Fi with FGF9 attenuated (P < 0.05) this inhibitory effect of FGF9. SM-follicle GC treated with E2Fi decreased (P < 0.05) E2F1 mRNA abundance by 70%. Collectively, our studies show that GC E2F1 mRNA is developmentally and hormonally regulated in cattle. Inhibition of E2F1 reduced FGF9-induced GC proliferation and attenuated FGF9-inhibited estradiol production, indicating that E2F1 may be involved in follicular development in cattle.

HIF1 driven transcriptional activity regulates steroidogenesis and proliferation of bovine granulosa cells

Hypoxia-inducible factor 1 (HIF1) is a heterodimeric transcription factor, consisting of a constitutively expressed β-subunit (HIF1B) and a regulated α-subunit (HIF1A). In the present study, we analyzed the HIF1 driven transcriptional activity in bovine granulosa cells (GC). Treatment of GC with FSH (follicle stimulating hormone) and IGF1 (insulin-like growth factor 1) resulted in the upregulation of HIF1A mRNA expression under normoxia. Immunohistochemistry of bovine ovarian sections showed distinct staining of HIF1A in the GC layer of different staged ovarian follicles. Suppression of HIF1 using echinomycin and gene knockdown procedures revealed that HIF1 transcriptionally regulates the genes associated with steroidogenesis (STAR, HSD3B and CYP19A1) and proliferation (CCND2 and PCNA) of GC. Further, our data suggest that CYP19A1, the key gene of estradiol production, is one of the plausible downstream targets of HIF1 in bovine GC as shown by gene expression, radioimmunoassay, and chromatin precipitation analysis. Based on these results, we propose that HIF1 driven transcriptional activity plays a crucial role in GC functionality, especially steroidogenesis and proliferation in developing bovine ovarian follicles.

Folliculogenesis and acquisition of oocyte competence in cows

IVF success depends on hundreds of factors and details but the oocyte quality remains the most important and problematic issue. All antral follicles contain oocytes and all of them have that have reached their full size, can be aspirated, can mature and can be fertilized in vitro. But only a few will make it to embryo unless harvested at a very specific time/status. The conditions impacting the oocyte competence are essentially dependant on the follicular status. Growing follicles contains oocytes that have not completed their preparation, as they are still writing information (RNA), later, dominant follicles or follicles at the plateau phase, stop transcription and become candidates for development. Once in transcriptional arrest, the oocytes, if not ovulated in a short amount of time, do not always make good embryos. This window is affected by time and follicle size and looks like a bell curve. The following review further explain the physiological and molecular evidences that we have to illustrate the competence window and provides clues on how to optimize ovarian stimulation to maximise oocyte quality.

Effect of heifer age on the granulosa cell transcriptome after ovarian stimulation

Genomic selection is accelerating genetic gain in dairy cattle. Decreasing generation time by using younger gamete donors would further accelerate breed improvement programs. Although ovarian stimulation of peripubertal animals is possible and embryos produced in vitro from the resulting oocytes are viable, developmental competence is lower than when sexually mature cows are used. The aim of the present study was to shed light on how oocyte developmental competence is acquired as a heifer ages. Ten peripubertal Bos taurus Holstein heifers underwent ovarian stimulation cycles at the ages of 8, 11 (mean 10.8) and 14 (mean 13.7) months. Collected oocytes were fertilised in vitro with spermatozoa from the same adult male. Each heifer served as its own control. The transcriptomes of granulosa cells recovered with the oocytes were analysed using microarrays. Differential expression of certain genes was measured using polymerase chain reaction. Principal component analysis of microarray data revealed that the younger the animal, the more distinctive the gene expression pattern. Using ingenuity pathway analysis (IPA) and NetworkAnalyst (www.networkanalyst.ca), the main biological functions affected in younger donors were identified. The results suggest that cell differentiation, inflammation and apoptosis signalling are less apparent in peripubertal donors. Such physiological traits have been associated with a lower basal concentration of LH.

KISS1 Suppresses Apoptosis and Stimulates the Synthesis of E2 in Porcine Ovarian Granulosa Cells

Previous studies have strongly recommended that KISS-1 metastasis suppressor (KISS1) plays an essential gatekeeper of the initiation of reproductive maturation in mammals. However, KISS1 has been recently reported to highly express in ovarian granulosa cells (GCs). But the biological functionalities of KISS1 on cell apoptosis, cell cycle, and synthesis of estradiol-17β (E2) have not been explored in GCs. In this study, using porcine GCs as a cellular model, the overexpression plasmid of KISS1 was built to explore the biological effects of KISS1 on the PI3K signaling pathway, estrogen signaling pathway, cell apoptosis, cell cycle, and E2 secretion. We found that mRNA of KISS1 highly expressed in the ovary and significantly increased from immature to mature follicles in gilts. Overexpression of KISS1 could significantly increase the mRNA expression of PIK3CG, PIK3C1, and PDK1, and significantly decreased the mRNA levels of FOXO3, TSC2, and BAD of PI3K signaling pathway. Furthermore, results of the flow cytometry showed that overexpression of KISS1 significantly inhibited the apoptosis of GCs and decreased the percentage of GCs at G0/G1 phase of the cell cycle. Additionally, overexpression of KISS1 could increase the mRNA levels of Star, CYP17, 3B-HSD, 17B-HSD of estrogen synthesis signaling pathway, significantly increase the concentration of E2 in the supernatant of the cultured GCs, and up-regulate the mRNA expression levels of ESR1 and ESR2. These results suggested that KISS1 might suppress cell apoptosis through activating the PI3K signaling pathway and stimulate synthesis of E2 via boosting the estrogen synthesis signaling pathway. This study would be of great interests for exploring the biological functionalities of KISS1 in the folliculogenesis and sex steroid production of the ovaries in mammals.

Transcriptome Landscape of Human Folliculogenesis Reveals Oocyte and Granulosa Cell Interactions

The dynamic transcriptional regulation and interactions of human germlines and surrounding somatic cells during folliculogenesis remain unknown. Using RNA sequencing (RNA-seq) analysis of human oocytes and corresponding granulosa cells (GCs) spanning five follicular stages, we revealed unique features in transcriptional machinery, transcription factor networks, and reciprocal interactions in human oocytes and GCs that displayed developmental-stage-specific expression patterns. Notably, we identified specific gene signatures of two cell types in particular developmental stage that may reflect developmental competency and ovarian reserve. Additionally, we uncovered key pathways that may concert germline-somatic interactions and drive the transition of primordial-to-primary follicle, which represents follicle activation. Thus, our work provides key insights into the crucial features of the transcriptional regulation in the stepwise folliculogenesis and offers important clues for improving follicle recruitment in vivo and restoring fully competent oocytes in vitro.

ASAS-SSR Triennial Reproduction Symposium: The use of natural cycle's follicular dynamic to improve oocyte quality in dairy cows and heifers

The selection of the best dairy heifers is mainly driven by the genetic value of their parents. The phenotype analysis of cows and of the daughters of bulls has been used to identify the best genetic value for decades before being replaced by genomic selection of individuals that are not yet parents. Because it is possible to predict the future value of an individual by its genetic makeup, it becomes feasible to do it as early as the blastocyst stage and to decide which should be transferred or not. Because we know the genotype of an animal at birth, or even before, it is becoming desirable to reproduce this animal as soon as possible to reduce generation interval and improve selection speed. Nature provides constraints that can be overcome: a single oocyte per cycle and age at puberty. Indeed, it is now possible to super-stimulate the ovary at any age and to start collecting oocytes at 6 mo by trans-vaginal ultrasonography. The challenge becomes the production of good eggs and embryos capable of implanting and developing into healthy calves. Our understanding of ovarian follicular physiology has been instrumental in designing stimulation protocols that may be adjusted to any physiological context including age, and even the individual animal, to obtain a good response. Therefore, the combination of procedures developed in cows to optimize oocyte quality, for example, FSH coasting, in association with in vitro fertilization and optimal culture conditions can now result in the production of several female embryos twice a month from animals 6 to12 mo of age. The transcriptomic and epigenetic analyses of embryos produced from the same females at different ages were compared and few differences were noted in particular in relation to embryo metabolism. These embryos are as good as the ones obtained from adult animals and can be produced with sexed sperm of bulls 12 mo of age. This combination of these technical optimizations with blastocyst genotyping allows the selection of a second generation within a year.

Cumulus cells surrounding oocytes with high developmental competence exhibit down-regulation of phosphoinositol 1,3 kinase/protein kinase B (PI3K/AKT) signalling genes involved in proliferation and survival

STUDY QUESTION

Is the phosphoinositol 1,3-kinase/protein kinase B (PI3K/AKT) pathway expression profile in cumulus cells (CCs) a potential marker of oocyte competence and predictive of pregnancy outcome?

SUMMARY ANSWER

Eleven genes (AKT1, ARHGEF7, BCL2L1, CCND1, E2F1, HRAS, KCNH2, PIK3C2A, SHC1, SOS1 and SPP1) in the PI3K/AKT pathway were significantly down-regulated in CCs from oocytes that went on to produce a pregnancy compared to CCs associated with a negative outcome.

WHAT IS KNOWN ALREADY

The PI3K/AKT pathway plays a pivotal role in the interdependence and continuous feedback between the oocyte and CCs.

STUDY DESIGN SIZE, DURATION

The expression analysis of 92 transcripts in the PI3K/AKT pathway in CCs from patients with negative or positive pregnancy outcome, after single embryo transfer, was performed. Mouse CCs target gene expression was conducted to associate the expression profile of PI3K/AKT pathway to oocyte developmental profile.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Fifty-five good prognosis IVF patients who had been referred to IVF or intracytoplasmic sperm injection treatment for male-factor infertility or tubal disease were enroled. CCs from single cumulus-oocyte complexes (COCs) from 16 patients who underwent a single embryo transfer were analyzed. Twenty-five CD-1 mice were used to assess gene expression in CCs associated with oocytes with different competence in relation to hCG priming. A total 220 human COCs were collected. The RNA extracted from CCs of 16 selected patients was used to analyze PI3K/AKT pathway gene expression employing a 96-well custom TaqMan Array. Expression data of CCs associated to positive IVF outcome were compared to data from negative outcome samples. Mice were sacrificed after 9, 12, 15, 21 and 24 h post-hCG administration to obtain CCs from MII oocytes with different developmental competence. Akt1, Bcl2l2 and Shc1 expression were tested in the collected mouse CCs. In addition, the expression of upstream regulator ESR1, the gene encoding for the oestrogen receptor ERβ, and the downstream effectors of the pathway FOXO1, FOXO3 and FOXO4 was evaluated in human and mouse samples.

MAIN RESULTS AND THE ROLE OF CHANCE

Transcripts involved in the PI3K Signaling Pathway were selectively modulated according to the IVF/ICSI outcome of the oocyte. Eleven transcripts in this pathway were significantly down-regulated in all samples of CCs from oocytes with positive when compared those with a negative outcome. These outcomes were confirmed in mouse CCs associated with oocytes at different maturation stages. Expression data revealed that the down-regulation of ESR1 could be related to oocyte competence and is likely to be the driver of expression changes highlighted in the PI3K/AKT pathway.

LIMITATIONS REASONS FOR CAUTION

Small sample size and retrospective design.

WIDER IMPLICATIONS OF THE FINDINGS

The CCs expression profile of PI3K/AKT signaling genes, disclosed a specific CCs gene signature related to oocyte competence. It could be speculated that CCs associated with competent oocytes have completed their role in sustaining oocyte development and are influencing their fate in response to metabolic and hormonal changes by de-activating anti-apoptotic signals.

STUDY FUNDING/COMPETING INTEREST(S)

Supported by Merck Serono an affiliate of Merck KGaA, Darmstadt, Germany (research grant for the laboratory session; Merck KGaA reviewed the manuscript for medical accuracy only before journal submission. The authors are fully responsible for the content of this manuscript, and the views and opinions described in the publication reflect solely those of the authors). The authors declare no conflict of interest.

Expression of atresia biomarkers in granulosa cells after ovarian stimulation in heifers

The use of younger gamete donors in dairy cattle genetic selection programs significantly accelerates genetic gains by decreasing the interval between generations. Ovarian stimulation (OS) and the practice of follicle-stimulating hormone (FSH) withdrawal, also known as coasting, are intensively used in pre-pubertal heifers without detrimental effects on subsequent reproductive performance but generally with lower embryo yields. However, recent data from embryo transfer programs showed similar embryo yields in younger and sexually mature animals but with a significant difference in the coasting period. The aim of the present study was to identify a set of granulosa cell biomarkers capable of distinguishing optimal follicle differentiation from late differentiation and atresia in order to assess the differences in coasting dynamics between pre- and post-pubertal donors. We integrated transcriptomic data sets from a public depository and used vote counting meta-analysis in order to elucidate the molecular changes occurring in granulosa cells during late follicle differentiation and atresia. The meta-analysis revealed the gene expression associated with follicle demise, and most importantly, identified potential biomarkers of that status in bovine granulosa cells. The comparison of the expression of six biomarkers between pre- and post-pubertal donors revealed that younger donors had more signs of atresia after the same period of coasting. We found different follicular dynamics following coasting in younger donors. It is possible that younger donors are less capable to sustain follicular survival most likely due to insufficient luteinizing hormone signaling. In summary, the pre-pubertal status influences follicular dynamics and reduces the oocyte developmental competence curve following OS and FSH withdrawal in heifers.

Initiation of follicular atresia: gene networks during early atresia in pig ovaries

In mammals, more than 99% of ovarian follicles undergo a degenerative process known as atresia. The molecular events involved in atresia initiation remain incompletely understood. The objective of this study was to analyze differential gene expression profiles of medium antral ovarian follicles during early atresia in pig. The transcriptome evaluation was performed on cDNA microarrays using healthy and early atretic follicle samples and was validated by quantitative PCR. Annotation analysis applying current database (Sus scrofa 11.1) revealed 450 significantly differential expressed genes between healthy and early atretic follicles. Among them, 142 were significantly upregulated in early atretic with respect to healthy group and 308 were downregulated. Similar expression trends were observed between microarray data and quantitative RT-PCR confirmation, which indicated the reliability of the microarray analysis. Further analysis of the differential expressed genes revealed the most significantly affected biological functions during early atresia including blood vessel development, regulation of DNA-templated transcription in response to stress and negative regulation of cell adhesion. The pathway and interaction analysis suggested that atresia initiation associates with (1) a crosstalk of cell apoptosis, autophagy and ferroptosis rather than change of typical apoptosis markers, (2) dramatic shift of steroidogenic enzymes, (3) deficient glutathione metabolism and (4) vascular degeneration. The novel gene candidates and pathways identified in the current study will lead to a comprehensive view of the molecular regulation of ovarian follicular atresia and a new understanding of atresia initiation.

MicroRNA-424/503 cluster members regulate bovine granulosa cell proliferation and cell cycle progression by targeting SMAD7 gene through activin signalling pathway

BACKGROUND

The granulosa cells are indispensable for follicular development and its function is orchestrated by several genes, which in turn posttranscriptionally regulated by microRNAs (miRNA). In our previous study, the miRRNA-424/503 cluster was found to be highly abundant in bovine granulosa cells (bGCs) of preovulatory dominant follicle compared to subordinate counterpart at day 19 of the bovine estrous cycle. Other study also indicated the involvement of miR-424/503 cluster in tumour cell resistance to apoptosis suggesting this miRNA cluster may involve in cell survival. However, the role of miR-424/503 cluster in granulosa cell function remains elusive Therefore, this study aimed to investigate the role of miRNA-424/503 cluster in bGCs function using microRNA gain- and loss-of-function approaches.

RESULTS

The role of miR-424/503 cluster members in granulosa cell function was investigated by overexpressing or inhibiting its activity in vitro cultured granulosa cells using miR-424/503 mimic or inhibitor, respectively. Luciferase reporter assay showed that SMAD7 and ACVR2A are the direct targets of the miRNA-424/503 cluster members. In line with this, overexpression of miRNA-424/503 cluster members using its mimic and inhibition of its activity by its inhibitor reduced and increased, respectively the expression of SMAD7 and ACVR2A. Furthermore, flow cytometric analysis indicated that overexpression of miRNA-424/503 cluster members enhanced bGCs proliferation by promoting G1- to S- phase cell cycle transition. Modulation of miRNA-424/503 cluster members tended to increase phosphorylation of SMAD2/3 in the Activin signalling pathway. Moreover, sequence specific knockdown of SMAD7, the target gene of miRNA-424/503 cluster members, using small interfering RNA also revealed similar phenotypic and molecular alterations observed when miRNA-424/503 cluster members were overexpressed. Similarly, to get more insight about the role of miRNA-424/503 cluster members in activin signalling pathway, granulosa cells were treated with activin A. Activin A treatment increased cell proliferation and downregulation of both miRNA-424/503 members and its target gene, indicated the presence of negative feedback loop between activin A and the expression of miRNA-424/503.

CONCLUSION

This study suggests that the miRNA-424/503 cluster members are involved in regulating bovine granulosa cell proliferation and cell cycle progression. Further, miRNA-424/503 cluster members target the SMAD7 and ACVR2A genes which are involved in the activin signalling pathway.

40 years of bovine IVF in the new genomic selection context

The development of a complex technology such as in vitro fertilization (IVF) requires years of experimentation, sometimes comparing several species to learn how to create the right in vitro environment for oocytes, spermatozoa and early embryos. At the same time, individual species characteristics such as gamete physiology and gamete interaction are recently evolved traits and must be analysed within the context of each species. In the last 40 years since the birth of Louise Brown, IVF techniques progressed and are now used in multiple domestic and non-domestic animal species around the world. This does not mean that the technology is completely matured or satisfactory; a number of problems remain to be solved and several procedures still need to be optimized. The development of IVF in cattle is particularly interesting since agriculture practices permitted the commercial development of the procedure and it is now used at a scale comparable to human IVF (millions of newborns). The genomic selection of young animals or even embryos combined with sexing and freezing technologies is driving a new era of IVF in the dairy sector. The time has come for a retrospective analysis of the success and pitfalls of the last 40 years of bovine IVF and for the description of the challenges to overcome in the years to come.

Dysregulated genes and their functional pathways in luteinized granulosa cells from PCOS patients after cabergoline treatment

Polycystic ovarian syndrome (PCOS) is a common reproductive disorder frequently associated with a substantial risk factor for ovarian hyperstimulation syndrome (OHSS). Dopamine receptor 2 (D2) agonists, like cabergoline (Cb2), have been used to reduce the OHSS risk. However, lutein granulosa cells (LGCs) from PCOS patients treated with Cb2 still show a deregulated dopaminergic tone (decreased D2 expression and low dopamine production) and increased vascularization compared to non-PCOS LGCs. Therefore, to understand the PCOS ovarian physiology, it is important to explore the mechanisms that underlie syndrome based on the therapeutic effects of Cb2. Here, LGCs from non-PCOS and PCOS patients were cultured with hCG in the absence/presence of Cb2 (n = 12). Subsequently, a transcriptomic-paired design that compared untreated vs treated LGCs within each patient was performed. After transcriptomic analysis, functions and genes were prioritized by systems biology approaches and validated by RT-qPCR. We identified that similar functions were altered in both PCOS and non-PCOS LGCs treated with Cb2; however, PCOS-treated LGCs exhibited more significant changes than non-PCOS. Among the prioritized functions, dopaminergic synapse, vascular endothelial growth factor (VEGF) signaling, apoptosis and ovarian steroidogenesis were highlighted. Finally, network modeling showed CASP9, VEGFA, AKT1, CREB, AIF, MAOA, MAPK14 and BMAL1 as key genes implicated in these pathways in Cb2 response, which might be potential biomarkers for further studies in PCOS.

Morphogenesis-related gene-expression profile in porcine oocytes before and after in vitro maturation

Mammalian oocyte maturation is achieved when oocytes reach metaphase II (MII) stage, and accumulate mRNA and proteins in the cytoplasm following fertilization. It has been shown that oocytes investigated before and after in vitro maturation (IVM) differ significantly in transcriptomic and proteomic profiles. Additionally, folliculogenesis and oogenesis is accompanied by morphogenetic changes, which significantly influence further zygote formation and embryo growth. This study aimed to determine new transcriptomic markers of porcine oocyte morphogenesis that are associated with cell maturation competence. An Affymetrix microarray assay was performed on an RNA template isolated from porcine oocytes before (n = 150) and after (n = 150) IVM. The brilliant cresyl blue (BCB) staining test was used for identification of cells with the highest developmental capacity. DAVID (Database for Annotation, Visualization, and Integrated Discovery) software was used for the extraction of the genes belonging to a cell morphogenesis Gene Ontology group. The control group consisted of freshly isolated oocytes. In total, 12,000 different transcripts were analysed, from which 379 genes were downregulated and 40 were upregulated in oocytes following IVM. We found five genes, SOX9, MAP1B, DAB2, FN1, and CXCL12, that were significantly upregulated in oocytes after IVM (in vitro group) compared with oocytes analysed before IVM (in vivo group). In conclusion, we found new transcriptomic markers of oocyte morphogenesis, which may be also recognized as significant mediators of cellular maturation capacity in pigs. Genes SOX9, MAP1B, DAB2, FN1, and CXCL12 may be involved in the regulation of the MII stage oocyte formation and several other processes that are crucial for porcine reproductive competence.

miRNA-mRNA network involved in folliculogenesis interactome: systems biology approach

At later phases of folliculogenesis, the mammalian ovarian follicle contains layers of granulosa cells surrounding an antral cavity. To better understand the molecular basis of follicular growth and granulosa cell maturation, we study transcriptome profiling of granulosa cells from small (<5 mm) and large (>10 mm) bovine follicles using simultaneous method of Affymetrix microarrays (24,128 probe sets) and RNA-Seq data sets. This study proposes a computational method to discover the functional miRNA-mRNA regulatory modules, that is, groups of miRNAs and their target mRNAs that are believed to take part cooperatively in post-transcriptional gene regulation under specific conditions. The reconstructed network was named Integrated miRNA-mRNA Bipartite Network. 277 genes and 6 key modules were disclosed through clustering for mRNA master list. The 66 genes are among the genes that belong to at least two modules. All these genes, being involved in at least one of the phenomena, namely cell survival, proliferation, metastasis and apoptosis, have an overexpression pattern (P < 0.01). For miRNA master list, a total of 172 sequences were differentially expressed (P < 0.01) between dominant (large) and each of subordinate (small) follicles. Within the follicle, these miRNAs were predominantly expressed in mural granulosa cells. Finally, predicted and validated targets of these miRNAs enriched in dominant (large) follicles were identified, which are mapped to signaling pathways involved in follicular cell proliferation, steroidogenesis, PI3K/AKT/mTOR and Ras/Raf/MEK/ERK. The identification of miRNAs and their target mRNAs and the construction of their regulatory networks may give new insights into biological procedures.

Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers

Background

The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation.

Results

The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean H_O and H_E was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures F_EH, F_VR, F_LEUT, F_ROH and F_PED was 0.129, −0.012, −0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis.

Conclusions

In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds’ development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed F_ST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3610-0) contains supplementary material, which is available to authorized users.

Differentially expressed genes and gene networks involved in pig ovarian follicular atresia

Ovarian folliculogenesis corresponds to the development of follicles leading to either ovulation or degeneration, this latter process being called atresia. Even if atresia involves apoptosis, its mechanism is not well understood. The objective of this study was to analyze global gene expression in pig granulosa cells of ovarian follicles during atresia. The transcriptome analysis was performed on a 9,216 cDNA microarray to identify gene networks and candidate genes involved in pig ovarian follicular atresia. We found 1,684 significantly regulated genes to be differentially regulated between small healthy follicles and small atretic follicles. Among them, 287 genes had a fold-change higher than two between the two follicle groups. Eleven genes ( DKK3, GADD45A, CAMTA2, CCDC80, DAPK2, ECSIT, MSMB, NUPR1, RUNX2, SAMD4A, and ZNF628) having a fold-change higher than five between groups could likely serve as markers of follicular atresia. Moreover, automatic confrontation of deregulated genes with literature data highlighted 93 genes as regulatory candidates of pig granulosa cell atresia. Among these genes known to be inhibitors of apoptosis, stimulators of apoptosis, or tumor suppressors INHBB, HNF4, CLU, different interleukins ( IL5, IL24), TNF-associated receptor ( TNFR1), and cytochrome-c oxidase ( COX) were suggested as playing an important role in porcine atresia. The present study also enlists key upstream regulators in follicle atresia based on our results and on a literature review. The novel gene candidates and gene networks identified in the current study lead to a better understanding of the molecular regulation of ovarian follicular atresia.

Comparative analysis of granulosa cell gene expression in association with oocyte competence in FSH-stimulated Holstein cows

Ovarian stimulation with exogenous FSH followed by FSH withdrawal or ‘coasting’ is an effective means of increasing the number of oocytes obtainable for the in vitro production of cattle embryos. However, the quality of the oocytes thus obtained varies considerably from one cow to the next. The aim of the present study was to gain a better understanding of the follicular conditions associated with low oocyte developmental competence. Granulosa cells from 94 Holstein cows in a commercial embryo production facility were collected following ovarian stimulation and coasting. Microarray analysis showed 120 genes expressed with a differential of at least 1.5 when comparing donors of mostly competent with donors of mostly incompetent oocytes. Using ingenuity pathway analysis, we revealed the main biological functions and potential upstream regulators that distinguish donors of mostly incompetent oocytes. These are involved in cell proliferation, apoptosis, lipid metabolism, retinol availability and insulin signalling. In summary, we demonstrated that differences in follicle maturity at collection could explain differences in oocyte competence associated with individual animals. We also revealed deficiencies in lipid metabolism and retinol signalling in granulosa cells from donors of mostly incompetent oocytes.

Transcriptome meta-analysis of three follicular compartments and its correlation with ovarian follicle maturity and oocyte developmental competence in cows

Oocyte developmental competence in superstimulated cows is dependent in part on the duration of the FSH coasting. FSH coasting refers to superstimulation with FSH (2 days of endogenous FSH following follicle ablation and 3 days of FSH injections) followed by no FSH for a specific duration. The optimal duration varies among individuals. FSH coasting appears to modulate the transcriptome of different follicular compartments, which cooperate as a single functional unit. However, the integrative effects of FSH coasting on different follicular compartments remain ambiguous. Meta-analysis of three independent transcriptome studies, each focused on a single cell type (granulosa, cumulus, and oocyte) during FSH coasting, allowed the identification of 12 gene clusters with similar time-course expression patterns in all three compartments. Network analysis identified HNF4A (involved in metabolic functions) and ELAVL1 (an RNA-binding protein) as hub genes regulated respectively upward and downward in the clusters enriched at the optimal coasting time, and APP (involved in mitochondrial functions) and COPS5 (a member of the COP9 signalosome) as hub genes regulated respectively upwards and downwards in the clusters enriched progressively throughout the coasting period. We confirmed the effects on HNF4A downstream targets (TTR, PPL) and other hub genes (ELAVL1, APP, MYC, and PGR) in 30 cows with RT-quantitative PCR. The correlation of hub gene expression levels with FSH coasting indicated that a combination of these genes could predict oocyte competence with 83% sensitivity, suggesting that they are potential biomarkers of follicle differentiation. These findings could be used to optimize FSH coasting on an individual basis.

MicroRNA-183-96-182 Cluster Regulates Bovine Granulosa Cell Proliferation and Cell Cycle Transition by Coordinately Targeting FOXO1

Large-scale expression profiling of micro-RNAs (miRNAs) in bovine granulosa cells from dominant and subordinate follicles on Day 19 of the estrous cycle revealed enriched micro-RNA-183-96-182 cluster miRNAs in preovulatory dominant follicles that coordinately regulate the forkhead box protein O1 (FOXO1) gene. However, little is known about the role of this cluster in bovine granulosa cell function. We used an in vitro granulosa cell culture model to investigate this role. Granulosa cells aspirated from small growing follicles (3-5 mm in diameter) were cultured in Dulbecco modified Eagle medium/F-12 medium supplemented with fetal bovine serum and transfected with locked nucleic acid-based miRNA mimics, inhibitors, and corresponding negative controls. Overexpression of the miRNA cluster resulted in suppression of FOXO1 mRNA and protein, whereas inhibition of the cluster increased expression of FOXO1 mRNA. Overexpression also increased the relative rate of cell proliferation, whereas inhibition slowed it down. Similarly, the proportion of cells under G0/G1 arrest declined, whereas the ratio of cells in S phase increased in response to miR-183-96-182 overexpression. Selective knockdown of FOXO1 mRNA using anti-FOXO1 small interfering RNA increased the rate of granulosa cell proliferation, decreased the proportion of cells under G0/G1 arrest, and increased the proportion of cells in the S phase of cell cycle. Our data suggest that miR-183-96-182 cluster miRNAs promote proliferation and G1/S transition of bovine granulosa cells by coordinately targeting FOXO1, suggesting a critical role in granulosa cell function. MicroRNA-183-96-182 cluster regulates bovine granulosa cell function by targeting FOXO1 gene.

Gene expression patterns in granulosa cells and oocytes at various stages of follicle development as well as in in vitro grown oocyte-and-granulosa cell complexes

Follicle development is accompanied by proliferation of granulosa cells and increasing oocyte size. To obtain high-quality oocytes in vitro, it is important to understand the processes that occur in oocytes and granulosa cells during follicle development and the differences between in vivo and in vitro follicle development. In the present study, oocytes and granulosa cells were collected from early antral follicles (EAFs, 0.5–0.7 mm in diameter), small antral follicles (SAFs, 1–3 mm in diameter), large antral follicles (LAFs, 3–7 mm in diameter), and in vitro grown oocyte-and-granulosa cell complexes (OGCs), which were cultured for 14 days after collection from EAFs. Gene expression was analyzed comprehensively using the next-generation sequencing technology. We found top upstream regulators during the in vivo follicle development and compared them with those in in vitro developed OGCs. The comparison revealed that HIF1 is among the top regulators during both in vivo and in vitro development of OGCs. In addition, we found that HIF1-mediated upregulation of glycolysis in granulosa cells is important for the growth of OGCs, but the cellular metabolism differs between in vitro and in vivo grown OGCs. Furthermore, on the basis of comparison of upstream regulators between in vivo and in vitro development of OGCs, we believe that low expression levels of FLT1 (VEGFA receptor), SPP1, and PCSK6 can be considered causal factors of the suboptimal development under in vitro culture conditions.

Dibutyl Phthalate Inhibits the Effects of Follicle-Stimulating Hormone on Rat Granulosa Cells Through Down-Regulation of Follicle-Stimulating Hormone Receptor

Dibutyl phthalate (DBP) is used worldwide in solvents and plasticizers. The cytotoxicity and potential tumorigenic effect of DBP have been reported. DBP has also been shown to impact reproductive function. In this study, to further evaluate the effects of DBP on granulosa cells (GCs), we treated rat GCs in vitro with DBP before evaluation of the biological alterations of these GCs. We found that DBP did not induce significant GC death at the tested concentrations. However, follicle-stimulating hormone (FSH)-induced KIT ligand (KITLG) expression in GCs was significantly reduced at both mRNA and protein levels by DBP treatment in a dose-dependent manner. The down-regulation of KITLG was due to the down-regulation of expression of FSH receptor (FSHR) in GCs. Down-regulation of FSHR impaired FSH-induced intracellular signaling in GCs, demonstrated by decreased phosphorylation of AKT and mechanistic target of rapamycin (mTOR). Furthermore, DBP treatment also reduced FSH-induced expression of hypoxia-inducible factor 1-alpha (HIF1A), which is an important signaling component for KITLG expression. Other FSH-induced biological effects, such as production of estradiol and progesterone, as well as GC proliferation, were also suppressed by DBP. Therefore, our study discovered a unique mechanism underlying the toxicity of DBP on GCs. These findings may initiate the development of novel therapeutic interventions for DBP-induced damage to GCs.

Global gene expression in granulosa cells of growing, plateau and atretic dominant follicles in cattle

BACKGROUND

The physiological state of the dominant follicle is important as it may be linked to the competence of the oocyte within. The objective of this study was to analyze, by transcriptomic analysis, the changes occurring in granulosa cells from dominant follicles at different phases of follicular growth.

METHODS

Granulosa cells were collected from slaughterhouse dairy cattle follicles with a diameter greater than 9 mm, and were classified at different phases of follicle growth based on flow cytometry profiles of DNA content after staining with propidium iodide. Three phases were identified based on the proportion of cells in -G1 (less than 2n DNA), G0-G1 (2n DNA) or S-M (more than 2n DNA) and follicles were thus allocated to the growing, plateau or atresia group. Between group analysis (BGA) showed clear segregation of the three groups, and the groups were contrasted against each other in a loop design to identify differently expressed genes. Ingenuity Pathway Analysis (IPA) was used to identify the functions and upstream regulators associated with the observed differently expressed genes.

RESULTS

Major differences were observed between the growth phases. Granulosa cells from follicles in the plateau phase had increased expression of TYRO3 and downregulation of JAM2 compared to growing follicles, supporting the idea of a shift from proliferation to differentiation. On the other hand, genes regulating the response to oxidative stress (VNN1) and angiogenesis (ANGPT2) were upregulated in granulosa cells from atretic follicles. While the predicted activated functions in cells at the plateau stage compared to cells at the growing stage included synthesis and transport of molecules, the predictions for atretic follicles relative to plateau ones included an increase in apoptosis and cell death.

CONCLUSION

Consistent with previous studies, these observations allowed us to match the presence of specific gene transcripts to a particular physiological status and consequently to classify follicles. The results also demonstrated that the plateau phase is not a simple 'in between' status between growth and atresia, as several characteristics are unique to this stage.