MicroRNA-376a Regulates 78-Kilodalton Glucose-Regulated Protein Expression in Rat Granulosa Cells

The 78-kilodalton glucose-regulated protein (GRP78) is a molecular chaperone that assists in protein assembly, folding, and translocation. Recently, our laboratory reported that GRP78 regulates the expression of luteinizing hormone-human chorionic gonadotropin receptor (LHR) in the early stage of corpus luteum formation. In this study, we investigated whether microRNAs (miRNAs), which post-transcriptionally regulate mRNA, are involved in the regulation mechanism of GRP78 in the ovary. A miRNA microarray was performed to analyze the overall miRNA expression profile, and the results indicated that 44 miRNAs were expressed highly after ovulation was induced. The results from a bio-informative database analysis and in vitro granulosa cell culture studies led us to focus on rno-miR-376a for further analysis. In both in vivo and in vitro studies, rno-miR-376a levels increased 12 h after human chorionic gonadotropin (hCG) administration. To elucidate whether rno-miR-376a induced mRNA destabilization or translational repression of GRP78, rno-miR-376a was transfected into cultured granulosa cells, resulting in decreased GPR78 protein levels without an alteration in GRP78 mRNA levels. To confirm that rno-miR-376a binds to GRP78 mRNA, we cloned the 3′-end of GRP78 mRNA (nucleotides 2439–2459) into a reporter vector that contained a Renilla luciferase coding region upstream of the cloning site. The luciferase assays revealed that rno-miR-376a bound to the 3′-end of GRP78 mRNA. From these data, we conclude that rno-miR-376a potentially negatively regulates GRP78 protein expression through translational repression at an early stage transition from the follicular phase to luteinization.


Introduction
In females, ovulation-the release of a mature and fertilizable oocyte-is an essential process in the establishment of pregnancy. Follicle stimulating hormone (FSH) is a key factor that induces luteinizing hormone-human chorionic gonadotropin receptor (LHR) expression in the granulosa cells of pre-ovulatory follicles, whereby an LH surge triggers ovulation followed by corpus luteum formation [1,2,3]. Ovulation provokes a dramatic transformation of ovulated follicle to form the corpus luteum, which in turn induces synthesis of progesterone in order to sustain pregnancy. This indicates that granulosa cells are under stress, which has led to ovulation being considered an inflammation-like phenomenon [4,5].
When a cell is exposed to a variety of environmental and physiological stresses, the function of the endoplasmic reticulum (ER) is disturbed, and unfolded proteins accumulated in the ER [6]. In the ER, molecular chaperones are responsible for proper protein folding and for preventing aggregation of unfoldedintermediates, which could lead to cell death [7]. Recently, our laboratory determined that the 78-kilodalton glucose-regulated protein (GRP78), an ER-associated molecular chaperone that assists in proper protein folding to execute primary protein maturation in the ER [8,9], is involved in the recovery of LHR after its down-regulation [10]. Although GRP78 is known be an important molecule for the upregulation of LHR, the precise mechanism underlying the regulation of GRP78 expression in the ovary has not been fully elucidated.
MicroRNAs (miRNAs) are non-coding RNAs (approximately 22 nucleotides) that regulate gene expression by binding to the 39untranslated regions of target mRNAs to induce the degradation of target mRNAs or to induce translational repression [11]. Within the past decade, miRNAs have been recognized as important regulators in many biological and cellular processes, such as cell proliferation, differentiation, apoptosis, and tumorigenesis [12,13,14]. Although the functions of miRNAs have not been elucidated fully, recent emerging evidence demonstrates that miRNAs are involved in ovarian follicular and luteal functions [15,16,17,18]. We have also reported that miR-136-3p targets LHR mRNA to induce the transient downregulation of LHR mRNA in granulosa cells after ovulation [19]. This result prompted us to search for miRNAs involved in the regulation of GRP78 expression.
In the following experiments, we profiled the miRNAs that were expressed in PMSG-primed rat ovaries in which an injection of hCG induced ovulation. From that data and data obtained from a bio-informatics database, we focused on rno-miR-376a, which potentially binds to GRP78 mRNA, and characterized the function of rno-miR-376a in cultured granulosa cells.

