Next Article in Journal
More Attention than Usual”: A Thematic Analysis of Dog Ownership Experiences in the UK during the First COVID-19 Lockdown
Next Article in Special Issue
Gene Expression Profile in Peripheral Blood Nuclear Cells of Small Ruminant Lentivirus-Seropositive and Seronegative Dairy Goats in Their First Lactation
Previous Article in Journal
Mineral Element Deposition and Gene Expression across Different Tissues of Cherry Valley Ducks
Previous Article in Special Issue
Exploring the Rumen and Cecum Microbial Community from Fetus to Adulthood in Goat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Roles of the miRNAome and Transcriptome in the Ovine Ovary Reveal Poor Efficiency in Juvenile Superovulation

Institute of Animal Husbandry and Veterinary Medicine, Tianjin Academy of Agricultural Sciences, Tianjin 300381, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this study.
Animals 2021, 11(1), 239; https://doi.org/10.3390/ani11010239
Submission received: 26 November 2020 / Revised: 9 January 2021 / Accepted: 12 January 2021 / Published: 19 January 2021

Abstract

:

Simple Summary

Using the technology of juvenile superovulation, more follicles can be acquired in juvenile animals than in adult animals. However, oocytes derived from the follicles of juvenile animals are usually of poor quality, meaning that they have lower levels of subsequent maturation and embryonic development. In the present study, we used an exogenous hormone treatment to stimulate Hu sheep in order to compare the differences in ovarian superovulation effects and serum hormone secretion in juvenile and adult sheep. Differentially expressed microRNA (miRNA) and messenger RNA (mRNA) from the ovaries of juvenile and adult Hu sheep were then investigated using high-throughput sequencing technology to reveal the formation mechanism of large numbers of follicles and poor oocyte quality in juvenile ovaries under superovulation treatment. We found that molecules of oar-miR-143 and follicle-stimulating hormone receptor (FSHR), among others, might regulate follicle formation, while oar-miR-485-3p, oar-miR-377-3p, and pentraxin 3 (PTX3), among others, may be associated with oocyte quality. The results will help us to identify miRNAs and mRNAs that could be used to predict ovarian superovulation potential and oocyte quality in the future.

Abstract

Juvenile superovulation can provide a wealth of oocyte material for embryo production, animal cloning, and genetic modification research, but embryos derived from juvenile oocytes show poor efficiency in subsequent developmental capacity. In order to reveal the formation mechanism of large numbers of follicles and poor oocyte quality in juvenile ovaries under superovulation treatment, differentially expressed microRNAs (miRNAs) and messenger RNAs (mRNAs) were characterized and investigated in the ovaries of lambs and adult sheep using high-throughput sequencing technology. The majority of differentially expressed miRNAs (337/358) were upregulated in lamb libraries. The expression levels of mRNAs related to hormone receptors (follicle-stimulating hormone receptor, FSHR; luteinizing hormone/choriogonadotropin receptor, LHCGR; estrogen receptor 1, ESR1), steroid hormone secretion (cytochrome P450 family 11 subfamily A member 1, CYP11A1; cytochrome P450 family 17 subfamily A member 1, CYP17A1; cytochrome P450 family 19 subfamily A member 1, CYP19A1), and oocyte quality (pentraxin 3, PTX3; BCL2 apoptosis regulator, BCL2; caspase 3, CASP3) were significantly different between the lamb and adult libraries. The miRNA aor-miR-143, which targets FSHR, was highly and differentially expressed, and PTX3 was predicted to be targeted by oar-miR-485-3p and oar-miR-377-3p in the ovine ovary. A considerable number of miRNAs were predicted to inhibit ESR1 expression in lamb ovaries. In conclusion, oar-miR-143 and FSHR molecules, among others, might regulate follicle formation, and oar-miR-485-3p, oar-miR-377-3p, and PTX3, among others, may be associated with oocyte quality. These identified miRNAs and mRNAs will be beneficial for the prediction of ovarian superovulation potential and screening of oocytes.

1. Introduction

In the breeding of domestic animals (sheep, cattle, pigs, etc.), limiting factors such as the reproductive cycle and generation interval affect their genetic progress. Exploring the reproductive potential of juvenile animals can provide an exciting shortcut to solve the above-mentioned limiting factors [1,2]. Although the gonadal systems of prepubertal females are not fully developed, primordial germ cells are already present in their ovaries [3]. Under the stimulation of superovulation hormones, prepubertal females of a suitable age can produce considerably more antral follicles than adult females [1]. However, for unknown physiological reasons, these antral follicles from prepubertal females cannot ovulate normally, and follicular aspiration on the juvenile ovary is necessary for oocyte collection [4]. The oocytes are then matured and fertilized in vitro, and the produced embryos are transplanted into the recipients. Compared with the technique of multiple ovulation and embryo transfer (MOET) applied in adult females, the above process, so-called juvenile in vitro embryo transfer (JIVET), is not smooth, and it has been found that embryos derived from juvenile oocytes show lower rates of development [5]. We suspected that the level of ovarian development in juvenile animals may be related to miRNA and mRNA expression levels, which contribute to limit its application on a large scale.
The miRNAs are a class of highly conserved small non-coding RNAs that can modulate mRNA and protein expression [6]. It has been proven that miRNAs participate in various physiological processes, including hormone secretion [7], organogenesis [8], cell proliferation [9], apoptosis [10], differentiation [11], metabolism [12], and reproduction control [13]. It has been estimated that miRNAs, which occupy a very small proportion of predicted genomes, could regulate up to one-third of all protein-coding genes in higher eukaryotes [14]. In recent years, many studies have further shown that miRNAs are involved in mammalian follicular, oocyte, and luteal development [15,16]. Otsuka et al. proved that miR-17-5p and let-7b contribute to the infertility of mutant mice [17]. Thereafter, the role of miRNAs in animal ovaries has been an area of constant interest for researchers. Yin et al. suggested that in mouse granulosa cells, miRNAs of miR-320 target E2F1 and SF-1, which can regulate proliferation and steroid production in follicular development, while miR-320 expression is regulated by follicle-stimulating hormone (FSH) secretion; moreover, miR-383 enhances miR-320-mediated suppression of granulosa cell proliferation [18]. With the advancement of high-throughput sequencing technology, more miRNAs in ovaries have been identified and correlated with reproduction-related traits [19,20,21].
Many studies have shown that oocytes derived from juvenile animals show poor subsequent development [22]. In 2015, our team also reported that the low global DNA methylation and hydroxymethylation in oocytes derived from lambs are associated with less subsequent developmental potential [23]. Recently, ovine mRNA expression in prepubertal and adult superstimulated follicle granulosa cells was studied using RNA sequencing technology, and more than 300 differentially expressed mRNAs were reported; in particular, the beta-estradiol upstream regulator in the prepubertal ovary was of interest [24]. In the present study, we used an exogenous hormone treatment to stimulate Hu sheep in order to compare the differences in ovarian superovulation effects and serum hormone secretion in juvenile and adult sheep; then, differentially expressed miRNAs and mRNAs from the ovaries of juvenile and adult Hu sheep were investigated using high-throughput sequencing technology to reveal the formation mechanism of large numbers of follicles and poor oocyte quality in juvenile ovaries under superovulation treatment. The results will help us to identify miRNAs and mRNAs that might be used to predict ovarian superovulation potential and oocyte quality in the future.

2. Materials and Methods

2.1. Animals and Superovulation Treatment

The two groups comprised six Hu sheep each. Animals were aged 1 month (lamb group) or 24 months (adult group), and were selected for superovulation treatment in the spring. The lambs were given intramuscular injections of FSH (Sansheng, Ningbo, China) every 12 h for 2 days (250 IU in total); pregnant mare serum gonadotropin (PMSG, Sansheng, Ningbo, China) was administered intramuscularly at the last FSH treatment (250 IU in total). An intravaginal controlled progesterone release device (CIDR; Pharmacia and Upjohn Co., Hartwell, Australia) insert was placed inside each adult sheep on Day 0, and removed on Day 12. FSH (300 IU in total) was injected intramuscularly every 12 h for the last 2 days. PMSG (360 IU in total) was injected intramuscularly 6 h after the CIDR was removed. All experimental procedures were approved by the Science Research Department of the Tianjin Academy of Agricultural Sciences (TAAS; Tianjin, China). Ethical approval for animal surgery was authorized by the animal welfare committee of TAAS (No. 2020009).

2.2. Sample Collection and RNA Extraction

