MiRNAs expression profiling of rat ovaries displaying PCOS with insulin resistance

Purpose The present study established microRNA (miRNA) expression profiles for rat ovaries displaying polycystic ovary syndrome (PCOS) with insulin resistance and explored the underlying biological functions of differentially expressed miRNAs. Methods A PCOS with insulin resistance rat model was created by administering letrozole and a high-fat diet. Total RNA was extracted from the ovaries of PCOS with insulin resistance rats and normal rats. Three ovaries from each group were used to identify differentially expressed miRNAs by deep sequencing. A hierarchical clustering heatmap and volcano plot were used to display the pattern of differentially expressed miRNAs. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to explore the potential target genes of the differentially expressed miRNAs and identify their putative biological function. Nine of the differentially expressed miRNAs were selected for validation by Real-time Quantitative PCR (qRT-PCR). Results A total of 58 differentially expressed miRNAs were identified in the rat ovaries exhibiting PCOS with insulin resistance compared with control ovaries, including 23 miRNAs that were upregulated and 35 miRNAs that were downregulated. GO and KEGG pathway analyses revealed that the predicted target genes were related to metabolic processes, cellular processes, and metabolic pathways. Furthermore, qRT-PCR confirmed that miR-3585-5p and miR-30-5p were significantly upregulated and miR-146-5p was downregulated in the ovaries of PCOS with insulin resistance rats compared with the controls. Conclusion These results indicate that differentially expressed miRNAs in rat ovaries may be involved in the pathophysiology of insulin resistance in PCOS. Our study may be beneficial in establishing miRNAs as novel diagnostic and therapeutic biomarkers for insulin resistance in PCOS.


Introduction
Polycystic ovary syndrome (PCOS), which affects 6-20% of women during reproductive age, is one of the most common endocrine and metabolic disorders [1][2][3]. According to the Rotterdam criteria, PCOS may be diagnosed with at least two of the three following characteristics: anovulation, clinical and/or biochemical signs of hyperandrogenism and polycystic ovaries [4,5]. PCOS is also associated with infertility, insulin resistance, obesity, glucose tolerance, dyslipidemia, type 2 diabetes and increased risk for cardiovascular disease [6]. Insulin resistance and hyperandrogenism are not only the two main characteristic features of PCOS, but also important pathogenic factors for PCOS [7]. Furthermore, there are approximately 40-70% of PCOS patients with insulin 1 3 resistance. Insulin may promote luteinizing hormone release, increase androgen levels by ovarian theca cells, and inhibit the production of sex hormone binding protein by the liver, resulting in elevated levels of circulating free testosterone. Therefore, insulin synergizing with androgen could result in deteriorating anovulation, hyperandrogenism, hyperinsulinemia and infertility in PCOS patients. However, whether insulin resistance is the cause or the result of PCOS remains to be determined. Genetic studies and familial clustering of phenotypic features have indicated that PCOS is heritable [8][9][10][11]. Genome-wide association studies have found over 20 PCOS susceptibility genes; however, they account for < 10% of all PCOS cases [12][13][14][15][16]. Therefore, other mechanisms must play a role in the development of PCOS as genetic factors are responsible for only a small portion of the cases. One potential mechanism is an epigenetic process which could yield the same phenotypic characteristics as heritable genetics [17,18].
MicroRNAs (miRNAs) are small, single-stranded noncoding RNAs consisting of 18-24 nucleotides in length. They regulate post-transcriptional gene expression by binding to the 3′ untranslated region of target RNAs [19]. MiR-NAs are involved in many physiological processes including cell growth, differentiation, proliferation and metabolism [20]. Furthermore, several studies have shown that dysfunctional miRNAs are also associated with the pathological mechanism of PCOS [21]. Sang et al. identified highly expressed miRNAs in human follicular fluids [22]. In addition, Chen et al. found that miR-93 is upregulated in the adipose tissue of PCOS patients and women with insulin resistance [23]. In our previous study, we found differentially expressed circular RNAs (circRNAs) in the granulosa cells of reproductive-aged PCOS women [24]. However, there is no literature regarding miRNA expression profiling in the ovaries of PCOS patients with insulin resistance. In the present study, we identified differentially expressed miRNAs from PCOS ovaries in an insulin-resistant rat model and determined their putative function.