Ethic Statement
All procedures involving animals followed ethical principles according the NIH Guide for Care and Use of the Laboratory Animals and were approved by the Gunma University Animal Care and Use Committee (Permit Number: 12-011).

Animals
Female 21-day-old Wistar rats (Japan SLC, Inc., Hamamatsu, Japan) were housed in a temperature-and light-controlled room (12 h light, 12 h dark cycle; lights on at 0600 h) with food and water provided ad libitum.
For the in vivo study, 21-day-old female rats received an injection of PMSG (30 IU, subcutaneous [sc]), followed by an injection of hCG (20 IU, sc) 48 h later. We then sacrificed the rats by decapitation to remove the ovaries at each time point after the injection of hCG. The ovaries were stored immediately in RNAlater Tissue Collection solution (Applied Biosystems, Foster City, CA, USA).
For the in vitro study, 21-day-old female rats received DES injections (2 mg, sc) for 4 d, and the ovaries were removed for use in primary granulosa cell culture, as described below.
Primary rat granulosa cell culture Granulosa cells were obtained from DES-primed rats. Both ovaries were removed from the rats, and the granulosa cells were released by puncturing the follicles with a 26-gauge needle. The granulosa cells were washed and collected by brief centrifugation, and cell viability was determined by trypan blue exclusion. The granulosa cells were then cultured in DMEM/Ham's nutrient mixture F-12 supplemented with 20 mg/L gentamicin sulfate, 500 mg/L Fungizone, and 1 g/L bovine serum albumin (BSA) on collagen-coated plates or dishes in a humidified atmosphere containing 5% CO2 and 95% air at 37uC.

miRNA microarray
Rat ovaries were removed from PMSG/hCG-treated rats to analyze miRNA expression on a microarray. The miRNA microarray expression profiling was performed by TaKaRa Bio Inc. using the Agilent Technologies Rat miRNA Microarray 8615K (Santa Clara, CA, USA), which contains 350 rat miRNAs listed in the Sanger miRBase v.10.1. Total RNA, including miRNAs, was isolated from the PMSG/hCG-treated rat ovaries using the FastPure RNA Kit (TaKaRa Bio Inc., Otsu, Japan) according to the manufacturer's protocol. A quality check of the total RNA was performed using an Agilent 2100 Bioanalyzer (Agilent Technologies). The labeling and hybridization of the total RNA samples were performed using an Agilent miRNA Microarray System (Agilent Technologies). The microarray results were extracted using Agilent Feature Extraction software, and the data were analyzed by TaKaRa Bio Inc.

