Skip to main content

DNA methylome profiling of granulosa cells reveals altered methylation in genes regulating vital ovarian functions in polycystic ovary syndrome



Women with polycystic ovary syndrome (PCOS) manifest a host of ovarian defects like impaired folliculogenesis, anovulation, and poor oocyte quality, which grossly affect their reproductive health. Addressing the putative epigenetic anomalies that tightly regulate these events is of foremost importance in this disorder. We therefore aimed to carry out DNA methylome profiling of cumulus granulosa cells and assess the methylation and transcript expression profiles of a few differentially methylated genes contributing to ovarian defects in PCOS. A total of 20 controls and 20 women with PCOS were selected from a larger cohort of women undergoing IVF, after carefully screening their sera and follicular fluids for hormonal and biochemical parameters. DNA extracted from cumulus granulosa cells of three women each, from control and PCOS groups was subjected to high-throughput, next generation bisulfite sequencing, using the Illumina HiSeq 2500® platform. Remaining samples were used for the validation of methylation status of some identified genes by pyrosequencing, and the transcript expression profiles of these genes were assessed by quantitative real-time PCR.


In all, 6486 CpG sites representing 3840 genes associated with Wnt signaling, G protein receptor, endothelin/integrin signaling, angiogenesis, chemokine/cytokine-mediated inflammation, etc., showed differential methylation in PCOS. Hypomethylation was noted in 2977 CpGs representing 2063 genes while 2509 CpGs within 1777 genes showed hypermethylation. Methylation differences were also noted in noncoding RNAs regulating several ovarian functions that are dysregulated in PCOS. Few differentially methylated genes such as aldo-keto reductase family 1 member C3, calcium-sensing receptor, resistin, mastermind-like domain 1, growth hormone-releasing hormone receptor and tumor necrosis factor, which predominantly contribute to hyperandrogenism, premature luteolysis, and oocyte development defects, were explored as novel epigenetic candidates in mediating ovarian dysfunction. Methylation profiles of these genes matched with our NGS findings, and their transcript expression patterns correlated with the gene hypo- or hypermethylation status.


Our findings suggest that the epigenetic dysregulation of genes involved in important processes associated with follicular development may contribute to ovarian defects observed in women with PCOS.


Polycystic ovary syndrome (PCOS), one of the leading causes of anovulatory infertility, is characterized by ovarian, neuroendocrine, and metabolic perturbations in women of reproductive age group. With a global prevalence of 6–15% [1], it generally features irregularity or absence of menses and presence of multiple ovarian cysts on ultrasonography, in addition to hyperandrogenemia and systemic insulin excess. Increased pulsatility of GnRH neurons at the hypothalamo-hypophyseal interface elevates gonadotropic secretion of LH over FSH, which alongside co-gonadotropic actions of insulin, promotes increased production of androgens from ovarian folliculo-thecal cells [2, 3]. Adding to the pool of free androgens in circulation is the insulin-mediated suppression of hepatic sex hormone binding globulin (SHBG), which increases the bioavailability of testosterone, thereby affecting its clearance [4]. The cumulative actions of all these events trigger a series of physiological defects including ovarian cyst formation, amenorrhoea, anovulation, infertility, hyperandrogenism, insulin resistance, hyperinsulinemia, obesity, glucose intolerance, lipid abnormalities, type 2 diabetes mellitus (T2DM), hypertension, and cardiovascular disease. In parallel, there is an intrinsic elevation of anti-Mullerian hormone (AMH) due to the presence of cystic follicles arrested in preantral to antral stages [5]. These key events dictate the principal dogma behind the pathophysiology of PCOS understood so far.

Although genetic factors impacting the development and progression of PCOS have been amply investigated, identification and experimental corroboration of cognate epigenetic factors that may contribute to the pathophysiology of this multifaceted disorder remain enigmatic. Both environmental and physiological factors serve as strong determinants for epigenetic alterations. Factors such as intrinsic hormonal aberrations, dysregulation of intrauterine milieu by endocrine disruptors (EDCs) during gestational periods, and lifestyle modifications in subsequent phases of growth and development have been recently implicated in epigenetic predisposition to this disease. Pioneering investigations on epigenetic alterations in PCOS began with studies on peripheral blood leukocytes (PBLs) [6]. However, tissue specificity of epigenetic modifications renders it difficult to extrapolate epigenetic data derived from circulating cells like PBLs, to organs such as ovaries or adipose tissues that are highly affected in PCOS [7]. This necessitates the undertaking of clinical epigenetic studies at a tissue-specific level. Ovary withstands most of the hormonal assaults triggered by systemic aberrations in the neuroendocrine-ovarian axis. It is therefore a primary hotspot for epigenetic perturbations, which may contribute to the multiple follicular and oocyte defects observed in women with PCOS. Defects related to steroidogenesis, follicular growth and dominance, ovulation, oocyte developmental competence, cumulus-oophorus complex (COC) expansion, luteal maintenance, etc., which are under stringent control of gonadotropins and other hormones, have been well documented so far. Therefore, identification of locus/gene-specific epigenetic alterations in the ovaries of these women is of prime significance to understand the pathophysiology of this disorder.

So far, few studies conducted to identify global DNA methylation differences in PBLs and mural as well as cumulus granulosa cells (CGCs) of women with PCOS have yielded ambiguous findings [6, 8, 9]. Also, promoter methylation profiles of a few established candidate genes of PCOS including yes-associated protein (YAP1), follistatin (FST), aromatase (CYP19A1) and  luteinizing hormone chorionic gonadotropin receptor (LHCGR) have been investigated in these cells, and ovarian tissues by some groups till date [10,11,12,13,14]. Among these genes, LHCGR has been consistently reported to be hypomethylated in women with PCOS and in animal models of PCO [14, 15]. Subsequently, a few high-throughput attempts were made to identify some differentially methylated genes (DMGs) in women with PCOS. These included the use of diverse approaches such as methylated DNA immunoprecipitation (MeDIP) [16] and Illumina platform-based methylation microarray [17] in peripheral blood leukocytes; MeDIP coupled with methyl promoter enrichment microarray [18], as well as methylation microarray combined with microarray-based transcriptome analysis [19] in ovarian tissue biopsies and adipose tissue samples [20]; total RNA sequencing (RNA-seq) coupled with methylation measured by base cleavage and mass spectrometry (EpiTYPER) [21]; and lastly, methylation microarray in mural granulosa cells (MGCs) collected from women undergoing controlled ovarian hyperstimulation (COH) [22]. However, next generation sequencing (NGS)-based methylome studies spanning individual CpG sites, specifically in CGCs which participate in extensive cross talk between the developing oocyte and surrounding follicular milieu while facilitating meiotic maturation and ovulation of competent oocytes, are yet to be reported in women with PCOS. In this study, we have carried out comparative genome-wide bisulfite sequencing of CGCs obtained from women with PCOS and normovulatory healthy controls, using a NGS-based, multiplexed Methyl-Capture Sequencing (MC-Seq) approach. MC-Seq is advantageous over reduced representative bisulfite sequencing (RRBS) and MeDIP-Seq, in that it avoids the over-representation of recurring reads and also over Infinium 450 K methyl-microarray wherein it enables the user to opt for much greater coverage of the epigenome (3.4 million CpG sites and 20.8 million non-CpG sites, as opposed to 4,50,000 CpG sites and 3091 non-CpG loci, respectively) [23]. Also, Infinium microarrays target only those genes that are known to be differentially methylated in some cancers while MC-Seq differentiates between any genomic region bearing altered methylation marks. Additionally, compared to whole-genome bisulfite sequencing (WGBS), MC-Seq is cost effective and can identify novel genomic loci while reducing the processing time associated with WGBS. This study provides insights on genome-wide data on CpG sites that show altered methylation at a single base resolution in CGCs of women with PCOS.


Hormonal characterization and oocyte quality parameters of the study participants

All study participants were subjected to stringent anthropometric, hormonal, and biochemical characterization prior to further investigation (Table 1). Baseline hormonal and biochemical profiling of 20 controls and 20 women with PCOS (days 3–5 of the follicular phase of menstrual cycle) showed that the baseline levels (day 3 follicular phase estimates in serum) of luteinizing hormone (LH), LH/follicle-stimulating hormone (FSH) ratio, and of AMH were high, while the levels of FSH were low in women with PCOS, compared to controls. Prolactin and thyroid-stimulating hormone (TSH) levels were similar in both groups. Serum levels of estradiol (E2) and progesterone (P4) measured before the administration of recombinant human chorionic gonadotropin (rhCG) were comparable between controls and PCOS. Serum E2 levels on the day of rhCG administration were high in women with PCOS while P4 levels were similar between these groups. Analysis of oocyte parameters revealed that the total numbers of follicles, oocytes (mature + immature), mature MII oocytes, %MII oocytes, and number of fertilized oocytes were unchanged between controls and women with PCOS. However, the rates of fertilization of MII oocytes were low in women with PCOS compared to controls. In serum and follicular fluid (FF) samples collected on the day of oocyte pick up (d-OPU), total testosterone (TT) levels were high and sex hormone binding globulin (SHBG) levels were low in PCOS. Androgen-excess indices such as free testosterone (Free-T), bioavailable testosterone (Bio-T), and free androgen index (FAI) were high in FF, while in serum, only Bio-T and FAI were high in PCOS.

