Skip to main content

Polycystic ovarian syndrome is accompanied by repression of gene signatures associated with biosynthesis and metabolism of steroids, cholesterol and lipids

Abstract

Background

Polycystic ovarian syndrome (PCOS) is a spectrum of heterogeneous disorders of reproduction and metabolism in women with potential systemic sequel such as diabetes and obesity. Although, PCOS is believed to be caused by genetic abnormalities, the genetic background that can be associated with PCOS phenotypes remains unclear due to the complexity of the trait. In this study, we used a rat model which exhibits reproductive and metabolic abnormalities similar to the human PCOS to unravel the molecular mechanisms underlining this complex syndrome.

Methods

Female Sprague–Dawley rats were randomly assigned to DHT and control (CTL) groups. Rats in the DHT group were implanted with a silicone capsule continuous-releasing 83 μg 5α-dihydrotestosterone (DHT) per day for 12 weeks to mimic the hyperandrogenic state in women with PCOS. The animals were euthanized at 15 weeks of age and the pairs of ovaries were excised and the ovarian cortex tissues were used for gene expression analysis. Total RNA was from the ovarian cortex was amplified, labeled and hybridized to the Affymetrix GeneChip® Rat Genome 230 2.0 Array. A linear model system for microarray data analysis was used to identify genes affected in DHT treated rat ovaries and the molecular pathway of those genes were analyzed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) analysis tool.

Results

A total of 573 gene transcripts, including CPA1, CDH1, INSL3, AMH, ALDH1B1, INHBA, CYP17A1, RBP4, GAS6, GAS7 and GATA4, were activated while 430 others including HSD17B7, HSD3B6, STAR, HMGCS1, HMGCR, CYP51, CYP11A1 and CYP19A1 were repressed in DHT-treated ovaries. Functional annotation of the dysregulated genes revealed that biosynthesis and metabolism of steroids, cholesterol and lipids to be the most top functions enriched by the repressed genes. However, cell differentiation/proliferation, transcriptional regulation, neurogenesis, cell adhesion and blood vessel development processes were enriched by activated genes.

Conclusion

The dysregulation of genes associated with biosynthesis and metabolism of steroids, cholesterol and lipids, cell differentiation/proliferation in DHT- treated ovaries could be a molecular clue for abnormal steroidogenesis, estrous cycle irregularity, abnormal folliculogenesis, anovulation and lipid metabolism in PCOS patients.

Introduction

The ovary is a key organ in the female reproductive system and its malfunction due to endocrine abnormalities could result in female infertility. Polycystic ovarian syndrome (PCOS) is one of the common hormonal and metabolic disorders in women of reproductive age [1]. However, due to its heterogeneity and complexity, universally accepted clinical definition of PCOS remains ambiguous [2,3]. Indeed, the presence of polycystic ovarian morphology is one of the common phenomena that can occur in the majority of PCOS patients. About 95% of women with PCOS at their early follicular phase could have polycystic ovaries and reduced level of follicle stimulating hormone [4] which may lead to antral follicle growth arrest and increased luteinizing hormone level [5]. In addition, PCOS is also associated with hyperandrogenism, menstrual dysfunction, oligo-ovulation and insulin resistance [6]. In this context, PCOS is considered as a complex androgen excess accompanied by different degrees of gonadotropic and metabolic dysregulation controlled by multiple gene interaction and environmental factors [7]. However, to what extent this trait is transmitted to the next generation and the intrinsic molecular factors underlining the occurrence of PCOS is unclear.

Although, the genetic basis of abnormal follicular development, anovulation, metabolic disorder and other heterogenous clinical abnormalities of PCOS patients seems to require detailed investigation, it is suggested that daughters from women exhibiting a characteristics of PCOS could have a higher chance of acquiring hyperandrogenism and other PCOS phenotypes [8]. Moreover, single nucleotide polymorphism in thyroid adenoma associated (THADA), DENN/MADD domain containing 1A (DENND1A), interleukin 6 (IL6) and adiponectin genes has been suggested to be the genetic causes of PCOS [9-11]. In addition, in vitro studies also showed altered expression of CYP11A and CYP17 genes in theca cell derived from PCOS woman [12]. Furthermore, changes in the granulosa and theca cell gene expression have been reported in women with PCOS [13-15]. Although these association studies were performed using the samples of PCOS patients, the majority of gene expression studies were based on the cell culture models which may not necessarily represent and describe the biological and molecular networks governing its complex phenotype. Indeed, this can in part be due to the availability and accessibility of the human sample or small sample size of the study populations [7] and the complexity of the trait between individuals. However, to addresse the clinical heterogeneity of PCOS, animal models have been described to be the best option to investigate the pathophysiologic mechanisms associated with the etiology of PCOS [16-21]. Therefore, to uncover the broad basis of molecular mechanisms associated with physiological and anatomical changes induced by PCOS, we generated a rat PCOS model that exhibit both polycystic ovaries (PCO) and metabolic abnormalities by implanting silastic capsules containing 5α-dihydrotestosterone (DHT) into their ovary in similar way as previously described by others [22]. Using this rat PCOS model that exhibits both polycystic ovaries (PCO) and metabolic abnormalities, we have previously demonstrated altered expression of 89 miRNAs following chronic androgen treatment [23]. However, the genes that are activated or repressed as well as their molecular functions, gene networks and molecular pathways associated with PCOS phenotypes, remained unclear. Therefore, this study was conducted to gain insight into the genes that are associated with follicular arrest, abnormal steroid and metabolite biosynthesis and metabolism, insulin resistance and ovarian dysfunction.

Materials and methods