Animals and treatment
All animal experiments were carried out based on the Guidelines for the Care and Use of Experimental Animals and approved by the Animal Ethics Committee of the Guangzhou University of Chinese Medicine (20180715002). Three-week-old female Sprague Dawley rats were purchased from the Animal Centre of Guangdong Province (Guangdong, China). Animals were fed adaptively for 1 week under conditions of 55-65% humidity, 21-22 °C temperature, and a 12 h light/12 h dark cycle. Twenty rats were randomly divided into two groups. Ten rats in the control group (CON) were fed a normal diet (Research Diet, D12450, 10% fat), while 10 rats in the treatment group were fed a high-fat diet (Research Diet, D12492, 60% fat) for 8 weeks (LEHF group). In addition, the con rats were gavaged with 1% carboxymethylcellulose (CMC) solution once daily, and LEHF rats were gavaged with 1 mg/kg letrozole in 1% CMC once daily for 21 days [25]. All rats were weighed each week.

Estrus cycle
The estrus cycle was examined daily. Depending on the results of microscopic analysis of a vaginal smear, the animals were flushed with physiological saline from the 12th day to the end of the experiment.

Blood and tissue sampling
All rats were anesthetized with isoflurane after overnight fasting. Blood samples were collected from the heart and the resulting serum samples were stored at -80 °C. Serum levels of follicle stimulating hormone (FSH), luteinizing hormone (LH), estradiol (E2), testosterone (T), progesterone (P), fasting blood glucose (FBG) and fasting insulin (FINS) were measured using ELISA kits (Cloud-Clone Corp., Wuhan, China) according to the manufacturer's protocol. The homeostasis model assessment for insulin resistance (HOMA-IR) was calculated as previously described using the following formula: HOMA-IR = (FBG × FINS) / 22.5 [26].
Ovary samples were quickly collected from the animals, and one ovary from each rat was fixed in 4% paraformaldehyde. The ovary tissue samples were then dehydrated, embedded in paraffin, and sectioned into 4 µm sections. Each section was stained with hematoxylin and eosin (H&E). The other rat ovary was stored at -80 °C and used to establish miRNA expression profiles.