Table 1 Clinical characteristics of study participants undergoing controlled ovarian hyperstimulation (COH) assessed before and after oocyte retrieval

Identification of differentially methylated targets and their gene ontology analysis

MC-Seq of CGCs of women with PCOS and controls identified a total of 6486 differentially methylated CpG sites associated with 3403 unique genes across the genome, of which 2977 CpG sites were hypomethylated and 2509 CpG sites were hypermethylated. Hypomethylated CpG sites were representative of 2063 (Additional file 1) genes in all, while the hypermethylated sites were linked to a total of 1777 genes (Additional file 2). Of the total DMGs, 438 genes harbored both hyper- and hypomethylated CpG sites. Additionally, many noncoding RNAs including 44 microRNAs (miRs) and 121 pseudogenes also showed differential methylation in women with PCOS (Additional files 1 and 2). Of the 44 differentially methylated miRs, several miRs (miR23A, miR127, miR10B, miR193A mir200B, miR182, and miR140) have been reported to be implicated in impaired follicle growth and steroidogenesis, anovulation, obesity, glucose metabolism, and so on [24]. In pathway enrichment and gene ontology (GO) analyses of the hypomethylated, hypermethylated, and combined gene lists, 227 (11%), 210 (11.82%), and 396 (11.64%) genes from the three respective categories could not be annotated. Among the identified pathways, those for Wnt signaling (Panther-GO-ID, P00057), integrin signaling (P00034), endothelin signaling (P00019), and cadherin signaling (P00012) were enriched in both the hypo- and hypermethylated gene sets (Fig. 1). Other prominent pathways included the platelet-derived growth factor (PDGF) signaling (P00047), inflammation mediated by chemokine and cytokine signaling (P00031), angiogenesis (P00005) and vascular endothelial growth factor (VEGF) signaling (P00056), fibroblast growth factor (FGF) signaling (P00021), G protein signaling (P00026, P00027), T cell activation (P00053), and nicotinic acetylcholine receptor signaling pathways (P00044).

Fig. 1
figure 1

Pathway analysis of differentially methylated genes identified by MC-Seq of cumulus granulosa cells (CGCs). Pie charts displaying gene ontology (GO) and pathway analyses for the lists of hypomethylated genes (a), hypermethylated genes (b), and all differentially methylated genes (DMGs) (c) between CGCs of controls (C-CC, n = 3) and women with PCOS (P-CC, n = 3) have been shown in the figure. Analysis was carried out using the GeneCodis3 web tool. Each chart represents the top ten pathways enriched in each of the three datasets

Validation of genes showing differential methylation in regions upstream to transcription start sites

Among the DMGs selected for validation, the upstream CpG sites of five genes, namely aldo-keto reductase 1 family C3 (AKR1C3), calcium-sensing receptor (CASR), growth hormone-releasing hormone receptor (GHRHR), resistin (RETN), and mastermind-like domain 1 (MAMLD1) were hypomethylated while those of transferrin (TF) and tumor necrosis factor (TNF) were hypermethylated in PCOS in our methylome analysis. We first investigated whether these seven genes are expressed in CGCs, and whether they are differentially expressed in women with PCOS using qPCR (Fig. 2). Apart from GHRHR, which was expressed only at baseline levels, the remaining six genes were abundantly expressed in CGCs. Transcripts of AKR1C3, CASR, GHRHR, RETN, and MAMLD1 were upregulated while those of TF and TNF were downregulated in CGCs of PCOS women (Fig. 2). To verify whether our NGS findings were replicative in a larger study cohort, we performed pyrosequencing to assess the average percent methylation of these genes in 17 controls and 17 women with PCOS. Hypomethylation of AKR1C3, CASR, GHRHR, RETN, and MAMLD1 genes at the indicated CpG sites was confirmed by pyrosequencing (Fig. 3). Decreased methylation in these genes was consistent with the upregulation of their respective transcripts in women with PCOS (Fig. 2), and the two variables showed an inverse correlation with one another (Table 2). Also, hypermethylation of TNF was correlated with downregulation of its transcript expression in women with PCOS (Table 2). However, TF which showed hypermethylation in NGS data was found to be hypomethylated upon pyrosequencing, though its transcript was downregulated in PCOS (Figs. 2 and 3). Also, there was no correlation between TF transcript levels and its methylation status (Table 2). As supporting evidence, we also evaluated the transcript expression of additional genes such as prostaglandin E receptor (PTGER1), leukemia inhibitory factor (LIF), and hyaluronan and proteoglycan link protein 1 (HAPLN1) which showed hypermethylation in NGS analysis. Transcripts of these genes were downregulated in PCOS (Fig. 2).

Fig. 2
figure 2

Transcript expression profiling of few differentially methylated genes identified by MC-Seq. Bar graphs representing transcript expression profiles of 10 differentially methylated genes (DMGs) between CGCs (C-CCs) of controls (n = 17) and women with PCOS (P-CCs, n = 17) are shown in the figure. Panel (a) represents transcript expression data of genes identified as hypomethylated in PCOS by NGS analysis while (b) shows transcript profiling of hypermethylated genes. Fold change was evaluated using the 2-ΔΔCt method, where the expression was normalized to 18s levels, using a CGC calibrator sample. Data are presented as “mean + SEM.” *P < 0.05, **P < 0.01, 'ns' denotes no significant change. Data are analyzed using the Mann-Whitney U test.

Fig. 3
figure 3

Validation of CpG methylation status of a few differentially methylated genes identified by MC-Seq. The figure depicts “box and whisker” plots for genes whose single CpG site was validated or “bar graphs” for genes in which validation was performed for CpG sites > 1. The plots compare percent (%) methylation for the stated CpG sites between cumulus granulosa cells (CGCs) of controls (C-CC, n = 17) and women with PCOS (P-CC, n = 17). Data for box plots are presented as whiskers ranging from minimum to maximum values and data for bar graphs are presented as “mean + SEM” using the Mann-Whitney U test. *P <  0.05, **P <  0.01, 'ns' denotes no significant change.

Table 2 Correlation analysis of CpG methylation levels of selected genes with their transcript expression profilesᅟ


Tissue-specific DNA methylation changes sired by alterations in the environmental or physiological milieu of an individual can bring about significant changes in gene and protein expression, and therefore predispose them to disease development. Alterations in both transcriptome and proteome profiles of ovarian cells/tissues and FF have been previously reported in PCOS [25,26,27,28]. With this background information, we had initially conducted a pilot study to screen for the tissue-specific global DNA methylation changes in PBLs and CGCs of controls and women with PCOS. Here, subtle alterations were detected in CpG methylation profiles of long interspersed nucleotide element 1 (LINE1) in PCOS, and these changes were found to be more prominent in CGCs of women with PCOS, compared to PBLs [8]. Since LINE1s are self-replicating transposons and occupy ~ 17% of the human genome, even slight changes in their methylation patterns can be reflective of genomic dysregulation. Therefore, the primary goal of this study was to identify the genome-wide methylation differences in CGCs of women with PCOS, at a single base resolution. GO analysis of the current methylome data revealed that genes regulating cell growth, adhesion, differentiation, proliferation, cell polarity and fate determination, apoptosis, signal transduction, transcription, post-translational modifications, protein binding, metal and nonmetal ion binding, ATP binding, vesicular transport, etc., were differentially methylated in PCOS (Additional file 3). The implications of differential methylation observed in few of these identified genes, which may contribute to hyperandrogenism, defects in COC expansion, oocyte maturation/and ovulation, premature luteolysis, and oxidative stress observed in PCOS, have been discussed here.

Androgen overproduction