Details of the materials and methods used in the present study have been described in our previous publication [23]. Briefly, female Sprague–Dawley rats were randomly assigned to DHT and control (CTL) groups. Rats in the DHT group were implanted with a silicone capsule continuous-releasing 83 μg 5α-dihydrotestosterone (DHT) per day for 12 weeks to mimic the hyperandrogenic state in women with PCOS, whose plasma DHT levels are approximately 1.7-fold higher than those of healthy control and those in CTL group received empty capsule [22]. The animals were euthanized at 15 weeks of age, ovaries were excised and extraneous tissues carefully removed. Corpus luteum (CL) was present in most of control rat ovaries while none or very few CL were observed in DHT-treated rat ovaries. Ovarian cortex tissues were snap-frozen in liquid nitrogen and stored at −80°C for further analysis. The PCOS phenotypic characteristics of DHT-treated rats have been described [23].

Gene expression analysis using GeneChip@rat genome array

Total RNAs were isolated from 3 independent DHT and CTL rat ovaries using miRNeasy mini kit (Qiagen, Hilden, Germany). Genomic DNA contamination was removed from the RNA samples using TURBO DNA-free™ kit (Ambion, Foster City, CA). The concentration of the RNA was analyzed using the Nanodrop 8000 Spectrophotometer (Thermo Fisher Scientific Inc, DE, USA). The RNA integrity and quality was evaluated using Agilent 2100 Bioanalyzer with RNA 6000 Nano LabChip® Kit (Agilent Technologies Inc, CA, USA).

RNA amplification

250 ng of total RNAs from DHT-treated or CTL rat groups in four replicates was amplified and labeled as per the GeneChip®3’ IVT Express Kit (Affymetrix, CA, USA). Eukaryotic poly-A RNA control kit (Affymetrix, CA, USA) was used as a SPIKE-IN control to monitor the entire target labeling process. Following amplification, the biotin labeled amplified RNA (aRNA) was purified and fragmented. The distribution of aRNA fragments were evaluated using Agilent 2100 bioanalyzer with RNA 6000 Nano LabChip® Kit (Agilent Technologies Inc, CA, USA).

Sample hybridization, array washing, staining and scanning

Prior to hybridization, the fragmented and biotin labelled cRNA from each rat ovary was mixed with control oligonucleotide B2 (3 nM), 20× eukaryotic hybridization controls (bioB, bioC, bioD, cre) (Affymetrix, CA, USA), 2× hybridization mix and DMSO. The hybridization cocktail were then incubated at 99°C (5 min) and subsequently at 45°C (5 min). Each sample was then transferred to independent GeneChip® Rat Genome 230 2.0 Array chip. Three biological replicates and one technical replicate (pool of three biological replicates) were hybridized for each rat group for 16 h. The array slides were washed and stained using the Fluidics Station 450/250 (Affymetrix, CA, USA), according to the GeneChip® expression user manual (P/N 702232 Rev. 3). Arrays were scanned with the GeneChip™ 3000 laser confocal slide scanner (Affymetrix, CA, USA) integrated with GeneChip® Operating System (GCOS).

Array data analysis and visualization