RNA-sequencing (RNA-seq) and bioinformatics analysis
RNA isolation and RNA sequencing were performed as previously described [24]. Briefly, total RNA was extracted from ovaries using TRIzol reagent according to the manufacturer's instructions. The concentration of RNA was measured with a NanoDrop 2000 spectrophotometer. An RNA library was generated from total RNA using the NEB Next Ultra Directional RNA Library Prep Kit for Illumina (NEB, MA, USA). The RNA library quality and quantity were assessed with an Agilent 2100 Bioanalyzer and an ABI Step One Plus Real-Time PCR System, respectively. The RNA library was run on a HiSeq 2000 platform (Illumina, CA, USA) for sequence analysis. The clean reads were filtered from raw reads by FastQC and used for further bioinformatics analysis. The differentially expressed miR-NAs were selected based on log 2 (LEHF/CON) > 1 and P value < 0.05.
Hierarchical clustering heatmaps and volcano plots generated by R package version 1.0.8 software (https ://cran.rproje ct.org/web/packa ges/pheat map/) were used to display the differentially expressed miRNA patterns between the two groups. Gene ontology (GO) analysis was performed to explore the molecular function, cellular components and biological processes of the differentially expressed miRNAs (https ://www.geneo ntolo gy.org). The biological pathways of the differentially expressed miRNAs were further analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG; https ://www.genom e.jp/kegg/) database.

Real-time quantitative PCR (qRT-PCR)
To verify the results of RNA sequencing, we measured relative miRNA expressions by qRT-PCR analysis using SYBR Premix Ex Taq II (Takara, China). The U6 gene was used as an internal control. The relative expression of miRNAs was calculated using the comparative Ct (2 −ΔΔCT ) method [27]. The sequences of primers used in this study are shown in Table 1.

Statistical analysis
All values are presented as the mean ± SD. Comparisons between experimental and control groups were conducted by a Student's t test. Data analyses were performed using SPSS 19.0 software (SPSS Inc., IL, USA). P value < 0.05 was considered statistically significant.

Animal model of PCOS with insulin resistance
We used letrozole and a high-fat diet to establish a rat model with a PCOS and insulin-resistance phenotype. The CON rats exhibited a normal estrus cycle, whereas the LEHF rats had a disrupted estrous cycle (Fig. 1a). H&E staining results showed that the ovary in the CON group had follicles and corpora lutea at various stages of development, a normal theca and granulosa cell layer and no ovarian cysts. In the LEHF group, the ovary exhibited atretic antral follicles, more large cystic follicles with a thinner granulosa cell layer and fewer corpora lutea (Fig. 1b). Compared with the CON group, the body weight, FINS and HOMA-IR were significantly increased in rats from the LEHF group (Fig. 2a, b). Moreover, rats in the LEHF group exhibited a dramatic increase in the LH, FSH, E2, T and P levels (Fig. 2c). Taken together, these results indicated that we successfully established a rat PCOS model with insulin resistance. Table 1 The primer sequences in this study

Differentially expressed miRNA profiles in the ovary of the LEHF and CON groups
We conducted high-throughput sequencing to establish miRNA expression profiles for three ovary samples from the LEHF and CON groups. A total of 767 miRNAs were detected in the ovary samples. Among these, 58 differentially expressed miRNAs were found between the two groups (fold change ≥ 2 and P value < 0.05). Among the differentially expressed miRNAs, 23 were upregulated, while 35 were significantly downregulated in the LEHF group. The list of the top 10 differentially expressed  .3 Differentially expressed miRNAs were identified by deep sequencing. A hierarchical clustering heatmap (a) and a volcano plot (b) were used to display the differentially expressed miRNA patterns between the two groups miRNAs identified by deep sequencing is shown in Table 2. A hierarchical clustering heatmap (Fig. 3a) indicated that the expression of miRNAs in the LEHF group was significantly different from that of the CON group. A volcano plot (Fig. 3b) also displayed the pattern of upregulated and downregulated miRNAs.

Functional analysis of differentially expressed miRNAs
Ovary dysfunction and abnormal follicular development are two characteristics of PCOS that are aggravated by insulin resistance. To gain insight into the underlying mechanism of PCOS with insulin resistance, GO analysis was conducted to identify the biological processes, cellular components and molecular functions of the differentially expressed miRNAs (Fig. 4a). We discovered that the most enriched GO term for the cellular component was cell apart (GO:0044464), the most enriched term for the molecular function was cellular process (GO:0009987) and the most enriched term for the biological process was metabolic process (GO:0008152). Moreover, the metabolic pathway (path: rno01100) was the most enriched pathway in the KEGG database (Fig. 4b).
The biological functions and pathways for the differentially expressed miRNAs are closely related to the development of insulin resistance in PCOS.

Discussion
PCOS is a complex heterogeneous disorder characterized by hyperandrogenism and insulin resistance [30]. Insulin resistance can cause infertility and increased androgen levels, which exacerbate PCOS traits [31,32]. In addition, the insulin-sensitizing drug metformin can improve menstrual disorder, promote ovulation, and elevate birth rate in PCOS, especially when accompanied by insulin resistance [33,34]. Therefore, there should be a strong relationship between PCOS and insulin resistance. However, the underlying mechanism of insulin resistance in PCOS remains unclear. In the present study, we successfully established a model for PCOS with insulin resistance in rats by administrating letrozole and a high-fat diet to the rats. Letrozole is an aromatase inhibitor that blocks the conversion of androgen Fig.4 Gene ontology (GO) and KEGG pathway analysis for predicted targets of the differentially expressed miRNAs to estrogen. This results in elevated androgen levels, which is the main pathogenetic mechanism occurring in PCOS. Pubertal female rats were continuously treated with letrozole, inducing increased LH and T levels, abnormal estrous cycle, anovulation, increased ovarian weight, cysts, atretic follicles and no CL in adulthood [35,36]. Furthermore, a high-fat diet in rats induces abnormal metabolism, including increased body weight, impaired glucose tolerance, and insulin resistance [37]. Thus, letrozole and a high-fat diet in rats contain both the reproductive and metabolic features associated with PCOS and insulin resistance. Several miRNA, long-non-coding-RNA and circRNA expression profiling studies of PCOS have been conducted using different human components including serum, adipose tissue, uterine tissue, and follicular fluid [24,38,39]. Studies have implicated miRNAs in the pathogenesis of PCOS, including increased androgen [38,39]. However, the exact role of miRNAs in the pathophysiology of PCOS with insulin resistance remains unclear.
The terms biological process, cellular component, and molecular function were predicted by a GO analysis of the differentially expressed miRNAs. Of these, the most significantly enriched term was metabolic process, indicating that these differentially expressed miRNAs participate in metabolism. Likewise, the most enriched pathway from the KEGG pathway analysis was also metabolism. It is known that insulin resistance, obesity, type 2 diabetes mellitus and dyslipidemia are common features of PCOS [1][2][3]40]. Moreover, insulin resistance is an important characteristic of PCOS and contributes to the development of PCOS. Thus, these aberrantly expressed miRNAs may be the cause of insulin resistance in PCOS, though elucidating the underlying mechanism awaits further study.
We established miRNA expression profiles from the ovaries of PCOS rats with insulin resistance along with control rats. A total of 58 differentially expressed miRNAs in the PCOS ovaries exhibited significant changes compared with the control group, including 23 that were upregulated and 35 that were downregulated. Of note, we found that miR-3585-5p and miR-30-5p were significantly upregulated, and miR-146-5p was downregulated in the LEHF group, which was consistent with the deep sequencing results. MiR-30 was significantly upregulated in ovaries from PCOS with insulin resistance rats which is consistent with previous studies [28,41,42]. However, miR-146 was significantly upregulated in serum from PCOS patients compared with that of healthy women, which is inconsistent with our findings [28,42]. This difference may be explained by the natural heterogeneity of PCOS as well as the different sequencing method employed. Interestingly, the expression of miR-146 was significantly decreased in peripheral blood mononuclear cells from type-2-diabetes patients compared with control subjects, and miR-146 was negatively correlated to insulin resistance [43]. This indicates that miR-146 may play a key role in the pathogenesis of insulin resistance in PCOS, and miR-146 may represent a novel diagnostic marker or therapeutic target for insulin resistance in PCOS. Therefore, miR-146 should be further investigated in this context. However, this study has some limitations. In view of the methods and results we have got, this paper should be described as a pilot study, and this study is not large enough to do justice to the heterogeneity of PCOS. But these results might provide some assistance for further study of insulin downregulated miRNAs (miR-132-3p, miR-146-5p, miR-21-5p, miR-124-3p and miR-122-5p) identified by deep sequencing were determined. *P < 0.05; **P < 0.01 resistance in PCOS. Therefore, we will make further experiments to study insulin resistance in PCOS and the comprehensive experiments will be obtained and reported in future.
In conclusion, we found 58 differentially expressed miR-NAs in a rat model of PCOS with insulin resistance by deep sequencing analysis and verified the differential expression of miR-146-5p, miR-3585-5p and miR-30-5p by qRT-PCR. Furthermore, miR-146 may represent a novel diagnostic marker and therapeutic target for insulin resistance in PCOS. Additional studies are needed to identify the underlying mechanism of miRNA function in the development of PCOS with insulin resistance.

Compliance with ethical standards
Conflict of interest The authors declare no potential conflicts of interest.
Ethical approval All animal experiments were conducted under the principles of the Guidelines to the Care and Use of Experimental Animals and approved by the animal Ethics Committee of Guangzhou University of Chinese Medicine.
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 http://creat iveco mmons .org/licen ses/by/4.0/.