CpG hypomethylation in genes such as AKR1C3, GHRHR, MAMLD1 and RETN, and hypermethylation in TNF, which can indirectly contribute to androgen excess, were consistent with increased and decreased levels of the respective gene transcripts (Figs. 2 and 3). AKR1C3 is a steroidogenic enzyme that converts androstenedione (A4) to biologically active testosterone in non-testicular tissues [29]. In PCOS, the increased expression of AKR1C3 and AKR1C3-mediated androgen production, have been reported in the adrenal cortex and visceral adipose tissues [29, 30]. This enzyme is also expressed by GCs (both mural and cumulus granulosa cells) of periovulatory follicles [31]. Since AKR1C3 expression was found to be high in CGCs of women with PCOS, it may contribute to the high androgen production in their ovaries. Next, GHRHR is a gene encoding the class B GPCR subfamily receptor, which regulates the release of somatotropin (GH) in the brain and other tissues including the ovary, via binding to the hypothalamic neuropeptide GHRH [32, 33]. Although the classical ovarian function of GH is to enhance sexual maturation at puberty via binding to its receptors (GHRs) [34], the GH-GHR interaction also stimulates the release of insulin-like growth factor (IGF1) via transcriptional activation [34]. IGF1, like insulin, can increase LH production from the pituitary and augment ovarian androgen synthesis in PCOS [35, 36]. Additionally, hyperinsulinemia in PCOS is known to increase the bioavailability of IGF1 via the downregulation of its carrier protein, i.e., IGFBP1 [25, 34]. High levels of GH and IGF1 also increase the sensitivity of developing follicles to gonadotropins [37, 38], and PCOS follicles have been reported to exhibit increased sensitivity and responsiveness to FSH [39]. As a result, the expression of LHCGR which is under the direct control of FSH is found to be elevated in the follicles of women with PCOS [13]. This can further augment LH-mediated ovarian androgen production in PCOS. Thus, hypomethylation and overexpression of GHRHR observed in CGCs can be an indirect mediator of androgen excess in PCOS (Fig. 4). Further, the pro-inflammatory cytokine TNF, which suppresses FSH-induced LHCGR promoter activation via NF-κB p65, was hypermethylated and low in PCOS [40], thus also making it an additional factor contributing to hyperandrogenemia.

Fig. 4
figure 4

Altered gene methylation can contribute to ovarian dysfuncion. The figure summarizes some of the processes that could be dysregulated in the follicular compartment of women, due to hypomethylation (indicated by green boxes) or hypermethylation (indicated by red boxes) of genes in cumulus granulosa cells (CGCs). AKR1C3, aldo-keto reductase family 1 member C3; AR, androgen receptor; CASR, calcium-sensing receptor; COC, cumulus-oophorus complex; FF, follicular fluid; FSH, follicle-stimulating hormone; FSHR, follicle-stimulating hormone receptor; GC, granulosa cells; IGF1; insulin-like growth factor 1; INS, insulin; MGC, mural granulosa cells; GHRHR, growth hormone-releasing hormone receptor; HAPLN1, hyaluronan and proteoglycan link protein 1; LH, luteinizing hormone; LHCGR, luteinizing hormone chorionic gonadotropin receptor; LIF, leukemia inhibitory factor; MAMLD1, mastermind-like domain containing 1; PTGER1, prostaglandin receptor E1; RETN, resistin; TC, theca cells; TF, transferrin; TNF, tumor necrosis factor alpha

Mutations in MAMLD1, a transcriptional coactivator, have been reported to result in compromised androgen synthesis during male fetal sexual development [41, 42]. However, information regarding a definite role of MAMLD1, its regulation and mechanisms of action in the context of other reproductive functions, is limited. In a murine study, MAMLD1 knockout male mice showed a reduced testosterone production in Leydig tumor cells while the activation of MAMLD1 promoter by the transcription factor, SF1, augmented testosterone production via transactivation of the hairy/enhancer of split 3 HES3 promoter [43]. We therefore propose that MAMLD1 hypomethylation and upregulation of its transcript in CGCs may be important in contributing to androgen excess in PCOS ovaries (Fig. 4). Next, increased circulatory levels of the adipokine, RETN, which modulates glucose tolerance and insulin action, have been linked to a higher incidence of insulin resistance and PCOS [44, 45]. However, there is some ambiguity regarding the role of RETN in PCOS, since its serum and FF levels have been either found to be high or unchanged in PCOS [46,47,48]. In theca cells, the dose-dependent increase in RETN showed augmented androgen production [49], while a RETN-like molecule β impaired the glucose tolerance and insulin actions in HEK293T cells and adipocytes [50]. Thus, the overexpression of RETN due to the hypomethylation of its promoter may be an important factor contributing to the androgen excess in PCOS (Fig. 4).

Oocyte development, ovulation, and COC matrix expansion defects

Enrichment of calcium (Ca2+) signaling pathway and pathways for regulation of cytoskeletal and focal adhesion elements in our NGS analysis indicated impairment of calcium homeostasis and cellular architecture in CGCs of women with PCOS (Fig. 1, Additional file 3). Ca2+ signaling pathways are crucial for the development and maturation of healthy oocytes and CASR, which is an important mediator of this pathway, responds to subtle changes in extracellular Ca2+ concentrations and activates or ameliorates the mobilization of intracellular Ca2+ stored in tissues [51]. Expression and localization of CASR have been reported in human oocytes and CGCs, wherein it supposedly facilitates bidirectional communication between these cells to either keep oocytes arrested in MI phase, or assist in the full resumption of their cytoplasmic and nuclear maturation upon entering the MII phase [52]. Therefore, alterations in CASR expression may affect oocyte maturation and yield poor quality oocytes as seen in PCOS (Fig. 4). In PCOS, so far, a single report exists on the association of a CASR polymorphism (Hin1I) with altered global calcium homeostasis [53]. Altered methylation has been previously reported in human CASR promoter in a few cancer conditions [54,55,56]. Our results indicate that altered CASR expression in CGCs of PCOS women can be also influenced by altered methylation.

TNF expression in CGCs has been reported to be either unchanged or reduced in CGCs of women with PCOS [57, 58]; however, its circulating levels have been found to be high in their serum and FF [59, 60]. Few studies demonstrated that treatment with high levels of TNF increased GC apoptosis, impaired P4 production from GCs, and caused other steroidogenic defects in these cells [61,62,63]. However, these studies utilized TNF at 10–20-fold higher doses relative to its physiological levels. In alternate studies, optimal TNF levels have been reported to impart a protective function in the maintenance of bovine GCs and oocytes [64], facilitate ovulation [65], and increase the GC proliferation in animal models [66]. Decreased levels of endogenous TNF in GCs has been attributed to diminished oocyte competence due to a reduction in its downstream effector, i.e. tumor necrosis factor-inducible gene 6 (TNFAIP6) [25], as well as compromised ovulation, and GC proliferation in ovarian follicles [67]. Thus, lowered TNF in CGCs of women with PCOS owing to hypermethylation, may hamper COC expansion and compromise ovulation. LIF, which was hypermethylated and downregulated in our study, has been reported to be low in FF and serum of women with PCOS [68, 69]. LIF has demonstrated embryotrophic effects in mice and humans [70, 71], and decrease in its levels has shown a positive association with low implantation rates and poor IVF outcome [68, 69]. Since dose-dependent administration of LIF also showed induction of COC expansion and improvement in oocyte competence in humans, it is imperative to investigate the role of epigenetics in the regulation of LIF expression in women with PCOS. Our earlier data on proteomics of FF in women with PCOS demonstrated the downregulation of COC matrix proteins including amphiregulin, TNFAIP6, and bikunin, whose diminished expression is implicated in COC matrix expansion defects [25]. Supporting these findings, we also observed hypermethylation and downregulation of HAPLN1, which is also a COC matrix-associated protein. HAPLN1 facilitates the expansion of COC matrix and imparts stability to the COC complex by binding to other matrix proteins and proteoglycans like hyaluronic acid, versican, aggrecan, and IαI [72]. Exogenous treatment with HAPLN1 increased the CGC viability in vitro and enabled their transformation into granulosa lutein cells, while knockdown of HAPLN1 decreased the cell viability [72]. Therefore, altered HAPLN1 methylation may contribute to changes in COC expansion dynamics in PCOS (Fig. 4).

Luteal insufficiency/premature luteolysis