The array data was normalized by integrating the bioconductor packages (http://bioconductor.org) in R environment (www.r-project.org), using GC robust multi-array average analysis [24]. Briefly, the cell intensity (CELL) files were imported into R software after loading bioconductor packages (http://bioconductor.org) that suit to the Rat GeneChip affymetrix array. The normalized data and the CELL files are stored in the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/, series entry number GSM1437398). The linear models for microarray data analysis system (LIMMA) [25] and the Benjamini–Hochberg procedure of false discovery rate adjustment [26] were employed to discriminate the gene expression profile between the samples. The differentially expressed genes were tested for their gene ontology (GO) terms for over- or under-representation, using a classical hypergeometric test [27]. The molecular pathway enriched by differentially expressed genes were obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Panther pathway data bases, using The Database for Annotation, Visualization and Integrated Discovery (DAVID) analysis tool [28]. The heatmaps and clustering of differentially expressed genes were constructed using Bioconductor (http://www.bioconductor.org) in R software environment http://www.r-project.org/ and/or PermutMatrix [29].

Validation of differentially expressed genes

Some of the differentially expressed genes were randomly selected for validation, using SYBR Green based quantitative real time polymerase chain reaction (qPCR) and with sequence specific primers designed online (http://frodo.wi.mit.edu/primer3/; Table 1). The specificity and identity of each primer pair was confirmed by sequencing with the GenomeLab™ GeXP Genetic Analysis System (Beckman Coulter). The mRNA level was subsequently quantified using the cDNAs obtained from reverse transcription of total RNA samples used for microarray study. The qPCR was then performed in 20 μl reaction volume containing iTaq SYBR Green Supermix with ROX (Bio-Rad laboratories, Munich, Germany), the cDNA samples of DHT or CTL samples and the specific forward and reverse primer in the StepOnePlus™ Real-Time PCR Systems (Applied Biosystems, Foster City, CA). The presence of specific amplification was monitored by evaluating the dissociation curve. The abundance of each transcript was determined using the comparative threshold cycle (ΔCT) method. Data were analyzed after normalizing the Ct value of the target genes against the housekeeping gene β actin (β-actin). The Student’s t-test or least significant difference test procedures was employed to detect differences in mRNA levels between samples. The level of activation or repression of a gene in DHT relative to CTL was determined using the formula 2-ΔΔCT. Differences with p < 0.05 were considered significant.

Table 1 Gene specific primers used for validation of differentially expressed genes

Results

Ovarian polycystic syndrome is associated with dysregulation of gene expression

This study was a part of our previously published results and thus a detailed description of the PCOS phenotypic data has been provided [23]. The DHT- treated rats had higher body weight, oestrous cycle irregularity and reduced insulin sensitivity. Ovaries from these rats had lower weight and exhibited absence of corpus luteum, higher percentage of follicle cysts, relatively thin membrane granulosa and theca hyperplasia. Thus, to unravel the genes that are abnormally activated or repressed due to PCOS, the ovarian transcriptome profile of DHT-treated and CTL rats were compared using the GeneChip @rat Genome array. The genes of which the pattern of expression was different between the DHT and CTL groups were identified by a linear model described in the LIMMA pipeline and bioconductor packages. 1774 of the 31046 gene transcripts available in GeneChip® Rat Genome 230 2.0 A array showed ≥ 2 fold change dysregulation in DHT-treated rat ovaries compared to CTL group (data not shown). However, using a stricter selection criteria (fold change ≥ 2, p ≤ 0.01 and false discovery rate (FDR) ≤ 10%), a total of 1003 transcripts were dysregulated in DHT-treated compared to CTL group (Figure 1A and B). Hierarchical clustering of the dysregulated genes exhibited a similar expression pattern within biological replicates and distinct differences between the groups (Figure 1C).

Figure 1
figure 1

Ovarian transcriptome profile differences between DHT-treated and CTL rat groups. (A) The number of annotated (with known gene symbol) and non-annotated (unknown gene symbol) dysregulated genes in DHT rat groups. (B) The fold change (FC) distribution of dysregulated genes. (C) The heatmap showing the expression pattern of differentially expressed genes within and between biological replicates. DHT represents biological replicates in DHT-treated rat groups while CTL represents biological replicates in rats received an empty silastic capsule. Heatmap was generated using normalized log2 transformed values and the normalized log2 transformed expression values are described by pseudo color scale with red indicating the activated transcript level while the green color shows repressed expression pattern of a specific gene.

Gene transcripts repressed in DHT-treated rats

Following identification of differentially expressed genes between DHT-treated and CTL groups, the genes of which their transcript level was significant reduced by ≥ 2 folds changes in DHT compared to CTL were considered as repressed genes. However, those genes which showed by ≥ 2 folds reduction but displayed a false discovery rate higher that 0.1 were excluded in the analysis. Based on these stringent criteria, including 123 express sequence tags (ESTs), a total of 430 gene transcripts were found to be repressed in the DHT-treated ovaries. Of all the repressed genes, 320 transcripts displayed fold changes between 2 and 5 whereas the expression level of 110 genes was dysregulated between 5 and 134 folds in the DHT group. The list and expression pattern of genes repressed by 11.9 or more folds are described in Figure 2. Among the top repressed genes, the expression pattern of relaxin 1 (RLN1), Phospholipase A2 (PLA2G1B), x-prolyl aminopeptidase (XPNPEP2), Alpha-2-macroglobulin (A2M) and fatty acid binding protein 6 (FABP6) was repressed by 133.9, 70.6, 61.6, 58 and 45.8 folds, respectively. In addition, genes involved in steroid and lipid biosynthesis and metabolism, including hydroxysteroid (17-beta) dehydrogenase 7 (HSD17B7), hydroxy-delta-5-steroid dehydrogenase 3 beta- and steroid delta-isomerase 6 (HSD3B6), steroidogenic acute regulatory protein (STAR), 3-hydroxy-3-methylglutaryl-coenzyme A synthase 1 (HMGCS1), 3-hydroxy-3-methylglutaryl-coenzyme reductase (HMGCR) and cytochrome P450 (CYP51, CYP11A1, CYP19A1) were repressed in DHT-treated ovaries.

Figure 2
figure 2

Significantly repressed genes by ≥ 11.9 folds in DHT-treated rat ovaries. The hierarchical clustering, the expression pattern, the average fold change (FC) and the false discovery rate of significantly highly repressed genes in DHT treated rat ovaries. DHT represents biological replicates in DHT-treated rat groups while CTL represents biological replicates in rats received an empty silastic capsule. The red color indicates the activated transcript level in CTL groups while the green color shows repressed expression pattern of a specific gene in DHT treated rat ovaries.

Functional classification of the genes repressed in ovaries of DHT-treated rats

To understand the biological processes and molecular function over- or under-represented in ovaries associated with repressed genes in the DHT-treated rats, the functional annotation of the repressed genes were interrogated using the DAVID bioinfomatic tool and the result revealed a total 31 crucial biological processes to be repressed in DHT-treated ovaries (Additional file 1: Table S1). Most importantly, closely interlinked biosynthetic and metabolic functions, namely the biosynthesis and metabolism of sterols/steroids, cholesterol, isoprenoids and lipids, were the top significant biological processes enriched by repressed genes (Figure 3). In addition, the repressed genes were also found to be associated with oxidation/reduction, reactive oxygen species, metabolism and immune processes. Moreover, categorization of the repressed genes into their corresponding activities (molecular functions) revealed that 28% of them were known to be involved in catalytic activity, whereas 15% of the repressed genes are associated with binding of cofactors, coenzymes, irons, carbohydrates, SH3 domains, FADs, heparins, glycosaminoglycans, NADPs/NADPHs, cholesterols, glucoses, immunoglobulins and sterols (Figure 4, Additional file 2: Table S2). Moreover, 34 repressed genes, including IDH1, ME1, NSDHL, FASN, GPD1, CYBB, FDXR, ALDH3B1 and ACAD9, are believed to regulate oxidoreductase activity (Figure 4, Additional file 2: Table S2).

Figure 3
figure 3

The top biological processes enriched by genes which were repressed in DHT-treated rat ovaries. The heatmaps on the right (A, B, C, D, E, F) describe clusters of genes involved in a particular function described by A, B, C, D, E, F on the left. The Heatmaps were generated using normalized log2 transformed values and the normalized log2 transformed expression values are described by pseudo color scale with red in CTL groups indicating the activated transcript level while the green color in DHT shows the repressed expression pattern of a specific gene. DHT represents biological replicates in DHT-treated ovaries while CTL represents biological replicates in rats which receive an empty silastic capsule.

Figure 4
figure 4

Molecular functions overrepresented by repressed genes in DHT-treated ovaries.

Gene transcripts activated in ovaries of DHT-treated rats

To identify the transcripts activated in rat ovaries by DHT treatment, genes of which the expression was increased by ≥ 2-fold (p ≤ 0.01, FDR < 0.1) were investigated, using bioinformatic tools and literature mining. A total of 573 gene transcripts (57.1% of the total dysregulated genes) were activated in ovaries of DHT-treated rats. The list of top activated genes (with ≥ 4 fold changes) and their expression pattern within and between biological replicates is provided in Figure 5. The most activated genes (between 8- and 35-fold change) include carboxypeptidase A1 (CPA1), chymotrypsinogen B1 (CTRB1), myosin heavy chain (MYH6), carboxypeptidase A2 (CPA2), dual specificity phosphatase 27 (DUSP27) and cadherin 1 (CDH1). In addition, increased expression of genes involved in follicular growth and function was also evident in DHT-treated ovaries, including, anti-Mullerial hormone (AMH), aldehyde dehydrogenase 1B1 (ALDH1B1), inhibins (INHBA, INHBB, INHA), CYP17A1, retinol binding protein (RBP4), growth arrest specific (GAS6, GAS7) and GATA binding protein 4 (GATA4).

Figure 5
figure 5

The hierarchical clustering and expression pattern of highly activated genes in DHT-treated ovaries relative to CTL group. DHT and CTL on the dendogram represent biological replicates in DHT-treated and CTL rat ovaries, respectively. The Heatmap was generated using normalized log2 transformed values and the normalized log2 transformed expression values are described by pseudo color scale with red in DHT-treated group indicating the activated transcript level while the green color in CTL group describes the repressed expression pattern of a specific gene.

Molecular functions activated in ovaries of DHT-treated rats

Gene set enrichment analysis showed 18 candidate biological processes, including cell differentiation/proliferation, transcriptional regulation/gene expression, neurogenesis, cell adhesion, RNA metabolism, macromolecule biosynthesis and blood vessel development processes to be affected in DHT-treated ovaries due to gene activation (Figure 6). Moreover, the activated genes known to be involved in selective and non-covalent binding of zinc ions, receptors, growth factors, protein phosphatases, DNA secondary structure and peptide antigens. In addition, some of the activated genes are known to be involved in growth factor activity, initiating for cell growth or proliferation and transcription corepressor activity (Table 2).

Figure 6
figure 6

The most significantly molecular functions enriched by genes activated in DHT-treated ovaries. The direction of the arrow shows the significant level of the gene ontology terms. The heatmaps describe the expression pattern of clusters of genes involving in selected functions. Heatmaps were generated using normalized log2 transformed values and the normalized log2 transformed expression values are described by pseudo color scale with red in DHT group indicating the activated transcript level while the green color in CTL shows repressed expression pattern of a specific gene. DHT and CTL on the bottom of the heatmaps represent biological replicates in DHT-treated rat ovaries and rats which receive an empty silastic capsule, respectively

Table 2 List of molecular functions containing genes with increased level of mRNA in DHT treated compared to CTL rat groups

Molecular pathways activated or repressed in DHT-treated rats

In addition to their biological or molecular functions, an important and significantly dysregulated molecular interactions and relations associated with the activated or repressed genes were identified using KEGG and Panther gene enrichment analysis. The studies suggest 19 molecular pathways being affected in DHT-treated group (Table 3). Among these, cholesterol and trepenoid biosynthesis pathways, citrate cycle, androgen and estrogen metabolism were enriched by repressed genes while several metabolic pathways including glycolysis/gluconeogenesis, fatty acid metabolism, pyruvate metabolism, butanoate metabolism, lysine, leucine, valine and isoleucine degradation and glutathione metabolic pathways were enriched by both activated and repressed genes (Table 3).

Table 3 List of molecular pathways containing genes with activated (↑) and repressed (↓) level of mRNA in DHT treated compared to CTL rat ovary groups

Validation of microarray data using real time quantitative PCR (qPCR)

To validate the expression data generated by microarray analysis, the expression level of 5 activated (CTRB1, MYH6, PLAU, ST3GAL6 and INHBA) and 4 repressed (STAR, HSD17B7, OSGIN1 and RLN1) genes were assessed by qPCR. The validation result shows that all the randomly selected genes displayed a similar trend to the microarray data, thus confirming the validity of the microarray results (Table 4).

Table 4 The array and qPCR results for selected differentially expressed genes

Discussion

Using a rat PCOS model, we have previously reported altered ovarian expression pattern of 83 miRNAs following DHT treatment [23]. In this study, we investigated the gene expression profile of the same ovarian samples and identified 573 activated and 430 repressed gene transcripts in DHT-treated rats, suggesting the presence of transcriptome profile dysregulation due to hyperandrogenism. In addition, the cellular localization of the products of the activated or repressed gene showed 180 of the dysregulated ones were present in the nucleus while the majority were localized in the cytoplasm (Additional file 3: Figure S1), suggesting a possible dysregulation of genes function associated with specific ovarian subcellular localization.

The number of dysregulated ovarian genes in the DHT-treated rats appeared consistent with an earlier report by human ovarian cDNA microarray [30] in which the number of up-regulated genes (n = 88) was relatively higher than the down-regulated ones (n = 31) in ovaries of PCOS women compared to the non-PCOS subjects. However in another study where the GeneChips HG_U133A and HG_U133B arrays from Affymetrix used, majority of the dysregulated genes in PCOS were found to be down-regulated [31]. Although the reasons for these apparent discrepancies are not clear, the possibility that this could be due to the differences in tissue sampling, the microarray platform used and the statistical analysis cannot be excluded. To evaluate whether altered ovarian genes in DHT-treated rat resemble those of PCOS women, we merged our data with the supplemental data previously reported [31] (http://press.endocrine.org/doi/suppl/10.1210/me.2004-0074/suppl_file/suppltable1part1me_04_0074.xls) and noted that 160 dysregulated genes were observed in ovaries of both DHT-rats and PCOS women (Additional file 4: Table S3).

In the current study, we have described the key biological processes, molecular functions and pathways affected by dysregulated genes in DHT-treated rats. While the text mining approach would shed light on the relevance of the dysregulated genes with respect to the PCOS phenotype, the present global approach involving pathway analysis and molecular functions provided a more comprehensive understanding of the potential genetic mechanism underling PCOS phenotypes. Our findings also demonstrated ovarian cell type-specific changes in expression of genes involved in granulosa cell proliferation and progesterone biosynthesis [AMH, RBP4 and cytochrome P450s (CYP51, CYP19A1, CYP11A1 and CYP17A1) and dysregulation of the genes associated with cholesterol biosynthesis and metabolism (acetyl-CoA acetyltransferases enzymes, ACTAs and HMGCR) in DHT-treated ovaries. These findings are consistent with the dependence of granulosa cell progesterone biosynthesis on the de novo cholesterol synthesis [32]. More importantly, the action of HMGCR is believed to be the rate limiting step in cholesterol biosynthesis [33,34].

Our present findings suggest that the dysregulation of the isoprene pathway in cholesterol biosynthesis, leading to abnormal accumulation of cholesterol precursors in hyperandrogenic state. Moreover, down-regulation of cholesterol biosynthesis may limit the amount of cholesterol influx and in turn affect ovarian cellular functions including fertility in DHT-treated rats. Here, a series of the enzymes involving in the conversion of mevalonate to arnesyl-diphosphate, isoprepene to saquale, and squalene in cholesterol were also repressed in ovaries of DHT-treated rats (Figure 7).

Figure 7
figure 7

The expression pattern of genes involving in cholesterol biosynthesis pathway in DHT-treated rat ovaries. (A) A modified cholesterol biosynthesis pathway outlined by [51,52] and the intermediate and final products performed by repressed genes in DHT-treated ovaries. The genes are indicated in green box while the intermediate products are described in white box. (B) The expression pattern of genes involving in cholesterol biosynthesis in DHT-treated ovaries compared to CTL. FC: The average fold change reduction in DHT compared to CTL. P value: the significant level. The red color indicates the activated transcript in CTL level while the green color shows the repressed expression pattern of a specific gene in DHT rat groups.

One of the major functions of the ovary is the production of steroid hormone [35]. Here, we identified dysregulation of several genes associated with steroidogenesis in ovaries of DHT-treated rats, including repression of cytochrome P450s, STAR and 3β-hydroxysteroid dehydrogenases (Figure 8). The repression of STAR in DHT-treated group suggests a reduced level of cholesterol influx for downstream steroidogenesis. In addition, the down-regulation of CYP11A1 and CYP19A1 in DHT-treated ovaries was consistent with our previous publication showing decreased aromatase expression and estradiol secretion in granulosa cells from DHT-treated rats [36]. Therefore, the dysregulation of genes associated with steroidogenesis in DHT-treated rats could result in abnormal sex hormone levels and ultimately PCOS phenotypes, including cycle irregularity, abnormal folliculogenesis and anovulation.

Figure 8
figure 8

The expression pattern of genes involving in sex steroid biosynthesis pathway in DHT-treated rat ovaries. (A) A modified steroid biosynthesis pathway from [53]. The intermediate and final products performed by dysregulated genes in DHT-treated ovaries are indicated in white box while the number in red or green circle corresponds the genes (enzymes) indicated in Figure B. (B) The heatmap displaying the expression pattern of repressed and activated genes involving in steroid biosynthesis in DHT compared to CTL. FC indicates the average fold change gene expression reduction or activation in DHT compared to CTL. P value: the significant levels.

Women with PCOS may have higher abdominal body fat distribution, due to hyperandrogenism and insulin resistance [37]. This phenomenon may be attributed by increased level of lactate, long-chain fatty acids and triglyceride [38-40]. We have previously demonstrated that the rats treated with DHT exhibited higher body weight compared to control [23]. Although data regarding the total fat content of the ovary is lacking, examination into the gene set enrichment analysis revealed 42 genes associated with lipid metabolism and biosynthesis were repressed in DHT-treated rats. Among those, 26 genes including acyl-CoA synthetase and fatty acid synthase are known to be involved in dual roles of lipid synthesis and metabolism, while 16 other genes were related only to lipid metabolism. The dysregulation of genes involving de novo lipid synthesis and metabolism in the DHT group may result in the accumulation of lipid precursors or lack of essential fatty acids which are required for normal ovarian physiology.

One of the characteristics of PCOS is the presence of atretic follicles or premature growth arrest without atresia [41,42]. These phenotypic manifestations could be due to defects in steroid biosynthesis and energy metabolism. In line with this notion, excess androgen, luteinizing hormone and insulin are associated with the recruitment of several but small preovulatory follicles [5]. Indeed, cell cycle progression and proliferation is thought to be controlled by several regulators [43]. In the current study, a total of 45 genes that associated with cell proliferation and differentiation, including AMH and BMP2, were activated in ovaries of DHT-treated rats (Figure 6). AMH inhibits primordial follicle recruitment and decreases the sensitivity of growing follicles to FSH [44]. AMH nulls and heterozygous mice exhibited early depletion of primordial follicles [45]. Similarly, it is possible that BMP2 gene activated in DHT-treated ovaries could participate in the regulation of folliculogenesis and luteinization by modulating gonadotropin receptor expression [46]. Activation of these genes may induce small follicle growth but dominant follicle growth arrest in DHT-treated rats, although this possibility needs further investigation.

In addition, genes involved in glycolysis/gluconeogenesis are also dysregulated by DHT treatment. It is known that hexokinases (HK1/2) convert glucose to glucose 6-phosphate [47] while pyruvate dehydrogenase complex component x (PDHX) catalyzes the conversion of pyruvate to acetyl coenzyme A [48,49]. In our study, ovarian HK1, HK2, PDHX and ACSS2 were repressed but FBP2 and ALDH1B1 were activated in DHT-treated rats (Table 3), further complicating the metabolic disorders in those groups. This finding is consistent with earlier report indicating the down-regulation of several genes regulating glucose synthesis and consumption in PCOS patients [50]. In addition, down-regulation of the oxidative reductase gene (Figure 3E) and those of citrate acid cycle pathway (Table 3) adds further evidence for dysregulated energy metabolism in DHT-treated ovaries.

In conclusion, we have provided detailed evidence for transcriptome profile changes in a chronically androgenized PCOS rat model. Our data suggest biosynthesis and metabolism of cholesterol, sterols/steroids, lipids and oxidation/reduction are key molecular functions associated with repressed gene expression in DHT-treated rats. On the other hand, cell differentiation/proliferation, transcriptional regulation, neurogenesis, cell adhesion and blood vessel development were enriched by activated genes in this animal model. It is therefore conceivable that these molecular functional alterations could be a molecular clue for abnormal steroidogenesis, estrous cycle irregularity, abnormal folliculogenesis, anovulation, and disorders in carbohydrate regulation and lipid metabolism occurring in PCOS patients. This study contributes significantly to our understanding of the ovarian transcriptome profile and associated molecular functional alterations in DHT-treated rats, and provides the basis for future in-depth functional and mechanistic studies that to shed light on the pathophysiologic significance of the current findings in PCOS.

References

  1. Azziz R, Woods KS, Reyna R, Key TJ, Knochenhauer ES, Yildiz BO. The prevalence and features of the polycystic ovary syndrome in an unselected population. J Clin Endocr Metab. 2004;89:2745–9.

    Article  CAS  PubMed  Google Scholar 

  2. Legro RS. Diagnostic criteria in polycystic ovary syndrome. Sem Reprod Med. 2003;21:267–75.

    Article  Google Scholar 

  3. Lujan ME, Chizen DR, Pierson RA. Diagnostic criteria for polycystic ovary syndrome: pitfalls and controversies. JOGC. 2008;30:671–9.

    PubMed Central  PubMed  Google Scholar 

  4. Legro RS, Chiu P, Kunselman AR, Bentley CM, Dodson WC, Dunaif A. Polycystic ovaries are common in women with hyperandrogenic chronic anovulation but do not predict metabolic or reproductive phenotype. J Clin Endocrinol Metab. 2005;90:2571–9.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  6. Azziz R, Carmina E, Dewailly D, Diamanti-Kandarakis E, Escobar-Morreale HF, Futterweit W, et al. The androgen excess and PCOS Society criteria for the polycystic ovary syndrome: the complete task force report. Fertil Steril. 2009;91:456–88.

    Article  PubMed  Google Scholar 

  7. Goodarzi MO. The genetic basis of the polycystic ovary syndrome. In Androgen excess disorders in women. Springer; 2007. p. 223-33.

  8. Welt CK, Carmina E. Lifecycle of Polycystic Ovary Syndrome (PCOS): From In Utero to Menopause. J Clin Endocrinol Metab. 2013;98:4629–38.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Cui L, Zhao H, Zhang B, Qu Z, Liu J, Liang X, et al. Genotype-phenotype correlations of PCOS susceptibility SNPs identified by GWAS in a large cohort of Han Chinese women. Hum Reprod. 2013;28:538–44.

    Article  CAS  PubMed  Google Scholar 

  10. Tumu VR, Govatati S, Guruvaiah P, Deenadayal M, Shivaji S, Bhanoori M. An interleukin-6 gene promoter polymorphism is associated with polycystic ovary syndrome in South Indian women. J Assist Reprod Genet. 2013;30:1541–6.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Radavelli-Bagatini S, de Oliveira IO, Ramos RB, Santos BR, Wagner MS, Lecke SB, et al. Haplotype TGTG from SNP 45 T/G and 276G/T of the adiponectin gene contributes to risk of polycystic ovary syndrome. J Endocrinol Invest. 2013;36:497–502.

    CAS  PubMed  Google Scholar 

  12. Nelson VL, Legro RS, Strauss 3rd JF, McAllister JM. Augmented androgen production is a stable steroidogenic phenotype of propagated theca cells from polycystic ovaries. Mol Endocrinol. 1999;13:946–57.

    Article  CAS  PubMed  Google Scholar 

  13. Catteau-Jonard S, Jamin SP, Leclerc A, Gonzales J, Dewailly D, di Clemente N. Anti-Mullerian hormone, its receptor, FSH receptor, and androgen receptor genes are overexpressed by granulosa cells from stimulated follicles in women with polycystic ovary syndrome. J Clin Endocrinol Metab. 2008;93:4456–61.

    Article  CAS  PubMed  Google Scholar 

  14. Wood JR, Nelson VL, Ho C, Jansen E, Wang CY, Urbanek M, et al. The molecular phenotype of polycystic ovary syndrome (PCOS) theca cells and new candidate PCOS genes defined by microarray analysis. J Biol Chem. 2003;278:26380–90.

    Article  CAS  PubMed  Google Scholar 

  15. Meng Y, Qian Y, Gao L, Cai LB, Cui YG, Liu JY. Downregulated expression of peroxiredoxin 4 in granulosa cells from polycystic ovary syndrome. PLoS One. 2013;8, e76460.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Abbott DH, Nicol LE, Levine JE, Xu N, Goodarzi MO, Dumesic DA. Nonhuman primate models of polycystic ovary syndrome. Mol Cell Endocrinol. 2013;373:21–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Maliqueo M, Benrick A, Stener-Victorin E. Rodent models of polycystic ovary syndrome: phenotypic presentation, pathophysiology, and the effects of different interventions. Sem Reprod Med. 2014;32:183–93.

    Article  Google Scholar 

  18. McNeilly AS, Duncan WC. Rodent models of polycystic ovary syndrome. Mol Cell Endocrinol. 2013;373:2–7.

    Article  CAS  PubMed  Google Scholar 

  19. Padmanabhan V, Veiga-Lopez A. Animal models of the polycystic ovary syndrome phenotype. Steroids. 2013;78:734–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Shi D, Vine DF. Animal models of polycystic ovary syndrome: a focused review of rodent models in relationship to clinical phenotypes and cardiometabolic risk. Fertil Steril. 2012;98:185–93.

    Article  PubMed  Google Scholar 

  21. Walters KA, Allan CM, Handelsman DJ. Rodent models for human polycystic ovary syndrome. Biol Reprod. 2012;86:149. 141-112.

    Article  PubMed  Google Scholar 

  22. Manneras L, Cajander S, Holmang A, Seleskovic Z, Lystig T, Lonn M, et al. A new rat model exhibiting both ovarian and metabolic characteristics of polycystic ovary syndrome. Endocrinology. 2007;148:3781–91.

    Article  PubMed  Google Scholar 

  23. Hossain MM, Cao M, Wang Q, Kim JY, Schellander K, Tesfaye D, et al. Altered expression of miRNAs in a dihydrotestosterone-induced rat PCOS model. J Ovarian Res. 2013;6:36.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Gharaibeh RZ, Fodor AA, Gibas CJ. Background correction using dinucleotide affinities improves the performance of GCRMA. BMC Bioinformatics. 2008;9:452.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer; 2005. p. 397–420.

    Chapter  Google Scholar 

  26. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Roy Statist Soc Ser B. 1995;57:289–300.

    Google Scholar 

  27. Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–8.

    Article  CAS  PubMed  Google Scholar 

  28. Da Wei Huang BTS, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Na Protoc. 2008;4:44–57.

    Article  Google Scholar 

  29. Caraux G, Pinloche S. PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005;21:1280–1.

    Article  CAS  PubMed  Google Scholar 

  30. Diao FY, Xu M, Hu Y, Li J, Xu Z, Lin M, et al. The molecular characteristics of polycystic ovary syndrome (PCOS) ovary defined by human ovary cDNA microarray. J Mol Endocrinol. 2004;33:59–72.

    Article  CAS  PubMed  Google Scholar 

  31. Jansen E, Laven JS, Dommerholt HB, Polman J, van Rijt C, van den Hurk C, et al. Abnormal gene expression profiles in human ovaries from polycystic ovary syndrome patients. Mol Endocrinol. 2004;18:3050–63.

    Article  CAS  PubMed  Google Scholar 

  32. Baranao JL, Hammond JM. FSH increases the synthesis and stores of cholesterol in porcine granulosa cells. Mo Cel Endocrinol. 1986;44:227–36.

    Article  CAS  Google Scholar 

  33. Tobert JA. Lovastatin and beyond: the history of the HMG-CoA reductase inhibitors. Nat Rev Drug Discov. 2003;2:517–26.

    Article  CAS  PubMed  Google Scholar 

  34. Olivier LM, Krisans SK. Peroxisomal protein targeting and identification of peroxisomal targeting signals in cholesterol biosynthetic enzymes. Biochim Biophys Acta. 2000;1529:89–102.

    Article  CAS  PubMed  Google Scholar 

  35. Sanderson JT. The steroid hormone biosynthesis pathway as a target for endocrine-disrupting chemicals. Toxico Sci. 2006;94:3–21.

    Article  CAS  Google Scholar 

  36. Wang Q, Kim JY, Xue K, Liu JY, Leader A, Tsang BK. Chemerin, a novel regulator of follicular steroidogenesis and its potential involvement in polycystic ovarian syndrome. Endocrinology. 2012;153:5600–11.

    Article  CAS  PubMed  Google Scholar 

  37. Diamanti-Kandarakis E. Role of obesity and adiposity in polycystic ovary syndrome. Int J Obes. 2007;31:S8–13.

    Article  Google Scholar 

  38. Cullberg G, Hamberger L, Mattsson LA, Mobacken H, Samsioe G. Lipid metabolic studies in women with a polycystic ovary syndrome during treatment with a low-dose desogestrel-ethinylestradiol combination. Acta Obstet Gynecol Scand. 1985;64:203–7.

    Article  CAS  PubMed  Google Scholar 

  39. Mattsson LA, Cullberg G, Hamberger L, Samsioe G, Silfverstolpe G. Lipid metabolism in women with polycystic ovary syndrome: possible implications for an increased risk of coronary heart disease. Fertli Steril. 1984;42:579–84.

    CAS  Google Scholar 

  40. Zhao Y, Fu L, Li R, Wang LN, Yang Y, Liu NN, et al. Metabolic profiles characterizing different phenotypes of polycystic ovary syndrome: plasma metabolomics analysis. BMC Med. 2012;10:153.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Tamura T, Kitawaki J, Yamamoto T, Osawa Y, Kominami S, Takemorit S, et al. Immunohistochemical localization of 17α-hydroxylase/C17-20 lyase and aromatase cytochrome P-450 in polycystic human ovaries. J Endocrinol. 1993;139:503–9.

    Article  CAS  PubMed  Google Scholar 

  42. Townson DH, Combelles CM. Ovarian follicular atresia. Basic Gynecology–Some Related Issues, Prof. Atef Darwish (Ed.), ISBN: 978-953-51-0166-6, InTech. doi:10.5772/32465. Available from: http://www.intechopen.com/books/basic-gynecology-some-related-issues/ovarian-follicular-atresia/.

  43. Salvetti NR, Stangaferro ML, Palomar MM, Alfaro NS, Rey F, Gimeno EJ, et al. Cell proliferation and survival mechanisms underlying the abnormal persistence of follicular cysts in bovines with cystic ovarian disease induced by ACTH. Anim Reprod Sci. 2010;122:98–110.

    Article  CAS  PubMed  Google Scholar 

  44. Durlinger AL, Visser JA, Themmen AP. Regulation of ovarian function: the role of anti-Mullerian hormone. Reproduction. 2002;124:601–9.

    Article  CAS  PubMed  Google Scholar 

  45. Durlinger AL, Kramer P, Karels B, de Jong FH, Uilenbroek JT, Grootegoed JA, et al. Control of primordial follicle recruitment by anti-Mullerian hormone in the mouse ovary. Endocrinology. 1999;140:5789–96.

    CAS  PubMed  Google Scholar 

  46. Shi J, Yoshino O, Osuga Y, Koga K, Hirota Y, Nose E, et al. Bone morphogenetic protein-2 (BMP-2) increases gene expression of FSH receptor and aromatase and decreases gene expression of LH receptor and StAR in human granulosa cells. Am J Reprod Immunol. 2011;65:421–7.

    Article  CAS  PubMed  Google Scholar 

  47. Walsh Jr CT, Spector LB. The glucose-glucose 6-phosphate exchange catalyzed by yeast hexokinase. Arch Biochem Biophys. 1971;145:1–5.

    Article  CAS  PubMed  Google Scholar 

  48. Harris RA, Bowker-Kinley MM, Huang B, Wu P. Regulation of the activity of the pyruvate dehydrogenase complex. Adv Enzyme Regul. 2002;42:249–59.

    Article  CAS  PubMed  Google Scholar 

  49. Liu S, Gong X, Yan X, Peng T, Baker JC, Li L, et al. Reaction mechanism for mammalian pyruvate dehydrogenase using natural lipoyl domain substrates. ArchBbioche Biophys. 2001;386:123–35.

    Article  CAS  Google Scholar 

  50. Kim JY, Song H, Kim H, Kang HJ, Jun JH, Hong SR, et al. Transcriptional profiling with a pathway-oriented analysis identifies dysregulated molecular phenotypes in the endometrium of patients with polycystic ovary syndrome. J Clin Endocrinol Metab. 2009;94:1416–26.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Liscurn L. Cholesterol biosynthesis. In: Biochemistry of lipids, lipoproteins and membranes. 4th ed. New York, N: Elsevier; 2002. p. 409–31.

    Chapter  Google Scholar 

  52. Waterham HR. Defects of cholesterol biosynthesis. FEBS Lett. 2006;580:5442–9.

    Article  CAS  PubMed  Google Scholar 

  53. Chen ZT, Wang IJ, Liao YT, Shih YF, Lin LL. Polymorphisms in steroidogenesis genes, sex steroid levels, and high myopia in the Taiwanese population. Mol Vis. 2011;17:2297–310.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the Canadian Institutes of Health Research (MOP-119381) and the World Class University (WCU) program through the Ministry of Education, Science and Technology and funded by the National Research Foundation of Korea (R31-10056) and research training awards [CIHR-REDIH Doctoral Scholarship and CIHR-QTNPR Doctoral Scholarship (QW)].

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Benjamin K Tsang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DSW performed the array hybridization, array data analysis, interpreting the results and drafting and writing the manuscript. QW prepared the PCOs rat model and revised the manuscript. DT assisted in the interpretation of the results, revised the manuscript and supervised the microarray data processing. KS supervised the gene expression analysis and revised the manuscript. MH revised the manuscript. MMH performed total RNA isolation, assisted in array hybridization and interpretation of the results and BKT conceived and oversaw the whole study, revised the manuscript and served as the corresponding author. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Significant biological processes enriched by genes repressed in DHT-treated ovaries.

Additional file 2: Table S2.

Significant molecular functions enriched by genes repressed in DHT-treated ovaries.

Additional file 3: Figure S1.

Cellular localization (cellular components) of dysregulated genes in DHT- treated ovaries.

Additional file 4: Table S3.

The expression pattern of some dysregulated genes in DHT rat ovaries in the current data set and their corresponding expression level in PCOS human ovary and long-term androgen-treated female-to-male transsexuals (TSX) compared to healthy ovaries [31].

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Salilew-Wondim, D., Wang, Q., Tesfaye, D. et al. Polycystic ovarian syndrome is accompanied by repression of gene signatures associated with biosynthesis and metabolism of steroids, cholesterol and lipids. J Ovarian Res 8, 24 (2015). https://doi.org/10.1186/s13048-015-0151-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13048-015-0151-5

Keywords