Jugular venous blood and ovaries of lambs were obtained 14 h after PMSG injection, and jugular venous blood and ovaries of adult sheep were obtained 36 h after CIDR removal. Samples of 10 mL jugular blood were collected in BD SST Tubes (BD, Franklin Lakes, NJ, USA) and held at room temperature for 30 min, and then centrifuged at 1500× g for 10 min. The product of serum was used for determination of FSH, luteinizing hormone (LH), progesterone (P4), and estradiol (E2) concentrations. Each hormone was measured using an iodine [125I] radioimmunoassay kit (BNIBT, Beijing, China) according to the instructions. Three sheep were randomly selected from each group for anesthesia; following anesthesia, the abdominal skin and muscle were cut open with a scalpel to a length of 5cm, and the left ovary was pulled out of the body. After tubal ligation, the left ovary was removed and frozen instantly in liquid nitrogen. In the liquid nitrogen freezing environment, each whole ovary was ground into powder. Finally, the powder was collected for RNA extraction using the Trizol method.

2.3. Library Construction and Sequencing

For miRNA library construction, small RNAs in the size range of 18–30 nt were isolated from total RNA based on 15% denaturing polyacrylamide gel electrophoresis (PAGE). Subsequently, the 3′ adapters were added and 36–44 nt RNAs were enriched. The 5′ adapters were then also ligated to the RNAs. The ligation products were reverse-transcribed by PCR amplification, and the PCR products with a size of 140–160bp were enriched to generate a cDNA library. The libraries were sequenced through an Illumina HiSeqTM 2500 at Gene Denovo Biotechnology Co. (Guangzhou, China).
For mRNA library construction, poly-T oligo-attached magnetic beads were used to purify the RNA samples, and mRNAs were subsequently fragmented. The fragments were then reversed for cDNA synthesis using random primers. With the action of polymerase I and RNase H, second-strand cDNAs were synthesized. Next, a QiaQuick PCR extraction kit (Qiagen, Venlo, Netherlands) was used to purify the cDNA fragments. The cDNA fragments underwent end repair with poly (A) added and were ligated to Illumina sequencing adapters. The second-strand cDNA was then digested and selected by agarose gel electrophoresis. After PCR amplification, the libraries were sequenced by Gene Denovo Biotechnology Co. (Guangzhou, China) using the Illumina HiSeqTM 4000 platform.

2.4. Read Filtering and Differential Expression Analysis