Apart from being an androgen synthesizing enzyme, AKR1C3 also acts as a prostaglandin F synthase (PGFS), which catalyzes the conversion of the luteotrophic prostaglandins, PGD and PGE2, to a luteolytic form, i.e., PGF2α to facilitate luteolysis [29]. Luteal insufficiency and premature luteolysis are frequent occurrences in PCOS [73], and these have been linked to aberrations in angiogenic mechanisms at the ovarian level [74]. Angiogenesis is largely under the control of prostaglandins [75, 76], and AKR1C family of enzymes are some major regulators of prostaglandins [77]. Disparities in levels of pro and anti-angiogenic factors in the ovary are largely responsible for defects in follicle development, premature degeneration of oocytes, and regression of CL due to an inefficient supply of oxygen and nutrients to the growing follicles [74]. Both, our present data and FF proteome study provide compelling evidence supporting angiogenic dysregulation in the ovarian compartment [25]. Since AKR1C3 has been implicated as a trigger for premature luteolysis [78] and was found to be hypomethylated and high in PCOS follicles, it may serve as an important epigenetic target to investigate this phenomenon in PCOS, Till now, no clear evidence exists on whether PGE2 or PGF2α levels are altered in PCOS. However, our data shows that the receptor for PGE2 gene, i.e., PTGER1, was hypermethylated in NGS and its transcript was low in CGCs of women with PCOS. Since activation of PTGER1 by specific agonists has demonstrated increased sprout formation in capillaries, both in vivo and in vitro [79, 80], it can be an important factor for the restoration of follicular angiogenesis. Also, PTGER1 has been shown to stimulate progesterone biosynthesis in human GCs [81]. Since the follicles of PCOS women lack optimum levels of progesterone required for maintenance of CL, low levels of PTGER1 in PCOS caused by promoter hypermethylation may explain this shortcoming.

Oxidative stress

Lastly, TF, which maintains the oxido-reductive homeostasis in proliferating cells and showed hypermethylation in our NGS data, was found to be hypomethylated by pyrosequencing, though its mRNA was downregulated in PCOS. Therefore, methylation changes at other CpG sites in TF promoter need to be analyzed. TF primarily transports iron released from hepatic, intestinal, and reticuloendothelial stores to the target tissues via its receptor endocytosis while also alleviating local oxidative stress, acting as a growth factor, and promoting follicle and oocyte maturation [82, 83]. Low levels of TF in follicular fluid of PCOS women have been previously reported [84], which may lead to a higher prevalence of unbound Fe2+; trigger oxidative damage to DNA, lipids, and proteins; and contribute to oxidative stress, which is reported in PCOS ovaries [85]. Figure 5 depicts a summary of biological processes which may be affected in PCOS due to the altered methylation of genes in follicles of women with PCOS.

Fig. 5
figure 5

Summary of the genes and their associated processes found to be affected in PCOS due to altered DNA methylation. The figure represents sa preovulatory follicle of women with PCOS, showing a few differentially methylated genes involved in the perpetuation of androgen excess, premature luteolysis, impaired calcium signaling, COC defects, oxidative stress and angiogenic defects in the follicles of women with PCOS as observed in our study. AKR1C3, aldo-keto reductase family 1 member C3; CASR, calcium-sensing receptor; COC, cumulus-oophorus complex; FF, follicular fluid; FSH, follicle-stimulating hormone; FSHR, follicle-stimulating hormone receptor; GC, granulosa cells; IGF1; insulin-like growth factor 1; INS, insulin; MGC, mural granulosa cells; GHRHR, growth hormone-releasing hormone receptor; HAPLN1, hyaluronan and proteoglycan link protein 1; LH, luteinizing hormone; LHCGR, luteinizing hormone chorionic gonadotropin receptor; LIF, leukemia inhibitory factor; MAMLD1, mastermind-like domain containing 1; PTGER1, prostaglandin receptor E1; RETN, resistin; TC, theca cells; TF, transferrin; TNF, tumor necrosis factor alpha


The primary limitation of this study was the low sample size. A number of controls and women with PCOS originally recruited for the study had to be excluded due to the presence of confounding factors such as hyperprolactinemia or thyroid dysfunction, and contamination of FF with blood, or due to recent treatment with metformin or thyroid medications. Further, women having low CGC counts were also excluded as the quantities of nucleic acids were insufficient for methylation and expression analyses. Also, since PCOS is a heterogeneous disorder, investigation of epigenetic changes in a large population based on phenotypic subgrouping of women with PCOS (A, B, C, and D phenotypes) as per the Rotterdam consensus may provide more accurate information on the effect of epigenetic components on PCOS development. Due to our limited sample size, it was not possible to carry out such subgrouping and subsequent analyses.

Our findings highlight that MC-Seq could identify several functionally important loci in CGCs, many of which were either known to be functionally dysregulated in PCOS ovary or have a compelling potential to be established as novel candidates influenced by epigenetic changes. The molecular sequelae of these alterations leading to ovarian dysfunction need to be further addressed via robust functional studies in women with PCOS.

Materials and methods

Study design, participants, sample collection, and estimated parameters

This study was carried out at the ICMR-National Institute for Research in Reproductive Health (NIRRH) as per ethical norms. All participants were recruited from the “Fertility Clinic and IVF Center” (Mumbai) after obtaining written informed consents and underwent IVF using a long, GnRH agonist protocol as reported earlier [8]. We initially recruited 35 women with PCOS and 38 age-BMI-matched, healthy, and regularly menstruating controls as per the Rotterdam consensus criteria [86]. Women showing normal ovarian morphology on ultrasound with no signs of hyperandrogenism or insulin resistance, and undergoing COH strictly owing to indications of male factor infertility in their spouses were recruited as controls. Women recruited as PCOS had at least 2 of the following 3 features, i.e., polycystic ovaries (PCO) on ultrasound, irregularity/absence of menses, and/or signs of hyperandrogenism during clinical screening, however the presence of PCO morphology was used as a mandate for all women recruited as PCOS. These women were carefully screened for their baseline hormonal estimates and biochemical characteristics using serum and follicular fluid samples (Table 1). Further, we had to exclude samples of women who were on metformin or thyroid medications, whose follicular fluid had blood contamination and whose CGC counts were low, since substantial amounts of DNA were required for high-throughput methylome sequencing and bisulfite PCRs. After this, 20 controls and 20 women with PCOS were selected for the study (Table 1). Baseline hormonal estimates (between days 3–7 of menstrual cycle) for LH, FSH, prolactin, TSH, and AMH could be obtained from IVF clinical records, while fasting serum and macroscopically clear FF collected on d-OPU were assayed for E2, P4, TT, and SHBG using commercial ELISA kits (Diagnostics Biochem Canada Inc., Dorchester, Ontario, Canada). Androgen excess indices were calculated using TT and SHBG values [8]. Upon follicle maturation, the levels of E2 and P4, which are routinely measured 1 day prior to and 1 day after rhCG administration (10,000 IU) to monitor ovarian response, were also recorded. COCs suspended in FF aspirates were separated from FF and manually stripped off to dissociate the cumulus granulosa cells (CGCs) from their oocytes. CGCs were washed, resuspended in ovum buffer, and transported from IVF center to the lab at 37 °C for further processing. The numbers of total retrieved oocytes and mature oocytes (in the MII phase of meiosis) were obtained from clinical records, and rates of fertilization of MII oocytes were calculated.

Methyl-capture sequencing (MC-Seq) by NGS approach and analysis of DMRs

From among the 20 women with PCOS and 20 controls selected for the study, individual DNA samples of CGCs from 3 PCOS women and 3 age-BMI-matched controls having a total yield of > 2 μg were subjected to MC-Seq using the NGS approach. Whole genomic DNA was extracted from CGCs of these women using QIAamp DNA mini kit (Qiagen, Hilden, Germany) and processed for library preparation. 1 μg of genomic DNA per sample was sheared in fragments of approximately 150 bp using the Covaris S220 ultrasonicator and further processed for DNA end repair, 3′-adenylation, methylated adapter ligation, and hybridization to the Methyl-Seq library probes using the SureSelect Methyl-Seq Library Prep Kit (Agilent Technologies, CA, USA) as per the manufacturer’s instructions. The hybrids were captured using streptavidin beads (Dynabeads) and subjected to bisulfite conversion using the EZ-DNA Methylation Gold Kit (Zymo Research, CA, USA). The library was PCR amplified, indexed using SureSelectXT Methyl-Seq indexing primers (Agilent Technologies, CA, USA), and pooled for 100 bp paired-end multiplexed sequencing on Illumina HiSeq 2500 platform.

Processing and alignment of reads

Whole-genome targeted methylation sequencing reads of the above six samples were aligned to human reference genome hg19 assembly. The hg19 reference genome was converted to a DNA methylation reference genome, and genome indexing was performed using Bismark genome preparation utility (v.14.3). The adapter trimmed sequencing reads were aligned to the converted methylation reference genome using Bismark tool with two allowed mismatches. The remaining parameters from Bismark were used as default. Further, deduplication of aligned reads was performed using “deduplicate_bismark” utility of Bismark tool.

Methylation analysis and functional annotation

Methylation extractor utility of Bismark tool was used to extract the methylation call for every methylated cytosine (C) in all three contexts: CpG, CHG, and CHH (where H is A or C or T). Differential methylation analysis was performed using “calculateDiffMeth” function of R-based methylKit package at q value cutoff threshold of ≤ 0.01 and methylation difference ≥ 25%. Hierarchical clustering of samples based on the similarity of their methylation profile was performed using the euclidean distance metric, and ward method clustering approaches of methylKit package. Further, enrichment and annotation of both hypo- and hypermethylated sites within up 10 kb of annotated transcription start site (TSS) to transcription end site (UCSC hg19) were performed using in-house perl scripts. GeneCodis3 web-based tool [87] was used for gene ontology (GO) and pathway enrichment analysis of the DMGs.

Selection of genes for validation studies

Genes showing differential methylation at CpG sites upstream up to 1000 bases relative to their TSS in our NGS data were selected for validation. These included a total of 354 hypermethylated and 397 hypomethylated genes (total n = 735 unique genes). The lists of these 735 genes, and the 4354 genes enlisted in Ovarian Kaleidoscope (OKdb), were compared using a Venn analysis, and a total of 132 genes that were common to both datasets could be identified. OKdb is an online search tool based on microarray-based transcriptome profiling and independent study reports on genes and proteins identified in ovarian tissues and cells [88]. Upon identification of genes common to our NGS dataset and OKdb, and a careful review of literature, 7 genes participating in the perpetration of androgen excess, impaired angiogenesis, luteal insufficiency, COC matrix defects, and oocyte defects, namely AKR1C3, CASR, GHRHR, MAMLD1, RETN, TF, and TNF, were selected for validation by both pyrosequencing and qPCR. Additionally, the transcript expression profiles of prostaglandin E receptor (PTGER1), leukemia inhibitory factor (LIF), and hyaluronan and proteoglycan link protein 1 (HAPLN1), were evaluated in CGCs of controls and PCOS women as supporting evidence for pathways associated with the above genes.

Extraction of nucleic acids from CGCs for qPCR and pyrosequencing

Total DNA and RNA were extracted from CGCs of 17 controls and 17 women with PCOS using the NucleoSpin TriPrep kit (Macherey-Nagel, Düren, Germany) for validation of the NGS data. The quality and yield of nucleic acids were assessed by agarose gel electrophoresis and by evaluating their spectrophotometric ratios at 260 and 280 nm. DNA samples were stored at − 20 °C while RNA was stored at − 80 °C until further use.

cDNA synthesis and quantitative real-time PCR

cDNA was synthesized from 500 ng of RNA (n = 34) using a first strand cDNA synthesis kit (Takara Bio USA Inc.). Transcript levels of DMGs selected for validation (Additional file 4) were assayed by TaqMan chemistry using the TaqMan™ Universal Master Mix II with UNG, and FAM-labeled probes (ThermoFisher Scientific, MA, USA). Assay containing VIC-labeled 18s rRNA probe was used as the housekeeping control. qPCR was carried out using cDNA dilutions ranging between neat to 1:100. Fold change in gene expression between controls and PCOS was evaluated using the 2-ΔΔCt method, where the expression was normalized to 18s levels, using a CGC calibrator sample.

Bisulfite primer design and pyrosequencing

Primers for the validation of selected genes were designed on the Pyromark Q96 ID machine using the PyroMark® Assay Design SW 2.0 software (Qiagen), the list of which has been provided in Additional file 4. Primers were procured from Sigma-Aldrich. The reverse pyrosequencing primers were tagged with biotin at the 5′-end and HPLC purified. Approximately 300–500 ng DNA from CGCs of 17 controls and 17 women with PCOS was bisulfite converted using the MethylCode bisulfite conversion kit (Invitrogen-ThermoFisher Scientific, MA, USA). The region of interest was amplified by two rounds of PCR. PCR product (1–3 μL from total volume of 20 μL) from the first round was used as template for the 2nd PCR round, which was scaled up to 45 μL. For both PCRs, initial denaturation was performed at 95 °C for 15 min followed by 40 rounds of amplification at 94 °C for 30s, annealing at the respective optimized temperatures for 10s, and 72 °C for 60s with a final extension at 72 °C for 10 min. Product amplification after both rounds was confirmed by agarose gel electrophoresis. For pyrosequencing, 40 μL of the PCR product was subjected to clean up on the Pyromark Q96ID sequencing workstation as per the manufacturer’s instructions and sequenced using 1.6 μL (16 picomoles) of each of the gene-specific sequencing primers.

Statistical analysis

Mann-Whitney U tests were employed for all univariate assessments of continuous variables including hormonal and biochemical parameters, DNA methylation, and transcript expression levels assessed between CGCs of controls and women with PCOS. Correlation between the CpG methylation status of genes and their respective transcripts was determined using a two-tailed Spearman’s correlation coefficient. FDR cutoffs were applied for genes showing multiple CpG sites as significant after Mann-Whitney testing.





Aldo-keto reductase family 1 member C3


Anti-Mullerian hormone


Anti-Mullerian hormone receptor 2


Androgen receptor


Calcium-sensing receptor


Cumulus granulosa cells of controls


Cumulus granulosa cell


Corpus luteum


Cumulus-oophorus complex


Controlled ovarian hyperstimulation


Differentially methylated gene

E2 :



Endocrine disruptor


Free androgen index


Follicular fluid


Fibroblast growth factor


Follicle-stimulating hormone


Follicle-stimulating hormone receptor




Mural and cumulus granulosa cells


Growth hormone


Growth hormone receptor


Growth hormone-releasing hormone receptor


Gonadotropin-releasing hormone


Gene ontology


Hexose-6-phosphate dehydrogenase/glucose 1-dehydrogenase


Hyaluronan and proteoglycan link protein 1


Human chorionic gonadotropin


Insulin-like growth factor receptor


Insulin-like growth factor binding proteins




Insulin receptor substrate 2


In vitro fertilization


Luteinizing hormone


Luteinizing hormone chorionic gonadotropin receptor


Leukemia inhibitory factor


Long interspersed nucleotide element 1


Mastermind-like domain containing 1


Methyl-capture sequencing


Methylated DNA immunoprecipitation


Mural granulosa cells


Next generation sequencing


Ovarian Kaleidoscope

P4 :



Peripheral blood leukocytes


Cumulus granulosa cells of PCOS


Polycystic ovary syndrome


Platelet-derived growth factor


Prostaglandin D


Prostaglandin E2


Prostaglandin F


Prostaglandin F synthase


Protein kinase C Beta


Prostaglandin receptor E1




Sex hormone binding globulin


Type 2 diabetes mellitus


Transcription end site




Tumor necrosis factor alpha


Tumor necrosis factor-inducible gene 6


Thyroid stimulating hormone (TSH)


Transcription start site


Total testosterone


Vascular endothelial growth factor


Whole genome bisulfite sequencing


Yes-associated protein 1


  1. Fauser BCJM, Tarlatzis BC, Rebar RW, Legro RS, Balen AH, Lobo R, et al. Consensus on women’s health aspects of polycystic ovary syndrome (PCOS): the Amsterdam ESHRE/ASRM-Sponsored 3rd PCOS Consensus Workshop Group. Fertil Steril. 2012;97(1):28–38 e25.

    Article  PubMed  Google Scholar 

  2. Norman RJ, Dewailly D, Legro RS, Hickey TE. Polycystic ovary syndrome. Lancet. 2007;370(9588):685–97.

    Article  CAS  PubMed  Google Scholar 

  3. Diamanti-Kandarakis E. Polycystic ovarian syndrome: pathophysiology, molecular aspects and clinical implications. Expert Rev Mol Med. 2008;10:e3.

    Article  PubMed  Google Scholar 

  4. Baptiste CG, Battista M-C, Trottier A, Baillargeon J-P. Insulin and hyperandrogenism in women with polycystic ovary syndrome. J Steroid Biochem Mol Biol. 2010;122(1–3):42–52.

    Article  CAS  PubMed  Google Scholar 

  5. Wiweko B, Maidarti M, Priangga MD, Shafira N, Fernando D, Sumapraja K, et al. Anti-mullerian hormone as a diagnostic and prognostic tool for PCOS patients. J Assist Reprod Genet. 2014;31(10):1311–6.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Xu N, Azziz R, Goodarzi MO. Epigenetics in polycystic ovary syndrome: a pilot study of global DNA methylation. Fertil Steril. 2010;94(2):781–3 e1.

    Article  CAS  PubMed  Google Scholar 

  7. Sandovici I, Naumova AK, Greenwood CM, editors. Epigenetics and establishment of tissue-specific epigenetic states during development traits. Epigenetics complex trait; 2013. p. 35–62.

    Google Scholar 

  8. Sagvekar P, Mangoli V, Desai S, Patil A, Mukherjee S. LINE1 CpG-DNA hypomethylation in granulosa cells and blood leukocytes is associated with PCOS and related traits. J Clin Endocrinol Metab. 2017;102(4):1396–405.

    Article  PubMed  Google Scholar 

  9. Pruksananonda K, Wasinarom A, Sereepapong W, Sirayapiwat P, Rattanatanyong P, Mutirangura A. Epigenetic modification of long interspersed elements in cumulus cells of mature and immature oocytes from patients with polycystic ovary syndrome. Clin Exp Reprod Med. 2016;43(2):82–9.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Le Jiang L, Xie JK, Cui JQ, Wei D, Yin BL, Zhang YN, et al. Promoter methylation of yes-associated protein (YAP1) gene in polycystic ovary syndrome. Med. 2017;96(2):e5768.

    Article  CAS  Google Scholar 

  11. Sang Q, Zhang S, Zou S, Wang H, Feng R, Li Q, et al. Quantitative analysis of follistatin (FST) promoter methylation in peripheral blood of patients with polycystic ovary syndrome. Reprod BioMed Online. 2013;26(2):157–63.

    Article  CAS  PubMed  Google Scholar 

  12. Yu Y-Y, Sun C-X, Liu Y-K, Li Y, Wang L, Zhang W. Promoter methylation of CYP19A1 gene in Chinese polycystic ovary syndrome patients. Gynecol Obstet Investig. 2013;76(4):209–13.

    Article  CAS  Google Scholar 

  13. Wang P, Zhao H, Li T, Zhang W, Wu K, Li M, et al. Hypomethylation of the LH/choriogonadotropin receptor promoter region is a potential mechanism underlying susceptibility to polycystic ovary syndrome. Endocrinology. 2014;155(4):1445–52.

    Article  PubMed  CAS  Google Scholar 

  14. Qu F, Wang F-F, Yin R, Ding G-L, El-Prince M, Gao Q, et al. A molecular mechanism underlying ovarian dysfunction of polycystic ovary syndrome: hyperandrogenism induces epigenetic alterations in the granulosa cells. J Mol Med. 2012;90(8):911–23.

    Article  CAS  PubMed  Google Scholar 

  15. Zhu JQ, Zhu L, Liang XW, Xing FQ, Schatten H, Sun QY. Demethylation of LHR in dehydroepiandrosterone-induced mouse model of polycystic ovary syndrome. Mol Hum Reprod. 2009;16(4):260–6.

    Article  PubMed  CAS  Google Scholar 

  16. Shen H-R, Qiu L-H, Zhang Z-Q, Qin Y-Y, Cao C, Di W. Genome-wide methylated DNA immunoprecipitation analysis of patients with polycystic ovary syndrome. PLoS One. 2013;8(5):e64801.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Li S, Zhu D, Duan H, Ren A, Glintborg D, Andersen M, et al. Differential DNA methylation patterns of polycystic ovarian syndrome in whole blood of Chinese women. Oncotarget. 2016;8(13):20656–66.

    PubMed Central  Google Scholar 

  18. Yu Y-Y, Sun C-X, Liu Y-K, Li Y, Wang L, Zhang W. Genome-wide screen of ovary-specific DNA methylation in polycystic ovary syndrome. Fertil Steril. 2015;104(1):145–53 e6.

    Article  CAS  PubMed  Google Scholar 

  19. Wang X, Wei J, Jiao J, Jiang S, Yu D, Li D. Genome-wide DNA methylation and gene expression patterns provide insight into polycystic ovary syndrome development. Oncotarget. 2014;5(16):6603–10.

  20. Kokosar M, Benrick A, Perfilyev A, Fornes R, Nilsson E, Maliqueo M, Johan Behre C, Sazonova A, Ohlsson C, Ling C, Stener-Victorina E. Epigenetic and transcriptional alterations in human adipose tissue of polycystic ovary syndrome. Sci Rep. 2016;6(22883).

  21. Pan JX, Tan YJ, Wang FF, Hou NN, Xiang YQ, Zhang JY, et al. Aberrant expression and DNA methylation of lipid metabolism genes in PCOS: a new insight into its pathogenesis. Clin Epigenetics. 2018;10(1):6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Xu J, Bao X, Peng Z, Wang L, Du L, Niu W, et al. Comprehensive analysis of genome-wide DNA methylation across human polycystic ovary syndrome ovary granulosa cell. Oncotarget. 2016;7(19):27899–909.

    PubMed  PubMed Central  Google Scholar 

  23. Ling Teh A, Pan H, Lin X, Lim YI, Pawan C, Patro K, et al. Comparison of methyl-capture sequencing vs. Infinium 450K methylation array for methylome analysis in clinical samples. Epigenetics. 2016;11(1):36–48.

    Article  Google Scholar 

  24. Sagvekar P, Dadachanji R, Patil K, Mukherjee S. Pathomechanisms of polycystic ovary syndrome: multidimensional approaches. Front Biosci (Elite Ed). 2018;10:384–422.

    Google Scholar 

  25. Ambekar AS, Kelkar DS, Pinto SM, Sharma R, Hinduja I, Zaveri K, et al. Proteomics of follicular fluid from women with polycystic ovary syndrome suggests molecular defects in follicular development. J Clin Endocrinol Metab. 2015;100(2):744–53.

    Article  CAS  PubMed  Google Scholar 

  26. Choi DH, Lee WS, Won M, Park M, Park HO, Kim E, et al. The apolipoprotein A-I level is downregulated in the granulosa cells of patients with polycystic ovary syndrome and affects steroidogenesis. J Proteome Res. 2010;9(9):4329–36.

    Article  CAS  PubMed  Google Scholar 

  27. Escobar-Morreale HF, Peral B, Corton M, Camafeita E, Lopez JA, San Millan JL, et al. Proteomic analysis of human omental adipose tissue in the polycystic ovary syndrome using two-dimensional difference gel electrophoresis and mass spectrometry. Hum Reprod. 2007;23(3):651–61.

    PubMed  Google Scholar 

  28. Ma X, Fan L, Meng Y, Hou Z, Mao YD, Wang W, et al. Proteomic analysis of human ovaries from normal and polycystic ovarian syndrome. Mol Hum Reprod. 2007;13(8):527–35.

    Article  CAS  PubMed  Google Scholar 

  29. Nakamura Y, Hornsby PJ, Casson P, Morimoto R, Satoh F, Xing Y, et al. Type 5 17beta-hydroxysteroid dehydrogenase (AKR1C3) contributes to testosterone production in the adrenal reticularis. J Clin Endocrinol Metab. The Endocrine Society. 2009;94(6):2192–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. O’Reilly MW, Kempegowda P, Walsh M, Taylor AE, Manolopoulos KN, Allwood JW, et al. AKR1C3-mediated adipose androgen generation drives lipotoxicity in women with polycystic ovary syndrome. J Clin Endocrinol Metab. The Endocrine Society. 2017;102(9):3327–39.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Dozier BL, Watanabe K, Duffy DM. Two pathways for prostaglandin F2 alpha synthesis by the primate periovulatory follicle. Reproduction. 2008;136(1):53–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Matsoukas M-T, Spyroulias GA. Dynamic properties of the growth hormone releasing hormone receptor (GHRHR) and molecular determinants of GHRH binding. Mol BioSyst. The Royal Society of Chemistry. 2017;13(7):1313–22.

    Article  CAS  PubMed  Google Scholar 

  33. Barabutis N, Schally AV. Growth hormone-releasing hormone: extrapituitary effects in physiology and pathology. Cell Cycle. 2010;9(20):4110–6.

    Article  CAS  PubMed  Google Scholar 

  34. Sharara FI, Giudice LC. Role of growth hormone in ovarian physiology and onset of puberty. J Soc Gynecol Investig. 1997;4(1):2–7.

    CAS  PubMed  Google Scholar 

  35. Bergh C, Carlsson B, Olsson J-H, Selleskog U, Hillensjö T. Regulation of androgen production in cultured human thecal cells by insulin-like growth factor I and insulin. Fertil Steril. 1993;59(2):323–31.

    Article  CAS  PubMed  Google Scholar 

  36. Poretsky L. The insulin-related ovarian regulatory system in health and disease. Endocr Rev Oxford University Press. 1999;20(4):535–82.

    Article  CAS  Google Scholar 

  37. van Dessel HJHMT, Lee PDK, Faessen G, Fauser BCJM, Giudice LC. Elevated serum levels of free insulin-like growth factor I in polycystic ovary syndrome 1. J Clin Endocrinol Metab. 1999;84(9):3030–5.

    Google Scholar 

  38. Bachelot A, Monget P, Imbert-BOLLORÉ P, Coshigano K, Kopchick JJ, Kelly PA, et al. Growth hormone is required for ovarian follicular growth. Endocrinology. 2002;143(10):4104–12.

    Article  CAS  PubMed  Google Scholar 

  39. Franks S, Stark J, Hardy K. Follicle dynamics and anovulation in polycystic ovary syndrome. Hum Reprod Update. 2008;14(4):367–78.

    Article  CAS  PubMed  Google Scholar 

  40. Nakao K, Kishi H, Imai F, Suwa H, Hirakawa T, Minegishi T. TNF-α suppressed FSH-induced LH receptor expression through transcriptional regulation in rat granulosa cells. Endocrinology. 2015;156(9):3192–202.

    Article  CAS  PubMed  Google Scholar 

  41. Kalfa N, Fukami M, Philibert P, Audran F, Pienkowski C, Weill J, et al. Screening of mamld1 mutations in 70 children with 46,xy dsd: identification and functional analysis of two new mutations. PLoS One. 2012;7(3):e32505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kalfa N, Cassorla F, Audran F, Oulad Abdennabi I, Philibert P, Béroud C, et al. Polymorphisms of MAMLD1 gene in hypospadias. J Pediatr Urol. 2011;7(6):585–91.

    Article  CAS  PubMed  Google Scholar 

  43. Fukami M, Wada Y, Okada M, Kato F, Katsumata N, Baba T, et al. Mastermind-like domain-containing 1 (MAMLD1 or CXorf6) transactivates the Hes3 promoter, augments testosterone production, and contains the SF1 target sequence. J Biol Chem. 2008;283(9):5525–32.

    Article  CAS  PubMed  Google Scholar 

  44. Çapoğlu İ, Erdem F, Uyanık A, Turhan H. Serum levels of resistin and hsCRP in women with PCOS. Open Med SP Versita. 2009;4(4):428–32.

    Google Scholar 

  45. Benomar Y, Gertler A, De Lacy P, Crépin D, Hamouda HO, Riffault L, et al. Central resistin overexposure induces insulin resistance through toll-like receptor 4. Diabetes. 2013;62(1):102–44.

    Article  CAS  PubMed  Google Scholar 

  46. Seow KM, Juan CC, Hsu YP, Ho LT, Wang YY, Hwang JL. Serum and follicular resistin levels in women with polycystic ovarian syndrome during IVF-stimulated cycles. Hum Reprod. 2005;20(1):117–21.

    Article  PubMed  Google Scholar 

  47. Panidis D, Koliakos G, Kourtis A, Farmakiotis D, Mouslech T, Rousso D. Serum resistin levels in women with polycystic ovary syndrome. Fertil Steril. 2004;81(2):361–6.

    Article  CAS  PubMed  Google Scholar 

  48. Dinka P, Baldani L, Skrgatic M, Kasum G, Zlopasa S, Kralik O, et al. Altered leptin, adiponectin, resistin and ghrelin secretion may represent an intrinsic polycystic ovary syndrome abnormality. Gynecol Endocrinol. 2019:1–5. In press

  49. Munir I, Yen HW, Baruth T, Tarkowski R, Azziz R, Magoffin DA, et al. Resistin stimulation of 17α-hydroxylase activity in ovarian theca cells in vitro: relevance to polycystic ovary syndrome. J Clin Endocrinol Metab. 2005;90(8):4852–7.

    Article  CAS  PubMed  Google Scholar 

  50. Kushiyama A, Shojima N, Ogihara T, Inukai K, Sakoda H, Fujishiro M, et al. Resistin-like molecule β activates MAPKs, suppresses insulin signaling in hepatocytes, and induces diabetes, hyperlipidemia, and fatty liver in transgenic mice on a high fat diet. J Biol Chem. 2005;280(51):42016–25.

    Article  CAS  PubMed  Google Scholar 

  51. Hofer AM, Brown EM. Extracellular calcium sensing and signalling. Nat Rev Mol Cell Biol. 2003;4(7):530–8.

    Article  CAS  PubMed  Google Scholar 

  52. Elena Dell’Aquila M, De Santis T, Cho YS, Reshkin SJ, Caroli AM, Maritato F, et al. Localization and quantitative expression of the calcium-sensing receptor protein in human oocytes. Fertil Steril. 2006;85:1240–7.

    Article  CAS  Google Scholar 

  53. Ranjzad F, Mahban A, Shemirani AI, Mahmoudi T, Vahedi M, Nikzamir A, et al. Influence of gene variants related to calcium homeostasis on biochemical parameters of women with polycystic ovary syndrome. J Assist Reprod Genet. 2011;28(3):225–32.

    Article  PubMed  Google Scholar 

  54. Hizaki K, Yamamoto H, Taniguchi H, Adachi Y, Nakazawa M, Tanuma T, et al. Epigenetic inactivation of calcium-sensing receptor in colorectal carcinogenesis. Mod Pathol. 2011;24(6):876–84.

    Article  CAS  PubMed  Google Scholar 

  55. Fetahu IS, Höbaus J, Aggarwal A, Hummel DM, Tennakoon S, Mesteri I, et al. Calcium-sensing receptor silencing in colorectal cancer is associated with promoter hypermethylation and loss of acetylation on histone 3. Int J Cancer. 2014;135(9):2014–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Casalà C, Gil-Guiñón E, Luis Ordóñez J, Miguel-Queralt S, Rodríguez E, Galván P, et al. The calcium-sensing receptor is silenced by genetic and epigenetic mechanisms in unfavorable neuroblastomas and its reactivation induces ERK1/2-dependent apoptosis. Carcinogenesis. 2013;34(2):268–76.

    Article  PubMed  CAS  Google Scholar 

  57. Lee JY, Tae JC, Kim CH, Hwang D, Kim KC, Suh CS, et al. Expression of the genes for peroxisome proliferator-activated receptor-γ, cyclooxygenase-2, and proinflammatory cytokines in granulosa cells from women with polycystic ovary syndrome. Clin Exp Reprod Med. 2017;44(3):146.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Schmidt J, Weijdegård B, Mikkelsen AL, Lindenberg S, Nilsson L, Brännström M. Differential expression of inflammation-related genes in the ovarian stroma and granulosa cells of PCOS women. Mol Hum Reprod. 2014;20(1):49–58.

    Article  CAS  PubMed  Google Scholar 

  59. Amato G, Conte M, Mazziotti G, Lalli E, Vitolo G, Tucker AT, et al. Serum and follicular fluid cytokines in polycystic ovary syndrome during stimulated cycles. Obstet Gynecol. 2003;101(6):1177–82.

    CAS  PubMed  Google Scholar 

  60. Gao L, Gu Y, Yin X. High serum tumor necrosis factor-alpha levels in women with polycystic ovary syndrome: a meta-analysis. Atkin SL, editor PLoS One. 2016;11(10):e0164021.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Yamamoto Y, Kuwahara A, Taniguchi Y, Yamasaki M, Tanaka Y, Mukai Y, et al. Tumor necrosis factor alpha inhibits ovulation and induces granulosa cell death in rat ovaries. Reprod Med Biol. 2015;14(3):107–15.

    Article  CAS  PubMed  Google Scholar 

  62. Sakumoto R, Shibaya M, Okuda K. Tumor necrosis factor-alpha (TNF alpha) inhibits progesterone and estradiol-17beta production from cultured granulosa cells: presence of TNFalpha receptors in bovine granulosa and theca cells. J Reprod Dev. 2003;49(6):441–9.

    Article  CAS  PubMed  Google Scholar 

  63. Morrison LJ, Marcinkiewicz JL. Tumor necrosis factor α enhances oocyte/follicle apoptosis in the neonatal rat ovary1. Biol Reprod. 2002;66(2):450–7.

    Article  CAS  PubMed  Google Scholar 

  64. Paulino LRFM, Cunha EV, Barbalho Silva AW, Souza GB, Lopes EPF, Donato MAM, et al. Effects of tumour necrosis factor-alpha and interleukin-1 beta on in vitro development of bovine secondary follicles. Reprod Domest Anim. 2018;53(4):997–1005.

    Article  CAS  PubMed  Google Scholar 

  65. Crespo D, Goetz FW, Planas JV. Luteinizing hormone induces ovulation via tumor necrosis factor α-dependent increases in prostaglandin F2α in a nonmammalian vertebrate. Sci Rep. 2015;5:14210.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Son D-S, Arai KY, Roby KF, Terranova PF. Tumor necrosis factor α (TNF) increases granulosa cell proliferation: dependence on c-Jun and TNF receptor type 1. Endocrinology. 2004;145(3):1218–26.

    Article  CAS  PubMed  Google Scholar 

  67. Yuan Y, Ida JM, Paczkowski M, Krisher RL. Identification of developmental competence-related genes in mature porcine oocytes. Mol Reprod Dev. 2011;78(8):565–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Li Z, Zhu Y, Li H, Jiang W, Liu H, Yan J, et al. Leukaemia inhibitory factor in serum and follicular fluid of women with polycystic ovary syndrome and its correlation with IVF outcome. Reprod BioMed Online. 2018;36(4):483–9.

    Article  CAS  PubMed  Google Scholar 

  69. Lédée-Bataille N, Laprée-Delage G, Taupin JL, Dubanchet S, Taieb J, Moreau JF, et al. Follicular fluid concentration of leukaemia inhibitory factor is decreased among women with polycystic ovarian syndrome during assisted reproduction cycles. Hum Reprod. 2001;16(10):2073–8.

    Article  PubMed  Google Scholar 

  70. Cheung LP, Leung HY, Bongso A. Effect of supplementation of leukemia inhibitory factor and epidermal growth factor on murine embryonic development in vitro, implantation, and outcome of offspring. Fertil Steril. 2003;80(SUPPL. 2):727–35.

    Article  PubMed  Google Scholar 

  71. De Matos DG, Miller K, Scott R, Tran CA, Kagan D, Nataraja SG, et al. Leukemia inhibitory factor induces cumulus expansion in immature human and mouse oocytes and improves mouse two-cell rate and delivery rates when it is present during mouse in vitro oocyte maturation. Fertil Steril. 2008;90(6):2367–75.

    Article  PubMed  Google Scholar 

  72. Liu J, Park E-S, Curry TE, Jo M, Jo M. Periovulatory expression of hyaluronan and proteoglycan link protein 1 (Hapln1) in the rat ovary: hormonal regulation and potential function. Mol Endocrinol The Endocrine Society. 2010;24(6):1203–17.

    Article  CAS  Google Scholar 

  73. Shah D, Patel S. Polycystic ovarian syndrome as a cause of recurrent pregnancy loss. J Obs Gynecol India. 2007;57(5):391–7.

    Google Scholar 

  74. Di Pietro M, Pascuali N, Parborell F, Abramovich D. Ovarian angiogenesis in polycystic ovary syndrome. Reproduction. 2018;155(5):R199–209.

    Article  CAS  PubMed  Google Scholar 

  75. Karnezis T, Shayan R, Fox S, Achen MG, Stacker SA. The connection between lymphangiogenic signalling and prostaglandin biology: a missing link in the metastatic pathway. Oncotarget. 2012;3(8):893–906.

    Article  PubMed  Google Scholar 

  76. Mantel A, Carpenter-Mendini AB, Vanbuskirk JB, De Benedetto A, Beck LA, Pentland AP. Aldo-keto reductase 1C3 is expressed in differentiated human epidermis, affects keratinocyte differentiation, and is upregulated in atopic dermatitis. J Invest Dermatol. 2012;132(4):1103–10.

    Article  CAS  PubMed  Google Scholar 

  77. Rižner TL, Penning TM. Role of aldo-keto reductase family 1 (AKR1) enzymes in human steroid metabolism. Steroids. 2014;79:49–63.

    Article  PubMed  CAS  Google Scholar 

  78. Nio-Kobayashi J, Kudo M, Sakuragi N, Iwanaga T, Duncan WC. Loss of luteotropic prostaglandin E plays an important role in the regulation of luteolysis in women. MHR Basic Sci Reprod Med. 2017;23(5):271–81.

    Article  CAS  Google Scholar 

  79. Trau HA, Davis JS, Duffy DM. Angiogenesis in the primate ovulatory follicle is stimulated by luteinizing hormone via prostaglandin E2. Biol Reprod. 2015;92(1):15.

    Article  PubMed  CAS  Google Scholar 

  80. Trau HA, Brännström M, Curry TE, Duffy DM, Duffy DM. Prostaglandin E2 and vascular endothelial growth factor A mediate angiogenesis of human ovarian follicular endothelial cells. Hum Reprod. 2016;31(2):436–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  81. Chandras C, Harris TE, López Bernal A, Abayasekara DRE, Michael AE. PTGER1 and PTGER2 receptors mediate regulation of progesterone synthesis and type 1 11β-hydroxysteroid dehydrogenase activity by prostaglandin E2in human granulosa-lutein cells. J Endocrinol. 2007;194(3):595–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Aleshire SL, Osteen KG, Maxson WS, Entman SS, Bradley CA, Parl FF. Localization of transferrin and its receptor in ovarian follicular cells: morphologic studies in relation to follicular development. Vol. 51. Fertil Steril. 1989;51(3):444–9.

    Article  CAS  PubMed  Google Scholar 

  83. Gomme PT, McCann KB. Transferrin: structure, function and potential therapeutic actions. Elsevier Current Trends. 2005;10(4):267–73.

    CAS  Google Scholar 

  84. Wu Y-T, Wu Y, Zhang J-Y, Hou N-N, Liu A-X, Pan J-X, et al. Preliminary proteomic analysis on the alterations in follicular fluid proteins from women undergoing natural cycles or controlled ovarian hyperstimulation. J Assist Reprod Genet. 2015;32(3):417–27.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Chattopadhayay R, Ganesh A, Samanta J, Jana SK, Chakravarty BN, Chaudhury K. Effect of follicular fluid oxidative stress on meiotic spindle formation in infertile women with polycystic ovarian syndrome. Gynecol Obstet Investig. 2010;69(3):197–202.

    Article  CAS  Google Scholar 

  86. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome. Fertil Steril. 2004;81(1):19–25.

  87. Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40(W1):W478–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Leo CP, Vitt UA, Hsueh AJW. The ovarian kaleidoscope database: an online resource for the ovarian research community. Endocrinology. 2000;141(9):3052–4.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors thank all the consenting study participants and collaborators who became a part of this study and Ms. Kavita Suri (Fertility Clinic and IVF Center, Mumbai) who collected the clinical samples of our study participants and provided us with their IVF data records.

Disclosure summary

The authors have nothing to disclose.


This work is partially supported by Department of Science and Technology (DST) (EMR/2015/001014), NIRRH (NIRRH/RA/623/04–2018) and Indian Council of Medical Research (ICMR), New Delhi, India. Financial assistance provided by DST and Department of Atomic Energy-Board of Research in Nuclear Sciences (DAE-BRNS), Government of India; to PS for pursuing her doctoral studies is acknowledged.

Availability of data and materials

Datasets analyzed during the current study will be made available from the corresponding author upon reasonable request. Supplementary data related to the methylome analysis of cumulus granulosa cells is available on Figshare, an online digital data repository

Author information

Authors and Affiliations



The study design was conceived by SM and PS. PS collected the patient samples, processed the samples, carried out all the experiments, performed data interpretation, and wrote the paper. PK carried out the bioinformatic analysis and interpretation of the NGS data. VM and SD are clinical collaborators on this manuscript who supervised the recruitment of study participants and collection of their samples at the Fertility Clinic and IVF Center, Mumbai. SM critically revised the manuscript for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Srabani Mukherjee.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all study participants who gave blood, follicular fluid and cumulus granulosa cell samples, in accordance with current edition of the NIRRH Ethics Committee for Clinical Studies (Ref: D/ICEC/Sci-54/56/2018 dated May 29, 2018). This committee is recognized by Strategic Initiative for Developing Capacity in Ethical Review (SIDCER) and Forum for Ethics Review Committees in Asia and the Western Pacific Region (FERCAP).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Hypomethylated CpG sites were representative of 2063 genes. (XLSX 408 kb)

Additional file 2:

Hypermethylated sites were linked to a total of 1777 genes (XLSX 354 kb)

Additional File 3:

Figure includes horizontal bar charts showing components such as A) molecular functions, B) biological processes, and C) cellular components that were most highly enriched in datasets obtained for all differentially methylated, hypomethylated, and hypermethylated genes identified in the NGS analysis. X-axis represents the number of genes present within each annotated category. (TIF 901 kb)

Additional File 4:

The table enlists pyrosequencing primer sets designed to evaluate upstream/promoter CpG methylation in differentially methylated genes (DMGs) selected for validation in controls (C-CC, n = 17) and women with PCOS (P-CC, n = 17) and commercial TaqMan Assay IDs used for validating the transcript expression levels of selected DMGs in controls and women with PCOS. (DOCX 13 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sagvekar, P., Kumar, P., Mangoli, V. et al. DNA methylome profiling of granulosa cells reveals altered methylation in genes regulating vital ovarian functions in polycystic ovary syndrome. Clin Epigenet 11, 61 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: