Altered expression of miRNAs in a dihydrotestosterone-induced rat PCOS model
© Hossain et al.; licensee BioMed Central Ltd. 2013
Received: 13 March 2013
Accepted: 6 May 2013
Published: 15 May 2013
The polycystic ovary syndrome (PCOS) is a complex and heterogeneous endocrine condition characterized by hyperandrogenism, hyperinsulinemia, insulin resistance and chronic anovulation. Regulation and interaction of a multitude of genes required for follicular development are found to be altered in PCOS. MicroRNAs (miRNAs) mediate posttranscriptional gene regulation by binding to the 3´ untranslated region of mRNAs to either inhibit or enhance translation. However, the extent and regulation of miRNA expression in PCOS is poorly understood and the current study is the first to describe altered expression of miRNAs in PCOS.
A chronically androgenized [5α-dihydrotestosterone (DHT)-treated] rat model which recapitulates many of the phenotypes of human PCOS, and miRNA PCR array was used to investigate the expression of 349 miRNAs in DHT treated rat ovaries. The ovarian expression of several selected miRNAs was also analyzed by in situ localization experiment.
DHT-treated rats exhibit increased body weight, disrupted estrus cyclicity, decreased insulin sensitivity and decreased ovarian weight, with the latter phenomenon readily rescued by gonadotropin treatment in vivo. In general, 24% of the 349 miRNAs investigated were found to be differentially expressed between DHT-treated and control rats. Most of the differentially expressed miRNAs were found to be predominantly localized in the theca cells of the follicles. In silico analysis of the potential target genes of dysregulated miRNAs revealed their possible involvement in various pathways in the regulation of ovarian function.
Our current findings suggest that miRNAs are differentially regulated in hyperandrogenism, a condition possibly involved in the dysregulation of steroid hormone receptors and intra-ovarian factors, and that miRNAs may be involved in the etiology of PCOS.
KeywordsPCOS miRNA Dihydrotestosterone
Polycystic ovary syndrome (PCOS), a multi-factorial endocrine condition, affects up to 5-10% of women in reproductive age and accounts for 75% of anovulatory infertility [1–3]. PCOS is a heterogeneous syndrome with complex pathologies with follicle growth arrest at the small antral stage, minimal granulosa cell proliferation, hyperthecosis and hyperandrogenemia, and chronic anovulation . It is associated not only with infertility, but also with increased risk of metabolic disorders, such as insulin resistance, diabetes, obesity, and cardiovascular diseases . Both genetic and environmental factors are known to contribute to the pathogenesis of PCOS; however, its etiology is unclear.
Differential gene expression is evidenced in human PCOS ovaries when compared to the normal ones . Genomic studies in PCOS have revealed a variety of altered genes that fall into many functional categories, including cell division, apoptosis and regulation of gene expression and metabolism. While genes involved in steroidogenesis, including those related to retinoic acid biosynthesis and luteinizing hormone-responsive gene, are generally up-regulated in theca cells of PCOS subjects , those associated with Wnt signalling and ovarian folliculogenesis appear to be down-regulated [7, 8]. Although ovarian apoptotic genes are altered in PCOS, the net effect of these changes is unclear . In addition, significant differences in mRNA abundance in oocytes from PCOS patients have also been reported and cluster analysis of 374 genes indicated an association with chromosomal alignment and segregation in mitosis and/or meiosis .
Despite ample evidences on gene dysregulation in PCOS ovary as reported above, precisely how these genes are transcriptional and post-transcriptional regulated is poorly understood. New advances in post-transcriptional gene regulation have resulted from the discovery of hundreds of miRNAs in different mammalian species. Diverse expression pattern of miRNAs and the high abundance of potential targets in the ovary have led to the notion of their importance in the regulation of ovarian function. A large number of miRNAs has now been identified in the ovary [10, 11] and found to be regulated by gonadotropins during follicular development [12, 13]. Recent studies indicate that miRNAs are involved in the regulation of steriodogenesis, cell proliferation and apoptosis in human granulosa cells [14, 15]. Differential expression of miR-23a, miR-23b,miR-542–3p, miR-211, and miR-17–5p in granulosa/cumulus cells from women undergoing assisted reproduction suggests aberrant miRNA expression may be an underlying etiology in female infertility [16, 17]. Furthermore, miRNAs have been implicated to play a role in obesity and metabolic syndrome, including insulin resistance , conditions often associated with PCOS. However, whether miRNAs contribute to the follicle arrest, dysregulated steroid production and insulin resistance in PCOS is not known. Understanding how miRNAs are regulated in the ovary and the identification of their specific targets and functions may offer novel insights into the etiology of PCOS and the development of target-specific gene regulation for its treatment.
Hyperandrogenism is one key feature in human PCOS patients and circulating levels of androgens, such as testosterone, androstenedione and dihydrotestosterone (DHT), are elevated in these patients [19, 20]. A number of animal models have been developed that allow one to focus on different aspects of PCOS pathology, including ovarian follicular and metabolic features [21–23]. A recent report on a rat PCOS model involving chronic DHT implant mimics the hyperandrogenic status and exhibits both polycystic ovaries and metabolic abnormalities similar to those in human PCOS , allowing one to examine not only ovarian characteristics but also metabolic features of PCOS.
We hypothesize that PCOS is associated with dysregulation of ovarian miRNAs which may be involved in post-transcriptional regulation of genes involved in the ovarian pathology. Taking advantage of DHT-induced rat PCOS model, we tested this hypothesis by comparing the expression of 349 miRNAs between DHT-treated vs untreated rat ovaries. Moreover, selected candidate miRNAs were investigated for their localization in follicular cells. Results of the present study revealed that PCOS is associated with differential expression of regulatory non-coding miRNAs in rat ovaries.
Materials and methods
DHT-filled silastic capsules preparation
The DHT-filled capsules were prepared as previously described . Silastic tubing (I.D. 1.98 mm x O.D. 3.18 mm; Dow Corning, Cat. No. 2415577, Midland, MI), cut to an appropriate length to achieve a desired surface area of 299.70 mm2, was filled with DHT (Sigma) from 1 ml syringe. The tubing was closed at each end (3 mm) with a sealant (Silicone Type A, Dow Corning Cat. No. 891, Midland, MI), ensuring that adhesive is in contact with tubing walls and that there are no air bubbles. Control capsules were empty with only sealant on both ends. After overnight dry, capsules were rinsed for 2 days in 3% bovine serum albumin (BSA) in phosphate buffer saline (PBS) with 0.1% NaN3 solution, thoroughly washed with PBS and sterilized by dipping briefly in 70% ethanol before use.
Sprague Dawley rats (Charles River, Montreal, Canada) were maintained on 12-h light, 12-h dark cycles and given food and water ad libitum. All procedures were carried out in accordance with the Guidelines for the Care and Use of Laboratory Animals, Canadian Council on Animal Care, and were approved by University of Ottawa and Ottawa Hospital Research Institute Animal Care Committee.
Animal surgery, DHT capsule implantation and animal care
Thirty-six female rats at 21 days of age were randomly divided into two groups (control vs. DHT group) and were implanted with a DHT-filled silicone capsule continuous-releasing (83 μg per day; empty capsule for controls) 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 . Animals were weighed weekly to monitor body weight gain and euthanized at 12 wks post-implantation. Ovaries were collected for analyses.
Determination of estrous cycle phases of rat
The estrous cycle of rats was monitored by microscopic analysis of vaginal smear obtained every morning for 12 days from 10–11 wks post-implantation. Predominant cells associated with each phase of the estrus cycle were: epithelial cells (proestrous; P), enucleated cornified cells (estrous; E), equal mix of leukocytes and cornified epithelial cells (metestrous; M) and leukocytes (diestrous; D) .
Determination of insulin sensitivity
Insulin sensitivity test (IST) was performed after 12 h overnight fasting at 10 wks of surgery. Human insulin (0.2 U/100 g body weight intravenously, Novo Nordisk Canada Inc. Mississauga, ON, Canada) was administered to control or DHT-treated rats through tail vein and blood samples were obtained from the tail at 0, 5, 10, 20, 40, 80, 160 min. after insulin injection. Glucose levels were determined by glucose test strips (ACCU-CHEK, Roche). The insulin sensitivity index K ITT was calculated from the slope of a linear portion of the regression line of natural logarithm of glucose versus time. Raw data were transformed by Log K ITT for analysis. K ITT = (0.693/t1/2) × 100, t1/2 represents the half-life of glucose decay after insulin injection .
Ovarian tissue collection and morphological evaluation
Pairs of ovaries were excised, and weighed. For morphological evaluation, whole ovaries were embedded in paraffin block and longitudinal sections (4 μm thick) were stained with hematoxylin and eosin. Since the considerable discrepancy in estimating follicle numbers exists between studies due to different section thickness and correction factors adopted , the percentage of number of cystic follicles (FC) and corpus lutea (CL) over total follicular structure (including all kinds of follicles and CL) per ovary were determined. A FC was defined as an collapsed anovulatory follicle exhibiting the absence of oocyte, reduced granulosa cell layer, but intact theca layer. The maximal longitudinal ovarian sections were investigated to count the numbers of FC and CL, mean of the percentage of FC and CL from 3–5 sections after sampling every fifth section were defined as the relative values per ovary from 5 individual animals. In an additional experiment, three rats from DHT group at each time point were administered with 20 IU of equine chorionic gonadotropin (eCG; Sigma; PBS as control) for designated time (0–30 h) to test their ovarian responsiveness to gonadotropin in vivo.
miRNA expression profiling
Ovarian cortex tissues were snap-frozen in liquid nitrogen and stored at −80°C until analysis. Total RNAs including small RNAs from 3 independent DHT (DHT1-3) treated and 3 non-treated control (CTL1-3) ovaries have been isolated using miRNeasy mini kit (Qiagen, Hilden, Germany) following manufacturer’s protocol. Immediately after isolation, contaminating DNA and divalent cations from the RNA samples were removed using TURBO DNA-free™ kit (Ambion, Foster City, CA) according to instruction of the manufacturer. The quality and the concentration of the RNA samples were assessed by NanoDrop 8000 spectrophotometer (NanoDrop, Wilmington, Delaware, USA). Subsequently, a total of 250 ng purified RNA from each samples were reverse transcribed for qPCR analysis, using the QuantiMir™ RT systems (System Biosciences, Mountain View, CA).
MiRNA expression analysis was performed using miRNA PCR array platform containing 349 well characterized Rat miRNAs including 3 endogenous controls, namely Rat U6 snRNA, RNU43 snoRNA & U1 snRNA (System Biosciences, Mountain View, CA, Cat # RA680A-1) according to the user manual of miRNome microRNA Profilers QuantiMir™ (System Biosciences, Mountain View, CA) with some modifications. Briefly, 20 μl of cDNA from each sample mixed with 1300 μl of Maxima 2x SYBR Green, (Fermentas Gmbh, Leon-Rot, Germany, Cat # K0222), 50 μl of universal reverse primer and 800 μl of nuclease-free water. Five microlitres of mixed cDNA cocktail were dispersed into each well of 384 well PCR plate containing 1 μl of pre-aliquoted miRNA specific primer. Multichannel laboratory automation workstation was used for pipetting the mix and mixing the reaction into 384 well plates (Biomek® NXP, Beckman Coulter, Krefeld, Germany). The assays were performed in ABI 7900 HT real time PCR system (Applied Biosystems, Foster City, CA, USA) with thermal cycling programme at 52°C for 2 min, 95°C for 10 min, 40 cycles of (95°C for 15 sec, 60°C for 1 min). Dissociation stage was included at the end of the programme to assess the Tm of the PCR amplicons and to verify the specificity of the amplification reactions. In addition, amplification was verified for randomly selected reactions through 2.5% agarose gel electrophoresis. Data were analysed by ΔΔCt method and normalized by the geometric mean of the three endogenous controls through System Bioscience’s free Software developed to analyse PCR Array data for Rat miRNome (System Biosciences, Mountain View, CA), equipped with t test. A fold regulation of ≥ 2. P ≤ 0.05 was considered as significant in expression differences.
Localization of selected miRNAs in ovarian cryo-section
Ovaries were fixed in 4% paraformaldehyde (pH 7.3; 4°C, overnight) followed by incubation in 0.5M Sucrose in DEPC-treated 1X PBS (4°C, overnight). They were embedded in Tissue-Tek Oct reagent (Shandon Cryomatrix, Thermo Cat# 6769006) in a tissue path base mold, underwent a slow freezing process (−1°C/min rate) in a freezing container (NALGENE™ Cryo, Cat# 5100–0001). Serial ovarian sections (10 μm) were cut at −23°C using rapid sectioning cryostat (Leica microsystem Nussloch GmbH, Heidelberger, Germany) and mounted on SuperFrost® Plus slides (Menzel-Glaeser, Germany) at −80°C until hybridization began. Sections were checked for quality and required cell content at a regular interval by staining with TB in a light microscope. Post-fixation, acetylation and proteinase K treatment was done as described previously . Two hours of pre-hybridization was performed at 52°C in hybridization solution (50% formamide, 5x sodium chloride/sodium citrate [SSC; pH 6.0], 0.1% Tween-20, 50 μg/ml heparin, and 500 mg/ml yeast tRNA). Ovarian sections were incubated overnight at 52°C in a hybridization buffer containing 3'-Digoxigenin (DIG) labeled LNA-modified oligonucleotide probes (1pM) for miRNA, U6 RNA and scrambled miRNA (Exiqon, Vedbaek, Denmark) in a humidified chamber. Blocking, incubation with anti-DIG-AP antibody, washing and color development with alkaline phosphatase substrates solution [5-Bromo-4-chloro-3-indolyl phosphate/Nitro blue tetrazolium (BCIP/NBT)] were performed as described previously . Vector® Methyl Green (Vector laboratories, Burlingame, CA) was used as nuclear counter-stain.
All data were analyzed using Prism 5.0 statistical software (GraphPad, San Diego, CA). The results were presented as means ± SEM as detailed in figure legends. Differences in body weight gains, follicle cysts and corpus luteum percentage were analyzed by two-way ANOVA, ovarian weight between groups were analyzed by one-way ANOVA and followed by Bonferroni post-hoc test. Insulin sensitivity in PCOS and control groups was compared by Student t-test. Statistical significant difference was defined at P < 0.05.
Characterisation of PCOS rat model
Body weight increased in PCOS rats
Irregularity in estrous cycle in DHT-treated rat
Estrous cycle was evaluated by microscopic analysis of the predominant cell type in vaginal smears obtained daily from 10–11 weeks post-implantation. Estrous cycle pattern in a representative rat from each group is shown (Figure 1B). During 12 days of evaluation, control rat # 22 had 3 regular estrous cycles, whereas DHT-treated rat # 30 started at metestrous and then showed prolong diestrous phase. 94% of animals in DHT group shown disrupted estrous cycle compared to 22% of animals in control group (Figure 1B, P < 0.001, N = 17).
DHT-treated rats showed insulin resistance
In order to examine the metabolic phenotype of rats which received DHT, insulin sensitivity test was performed. Rat plasma glucose levels were determined by glucose test strips after insulin challenge. The insulin sensitivity index K ITT was calculated and a lower K ITT suggests decreased insulin sensitivity or insulin resistance . Our data shows reduced insulin sensitivity in DHT-treated rats compared to control (Control rats Log K ITT = 0.34 ± 0.03 vs. DHT-treated rats 0.23 ± 0.05, N = 17, P < 0.05) (Figure 1C). As indicated as dotted line in Figure 1C, it took 10 minutes for the Control rat # 12 to change the basal glucose level from 6.9 mmol/L to 3 mmol/L, whereas it took 40 minutes for the glucose levels of DHT-treated rat # 5 to decrease from 6.3 to 3 mmol/L. A glucose curve shift to the right was evident in DHT-treated rats when compared to control animals.
Ovary weights were lower in DHT-treated rats
Pair ovary weight in DHT-treated rats was significantly lower than their control counterparts (0 h: 18.93 ± 4.85 mg vs. 91.90 ± 4.48 mg. P < 0.001) (Figure 1D). To test if DHT-treated ovaries are responsive to gonadotropin, three DHT-treated rats per group were administered with 20 IU of eCG (PBS as control) for designated time (0–30 h) and the ovaries were weighed in an additional experiment. Interestingly, the reduction of ovary weight in DHT-treated rats was markedly rescued with longer eCG treatment, indicating that this PCOS phenotype was responsive to gonadotropin (Figure 1D).
DHT-treated rats exhibited anovulation
To assess ovarian feature of PCOS rat, ovarian morphology in sections from CTL or DHT-treated rats was examined by H&E staining. Percentage of cystic follicle (FC) and corpus luteum (CL) over the total follicular structure were estimated in the longitudinal ovarian sections, as described in Materials and Methods. Estimation of relative value of FC and CL number, regardless of the total follicle numbers and their stages, is sufficient to evaluate disruption of ovarian follicular development as seen in PCOS. In control rat, CL accounted for 35% of total follicular structure per section, whereas those in DHT-treated rats were 4.3% (Figure 1E, P < 0.05). In contrast, follicle cysts in DHT-treated rats accounted for 43% of total follicle structure, which was significantly higher than in control rat (7%) (P < 0.05). Gonadotropin treatment reactivated follicular development in the cystic ovaries and induced ovulation, as the presence of CL increased to 43.6% in eCG-treated rats, while FC percentage declined to 6.3% over the total follicle structure (Figure 1E).
Differential expression of miRNAs in the rat PCOS
In situ localization of differentially expressed miRNAs in rat PCOS
Functional classification of differentially expressed miRNA
Molecular pathways regulated by miRNAs in DHT-treated ovaries
List of miRNA
Ovarian cancer signalling
rno-miR-96↓, rno-miR-25↑, rno-miR-24↓, rno-miR-221↓, rno-miR-21↓, rno-miR-183↓, rno-miR-182↓, rno-miR-132↓, rno-let-7d↑
rno-let-7d↑, rno-miR-132↓, rno-miR-182↓, rno-miR-183↓, rno-miR-184↓, rno-miR-21↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↓, rno-miR-26b↑, rno-miR-31↓, rno-miR-335↑, rno-miR-384-3p↓
Androgen & estrogen metabolism
rno-miR-96↓, rno-let-7d↑, rno-miR-24↓, rno-miR-26b↑, rno-miR-21↓, rno-miR-26b↑, rno-miR-31↓
Aldosterone signalling (epithelial cells)
rno-miR-182↓, rno-miR-183↓, rno-let-7d↑, rno-miR-132↓, rno-miR-182↓, rno-miR-21↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-335↑, rno-miR-384-3p↓, rno-miR-96↓
rno-miR-96↓, rno-miR-182↓, rno-miR-31↓
rno-miR-132↓, rno-miR-182v, rno-miR-183↓, rno-miR-384-3p↓, rno-miR-96↓
rno-miR-96↓, rno-miR-384-3p↓, rno-miR-335↑, rno-miR-31↓, rno-miR-26b↑, rno-miR-25↑, rno-miR-24↓, rno-miR-221↓, rno-miR-21↓, rno-miR-184↓, rno-miR-183↓, rno-miR-182↓, rno-miR-132↓, rno-let-7d↑
rno-miR-25↑, rno-let-7d↑, rno-miR-182↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-31↓, rno-miR-335↑, rno-miR-384-3p↓, rno-miR-96↓
rno-miR-384-3p↓, rno-let-7d↓, rno-miR-132v, rno-miR-183↓, rno-miR-21↓, rno-miR-221↓, rno-miR-2↓5, rno-miR-26b↑, rno-miR-384-3p↓, rno-miR-96↓
Estrogen receptor signalling
rno-miR-182↓, rno-miR-96↓, rno-let-7d↑, rno-miR-132↓, rno-miR-182↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-335↑, rno-miR-384-3p↓
rno-let-7d↑, rno-miR-132↓, rno-miR-182↓, rno-miR-183↓, rno-miR-21↓, rno-miR-221↓, rno-miR-384-3p↓, rno-miR-96
G-Protein coupled receptor signalling
rno-let-7d↑, rno-miR-182↓, rno-miR-183↓, rno-miR-21↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-31↓, rno-miR-335↑, rno-miR-384-3p↓, rno-miR-96↓
Focal adhesion kinase-1 (FAKn)
rno-miR-25↑, rno-miR-384-3p↓, rno-miR-96↓, rno-miR-183↓, rno-miR-26b↑, rno-miR-182↓, rno-miR-21↓, rno-miR-24↓, rno-miR-31↓
Glucocorticoid receptor Signalling
rno-let-7d↑, rno-miR-182↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑
Growth hormone Signalling
rno-let-7d↑, rno-miR-182↓, rno-miR-183↓, rno-miR-21↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑
Inositol phosphate metabolism
rno-miR-335↑, rno-let-7d↑, rno-miR-132↓, rno-miR-182↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-335↑, rno-miR-384-3p↓, rno-miR-96↓
rno-miR-132↓, rno-miR-183↓, rno-miR-21↓, rno-miR-221↓, rno-miR-24↓, rno-miR-25↑, rno-miR-26b↑, rno-miR-31v↓, rno-miR-96↓
PCOS is a complex and heterogeneous endocrine condition characterized by hyperandrogenemia, hyperinsulinemia, insulin resistance, and chronic anovulation. Using a rodent model for PCOS which recapitulates the reproductive and metabolic phenotypes of the human conditions , we have begun to examine molecular and cellular mechanisms underlying the dysregulated reproductive and metabolic functions of PCOS. We have demonstrated that chronically androgenised rats, exhibited increased body weight, disrupted estrus cyclicity, decreased insulin sensitivity and decreased ovarian weight. In the present study we have for the first time that these phenotypes are associated with dysregulation of gene regulatory noncoding miRNA.
Ovarian growth is the result of a series of complex and coordinated processes, which include morphological and functional changes in different follicular cells and their interactions. Among the small non-coding RNAs associated to various biological and physiological disorders (sncRNAs), miRNAs are the best characterized. Studies on the bovine and mouse ovary have demonstated the presence of miRNAs and suggested their potential role in ovarian functions [10, 29]. Several intra-ovarian regulators and their receptors have been predicted as targets of this ovarian miRNAs, which mediate important intracellular events necessary for normal follicular development and functions.
Although miRNAs are expressed in the human, mouse and cow ovaries [10, 14, 29–34], little is known of their role and dysregulation in PCOS. The present study represents the first report on the dysregulation of miRNAs in a rat ovary chronically exposed to DHT, a condition mimicking the hyperandrogenic condition in human PCOS. Among the miRNAs examined, 79 miRNAs (24%) responded to the hyperandrogenic condition and interestingly, 80% of which were upregulated compared to the control group supporting the notion that hyperandrogenic condition down-regulates androgen receptors in the granulosa cells  which could be mediated by these upregulated miRNAs (rno-miR-379*, rno-let-7d, rno-miR-24, rno-miR-673, rno-miR-26b, rno-miR-335, rno-miR-382*, rno-miR-412, rno-miR-99a*, rno-miR-543, rno-miR-674-3p, rno-miR-409-3p). It has been reported that the expression of miR-143, let-7a, miR-15b is under negative control of follicle stimulating hormone (FSH) during follicular development  and may be involved in FSH-induced rat granulosa cell progesterone production . Thus, it is possible that the down-regulation of miRNAs (rno-miR-770, rno-miR-466c, rno-miR-31, rno-miR-183, rno-miR-96, rno-miR-132, rno-miR-182, rno-miR-384-3p and rno-miR-184) observed in this study could be associated with promoted thecal hyperandrogenesis [37, 38]. Similarly, stepwise artificial neural networks analysis revealed predictive miRNA signatures (miR-342, miR-299, miR-217, miR-190, miR-135b, miR-218) corresponding to oestrogen receptor status in breast cancer . These findings raises the possibility that the modulated miRNA expression by exogenous DHT treatment activates estrogen receptors (ERs) to exert indirect steroidal actions on the follicular cells or directly on androgen receptor pathway for the regulation of androgen metabolism .
Recent evidences suggest that miRNAs are involved in the regulation of steroid hormone receptors (SHRs) in ovarian follicle growth and hormone-responsive cancers. Follicular development is the result of tightly regulated activities of FSH, estrogen, androgen and luteinizing hormone receptor (LHR) in granulosa and theca cells, as well as oocyte in response to hormonal stimulation and actions of intra-ovarian regulators . MicroRNAs expression profile and analysis of their target genes revealed that microRNAs exhibits characteristics important in regulating downstream LHR signaling in ovarian cancer cells . MiR-132, 5-fold down-regulation in DHT-treated group compared to control in this study, has previously been shown to be related to luteinizing hormones , raising the possibility that elevated LH in PCOS could involve the expression and action of miR-132.
Androgen Receptors (ARs) mRNA is predominantly expressed in rat granulosa cells and its expression decreases with the stage of follicle development . Hyperandrogenemia induced by DHT is characterized by down-regulation of granulosa cell AR in rat preantral which could be restored by FSH administration . Ovarian follicle expression of miR-221 and miR-222 has been reported to be repressed by androgens [42–44] which regulates cell proliferation by targeting p27/kip1 [44, 45]. Here we observed higher expression of miR-222 in the DHT-induced PCOS rats, a response most evidenced in both theca and granulosa of early stage follicles. Interestingly, its expression pattern changes with increased follicular growth and antrum formation. In the cystic follicle, expression of miR-222 becomes restricted to theca only, suggesting that miR-222 modulates AR expression and hereby paracrine regulation.
Clinical studies have shown that follicular ERα and ERβ expression is significant altered in PCOS . In breast cancer cells, over-expression of miR-221, miR-222, let-7 and miR-20b is associated with reduced of ERα protein content, signaling and expression of ERα target genes [47–49]. Our present study also showed higher expression of these miRNAs in the induced PCOS, consistent with a similar mechanism found in the breast cancer cells and their possible role in the regulation of ERα signaling. In situ localization shows specific expression of miR-222 in the theca cells (Figure 5) where it may be involved in regulation of the ERα. In addition to miR-221/222, several studies also highlighted the differential regulation of let-7d, let-7f, miR-25 and miR-26b in prostate and breast cancer, as well as in leukemia by the estrogen receptor pathways and that their expression was up-regulated in ERα-positive cells [50–52]. Expression analysis of these miRNAs revealed their higher expression in PCOS, raising the possibility that estrogen-mediated regulation of these miRNAs might be involved in the etiology of PCOS. Unlike other miRNAs in DHT-treated PCOS rats, miR-183 was found to be highly expressed in granulosa cells of control ovaries but was down-regulated by 10-fold in the PCOS counterpart. Interestingly, this miRNA was found to be correlated with estrogen receptor in breast cancer .
Hyperandrogenism in PCOS is one of the main causes for early follicle excess, which induces nuclear forkhead box transcription factors 3a (FOXO3a) exclusion and follicular arrest . FOXO3a, a target of miR-96, promotes cancer cell proliferation . Accordingly, expression of miR-96 in PCOS shows a 6-fold down-regulation compared to control ovary and is localized in the theca and cumulus-granulosa of cystic follicles. Thus, the cell proliferation in the preantral follicles in hyperandrogenic condition could be induced through a mechanism utilizing miR-96-mediated regulation.
PCOS is a complex syndrome with both reproductive and metabolic disorders. Although the rat PCOS model (DHT-treatment) simulates the hyperandrogenic state found in human PCOS and recapitulates some of the phenotypes of the human syndrome (obesity, disruption of menstrual cycle, insulin resistance), an inherent problem of any animal model for the complex syndrome is that it does not address all the phenotypes in PCOS. While the rat PCOS model may provide some mechanistic insight in the pathophysiology of the human syndrome, the significant findings from the animal studies need to be validated in the human context. In addition, although miRNAs are known to regulate multiple gene targets in different system [56, 57], they can also be regulated by other paracrine/autocrine factors [44, 58]. The precise target (s) for the reported dysregulated miRNA and their immediate down-stream processes need to be defined.
In summary, using an androgenized rat model which recapitulates many of the reproductive and metabolic phenotypes of human PCOS, we have shown that miRNAs are differentially regulated in a follicular cell-specific manner. Since one given miRNA could target multiple mRNAs of different genes and a given target may be regulated by multiple miRNAs , whether these miRNAs involve in the dysregulation of the hormone receptors (AR, FSHR and ER) and intra-ovarian regulators of ovarian follicle growth and function in PCOS is unknown. The expression of predicted target genes and their regulation by these miRNAs needs to be clarified by further experiments. Our ongoing studies not only include target gene expression analysis but also their correlation with the miRNAs, and functional analysis. Taken together, the present findings support the notion that the present DHT-treated rat model is useful for investigating the role and regulation of miRNAs in the molecular and cellular mechanism of PCOS and may offer new insights for the treatment of reproductive and metabolic disorders associated to PCOS.
Polycystic ovarian syndrome
Equine chorionic gonadotropin
Insulin sensitivity test
Sodium chloride/sodium citrate.
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 [IHDCYH-STIRRHS Postdoctoral Fellowship (MJC), CIHR-QTNPR Doctoral Scholarship (QW) and CIHR-QTNPR Postdoctoral Fellowship (JYK)].
- Dunaif A: Insulin resistance and the polycystic ovary syndrome: mechanism and implications for pathogenesis. Endocr Rev 1997, 18: 774–800. 10.1210/er.18.6.774PubMedGoogle Scholar
- Franks S: Polycystic ovary syndrome. N Engl J Med 1995, 333: 853–861. 10.1056/NEJM199509283331307PubMedView ArticleGoogle Scholar
- Legro RS, Barnhart HX, Schlaff WD, Carr BR, Diamond MP, Carson SA, Steinkampf MP, Coutifaris C, McGovern PG, Cataldo NA, Gosman GG, Nestler JE, Giudice LC, Leppert PC, Myers ER: Clomiphene, metformin, or both for infertility in the polycystic ovary syndrome. N Engl J Med 2007, 356: 551–566. 10.1056/NEJMoa063971PubMedView ArticleGoogle Scholar
- Guzick DS: Polycystic ovary syndrome. Obstet Gynecol 2004, 103: 181–193. 10.1097/01.AOG.0000104485.44999.C6PubMedView ArticleGoogle Scholar
- Diao FY, Xu M, Hu Y, Li J, Xu Z, Lin M, Wang L, Zhou Y, Zhou Z, Liu J, Sha J: The molecular characteristics of polycystic ovary syndrome (PCOS) ovary defined by human ovary cDNA microarray. J Mol Endocrinol 2004, 33: 59–72. 10.1677/jme.0.0330059PubMedView ArticleGoogle Scholar
- Wood JR, Nelson VL, Ho C, Jansen E, Wang CY, Urbanek M, McAllister JM, Mosselman S, Strauss JF III: 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–26390. 10.1074/jbc.M300688200PubMedView ArticleGoogle Scholar
- Hughes C, Elgasim M, Layfield R, Atiomo W: Genomic and post-genomic approaches to polycystic ovary syndrome–progress so far: Mini Review. Hum Reprod 2006, 21: 2766–2775. 10.1093/humrep/del222PubMedView ArticleGoogle Scholar
- Wood JR, Ho CK, Nelson-Degrave VL, McAllister JM, Strauss JF III: The molecular signature of polycystic ovary syndrome (PCOS) theca cells defined by gene expression profiling. J Reprod Immunol 2004, 63: 51–60. 10.1016/j.jri.2004.01.010PubMedView ArticleGoogle Scholar
- Wood JR, Dumesic DA, Abbott DH, Strauss JF III: Molecular abnormalities in oocytes from women with polycystic ovary syndrome revealed by microarray analysis. J Clin Endocrinol Metab 2007, 92: 705–713.PubMedView ArticleGoogle Scholar
- Ro S, Song R, Park C, Zheng H, Sanders KM, Yan W: Cloning and expression profiling of small RNAs expressed in the mouse ovary. RNA 2007, 13: 2366–2380. 10.1261/rna.754207PubMed CentralPubMedView ArticleGoogle Scholar
- Tang F, Kaneda M, O'Carroll D, Hajkova P, Barton SC, Sun YA, Lee C, Tarakhovsky A, Lao K, Surani MA: Maternal microRNAs are essential for mouse zygotic development. Genes Dev 2007, 21: 644–648. 10.1101/gad.418707PubMed CentralPubMedView ArticleGoogle Scholar
- Fiedler SD, Carletti MZ, Hong X, Christenson LK: Hormonal regulation of MicroRNA expression in periovulatory mouse mural granulosa cells. Biol Reprod 2008, 79: 1030–1037. 10.1095/biolreprod.108.069690PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao H, Rajkovic A: MicroRNAs and mammalian ovarian development. Semin Reprod Med 2008, 26: 461–468. 10.1055/s-0028-1096126PubMedView ArticleGoogle Scholar
- Sirotkin AV, Laukova M, Ovcharenko D, Brenaut P, Mlyncek M: Identification of microRNAs controlling human ovarian cell proliferation and apoptosis. J Cell Physiol 2009, 223: 49–56.Google Scholar
- Sirotkin AV, Ovcharenko D, Grossmann R, Laukova M, Mlyncek M: Identification of microRNAs controlling human ovarian cell steroidogenesis via a genome-scale screen. J Cell Physiol 2009, 219: 415–420. 10.1002/jcp.21689PubMedView ArticleGoogle Scholar
- Alford C, Toloubeydokhti T, Al-Katanani Y, Drury KC, Williams R, Chenini N: The expression of microRNA (miRNA) mir-23a and 23b and their target gene, CYP19A1 (aromatase) in follicular cells obtained from women undergoing ART. Fertil Steril 2007, 88: S166-S167.View ArticleGoogle Scholar
- Toloubeydokhti T, Alford C, Al-Katanani Y, Luo X, Williams RS, Nasser C: The expression of microRNA (miRNA), mir-17, mir-211 and mir-542 and their target genes, StAR, IL-1b and Cox2 in follicular cells derived from women undergoing ART. Fertil Steril 2007, 88: S165-S166.View ArticleGoogle Scholar
- Heneghan HM, Miller N, Kerin MJ: Role of microRNAs in obesity and the metabolic syndrome. Obes Rev 2010, 11: 354–361. 10.1111/j.1467-789X.2009.00659.xPubMedView ArticleGoogle Scholar
- Balen AH, Conway GS, Kaltsas G, Techatrasak K, Manning PJ, West C, Jacobs HS: Polycystic ovary syndrome: the spectrum of the disorder in 1741 patients. Hum Reprod 1995, 10: 2107–2111.PubMedView ArticleGoogle Scholar
- Wachs DS, Coffler MS, Malcom PJ, Shimasaki S, Chang RJ: Increased androgen response to follicle-stimulating hormone administration in women with polycystic ovary syndrome. J Clin Endocrinol Metab 2008, 93: 1827–1833. 10.1210/jc.2007-2664PubMed CentralPubMedView ArticleGoogle Scholar
- Motta AB: Dehydroepiandrosterone to induce murine models for the study of polycystic ovary syndrome. J Steroid Biochem Mol Biol 2010, 119: 105–111. 10.1016/j.jsbmb.2010.02.015PubMedView ArticleGoogle Scholar
- Singh KB: Rat models of polycystic ovary syndrome. In Sourcebook of models for biomedical research. Edited by: Conn PM. Totowa, New Jersey: Humana Press Inc; 2008:405–410.View ArticleGoogle Scholar
- Manneras L, Cajander S, Holmang A, Seleskovic Z, Lystig T, Lonn M, Stener-Victorin E: A new rat model exhibiting both ovarian and metabolic characteristics of polycystic ovary syndrome. Endocrinology 2007, 148: 3781–3791. 10.1210/en.2007-0168PubMedView ArticleGoogle Scholar
- Brawer J, Schipper H, Robaire B: Effects of long term androgen and estradiol exposure on the hypothalamus. Endocrinology 1983, 112: 194–199. 10.1210/endo-112-1-194PubMedView ArticleGoogle Scholar
- Marcondes FK, Bianchi FJ, Tanno AP: Determination of the estrous cycle phases of rats: some helpful considerations. Braz J Biol 2002, 62: 609–614. 10.1590/S1519-69842002000400008PubMedView ArticleGoogle Scholar
- Anwer T, Sharma M, Pillai KK, Iqbal M: Effect of Withania somnifera on insulin sensitivity in non-insulin-dependent diabetes mellitus rats. Basic Clin Pharmacol Toxicol 2008, 102: 498–503. 10.1111/j.1742-7843.2008.00223.xPubMedView ArticleGoogle Scholar
- Tilly JL: Ovarian follicle counts–not as simple as 1, 2, 3. Reprod Biol Endocrinol 2003, 1: 11. 10.1186/1477-7827-1-11PubMed CentralPubMedView ArticleGoogle Scholar
- Obernosterer G, Martinez J, Alenius M: Locked nucleic acid-based in situ detection of microRNAs in mouse tissue sections. Nat Protoc 2007, 2: 1508–1514. 10.1038/nprot.2007.153PubMedView ArticleGoogle Scholar
- Hossain MM, Ghanem N, Hoelker M, Rings F, Phatsara C, Tholen E, Schellander K, Tesfaye D: Identification and characterization of miRNAs expressed in the bovine ovary. BMC Genomics 2009, 10: 443. 10.1186/1471-2164-10-443PubMed CentralPubMedView ArticleGoogle Scholar
- Dahiya N, Sherman-Baust CA, Wang TL, Davidson B, Shih Ie M, Zhang Y, Wood W III, Becker KG, Morin PJ: MicroRNA expression and identification of putative miRNA targets in ovarian cancer. PLoS One 2008, 3: e2436. 10.1371/journal.pone.0002436PubMed CentralPubMedView ArticleGoogle Scholar
- Giannakakis A, Sandaltzopoulos R, Greshock J, Liang S, Huang J, Hasegawa K, Li C, O'Brien-Jenkins A, Katsaros D, Weber BL, Simon C, Coukos G, Zhang L: miR-210 links hypoxia with cell cycle regulation and is deleted in human epithelial ovarian cancer. Cancer Biol Ther 2008, 7: 255–264. 10.4161/cbt.7.2.5297PubMed CentralPubMedView ArticleGoogle Scholar
- Laios A, O'Toole S, Flavin R, Martin C, Kelly L, Ring M, Finn SP, Barrett C, Loda M, Gleeson N, D'Arcy T, McGuinness E, Sheils O, Sheppard B, O’Leary J: Potential role of miR-9 and miR-223 in recurrent ovarian cancer. Mol Cancer 2008, 7: 35. 10.1186/1476-4598-7-35PubMed CentralPubMedView ArticleGoogle Scholar
- Nam EJ, Yoon H, Kim SW, Kim H, Kim YT, Kim JH, Kim JW, Kim S: MicroRNA expression profiles in serous ovarian carcinoma. Clin Cancer Res 2008, 14: 2690–2695. 10.1158/1078-0432.CCR-07-1731PubMedView ArticleGoogle Scholar
- Zhang L, Volinia S, Bonome T, Calin GA, Greshock J, Yang N, Liu CG, Giannakakis A, Alexiou P, Hasegawa K, Johnstone CN, Megraw MS, Adams S, Lassus H, Huang J, Kaur S, Liang S, Sethupathy P, Leminen A, Simossis VA, Sandaltzopoulos R, Naomoto Y, Katsaros D, Gimotty PA, DeMichele A, Huang Q, Butzow R, Rustgi AK, Weber BL, Birrer MJ, Hatzigeorgiou AG, Croce CM, Coukos G: Genomic and epigenetic alterations deregulate microRNA expression in human epithelial ovarian cancer. Proc Natl Acad Sci U S A 2008, 105: 7004–7009. 10.1073/pnas.0801615105PubMed CentralPubMedView ArticleGoogle Scholar
- Tetsuka M, Hillier SG: Androgen receptor gene expression in rat granulosa cells: the role of follicle-stimulating hormone and steroid hormones. Endocrinology 1996, 137: 4392–4397. 10.1210/en.137.10.4392PubMedGoogle Scholar
- Yao N, Lu CL, Zhao JJ, Xia HF, Sun DG, Shi XQ, Wang C, Li D, Cui Y, Ma X: A network of miRNAs expressed in the ovary are regulated by FSH. Front Biosci 2009, 14: 3239–3245.View ArticleGoogle Scholar
- Yao N, Yang BQ, Liu Y, Tan XY, Lu CL, Yuan XH, Ma X: Follicle-stimulating hormone regulation of microRNA expression on progesterone production in cultured rat granulosa cells. Endocrine 2010, 38: 158–166. 10.1007/s12020-010-9345-1PubMedView ArticleGoogle Scholar
- Cardenas H, Herrick JR, Pope WF: Increased ovulation rate in gilts treated with dihydrotestosterone. Reproduction 2002, 123: 527–533. 10.1530/rep.0.1230527PubMedView ArticleGoogle Scholar
- Lowery AJ, Miller N, Devaney A, McNeill RE, Davoren PA, Lemetre C, Benes V, Schmidt S, Blake J, Ball G, Kerin MJ: MicroRNA signatures predict oestrogen receptor, progesterone receptor and HER2/neu receptor status in breast cancer. Breast Cancer Res 2009, 11: R27. 10.1186/bcr2257PubMed CentralPubMedView ArticleGoogle Scholar
- Mercier L, le Guellec C, Thieulant ML, Samperez S, Jouan P: Androgen and estrogen receptors in the cytosol from male rat anterior hypophysis: further characteristics and differentiation between androgen and estrogen receptors. J Steroid Biochem 1976, 7: 779–785. 10.1016/0022-4731(76)90179-5PubMedView ArticleGoogle Scholar
- Cui J, Eldredge JB, Xu Y, Puett D: MicroRNA expression and regulation in human ovarian carcinoma cells by luteinizing hormone. PLoS One 2011, 6: e21730. 10.1371/journal.pone.0021730PubMed CentralPubMedView ArticleGoogle Scholar
- Chen Y, Robles AI, Martinez LA, Liu F, Gimenez-Conti IB, Conti CJ: Expression of G1 cyclins, cyclin-dependent kinases, and cyclin-dependent kinase inhibitors in androgen-induced prostate proliferation in castrated rats. Cell Growth Differ 1996, 7: 1571–1578.PubMedGoogle Scholar
- Mercatelli N, Coppola V, Bonci D, Miele F, Costantini A, Guadagnoli M, Bonanno E, Muto G, Frajese GV, De Maria R, Spagnoli LG, Farace MG, Ciafre SA: The inhibition of the highly expressed miR-221 and miR-222 impairs the growth of prostate carcinoma xenografts in mice. PLoS One 2008, 3: e4029. 10.1371/journal.pone.0004029PubMed CentralPubMedView ArticleGoogle Scholar
- Sun T, Wang Q, Balk S, Brown M, Lee GS, Kantoff P: The role of microRNA-221 and microRNA-222 in androgen-independent prostate cancer cell lines. Cancer Res 2009, 69: 3356–3363. 10.1158/0008-5472.CAN-08-4112PubMed CentralPubMedView ArticleGoogle Scholar
- Galardi S, Mercatelli N, Giorda E, Massalini S, Frajese GV, Ciafre SA, Farace MG: miR-221 and miR-222 expression affects the proliferation potential of human prostate carcinoma cell lines by targeting p27Kip1. J Biol Chem 2007, 282: 23716–23724. 10.1074/jbc.M701805200PubMedView ArticleGoogle Scholar
- Jakimiuk AJ, Weitsman SR, Yen HW, Bogusiewicz M, Magoffin DA: Estrogen receptor alpha and beta expression in theca and granulosa cells from women with polycystic ovary syndrome. J Clin Endocrinol Metab 2002, 87: 5532–5538. 10.1210/jc.2002-020323PubMedView ArticleGoogle Scholar
- Castellano L, Giamas G, Jacob J, Coombes RC, Lucchesi W, Thiruchelvam P, Barton G, Jiao LR, Wait R, Waxman J, Hannon GJ, Stebbing J: The estrogen receptor-alpha-induced microRNA signature regulates itself and its transcriptional response. Proc Natl Acad Sci U S A 2009, 106: 15732–15737. 10.1073/pnas.0906947106PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao JJ, Lin J, Yang H, Kong W, He L, Ma X, Coppola D, Cheng JQ: MicroRNA-221/222 negatively regulates estrogen receptor alpha and is associated with tamoxifen resistance in breast cancer. J Biol Chem 2008, 283: 31079–31086. 10.1074/jbc.M806041200PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao Y, Deng C, Wang J, Xiao J, Gatalica Z, Recker RR, Xiao GG: Let-7 family miRNAs regulate estrogen receptor alpha signaling in estrogen receptor positive breast cancer. Breast Cancer Res Treat 2010, 127: 69–80.PubMedView ArticleGoogle Scholar
- Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin SF, Dunning MJ, Barbosa-Morais NL, Teschendorff AE, Green AR, Ellis IO, Tavare S, Caldas C, Miska EA: MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biol 2007, 8: R214. 10.1186/gb-2007-8-10-r214PubMed CentralPubMedView ArticleGoogle Scholar
- Kovalchuk O, Tryndyak VP, Montgomery B, Boyko A, Kutanzi K, Zemp F, Warbritton AR, Latendresse JR, Kovalchuk I, Beland FA, Pogribny IP: Estrogen-induced rat breast carcinogenesis is characterized by alterations in DNA methylation, histone modifications and aberrant microRNA expression. Cell Cycle 2007, 6: 2010–2018. 10.4161/cc.6.16.4549PubMedView ArticleGoogle Scholar
- Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, Fedele V, Ginzinger D, Getts R, Haqq C: Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer 2006, 5: 24. 10.1186/1476-4598-5-24PubMed CentralPubMedView ArticleGoogle Scholar
- Lowery AJ, Miller N, Dwyer RM, Kerin MJ: Dysregulated miR-183 inhibits migration in breast cancer cells. BMC Cancer 2010, 10: 502. 10.1186/1471-2407-10-502PubMed CentralPubMedView ArticleGoogle Scholar
- Yang JL, Zhang CP, Li L, Huang L, Ji SY, Lu CL, Fan CH, Cai H, Ren Y, Hu ZY, Gao F, Liu YX: Testosterone induces redistribution of forkhead box-3a and down-regulation of growth and differentiation factor 9 messenger ribonucleic acid expression at early stage of mouse folliculogenesis. Endocrinology 2010, 151: 774–782. 10.1210/en.2009-0751PubMedView ArticleGoogle Scholar
- Lin H, Dai T, Xiong H, Zhao X, Chen X, Yu C, Li J, Wang X, Song L: Unregulated miR-96 induces cell proliferation in human breast cancer by downregulating transcriptional factor FOXO3a. PLoS One 2010, 5: e15797. 10.1371/journal.pone.0015797PubMed CentralPubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004, 116: 281–297. 10.1016/S0092-8674(04)00045-5PubMedView ArticleGoogle Scholar
- Toloubeydokhti T, Bukulmez O, Chegini N: Potential regulatory functions of microRNAs in the ovary. Semin Reprod Med 2008, 26: 469–478. 10.1055/s-0028-1096127PubMed CentralPubMedView ArticleGoogle Scholar
- Sun Q, Zhang Y, Yang G, Chen X, Cao G, Wang J, Sun Y, Zhang P, Fan M, Shao N, Yang X: Transforming growth factor-beta-regulated miR-24 promotes skeletal muscle differentiation. Nucleic Acids Res 2008, 36: 2690–2699. 10.1093/nar/gkn032PubMed CentralPubMedView ArticleGoogle Scholar
- Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet 2005, 37: 495–500. 10.1038/ng1536PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.