For miRNA libraries, contaminants with an adaptor, poly A, and low-quality reads, or those shorter than 18nt, were filtered to generate clean reads. In order to remove rRNAs, tRNAs, snRNAs, and snoRNAs, we aligned the clean reads in GenBank and Rfam11.0 [25]. Clean reads were also aligned with the miRNA precursors/mature miRNAs of all animals in miRBase 22.0 (http://www.mirbase.org/), and the sequence and count of miRNAs in each library were displayed. The prefix “oar-miR” indicates that the miRNA has been included in the miRBase database for the species Ovis aries; no prefix indicates that the miRNA has not been included in the species database. Moreover, “miR-number-x” and “miR-number-y” indicate that the sequence can be aligned to -5 p and -3 p for the corresponding miRNA in miRBase for other animal species. Principal component analysis (PCA) was performed based on the sample correlation matrix. Differentially expressed miRNAs in lamb and adult groups were screened with a threshold of false discovery rate (FDR) < 0.05, log2 (fc) > 1 or log2 (fc) < –1.
For mRNA libraries, contaminants with an adapter, poly-N, and low-quality reads were removed to generate clean reads. Based on Hisat2 v2.1.0 (https://daehwankimlab.github.io/hisat2/) [26], the clean reads were mapped into the ovine genome (Oar_v4.0). The read number of each mRNA was calculated using featureCounts v1.5.0-p3 (https://sourceforge.net/projects/subread/files/subread-1.5.0-p3/) [27]. PCA was again performed based on the sample correlation matrix. Differentially expressed mRNAs were screened with a threshold of FDR < 0.05, log2 (fc) > 1 or log2(fc) < –1.

2.5. Targeted Relationship between Differentially Expressed miRNAs and mRNAs

The 2–8 nt sequences, starting from the 5′ end of miRNA, were chosen as seed sequences for prediction with the 3′-UTR of target mRNAs. RNAhybrid (version 2.1.2) [28], svm_light (version 6.01) [29], Miranda (version 3.3a) [30], and TargetScan (version 7.0) [31] were used to predict miRNA targets. The intersections of these results were more likely to be chosen as predicted miRNA target genes. Upon combining the mRNAs targeted by differentially expressed miRNAs and the differentially expressed mRNAs in the present study, an intersection was generated. We then constructed an miRNA–mRNA expression regulatory network based on the negative correlation between the expression of miRNAs and mRNAs using Cytoscape_v3.7.1 [32].

2.6. Enrichment Analysis by GO and KEGG

Predicted target mRNAs of differentially expressed miRNAs and differentially expressed mRNAs were used separately for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. These mRNAs were blasted to all known animals in the GO and KEGG databases. Gene numbers for terms in GO and KEGG were calculated, and a hyper-geometric test was used to find significantly enriched GO terms and KEGG pathways. A term or pathway with a corrected p-value (Q-value) < 0.05 was defined as significantly enriched.

2.7. Real-Time PCR Validation

Real-time PCR was performed on the miRNAs of oar-miR-143, miR-224-x, oar-miR-377-3p, oar-miR-485-3p, and oar-miR-487a-3p, as well as mRNAs of the follicle stimulating hormone receptor (FSHR), luteinizing hormone/choriogonadotropin receptor (LHCGR), estrogen receptor 1 (ESR1), pentraxin 3 (PTX3), and the BCL2 apoptosis regulator (BCL2). The primers of these molecules are shown in Supplementary File 12. Total RNA from sequencing samples was separately reversed into cDNA for miRNA and mRNA validation assay using the miScript Reverse Transcription Kit (Qiagen, Dusseldorf, Germany) and Prime Script® RT reagent Kit (Takara Bio Inc., Dalian, China). All cDNAs from miRNAs and mRNAs were separately mixed to generate standard samples for miRNA and mRNA validation assay. U6 and RPL19 were separately used as reference miRNA and mRNA to calculate the relative expression level with the 2−ΔΔCT method. The real-time PCR procedure for miRNA was as follows: incubation at 95 °C for 10 min, followed by 40 cycles at 94 °C for 15 s and 60 °C for 45 s. The real-time PCR procedure for mRNAs was as follows: incubation at 95 °C for 5 min, followed by 40 cycles at 95 °C for 10 s and 60 °C for 30 s. All reactions were performed in triplicate.

2.8. Statistical Analysis

The chi-square test was used to analyze categorical variables of follicle numbers. Duncan’s multiple range test program in ANOVA was adopted to analyze continuous variables of hormone concentration and mRNA/miRNA expression based on the SAS 8.0 software (SAS Institute Inc., Cary, NC, USA). All of the results are presented as mean ± SE. To assess the relationship between high-throughput sequencing and real-time PCR in mRNA/miRNA expression assay, the Spearman correlation was calculated based on SPSS version 20.0 (SPSS Inc., Chicago, IL, USA)).

3. Results

3.1. Measurement of Follicle Numbers and Serum Hormone in Lamb and Adult Sheep under Superovulation Treatment

Under stimulation with the exogenous hormone, superovulation was successfully performed in lambs and adult Hu sheep. The follicle numbers and serum reproduction hormone levels are presented in Table 1. The mean super-stimulated follicle number (total of two ovaries) in lambs was 70.17 ± 5.14, which was much higher than that in adult sheep (20.17 ± 2.65) with a high significance (p < 0.01). However, the concentrations of FSH, LH, progesterone (P4), and estradiol (E2) in jugular vein serum all showed no significant difference between the lamb and adult groups after superovulation treatment.

3.2. Overview of Sequencing Data

Three biological repetitions in each group (lamb and adult) were selected to construct small RNA and mRNA libraries for sequencing analysis. For the six small RNA libraries, a total of 98,479,314 high-quality reads were obtained, and a total of 94,250,039 clean tags (Supplementary File 1), with lengths ranging from 18 to 30 nt (Figure 1A), were generated after discarding contaminants. As shown in Supplementary File 2, average conserved miRNAs (including exist_mirna, exist_mirna_edit and known_mirna) accounted for 89.59% ± 1.60% of the total clean tags in the lamb libraries, which was significantly lower than 96.29 ± 0.40% in the adult libraries (p < 0.05). This means that more miRNAs may be found in prepubertal ovine ovaries in the future.
For the six mRNA libraries, an average of 77,998,161 clean reads in each library were generated after discarding adapter–adapter and low-quality reads (Supplementary File 3). The Q20 values of the mRNA libraries were all higher than 97%. Supplementary File 4 shows the results of clean reads mapping the ovine genome. More than 96% and 93% of clean reads were totally mapped and uniquely mapped, respectively, in the ovine genome for all six mRNA libraries. In summary, the data quality of the miRNA and mRNA libraries in the present study was high enough for them to be used for subsequent statistical analyses.

3.3. Differentially Expressed miRNA and mRNA in the Lamb and Adult Libraries

A total of 1819 miRNAs were detected in all six miRNA libraries, and random clustering of the replicates in the PCA of all detected miRNAs showed that the first principal component explained 70.6% of the total variation between lamb and adult groups (Figure 1B). As shown in Figure 1C, 358 miRNAs were differentially expressed between the lamb and adult libraries (FDR < 0.05, log2 (fc) > 1 or log2 (fc) < –1) (Supplementary File 5). Compared to the adult libraries, 337 miRNAs were upregulated in the lamb libraries, which accounted for the majority of differentially expressed miRNAs, and only 21 miRNAs were down regulated. Heat-map analysis of differentially expressed miRNAs showed a separation between the individuals of lamb and adult sheep (Figure 1D). The top 10 highly and differentially expressed miRNAs that might play important roles in sexual maturation of the ovary are listed in Table 2.
A read count greater than 1 defined detectable genes in mRNA expression level, and a total of 17,842 mRNAs were detected in the six mRNA libraries (Supplementary File 6). Correlation coefficients within the lamb group were all higher than 0.979, and correlation coefficients within the adult group were all higher than 0.969 (Figure 2A). PCA revealed a clear distinction in mRNA expression of ovaries between lamb and adult individuals (Figure 2B). Finally, 3150 mRNAs were found to be differentially expressed in the two groups (Figure 2C, Supplementary File 6). Compared to the adult libraries, 1359 and 1791 mRNAs were upregulated and downregulated, respectively, in the lamb libraries. In order to examine the differences in the sensitivity of juvenile and adult sheep to superovulation hormones, as well as expression levels of key mRNAs related to oocyte quality, we selected some mRNAs to emphasize their expression differences between the two groups. As shown in Table 3, the expression levels of hormone receptor genes (follicle-stimulating hormone receptor, FSHR; luteinizing hormone/choriogonadotropin receptor, LHCGR; estrogen receptor 1, ESR1), steroid hormone secretion genes (cytochrome P450 family 11 subfamily A member 1, CYP11A1; cytochrome P450 family 17 subfamily A member 1, CYP17A1; cytochrome P450 family 19 subfamily A member 1, CYP19A1; steroidogenic acute regulatory protein, STAR; 3 beta-hydroxysteroid dehydrogenase, HSD3B; 17-beta-hydroxysteroid dehydrogenase type 1, HSD17B1), and genes related to oocyte quality (pentraxin 3, PTX3; growth differentiation factor 9, GDF9; DNA methyltransferase 1, DNMT1, caspase 3, CASP3; insulin like growth factor binding protein 4, IGFBP4; BCL2 apoptosis regulator, BCL2; apoptosis regulator BAX, BAX) were significantly different in the ovaries of lambs and adult sheep (FDR < 0.05). Heat map analysis of differentially expressed mRNAs also showed a separation between the individual lambs and adult sheep (Figure 2D).

3.4. Analysis of the Targeted Relationship between Differentially Expressed miRNAs and mRNAs

Based on pairing between the 5′ end of miRNAs and the 3′ untranslated region of target mRNAs, miRNAs can modulate mRNA expression at the post-transcriptional level [6]. Taking the intersection of the differentially expressed miRNA target mRNAs with differentially expressed mRNAs, and screening with a negative correlation between the expression of the miRNA and mRNA, we constructed an miRNA–mRNA expression regulatory network. As shown in Supplementary File 7, 2194 mRNAs were predicted to be targeted by the 358 differentially expressed miRNAs. The differentially expressed miRNAs that targeted key genes (FSHR, LHCGR, CYP11A1, STAR, PTX3, and ESR1) are listed in Table 4 and Figure 3. The miRNA oar-miR-143, which targeted FSHR at a high level and with differential expression in ovaries between lamb and adult sheep [33,34], and the PTX3 gene, which was related to oocyte quality and targeted by oar-miR-377-3p, oar-miR-485-3p, oar-miR-487a-3p, miR-496-x, and miR-767-x, are particularly noteworthy. The most notable feature was the part of the miRNA–mRNA network related to the ESR1 and PTX3 gene (Figure 3).

3.5. GO and KEGG Enrichment Analysis for Target Genes of Differentially Expressed miRNAs

For target genes of miRNAs that were differentially expressed in the lamb and adult groups, a total of 21,194 target genes were mapped to GO in terms of the cellular component (CC): 19,819 to molecular function (MF) and 21,888 to biological processes (BP) (Supplementary File 8). Comparing the reference gene background in CC, MF, and BP, a total of 105, 99, and 458 GO terms were significantly enriched (Q-value < 0.05), respectively. GO terms of the hormone-mediated signaling pathway (GO: 0009755), cellular response to hormone stimulus (GO: 0032870), response to hormone (GO: 0009725), histone methylation (GO: 0016571), protein methylation (GO: 0006479), and histone lysine methylation (GO: 0034968) were significantly enriched in the target genes (Q-value < 0.05). This suggests that the ovaries of lambs and adult sheep differ in their response to hormones and in their methylation levels. In addition, GO terms of epithelial cell proliferation (GO: 0050673), cell differentiation (GO: 0030154), programmed cell death (GO: 0012501), and regulation of MAPK cascade (GO: 0043408) were significantly enriched (Q-value < 0.05), which suggests that cell proliferation, differentiation, and apoptosis behaviors in the ovaries of lambs and adult sheep are also different.
KEGG pathway enrichment showed that a total of 16,627 background genes were annotated in 344 biological functions. In the present study, we found that 13,290 target genes for differentially expressed miRNAs were annotated in KEGG pathways. A total of 180 pathways were significantly enriched (Q-value < 0.05) (Supplementary File 9). It is worth mentioning that progesterone-mediated oocyte maturation (ko04914), oocyte meiosis (ko04114), the estrogen signaling pathway (ko04915), the TGF-beta signaling pathway (ko04350), and the MAPK signaling pathway (ko04010) were all significantly enriched in target genes for miRNAs differentially expressed in the lamb and adult groups (Q-value < 0.05). These results suggest that oocytes in lambs and adults have different maturation levels.

3.6. GO and KEGG Enrichment Analysis of Differentially Expressed mRNAs

A total of 1944, 1763, and 2004 differentially expressed mRNAs were annotated in the GO terms of BP, CC, and MF, respectively. Compared to the reference gene background in CC, MF, and BP, a total of 39, 11, and 309 GO terms were significantly enriched (Q-value < 0.05), respectively (Supplementary File 10). In particular, the GO terms of cell proliferation (GO: 0008283), steroid metabolic process (GO: 0008202), sterol metabolic process (GO: 0016125), steroid biosynthetic process (GO: 0006694), regulation of secretion (GO: 0051046), and MAPK cascade (GO: 0000165) were significantly enriched (Q-value < 0.05). This means that the significantly enriched GO terms in mRNA and miRNA analysis had similar results in terms of cell proliferation and hormone secretion regulation.
Subsequently, enrichment analysis of the KEGG pathway was performed based on the differentially expressed mRNAs. A total of 1442 mRNAs were mapped to the relevant KEGG pathways, and 48 KEGG pathways were shown to be significantly enriched (Q-value < 0.05) (Supplementary File 11). In the analysis of mRNA expression level, the pathways of glutathione metabolism (ko00480) and ovarian steroidogenesis (ko04913) attracted our research attention. However, neither of these two significant pathways presented as enriched in the miRNA KEGG analysis. In addition, many significant enrichment pathways in the miRNA analysis were also absent from the mRNA KEGG analysis.

3.7. Real-Time PCR Validation of mRNA and miRNA Associated with Hormone Receptor and Oocyte Quality

The expression levels of five miRNAs (oar-miR-143, miR-224-x, oar-miR-377-3p, oar-miR-485-3p, and oar-miR-487a-3p) and five mRNAs (FSHR, LHCGR, PTX3, BCL2, and CASP3) were validated in the ovaries of lambs and adult sheep using real-time PCR (Figure 4). In these selected molecules, FSHR was targeted by oar-miR-143 [34,36], and PTX3 was targeted by miR-224-x, oar-miR-377-3p, oar- oar-miR-485-3p, and miR-487a-3p [37]. mRNA and miRNA expression values from each sample, detected by high-throughput sequencing and real-time PCR were used for Pearson correlation analysis (Supplementary File 12). The correlation coefficients of the two methods in expression detection of the 10 selected molecules ranged from 0.635 to 0.978, which demonstrates that the expression levels detected by high-throughput sequencing were reliable.

4. Discussion

JIVET technology can be used to shorten the generation interval and accelerate the progress of genetic selection, and it can also provide a wealth of materials for research on embryonic stem cells, animal cloning, and genetic modification [3,22,38]. Nevertheless, embryos derived from JIVET exhibit a significantly lower survival rate [38]. The maturation status of the cumulus cell [39], oocyte cytoplasm [40], oocyte nuclear [41], and follicular fluid microenvironment in juvenile ovaries [42,43,44] have been suggested to explain the poor efficiency of JIVET. In the present study, we focused on the initial stage of JIVET: ovaries under the treatment of superovulation; therefore, the miRNAome and transcriptome were compared between juvenile and adult ovaries to identify the mechanism of inferior oocyte capture by JIVET.
We initially identified 358 miRNAs that were differentially expressed in the lamb and adult libraries, and the majority of these miRNAs were upregulated in the lamb libraries. Considering that miRNAs can negatively regulate gene expression at the post-transcriptional level, we speculate that much more gene expression in lamb ovaries was suppressed. Under superstimulatory treatment of prepubertal and adult sheep, Wu et al. found 175 downregulated genes vs. 136 upregulated genes in prepubertal granulosa cells using RNA sequencing [24], and they also found 155 downregulated proteins vs. 88 upregulated proteins in prepubertal follicular fluid using iTRAQ [43]. In our results, we also found 1791 downregulated mRNAs vs. 1359 upregulated mRNAs in the ovaries of lamb. The results from the research of Wu et al. and our mRNA results verify the findings of the miRNAs comparisons. Based on an enzyme immunoassay evaluating the effects of 80 miRNAs on human granulosa cells, Sirotkin et al. showed that miR-145 enhanced progesterone production up to 5-fold, and miR-144 reduced estradiol secretion up to 4.5-fold [45]. In this study, miR-145-x and miR-144-y were also highly expressed and were down- and upregulated, respectively, in the lamb libraries. In addition, downregulation of the miRNAs oar-miR-10a and miR-145-x in lambs was speculated in the present study, implicating suppression of granulosa cell proliferation by the targeting of brain-derived neurotropic factor (BDNF) and activin receptor IB (ACVRIB) [7,46,47].
The cyclic activity of the ovary is mainly influenced by reproductive hormones [48,49]; under the stimulation of FSH and LH, the follicles undergo primary, secondary, and antral stages, and only a few of them reach the pre-ovulatory stage [50]. After superovulation treatment in the present study, we found that there was no significant difference in the concentrations of FSH and LH in jugular vein serum between lambs and adult sheep. However, the exciting finding was that FSHR and LHCGR expression levels in lamb ovaries were significantly higher than those in adult sheep ovaries. Specifically, at the same hormone level, lambs have an increased sensitivity to hormones due to the increased expression of hormone receptors in ovarian tissue. Based on RT-PCR, immunohistochemistry, and Western blot results, Scarlet et al. also showed that FSHR was present in the ovaries of prepubertal fillies, and found that FSHR was abundant in the ovarian stroma cells of neonate fillies but not of adults [51]. LHCGR is also considered to be an important target gene induced by FSH and is closely related to steroid hormone metabolism and ovarian follicular development [52]. Recently, Michalovic et al. reported that cumulus–oocyte complexes (CoCs) and the granulosa cells (GCs) collected from Holstein heifers that matured under FSH treatment presented increased mRNA levels of genes promoting follicular growth (FSHR) and estradiol synthesis (CYP19A1) [53]. In this study, the mRNA expression levels of the steroid hormone or estradiol synthesis genes (CYP11A1, CYP17A1, CYP19A1, STAR, HSD3B, and HSD17B1) in lamb ovaries were also significantly higher than those in adult sheep, which means that estradiol synthesis is occurring. However, we found that mRNA expression of the ESR1 gene in lambs ovaries was significantly inhibited, which may be related to the physiological state of prepuberty and detrimental to follicular development. In mice, females lacking the ESR1 gene are infertile and non-receptive to males [54]. In the present study, a considerable number of miRNAs were predicted to target the ESR1 mRNA in lamb ovaries, and these miRNAs were significantly highly expressed in lamb. Couse and Korach demonstrated that ESR1 plays an important role in follicular growth and development and subsequent ovulation [55]. Liu et al. showed that the estradiol (E2) estrogen receptors (ERs) system in follicular granulosa cells plays a dominant role in controlling oocyte meiotic resumption in mammals [56]. In human and mouse cumulus cells, downregulation of ESR1 could be related to oocyte competence and is likely to be the driver of the expression changes highlighted in the PI3K/AKT pathway [57].
The poor efficiency of juvenile superovulation is mainly reflected in the poor quality of oocytes collected from prepubertal female animals [58]. Zhang et al. reported that the PTX3 gene in human cumulus cells is highly associated with oocyte development [59]. In the present study, the downregulated PTX3 mRNA levels in lamb ovaries suggested that oocytes derived from lambs may be of poor quality. However, growth differentiation factor-9 (GDF-9), a gene upstream of PTX3, was upregulated in lambs. Kona et al. reported that the expression of GDF9 was the highest in ovine primordial follicles compared with the preantral, early antral, antral, and large antral stages [60]. Upregulated GDF9 in the lamb group indicated that ovaries from lambs contain more primordial follicles than adult sheep. However, we found that mRNA expression levels of the BCL2 gene in lamb ovaries were downregulated, but we did not find a difference in BAX. It is well known that the expression ratio of Bax and Bcl-2 determines the propensity for cell apoptosis or survival [61]. In the present study, because the ratio of Bax and Bcl-2 was significantly lower in the lamb group than the adult group, we speculated that the proliferative capacity of granulosa cells in lamb ovaries is not good. DNMT1 is responsible for maintaining methylation patterns throughout DNA replication, and we found that its mRNA expression is higher in lamb ovaries. However, Fang et al. found that CoCs derived from lamb matured in vitro with a lower mRNA expression of DNMT1 [23]. This suggests that CoCs and ovarian tissue of lamb may be inconsistent in methylation levels. CASP3 involved in apoptosis was proven to be responsible for follicular atresia [62]; the higher level of mRNA expression of CASP3 in lamb ovaries in this study may reflect that a more follicular atresia events occur. IGFBP4 mRNA transcripts were detected in human follicular fluid, and this has been demonstrated to be associated with oocyte maturation and embryo development [63]. This suggests that the lower mRNA expression level of IGFBP4 detected in lamb ovaries may lead to the emergence of low-quality oocytes.
miRNA and mRNA together have been proven to play important roles in the process of follicular development and oocyte maturation in the ovary [50,64,65,66,67,68]. Using the combined miRNAome and transcriptome results, we constructed miRNA regulatory networks for key genes in the ovine ovary. What needs to be emphasized is that downregulation of ESR1 and PTX3 may be attributed to the targeted inhibition of a considerable number of miRNAs in the lamb ovary. Unfortunately, we found that miR-224, which targets the PTX3 gene in sow ovaries [37], was not negatively correlated with the expression level of PTX3 in this study. The miRNAs oar-miR-485-3p, oar-miR-487a-3p, and oar-miR-377-3p, which are species-specific to sheep and, target PTX3 in the ovine ovary, should be further studied in the future. Last year, Zhang et al. demonstrated using bovine granulosa cells that overexpression of miR-143 led to decreased secretion of estradiol, and silencing of miR-143 was associated with increased estradiol production [34]. This inhibitory action of miR-143 on estradiol synthesis is mediated via its target gene, FSHR [7,34]. A previous study also reported high expression of miR-143 in granulosa cells of primary, secondary, and antral follicles associated with a negative regulation of FSH-induced estradiol production [36]. In the present study, the downregulation of oar-miR-143 in lamb ovaries may have directly contributed to the upregulated of FSHR expression. The expression relationship between miR-513-x and LHCGR confirmed the results in human granulosa cells [35].
Each miRNA might target hundreds of genes, and in turn, each of these genes might be regulated by several miRNAs [69]. Therefore, placing too much emphasis on the expression levels of single miRNAs or mRNAs is not advisable [50]. In order to comprehensively analyze the differentially expressed miRNAs and mRNAs in lamb and adult sheep ovaries, GO term and KEGG pathway enrichment analysis was performed. Compared with the enrichment analysis of differentially expressed mRNAs, many more GO terms and KEGG pathways were enriched in the target genes of differentially expressed miRNAs. However, mRNA and miRNA target genes were all enriched in GO terms of cell proliferation, which is consistent with the results reported on the protein level [43,70]. However, we found little difference in KEGG enrichment between differentially expressed mRNA and target genes of miRNA. This result may be attributed to the difference in spatiotemporal expression between miRNA and mRNA, as miRNA needs to find the 3′ untranslated regions of target mRNA before it can play a regulatory role [6]. On the basis of our observations of combination of the mRNAs and target genes of miRNAs in KEGG analysis, we firmly believe that there are some significant differences between the ovaries of lambs and adult sheep in terms of cell proliferation, hormones response, and oocyte development level.

5. Conclusions

In conclusion, the majority of differentially expressed miRNAs (337/358) were upregulated in the lamb libraries. The expression levels of mRNAs related to hormone receptors (FSHR, LHCGR, and ESR1), steroid hormone secretion (CYP11A1, CYP17A1, and CYP19A1), and oocyte quality (PTX3, BCL2, and CASP3) showed a significant difference between the lamb and adult sheep libraries. The miRNA oar-miR-143, which targets FSHR, was highly differentially expressed, and PTX3 was predicted to be targeted by oar-miR-485-3p and oar-miR-377-3p in the ovine ovary in this study. A considerable number of miRNAs were predicted to inhibit ESR1 expression in lamb ovaries. We suppose that molecules of oar-miR-143 and FSHR, among others, might regulate follicle formation, and oar-miR-485-3p, oar-miR-377-3p, and PTX3, among others, may be associated with oocyte quality. These identified miRNAs and mRNAs will be beneficial for potential ovarian superovulation prediction and oocyte screening.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/2076-2615/11/1/239/s1. Supplementary File 1: Data filtration of the miRNA libraries; Supplementary File 2: Sequence distribution of small RNA; Supplementary File 3: Statistics of sequencing quality for the mRNA libraries; Supplementary File 4: Sequence distribution of mRNA; Supplementary File 5: miRNA expression profiles of ovaries in lamb and adult sheep, Supplementary File 6: mRNA expression profiles of ovaries in lamb and adult sheep; Supplementary File 7: Target gene prediction of differentially expressed miRNAs; Supplementary File 8: GO analysis of target genes for differentially expressed miRNAs in the lamb and adult sheep libraries; Additional File 9: KEGG analysis of target genes for differentially expressed miRNAs in the lamb and adult sheep libraries; Supplementary File 10: GO analysis of differentially expressed mRNAs in the lamb and adult sheep libraries; Supplementary File 11: KEGG analysis of differentially expressed mRNAs in the lamb and adult sheep libraries; Supplementary File 12: Pearson correlation analysis of high-throughput sequencing and real-time PCR in expression detection of five selected miRNAs and five selected mRNAs.

Author Contributions

X.Z., J.Z., and X.G. conceived and designed this study; X.Z., C.D., J.Y., Y.L., J.F., B.W. and J.Z. conducted the experiments; X.Z., X.G., and C.D. analyzed the data; X.Z., C.D., X.G., and J.Z. prepared the manuscript; X.Z., C.D. and X.G. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 31902150), the Youth Innovative Research and Experimental Project of Tianjin Academy of Agricultural Sciences (No. 201902 and 2020013), Tianjin Science and Technology Plan Project (No. 17ZXZYNC00040, 18ZXZYNC00180, and 19ZXZYSN00030), the Natural Science Foundation of Tianjin (No. 20JCQNJC00630), and Tianjin ‘131’ Innovative Talents Training Project.