miRNA target prediction
Target miRNAs that bound to rat GRP78 mRNA were predicted using MicroCosm Targets (http://www.ebi.ac.uk/ enright-srv/microcosm/htdocs/targets/v5/; hosted by the European Bioinformatics Institute) based on the miRNA microarray results.
miRNA expression profiling of rat ovaries and rat granulosa cells using TaqMan miRNA assays Total RNA, including miRNAs, was isolated from rat ovaries or rat granulosa cells using the mirVana miRNA Isolation Kit (Ambion Inc., Austin, TX, USA) according to the manufacturer's protocol. The concentration of the total RNA was measured using a Nanodrop-1000 spectrophotometer (Thermo scientific, Wilmington, DE, USA). The total RNA was then reverse transcribed using specific RT primers from a TaqMan MicroRNA Assay Kit (rno-miR-144: 197375, rno-miR-376a: 001069, rno-miR-451: 001141, and 4.5S RNA(H): 001716 as an internal control; Ambion Inc.). Single-stranded cDNA was synthesized from 10 ng of total RNA in a 15-mL reaction volume using the MicroRNA Reverse Transcription Kit (Applied Biosystems), according to the manufacturer's protocol. The reactions were incubated for 30 min at 16uC, 30 min at 30uC, and 5 min at 85uC in a thermal cycler. Real-time RT-PCR was performed using sequence-specific primers from the TaqMan MicroRNA Assay Kit according to the manufacturer's instructions. The reactions were performed using an ABI PRISM 7000 sequence detection system (Applied Biosystems) in a 20-mL reaction volume at 95uC for 10 min, followed by 40 cycles of 95uC for 15 s and 60uC for 60 s. All of the reactions were tested in triplicate. We determined the threshold cycle (Ct) values for each reaction, and their means were used to determine the relative miRNA expression levels using the DDCt method. Additionally, we calculated the percentage of rat GRP78 mRNA or miRNA expression using the DDCt values for each reaction as follows: These values were used to calculate, and the relative expression of GRP78 or each miRNA was normalized to the expression of hCG 0 h in the control, which was set at 1.

Expression of GRP78 mRNA expression in rat ovaries or rat granulosa cells using TaqMan Gene Expression Assays
Total RNA was isolated from rat ovaries or rat granulosa cells, and the concentrations were measured as described above. The isolated RNAs (2 mg of each sample) were treated with DNaseI (Invitrogen) to eliminate residual genomic DNA. These RNAs were reverse transcribed using random primers, 10 mM deoxynucleoside triphosphate mix, and SuperScript III reverse transcriptase (Invitrogen) according to the manufacturer's protocol. The reactions were incubated for 5 min at 25uC, 60 min at 50uC, and 15 min at 70uC in a thermal cycler. To remove complementary RNA, RNase H was added to the cDNAs, and the reactions were incubated for 20 min at 37uC. Real-time RT-PCR was performed using TaqMan Gene Expression Assays (GRP78-Hspa5: Rn00565250_m1, Eukaryotic 18S rRNA: Hs99999901_s1 as an internal control; Applied Biosystems) according to the manufacturer's instructions. Rat GRP78 mRNA expression levels were measured in a similar manner as in the miRNA assays.
Transfection with Pre-miR-376a miRNA (precursor) and Anti-miR-376a miRNA (inhibitor) Primary rat granulosa cells (8.0-10610 4 cells obtained from DES-primed rats) were cultured in 24-well tissue culture plates with serum-free medium. Twenty-four hours after seeding, the cells were incubated with FSH (30 ng/mL) and estradiol (10 nM) for 48 h. The cells were then transfected with Pre-miR-376a (precursor) or Anti-miR-376a (inhibitor) (50 nM each) purchased from Ambion (product IDs: PM10504 and AM10504, respectively) using siPORT NeoFX Transfection Agent according to the manufacturer's protocol. Twelve hours after transfection, hCG (30 ng/mL) was added to the culture medium to induce the downregulation or upregulation of GRP78 mRNA.
To evaluate the effects of precursor and inhibitor on GRP78 mRNA expression, we performed real-time RT-PCR using TaqMan gene expression assays 12 h after hCG treatment and compared the results with the mock-transfected cells. All of the reactions were tested in triplicate. Rat GRP78 mRNA expression was measured in a similar manner as in the miRNA assays.

Western blot analysis
Primary rat granulosa cells (150610 4 cells obtained from DESprimed rats) were cultured in 3.5-cm tissue culture dishes with serum-free medium. Twenty-four hours after seeding, the cells were incubated with FSH (30 ng/mL) and estradiol (10 nM) for 48 h. The cells were then transfected with Pre-miR Negative control #1 (product IDs: AM17110), precursor or inhibitor (50 nM each) purchased from Ambion, using Lipofectamine 2000 Transfection Reagent according to the manufacturer's protocol. Twelve hours after transfection, hCG (30 ng/mL) was added to the culture medium.

Reporter vectors and DNA constructs
To identify the rno-miR-376a-binding site at the 39-end of GRP78 mRNA, we generated a direct-match miRNA target site and cloned the insert into the multiple cloning site in the luciferase reporter vector from the pMIR-REPORT miRNA Expression Reporter Vector System (Applied Biosystems). The sense and antisense strands of the oligonucleotides were annealed by adding 2 mg of each oligonucleotide to 46 mL of annealing solution (100 mM potassium acetate, 30 mM HEPES-KOH, pH 7.4 and 2 mM magnesium acetate) and incubating at 90uC for 5 min and then at 37uC for 1 h. The annealed oligonucleotides were digested with HindIII and SpeI and ligated into the HindIII and SpeI restriction sites of pMIR-REPORT vector. The sequences of the inserts were confirmed by sequence analysis using a PRISM 3100 genetic analyzer (Applied Biosystems).

Luciferase assay
We transfected HEK293 cells with 200 ng of each vector (pMIR-REPORT luciferase vectors, as described above, and the pMIR-REPORT bgal vector as a control for transfection normalization) and 50 nM precursor or inhibitor using Lipofectamine 2000 Transfection Reagent according to the manufacturer's protocol. To measure the luciferase activity, we harvested the cells

Data analysis
The microarray data were first filtered by subtracting the control probe data. To identify miRNAs, a one-way ANOVA was performed for differentially expressed miRNAs between four groups (hCG 0 h, hCG 6 h, hCG 12 h, hCG 24 h). False discovery rates (FDR) were assessed using the Benjamini and Hochberg method. Hierarchical clustering was then performed using the complete linkage method.
Comparisons between groups were performed using one-way ANOVA followed by Dunnett's multiple comparison test. The data represent the mean 6SE from at least three independent experiments. A value of P,0.05 was considered significant.

miRNA microarray and target prediction
We utilized a miRNA microarray analysis to identify miRNAs that were expressed in rat ovaries. Based on the differentially expressed miRNAs, a cluster analysis was executed to generate subgroups characterized by the expression patterns in four groups (hCG 0 h, hCG 6 h, hCG 12 h, and hCG 24 h) (Fig. S1). A total of 44 miRNAs were found to differ significantly in response to hCG treatment, as measured using the Benjamini and Hochberg method, suggesting that these miRNAs potentially affected the expression of GRP78. From these 44 miRNAs, we focused on 23 miRNAs that increased 6 h after hCG treatment and then decreased by 24 h (Table 1), and 19 miRNAs that decreased 6 h after hCG treatment and then increased by 24 h ( Table 2). The results of the bioinformatics database inquiry (MicroCosm, http:// ebi.ac.uk/) revealed that rno-miR-144, rno-miR-376a, and rno-miR-451 can bind to the 39-UTR of GRP78 mRNA (from bp 2439-2459) and negatively regulate GRP78 expression. Figure 2. Rat GRP78 mRNA, rno-miR-144, rno-miR-376a, and rno-miR-451 expression in primary rat granulosa cells induced by FSH and hCG. Primary rat granulosa cells were prepared, and the indicated reagents were added to the medium after 24 h of culture. Cells were then incubated with FSH (30 ng/mL) and estradiol (10 nM) for 48 h. Subsequently, hCG (30 ng/mL) was added to the culture medium, as described in the Materials and Methods. Total RNA was isolated, and GRP78 mRNA (A), rno-miR-144 (B), rno-miR-376a (C), and rno-miR-451 (D) expression levels were measured using real-time RT-PCR as described in the Materials and Methods. The amounts of GRP78 mRNA, rno-miR-144, rno-miR-376a, and rno-miR-451 in the hCG 0 h group were set at 1. Data were normalized for 18S rRNA (for GRP78 mRNA) and 4.5S RNA(H) (for rno-miR-144, rno-miR-376a, and rno-miR-451) levels in each sample and represent the mean 6SE of 3 independent experiments. *, significantly different from the control value at hCG 0 h, P,0.05. doi:10.1371/journal.pone.0108997.g002 MicroRNA-376a Regulates GRP78 Protein Expression PLOS ONE | www.plosone.org miRNA expression in rat ovaries In both in vivo and in vitro experiments, we previously demonstrated that GRP78 mRNA levels peaked while LHR mRNA was downregulated by hCG and that the increase of LHR protein levels was dependent on the increment of GRP78 protein [10]. Therefore, we measured the expression of GRP78 mRNA, rno-miR-144, rno-miR-376a, and rno-miR-451 using real-time RT-PCR to examine their expression patterns in the context of LHR downregulation (Fig. 1). Consistent with previous studies, GRP78 mRNA upregulation peaked 12 h after the administration of an ovulatory dose of hCG and subsequently decreased. In a similar manner, rno-miR-144 and rno-miR-376a expression peaked 12 h after the hCG treatment, and rno-miR-451 expression peaked 24 h after the hCG treatment.

miRNA expression in rat granulosa cells
Next, we investigated the effects of rno-miR-144, rno-miR-376a, and rno-miR-451 on GRP78 mRNA expression in granulosa cells isolated from DES-treated immature rats (Fig. 2). GRP78 mRNA and rno-miR-376a expression increased significantly 12 h after the addition of hCG into the culture medium and subsequently decreased, which was consistent with in vivo study (Fig. 1). In contrast, rno-miR-144 and rno-miR-451 expression  . Effects of Pre-miR-376a (precursor) and Anti-miR-376a (inhibitor) transfection on rat GRP78 mRNA expression in primary rat granulosa cells. Primary rat granulosa cells were prepared, and the indicated reagents were added to the medium after 24 h of culture. The cells were incubated with FSH (30 ng/mL) and estradiol (10 nM) for 48 h. Pre-miR-376a (precursor) or Anti-miR-376a (inhibitor) was transfected into the cells, and 30 ng/mL hCG was added 12 h later. The effects of precursor and inhibitor on the expression of GRP78 mRNA were measured using real-time RT-PCR, as described in the Materials and Methods. The expression of GRP78 mRNA at hCG 0 h in the control was set at 1. Each value represents the mean 6SE of three independent experiments. doi:10.1371/journal.pone.0108997.g004 decreased significantly after the hCG treatment. Furthermore, miR-376a did not show a substantial change in the absence of hCG treatment (Fig. 3), confirming that hCG has an important role for the induction of miR-376a.

miR-376a does not affect GRP78 mRNA expression
We then examined the effects of rno-miR-376a on GRP78 mRNA levels following the transfection of granulosa cells with either precursor or inhibitor (Fig. 4). In this experiment, the granulosa cells were cultured in a similar manner to the cells utilized in Fig. 2 except that the cells were transfected with precursor or inhibitor 12 h prior to hCG treatment. The results indicated that the transfection with precursor or inhibitor did not affect GRP78 mRNA expression.

miR-376a represses GRP78 translation
Furthermore, to ascertain the expression of GRP78 protein, western blot analysis (Fig. 5) was performed using the same animal model designed for the real-time RT-PCR analysis. The results indicated that the transfection of the granulosa cells with precursor Figure 5. Effects of rno-miR-376a on GRP78 protein in granulosa cells. Primary rat granulosa cells were prepared, and the indicated reagents were added to the medium after 24 h of culture. The cells were incubated with FSH (30 ng/mL) and estradiol (10 nM) for 48 h. Negative control (A), precursor (B) or inhibitor (C) was transfected into the cells, and 30 ng/mL hCG was added 12 h later. Cells were harvested 24 h after the addition of hCG, and GRP78 protein levels were quantified using western blot analysis. (D) Levels of GRP78 protein were quantified by densitometric scanning. The expression of GRP78 protein in the control (hCG 0 h) was set at 1. Values represent the mean 6SE of three independent experiments. *, significantly different from the control value, P,0.05. doi:10.1371/journal.pone.0108997.g005 significantly decreased GRP78 protein (Fig. 5B). However, the transfection of granulosa cells with inhibitor upregulated GRP78 protein expression 24 h after the hCG treatment (Fig. 5C).

miR-376a binds to predicted site of GRP78 39-UTR
To confirm the presence of an rno-miR-376a binding site on GRP78 mRNA, we constructed a reporter vector that contained a putative rno-miR-376a binding site sequence (bp 2439-2459) in the 39-UTR downstream of a Renilla luciferase coding region. To avoid the influence of endogenous rno-miR-376a, we used HEK293 cells rather than rat granulosa cells (Fig. 6). To verify the effect of the miRNA transfection, we generated a reporter vector that included the whole miR-376a sequence. The results indicated that precursor reduced the luciferase activity to approximately 0%, while inhibitor did not affect the luciferase activity. These data demonstrate that the translation of GRP78 was blocked by precursor via a binding site that was the complementary sequence to rno-miR-376a.
The luciferase activity in the cells transfected with the reporter vector containing the putative rno-miR-376a binding site was reduced by approximately 50% 24 h after the cells were transfected with precursor. In contrast, the transfection with inhibitor did not affect the luciferase activity. This result indicates that rno-miR-376a binds to the 39-end of GRP78 mRNA from bp 2439-2459 in the 39-UTR. To identify the rno-miR-376a-binding site in the 39-UTR of GRP78 mRNA, luciferase reporter vectors were generated as described in the Materials and Methods. (C) Luciferase activity was measured to identify the rno-miR-376a-binding site in the 39-UTR of GRP78 mRNA. HEK293 cells were prepared, and the cells were transfected with 200 ng of each reporter vector with 50 nM Pre-miR-376a (precursor) or Anti-miR-376a (inhibitor) as described in the Materials and Methods. For transfection normalization, the cells were also transfected with the pMIR-REPORT bgal vector. Luciferase activity was measured 24 h after transfection. The activity of the control (empty vector) was set at 1. Each value represents the mean 6SE of three independent experiments. *, significantly different from the control value, P,0.05. doi:10.1371/journal.pone.0108997.g006

Discussion
The transition from the ovulatory follicle to corpus luteum formation is a critical phase for maintaining female reproductive functions. Within a follicle after ovulation, many different cell types are subjected to dramatic changes, including proliferation, apoptosis, and differentiation. Thus, the ovulatory stimuli provided by the LH surge trigger the synthesis of numerous proteins inside cells, indicating that ER-associated chaperones facilitate proper maturation of newly synthesized proteins. In fact, our recent finding suggests that GRP78 is involved in restoring LHR expression following the down-regulation induced by ovulation [10]. Therefore, to better understand the regulation of GRP78 expression, we focused on miRNAs that are induced in the ovary after ovulation because increasing evidence suggests that miRNAs play pivotal roles in the regulation of a wide range of biological processes [18,20,21].
To identify miRNAs that potentially regulate GRP78 expression, microRNA microarrays were performed using ovaries from PMSG-primed rats injected with an ovulatory dose of hCG. The array data along with the bioinformatic analysis provided by MicroCosm Targets, which indicated that several miRNAs bind to the GRP78 mRNA 39-UTR, led us to focus on rno-miR-144, rno-miR-376a, and rno-miR-451. From these, we narrowed the focus to rno-miR-376a based on the results of the in vitro experiments (Fig. 2) since rno-miR-144 and rno-miR-451 was not induced expression by hCG treatment. We believe that the discrepancy regarding in the expression of miR-144 and 451 between in the in vivo and in vitro experiments ( Fig. 1 and 2) can be explained by the presence of blood cells in the in vivo study which express both rno-miR-144 and rno-miR-451 [22,23].
Initially, miRNAs were thought to bind within the 39-UTR of their target mRNAs, but recent evidence has demonstrated that the 59-UTR [24] or coding region [25] can also be the target sites for miRNAs. As more than 60% of human protein coding genes are predicted to be under the control of miRNAs by binding within their 39-UTR [26,27], GRP78 mRNA in the rat ovary is regulated post-transcriptionally by miRNAs in a same manner as the majority of human mRNAs.
Predicting miRNA targets using computational analysis, which was employed for this study, is efficient but can result in tens or hundreds of targets with high false-positive rates [28]. To avoid this issue, we executed a luciferase assay to identify an rno-miR-376a-binding site in the 39-end of GRP78 mRNA, which confirmed that the 39-end of GRP78 mRNA from bp 2439 to 2459 is the binding site for rno-miR-376a. As mentioned above, miRNAs mainly control their target genes through base pairing to the 39-UTRs of target mRNAs, resulting in the repression of translation of target proteins or the degradation of target mRNAs. In initial reports, the extent of complementarity between miRNA and its target mRNA is thought to govern either translational mRNA cleavage or repression [29]: in plants, nearly perfect complementarity results in the degradation of target miRNAs [30], whereas in animals, partial complementarity results in a translational block [31]. However, recent studies have reported that miRNAs can also induce mRNA degradation in animals and, conversely, translational repression in plants, although the topic remains controversial due to conflicting data [32,33]. Moreover, the state of the complementarity between target mRNA and miRNA has been suggested to affect mRNA degradation [34], which is applicable to rno-miR-376a because the introduction of the complete identical sequence of rno-miR-376a into a reporter vector abolished the luciferase activity, whereas the luciferase activity was reduced to 50% using the reporter vector containing the 39-UTR rno-miR-376a binding site of GRP78 mRNA (Fig. 6).
The results presented in this study ( Fig. 4 and Fig. 5) demonstrate that rno-miR-376a decreases GRP78 protein production by translational repression without altering GRP78 mRNA levels and that the transfection of rno-miR-376a into granulosa cells repressed protein expression by approximately three-fold (Fig. 5). One can argue whether this reduction of GRP78 protein by rno-miR-376a has significant meaning for physiological functions in the ovary. In our previous report [10], three-fold reduction in the level of GRP78 protein significantly impaired restoration of LHR expression in the granulosa cells after the LHR down-regulation induced by ovulation. In addition, previous reports have concluded that the magnitude of repression by miRNA rarely exceeds two-fold [26,35]. Thus, we believe that rno-miR-376a is involved in the regulation of GRP78 expression in the ovary.
It has been accepted that a single miRNA can downregulate hundreds of mRNAs and the production of hundreds of proteins [26,36]. Combined with the MicroCosm analysis of miRNA targets, which revealed that 44 miRNAs can bind to the 39-UTR of GRP78 mRNA, we speculate that multiple miRNAs may constitute a network that is involved in the regulation of GRP78 mRNA, whereas rno-miR-376a is solely identified to be evoked by the ovulatory signal. In the regulation of gene expression, many intermediate steps occur during the processing of RNA to protein.
GRP78, which is regarded as a stress-inducible house-keeping gene, is transcribed constitutively. Therefore, the post-transcriptional control of GRP78 by rno-miR-376a has a significant impact on the regulation of GRP78 protein expression and not gene expression, which is evident by the lack of GRP78 transcriptional regulation.
In conclusion, the present results demonstrate that the induction of rno-miR-376a production by hCG leads to repression of GRP78 translation. Although the precise mechanism regulating GRP78 has not been elucidated, our finding may provide clues for those investigating miRNA-mediated regulation of ovarian functions. Figure S1 Hierarchical clustering analysis of miRNA expression in rat ovaries induced PMSG and hCG. The clustering analysis was conducted by TaKaRa Bio Inc. and was carried out using the complete linkage method. Forty-four differentially expressed miRNAs chosen with an adjusted false discovery rate ,0.05. Each row represents an individual miRNA and each column represents the time point after the injection of hCG. The miRNA clustering tree is shown on the left. The color scale illustrates the relative expression level of miRNAs. (TIF)