Institutional Review Board Statement

The study was conducted according to the regulations and ethical guidelines set by the Science Research Department of the Tianjin Academy of Agricultural Sciences (TAAS; Tianjin, China). Ethical approval for animal surgery was authorized by the animal welfare committee of TAAS (No. 2020009).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available to preserve privacy of the data.

Acknowledgments

We are grateful to Gene Denovo Biotechnology Co. (Yang Liu and Zhenyue Yuan) for their assistance in the analysis of the targeted relationship between differentially expressed miRNA and mRNA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morton, K.M.; Catt, S.L.; Maxwell, W.M.C.; Evans, G. An efficient method of ovarian stimulation and in vitro embryo production from prepubertal lambs. Reprod. Fertil. Dev. 2005, 17, 701–706. [Google Scholar] [CrossRef] [PubMed]
  2. Daly, J.; Smith, H.; McGrice, H.A.; Kind, K.L.; van Wettere, W. Towards Improving the Outcomes of Assisted Reproductive Technologies of Cattle and Sheep, with Particular Focus on Recipient Management. Animals (Basel) 2020, 10, 293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Armstrong, D.T.; Kotaras, P.J.; Earl, C.R. Advances in production of embryos in vitro from juvenile and prepubertal oocytes from the calf and lamb. Reprod. Fertil. Dev. 1997, 9, 333–340. [Google Scholar] [CrossRef] [PubMed]
  4. Valasi, I.; Menegatos, I.; Papanikolaou, T.; Goulas, P.; Amiridis, G.S. Oocyte pick-up in juvenile lambs affects neither onset of puberty nor their future fertility. Theriogenology 2006, 66, 2144–2151. [Google Scholar] [CrossRef] [PubMed]
  5. Ptak, G.; Loi, P.; Dattena, M.; Tischner, M.; Cappai, P. Offspring from one-month-old lambs: Studies on the developmental capability of prepubertal oocytes. Biol. Reprod. 1999, 61, 1568–1574. [Google Scholar] [CrossRef] [Green Version]
  6. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function; Elsevier Inc.: Philadelphia, PA, USA, 2004; Volume 116, pp. 281–297. [Google Scholar]
  7. Azhar, S.; Dong, D.; Shen, W.J.; Hu, Z.; Kraemer, F.B. The role of miRNAs in regulating adrenal and gonadal steroidogenesis. J. Mol. Endocrinol. 2020, 64, R21–R43. [Google Scholar] [CrossRef]
  8. Hayashi, T.; Hoffman, M.P. Exosomal microRNA communication between tissues during organogenesis. RNA Biol. 2017, 14, 1683–1689. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, G.; Guo, X.; Cheng, L.; Chu, P.; Chen, M.; Chen, Y.; Chang, C. An integrated analysis of the circRNA-miRNA-mRNA network reveals novel insights into potential mechanisms of cell proliferation during liver regeneration. Artif. Cells Nanomed. Biotechnol. 2019, 47, 3873–3884. [Google Scholar] [CrossRef] [Green Version]
  10. Xu, P.; Wu, Q.; Lu, D.; Yu, J.; Rao, Y.; Kou, Z.; Fang, G.; Liu, W.; Han, H. A systematic study of critical miRNAs on cells proliferation and apoptosis by the shortest path. BMC Bioinform. 2020, 21, 1–14. [Google Scholar] [CrossRef]
  11. Martin, E.C.; Qureshi, A.T.; Dasa, V.; Freitas, M.A.; Gimble, J.M.; Davis, T.A. MicroRNA regulation of stem cell differentiation and diseases of the bone and adipose tissue: Perspectives on miRNA biogenesis and cellular transcriptome. Biochimie 2016, 124, 98–111. [Google Scholar] [CrossRef]
  12. Rottiers, V.; Naar, A.M. MicroRNAs in metabolism and metabolic disorders. Nat. Rev. Mol. Cell Biol. 2012, 13, 239–250. [Google Scholar] [CrossRef] [PubMed]
  13. Dahiya, N.; Sherman-Baust, C.A.; Wang, T.L.; Davidson, B.; Shih, I.; Zhang, Y.; Wood, W.R.; Becker, K.G.; Morin, P.J. MicroRNA expression and identification of putative miRNA targets in ovarian cancer. PLoS ONE 2008, 3, e2436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hu, S.; Ren, G.; Liu, J.; Zhao, Z.; Yu, Y.; Su, R.; Ma, X.; Ni, H.; Lei, W.; Yang, Z. MicroRNA Expression and Regulation in Mouse Uterus during Embryo Implantation. J. Biol. Chem. 2008, 283, 23473–23484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Donadeu, F.X.; Schauer, S.N.; Sontakke, S.D. Involvement of miRNAs in ovarian follicular and luteal development. J. Endocrinol. 2012, 215, 323–334. [Google Scholar] [CrossRef]
  16. Tesfaye, D.; Gebremedhn, S.; Salilew-Wondim, D.; Hailay, T.; Hoelker, M.; Grosse-Brinkhaus, C.; Schellander, K. MicroRNAs: Tiny molecules with a significant role in mammalian follicular and oocyte development. Reproduction 2018, 155, R121–R135. [Google Scholar] [CrossRef] [Green Version]
  17. Otsuka, M.; Zheng, M.; Hayashi, M.; Lee, J.D.; Yoshino, O.; Lin, S.; Han, J. Impaired microRNA processing causes corpus luteum insufficiency and infertility in mice. J. Clin. Investig. 2008, 118, 1944–1954. [Google Scholar] [CrossRef] [Green Version]
  18. Yin, M.; Wang, X.; Yao, G.; Lu, M.; Liang, M.; Sun, Y.; Sun, F. Transactivation of microRNA-320 by microRNA-383 regulates granulosa cell functions by targeting E2F1 and SF-1 proteins. J. Biol. Chem. 2014, 289, 18239–18257. [Google Scholar] [CrossRef] [Green Version]
  19. Zhai, M.; Xie, Y.; Liang, H.; Lei, X.; Zhao, Z. Comparative profiling of differentially expressed microRNAs in estrous ovaries of Kazakh sheep in different seasons. Gene 2018, 664, 181–191. [Google Scholar] [CrossRef]
  20. Pokharel, K.; Peippo, J.; Honkatukia, M.; Seppala, A.; Rautiainen, J.; Ghanem, N.; Hamama, T.M.; Crowe, M.A.; Andersson, M.; Li, M.H.; et al. Integrated ovarian mRNA and miRNA transcriptome profiling characterizes the genetic basis of prolificacy traits in sheep (Ovis aries). BMC Genom. 2018, 19, 104. [Google Scholar] [CrossRef] [Green Version]
  21. Miao, X.; Luo, Q.; Zhao, H.; Qin, X. Ovarian transcriptomic study reveals the differential regulation of miRNAs and lncRNAs related to fecundity in different sheep. Sci. Rep.-UK 2016, 6, 35299. [Google Scholar] [CrossRef]
  22. Morton, K.M. Developmental capabilities of embryos produced in vitro from prepubertal lamb oocytes. Reprod. Domest. Anim. 2008, 43 (Suppl. 2), 137–143. [Google Scholar] [CrossRef]
  23. Fang, Y.; Zhang, X.; Zhang, J.; Zhong, R.; Zhou, D. Global DNA methylation and related mRNA profiles in sheep oocytes and early embryos derived from pre-pubertal and adult donors. Anim. Reprod. Sci. 2016, 164, 144–151. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, Y.; Lin, J.; Li, X.; Han, B.; Wang, L.; Liu, M.; Huang, J. Transcriptome profile of one-month-old lambs’ granulosa cells after superstimulation. Asian-Australas. J. Anim. Sci. 2017, 30, 20–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Griffiths-Jones, S.; Bateman, A.; Marshall, M.; Khanna, A.; Eddy, S.R. Rfam: An RNA family database. Nucleic Acids Res. 2003, 31, 439–441. [Google Scholar] [CrossRef] [Green Version]
  26. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. [Google Scholar] [CrossRef]
  27. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kruger, J.; Rehmsmeier, M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34, W451–W454. [Google Scholar] [CrossRef] [PubMed]
  29. Joachims, T. Svmlight: Support Vector Machine. Available online: http://svmlight.joachims.org/ (accessed on 2 January 2021).
  30. Turner, D.A. Miranda: A Non-Strict Functional Language with Polymorphic Types. Conference on Functional Programming Languages and Computer Architecture; Springer: Berlin/Heidelberg, Germany, 1985; Volume 201, pp. 1–16. [Google Scholar]
  31. Lewis, B.P.; Burge, C.B.; Bartel, D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20. [Google Scholar] [CrossRef] [Green Version]
  32. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  33. Du, X.; Zhang, L.; Li, X.; Pan, Z.; Liu, H.; Li, Q. TGF-beta signaling controls FSHR signaling-reduced ovarian granulosa cell apoptosis through the SMAD4/miR-143 axis. Cell Death Dis. 2016, 7, e2476. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, Z.; Chen, C.Z.; Xu, M.Q.; Zhang, L.Q.; Liu, J.B.; Gao, Y.; Jiang, H.; Yuan, B.; Zhang, J.B. MiR-31 and miR-143 affect steroid hormone synthesis and inhibit cell apoptosis in bovine granulosa cells through FSHR. Theriogenology 2019, 123, 45–53. [Google Scholar] [CrossRef]
  35. Troppmann, B.; Kossack, N.; Nordhoff, V.; Schuring, A.N.; Gromoll, J. MicroRNA miR-513a-3p acts as a co-regulator of luteinizing hormone/chorionic gonadotropin receptor gene expression in human granulosa cells. Mol. Cell. Endocrinol. 2014, 390, 65–72. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, L.; Zhang, X.; Zhang, X.; Lu, Y.; Li, L.; Cui, S. miRNA-143 mediates the proliferative signaling pathway of FSH and regulates estradiol production. J. Endocrinol. 2017, 234, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Li, X.; Wang, H.; Sheng, Y.; Wang, Z. MicroRNA-224 delays oocyte maturation through targeting Ptx3 in cumulus cells. Mech. Dev. 2017, 143, 20–25. [Google Scholar] [CrossRef] [PubMed]
  38. Paramio, M.T.; Izquierdo, D. Current status of in vitro embryo production in sheep and goats. Reprod. Domest. Anim. 2014, 49 (Suppl. 4), 37–48. [Google Scholar] [CrossRef] [Green Version]
  39. Fang, Y.; Zhang, J.; Li, Y.; Guo, X.; Li, J.; Zhong, R.; Zhang, X. Melatonin-induced demethylation of antioxidant genes increases antioxidant capacity through RORalpha in cumulus cells of prepubertal lambs. Free Radical Biol. Med. 2019, 131, 173–183. [Google Scholar] [CrossRef]
  40. Van Blerkom, J. Mitochondrial function in the human oocyte and embryo and their role in developmental competence. Mitochondrion 2011, 11, 797–813. [Google Scholar] [CrossRef]
  41. Ptak, G.; Matsukawa, K.; Palmieri, C.; Della, S.L.; Scapolo, P.A.; Loi, P. Developmental and functional evidence of nuclear immaturity in prepubertal oocytes. Hum. Reprod. 2006, 21, 2228–2237. [Google Scholar] [CrossRef] [Green Version]
  42. Warzych, E.; Pawlak, P.; Pszczola, M.; Cieslak, A.; Lechniak, D. Prepubertal heifers versus cows-The differences in the follicular environment. Theriogenology 2017, 87, 36–47. [Google Scholar] [CrossRef]
  43. Wu, Y.; Lin, J.; Han, B.; Wang, L.; Chen, Y.; Liu, M.; Huang, J. Proteomic profiling of follicle fluids after superstimulation in one-month-old lambs. Reprod. Domest. Anim. 2018, 53, 186–194. [Google Scholar] [CrossRef] [Green Version]
  44. Tian, H.; Liu, K.; Zhang, Y.; Qi, Q.; Wang, C.; Guan, H.; Yan, F.; Hou, J. Adult follicular fluid supplementation during in vitro maturation improves the developmental competence of prepubertal lamb oocytes. Theriogenology 2019, 130, 157–162. [Google Scholar] [CrossRef] [PubMed]
  45. Sirotkin, A.V.; 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. [Google Scholar] [CrossRef] [PubMed]
  46. Jiajie, T.; Yanzhou, Y.; Hoi-Hung, A.C.; Zi-Jiang, C.; Wai-Yee, C. Conserved miR-10 family represses proliferation and induces apoptosis in ovarian granulosa cells. Sci. Rep.-UK 2017, 7, 41304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yan, G.; Zhang, L.; Fang, T.; Zhang, Q.; Wu, S.; Jiang, Y.; Sun, H.; Hu, Y. MicroRNA-145 suppresses mouse granulosa cell proliferation by targeting activin receptor IB. FEBS Lett. 2012, 586, 3263–3270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ascari, I.J.; Alves, N.G.; Alves, A.C.; Garcia, I.F.F.; Junqueira, F.B. Resumption of cyclic ovarian activity in postpartum ewes: A review. Semin. Ciên. Agrár. 2016, 37, 1101–1116. [Google Scholar] [CrossRef] [Green Version]
  49. Satue, K.; Fazio, E.; Medica, P. Can the Presence of Ovarian Corpus Luteum Modify the Hormonal Composition of Follicular Fluid in Mares? Animals (Basel) 2020, 10, 646. [Google Scholar] [CrossRef] [Green Version]
  50. Ling, Y.H.; Guo, X.F.; Chen, T.; Ding, J.P.; Ma, Y.H.; Chu, M.X.; Di, R.; Zhang, Y.H.; Zhang, X.R. Characterization and analysis of differentially expressed microRNAs in hircine ovaries during the follicular and luteal phases. Anim. Reprod. Sci. 2016, 166, 47–57. [Google Scholar] [CrossRef]
  51. Scarlet, D.; Walter, I.; Hlavaty, J.; Aurich, C. Expression and immunolocalisation of follicle-stimulating hormone receptors in gonads of newborn and adult female horses. Reprod. Fertil. Dev. 2016, 28, 1340–1348. [Google Scholar] [CrossRef]
  52. Dierich, A.; Sairam, M.R.; Monaco, L.; Fimia, G.M.; Gansmuller, A.; LeMeur, M.; Sassone-Corsi, P. Impairing follicle-stimulating hormone (FSH) signaling in vivo: Targeted disruption of the FSH receptor leads to aberrant gametogenesis and hormonal imbalance. Proc. Natl. Acad. Sci. USA 1998, 95, 13612–13617. [Google Scholar] [CrossRef] [Green Version]
  53. Michalovic, L.; Currin, L.; Gutierrez, K.; Bellefleur, A.M.; Glanzner, W.G.; Schuermann, Y.; de Macedo, M.P.; Bohrer, R.C.; Dicks, N.; Lopez, R.; et al. Granulosa cells of prepubertal cattle respond to gonadotropin signaling and upregulate genes that promote follicular growth and prevent cell apoptosis. Mol. Reprod. Dev. 2018, 85, 909–920. [Google Scholar] [CrossRef]
  54. Lubahn, D.B.; Moyer, J.S.; Golding, T.S.; Couse, J.F.; Korach, K.S.; Smithies, O. Alteration of reproductive function but not prenatal sexual development after insertional disruption of the mouse estrogen receptor gene. Proc. Natl. Acad. Sci. USA 1993, 90, 11162–11166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Couse, J.F.; Korach, K.S. Estrogen receptor null mice: What have we learned and where will they lead us? Endocr. Rev. 1999, 20, 358–417. [Google Scholar] [CrossRef] [PubMed]
  56. Liu, W.; Xin, Q.; Wang, X.; Wang, S.; Wang, H.; Zhang, W.; Yang, Y.; Zhang, Y.; Zhang, Z.; Wang, C.; et al. Estrogen receptors in granulosa cells govern meiotic resumption of pre-ovulatory oocytes in mammals. Cell Death Dis. 2017, 8, e2662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Artini, P.G.; Tatone, C.; Sperduti, S.; D’Aurora, M.; Franchi, S.; Di Emidio, G.; Ciriminna, R.; Vento, M.; Di Pietro, C.; Stuppia, L.; et al. Cumulus cells surrounding oocytes with high developmental competence exhibit down-regulation of phosphoinositol 1,3 kinase/protein kinase B (PI3K/AKT) signalling genes involved in proliferation and survival. Hum. Reprod. 2017, 32, 2474–2484. [Google Scholar] [CrossRef] [PubMed]
  58. Palma, G.A.; Tortonese, D.J.; Sinowatz, F. Developmental capacity in vitro of prepubertal oocytes. Anat. Histol. Embryol. 2001, 30, 295–300. [Google Scholar] [CrossRef]
  59. Zhang, X.; Jafari, N.; Barnes, R.B.; Confino, E.; Milad, M.; Kazer, R.R. Studies of gene expression in human cumulus cells indicate pentraxin 3 as a possible marker for oocyte quality. Fertil. Steril. 2005, 83 (Suppl. 1), 1169–1179. [Google Scholar] [CrossRef]
  60. Kona, S.S.; Praveen, C.V.; Siva, K.A.; Srividya, D.; Padmaja, K.; Rao, V.H. Quantitative expression patterns of GDF9 and BMP15 genes in sheep ovarian follicles grown in vivo or cultured in vitro. Theriogenology 2016, 85, 315–322. [Google Scholar] [CrossRef]
  61. Wong, A.C.; Ryan, A.F. Mechanisms of sensorineural cell damage, death and survival in the cochlea. Front. Aging Neurosci. 2015, 7, 58. [Google Scholar] [CrossRef] [PubMed]
  62. Boone, D.L.; Tsang, B.K. Caspase-3 in the rat ovary: Localization and possible role in follicular atresia and luteal regression. Biol. Reprod. 1998, 58, 1533–1539. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, T.H.; Chang, C.L.; Wu, H.M.; Chiu, Y.M.; Chen, C.K.; Wang, H.S. Insulin-like growth factor-II (IGF-II), IGF-binding protein-3 (IGFBP-3), and IGFBP-4 in follicular fluid are associated with oocyte maturation and embryo development. Fertil. Steril. 2006, 86, 1392–1401. [Google Scholar] [CrossRef] [PubMed]
  64. Gu, B.; Liu, H.; Han, Y.; Chen, Y.; Jiang, H. Integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old Small Tail Han Sheep ovaries reveals that oar-miR-432 downregulates RPS6KA1 expression. Gene 2019, 710, 76–90. [Google Scholar] [CrossRef] [PubMed]
  65. Christenson, L.K. MicroRNA control of ovarian function. Anim. Reprod. 2010, 7, 129–133. [Google Scholar]
  66. Ling, Y.H.; Ren, C.H.; Guo, X.F.; Xu, L.N.; Huang, Y.F.; Luo, J.C.; Zhang, Y.H.; Zhang, X.R.; Zhang, Z.J. Identification and characterization of microRNAs in the ovaries of multiple and uniparous goats (Capra hircus) during follicular phase. BMC Genom. 2014, 15, 339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Liu, Y.; Li, M.; Bo, X.; Li, T.; Ma, L.; Zhai, T.; Huang, T. Systematic Analysis of Long Non-Coding RNAs and mRNAs in the Ovaries of Duroc Pigs During Different Follicular Stages Using RNA Sequencing. Int. J. Mol. Sci. 2018, 19, 1722. [Google Scholar] [CrossRef] [Green Version]
  68. Li, J.; Li, C.; Li, Q.; Li, W.T.; Li, H.; Li, G.X.; Kang, X.T.; Liu, X.J.; Tian, Y.D. Identification of the Key microRNAs and miRNA-mRNA Interaction Networks during the Ovarian Development of Hens. Animals (Basel) 2020, 10, 1680. [Google Scholar] [CrossRef]
  69. Friedman, R.C.; Farh, K.K.H.; Burge, C.B.; Bartel, D.P. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2008, 19, 92–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Lin, J.; Wu, Y.; Han, B.; Chen, Y.; Wang, L.; Li, X.; Liu, M.; Huang, J. iTRAQ-based proteomic profiling of granulosa cells from lamb and ewe after superstimulation. Theriogenology 2017, 101, 99–108. [Google Scholar] [CrossRef]
Figure 1. miRNA expression analysis of the ovaries of lamb and adult sheep under superovulation treatment. (A) Sequence length distribution for miRNA sequencing. (B) Principal component analysis (PCA) of expression of miRNA values in the lamb and adult sheep libraries. (C) Screen of differentially expressed miRNAs. (D) Heat map analysis of differentially expressed miRNAs.
Figure 1. miRNA expression analysis of the ovaries of lamb and adult sheep under superovulation treatment. (A) Sequence length distribution for miRNA sequencing. (B) Principal component analysis (PCA) of expression of miRNA values in the lamb and adult sheep libraries. (C) Screen of differentially expressed miRNAs. (D) Heat map analysis of differentially expressed miRNAs.
Animals 11 00239 g001
Figure 2. mRNA expression analysis for the ovaries of lambs and adult sheep under superovulation treatment. (A) Pearson correlation for mRNA expression values across samples. (B) PCA of mRNA expression values of the lamb and adult libraries. (C) Volcano plot of distribution for differentially expressed mRNAs. (D) Heat map analysis of differentially expressed mRNAs.
Figure 2. mRNA expression analysis for the ovaries of lambs and adult sheep under superovulation treatment. (A) Pearson correlation for mRNA expression values across samples. (B) PCA of mRNA expression values of the lamb and adult libraries. (C) Volcano plot of distribution for differentially expressed mRNAs. (D) Heat map analysis of differentially expressed mRNAs.
Animals 11 00239 g002
Figure 3. miRNA–mRNA expression regulatory network related to genes associated with the hormone receptor, steroid hormone secretion, and oocyte quality. Molecules in red and green denote upregulation and downregulation in the lamb libraries, respectively.
Figure 3. miRNA–mRNA expression regulatory network related to genes associated with the hormone receptor, steroid hormone secretion, and oocyte quality. Molecules in red and green denote upregulation and downregulation in the lamb libraries, respectively.
Animals 11 00239 g003
Figure 4. Real-time PCR validation of mRNA associated with hormone receptor and oocyte quality. * Significant difference at p < 0.05, ** significant difference at p < 0.01.
Figure 4. Real-time PCR validation of mRNA associated with hormone receptor and oocyte quality. * Significant difference at p < 0.05, ** significant difference at p < 0.01.
Animals 11 00239 g004
Table 1. Measurement of follicle numbers and serum hormone levels in lamb and adult sheep.
Table 1. Measurement of follicle numbers and serum hormone levels in lamb and adult sheep.
GroupNumber of AnimalsAverage Number of FolliclesFSH (mIU/mL)LH (mIU/mL)P4 (10−2 ng/mL)E2 (pg/mL)
Lamb670.17 ± 5.14 a3.65 ± 0.353.66 ± 0.371.36 ± 0.244.43 ± 0.37
Adult620.17 ± 2.65 b3.61 ± 0.323.88 ± 0.101.33 ± 0.114.67 ± 0.30
Note: Values with different superscript letters in the same column represent highly significant difference (p < 0.01). Values with no superscript letters in the same column have no significant difference (p > 0.05).
Table 2. Top 10 highly and differentially expressed miRNAs in the lamb and adult sheep libraries.
Table 2. Top 10 highly and differentially expressed miRNAs in the lamb and adult sheep libraries.
miRNALamb_MeanAdult_Meanlog2(fc)p-ValueFDR
miR-145-x7649.9223730.58–1.633.43 × 10−251.04 × 10−22
miR-144-y6678.34286.184.543.79 × 10−99.32 × 10−8
miR-451-x4928.68578.773.094.28 × 10−43.33 × 10−3
miR-450-x3914.41721.042.446.42 × 10−69.05 × 10−5
miR-486-x3379.14327.353.371.22 × 10−41.11 × 10−3
miR-424-x2801.34448.962.643.14 × 10−75.34 × 10−6
miR-202-x2696.87316.413.091.93 × 10−117.79 × 10−10
oar-miR-10a1594.464637.33–1.542.30 × 10−119.09 × 10−10
oar-miR-3958-3p1497.39392.701.934.71 × 10−21.65 × 10−1
oar-miR-409-3p1485.38197.242.915.88 × 10−112.02 × 10−9
Note: miRNAs with a red background represent upregulation in the lamb libraries compared to the adult sheep libraries (false discovery rate (FDR) < 0.05, log2(fc) > 1); miRNAs with a green background represent downregulation in the lamb libraries (FDR < 0.05, log2(fc) < –1).
Table 3. mRNAs associated with hormone receptors, steroid hormone secretion and oocyte quality were differentially expressed in the lamb and adult sheep libraries.
Table 3. mRNAs associated with hormone receptors, steroid hormone secretion and oocyte quality were differentially expressed in the lamb and adult sheep libraries.
GeneLamb_MeanAdult_Meanlog2(fc)FDRGene Classification
FSHR5.931.491.994.63 × 10−7HC
LHCGR84.1016.312.371.77 × 10−29HC
ESR18.8922.16−1.323.77 × 10−13HC
CYP11A1251.2941.332.603.66 × 10−42SHE
CYP17A174.874.644.011.71 × 10−15SHE
CYP19A150.231.375.202.46 × 10−14SHE
STAR32.5611.311.535.09 × 10−13SHE
HSD3B112.8439.141.536.76 × 10−3SHE
HSD17B187.116.593.731.15 × 10−7SHE
PTX34.0220.19−2.332.75 × 10−40ROQ
GDF94.071.331.616.91 × 10−4ROQ
DNMT128.3312.301.202.21 × 10−39ROQ
CASP37.191.582.182.75 × 10−8ROQ
IGFBP444.0174.65−0.762.71 × 10−6ROQ
BCL23.285.38−0.718.84 × 10−3ROQ
BAX16.8413.980.272.60 × 10−1ROQ
BAX/BCL25.132.600.987.75 × 10−3NA
Note: mRNAs with a red background represent upregulation in the lamb libraries compared to the adult sheep libraries (FDR < 0.05, log2 (fc) > 1); mRNA with green background represent downregulation in lamb libraries (FDR < 0.05, log2 (fc) < – 1). HC—represents the gene classified into the hormone receptor gene; SHE—represents the gene classified into the steroid hormone secretion gene; ROQ—represents the gene is classified into genes related to oocyte quality; NA—represents not applicable.
Table 4. List of differentially expressed miRNAs targeting key genes.
Table 4. List of differentially expressed miRNAs targeting key genes.
idLamb_Mean Adult_Meanlog2(fc)FDRTarget GeneTarget Gene Classification
oar-miR-143 [34]70,716.65 122,228.40 −0.791.98 × 10−11FSHRHC
miR-513-x [35]1.09 2.41 −1.142.96 × 10−4LHCGRHC
miR-4324-y2.84 10.16 −1.841.20 × 10−16CYP11A1SHE
miR-57-x0.050.42−3.131.47 × 10−3STARSHE
miR-993-x0.340.73−1.117.74 × 10−2STARSHE
oar-miR-485-3p60.89 5.71 3.41 4.33 × 10−10PTX3ROQ
oar-miR-487a-3p55.39 4.55 3.61 6.36 × 10−11PTX3ROQ
oar-miR-377-3p9.39 1.46 2.69 6.53 × 10−3PTX3ROQ
miR-767-x3.810.433.161.09 × 10−2PTX3ROQ
miR-496-x0.51 0.01 5.67 2.67 × 10−2PTX3ROQ
Note: Molecules with a red background represent upregulation in the lamb libraries compared to the adult sheep libraries (FDR < 0.05, log2 (fc) > 1); molecules with a green background represent downregulation in the lamb libraries (FDR < 0.05, log2 (fc) < –1). HC—the gene is classified into the hormone receptor gene; SHE—the gene is classified into the steroid hormone secretion gene; ROQ—the gene is classified into genes related to oocyte quality.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, X.; Dong, C.; Yang, J.; Li, Y.; Feng, J.; Wang, B.; Zhang, J.; Guo, X. The Roles of the miRNAome and Transcriptome in the Ovine Ovary Reveal Poor Efficiency in Juvenile Superovulation. Animals 2021, 11, 239. https://doi.org/10.3390/ani11010239

AMA Style

Zhang X, Dong C, Yang J, Li Y, Feng J, Wang B, Zhang J, Guo X. The Roles of the miRNAome and Transcriptome in the Ovine Ovary Reveal Poor Efficiency in Juvenile Superovulation. Animals. 2021; 11(1):239. https://doi.org/10.3390/ani11010239

Chicago/Turabian Style

Zhang, Xiaosheng, Chunxiao Dong, Jing Yang, Yihai Li, Jing Feng, Biao Wang, Jinlong Zhang, and Xiaofei Guo. 2021. "The Roles of the miRNAome and Transcriptome in the Ovine Ovary Reveal Poor Efficiency in Juvenile Superovulation" Animals 11, no. 1: 239. https://doi.org/10.3390/ani11010239

APA Style

Zhang, X., Dong, C., Yang, J., Li, Y., Feng, J., Wang, B., Zhang, J., & Guo, X. (2021). The Roles of the miRNAome and Transcriptome in the Ovine Ovary Reveal Poor Efficiency in Juvenile Superovulation. Animals, 11(1), 239. https://doi.org/10.3390/ani11010239

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop