- Research article
- Open Access
Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene
BMC Genomics volume 22, Article number: 217 (2021)
Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays.
A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3’UTR at a unique binding site.
Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.
MicroRNAs (miRNAs) belong to a large family of endogenous noncoding RNAs (ncRNAs). MiRNAs may regulate the expression of target genes by binding to complementary regions in their 3′ untranslated regions (3’UTRs) [1, 2]. Many studies have shown that miRNAs play important regulatory roles in animal reproduction [3,4,5,6,7,8,9]. In recent years, miRNAs have also been found to be involved in the regulation of animal seasonal reproduction [10,11,12,13]. The pineal gland is a key organ known to transduce the photoperiod signal to the reproductive axis and is therefore crucial for seasonal reproduction [14, 15]. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (nonbreeding and breeding seasons). Sheep are regarded as a typical animal model for seasonal reproduction studies [16,17,18,19,20]. Therefore, one objective of this study was to reveal the expression characteristics of pineal miRNAs in sheep during different reproductive stages. The second purpose of this study was to screen differentially expressed miRNAs among the different reproductive stages and predict their probable functions in the pineal gland.
The joint analysis of miRNAs and gene expression profiles has become feasible in recent years [21,22,23,24], and miRNA-target gene pairs can be predicted. Furthermore, the specific target sites at which miRNAs regulate gene expression can also be revealed by a double luciferase reporter assay [25, 26]. The pineal gland participates in the regulation of seasonal reproduction mainly through melatonin [27,28,29,30]. Known rate-limiting enzymes for melatonin synthesis include aralkylamine N-acetyltransferase (AANAT) and hydroxyindole O-methyltransferase (HIOMT). Thus, the third objective of this study was to identify the miRNAs targeting the genes encoding the above-mentioned rate-limiting enzyme that were differentially expressed during the different reproductive stages. Together, this study provides new insights into the expression characteristics and roles of pineal miRNAs in sheep during different reproductive stages.
Expression characteristics of pineal small RNAs in sheep at different reproductive stages
Pineal tissues were collected at different reproductive stages, and the accuracy of the sampling stage determination was confirmed by the results obtained from ovarian sections (Fig. 1). To reveal the expression patterns of pineal small RNAs in different ovine reproductive stages, high-throughput deep sequencing was used to obtain the pineal sRNA expression profiles at each stage. Generally, the clean reads in the pineal gland showed a three-peaked length distribution, with peaks at 18, 22 and 32 nucleotides (nt) (Fig. 2a). The majority of expressed reads from the anestrus sample were 18 nt long, whereas reads from samples during the breeding season were predominantly 22 or 32 nt long (Additional file 1). To further assess the efficiency of high-throughput sequencing in detecting miRNA, all the mapped clean reads (clean read counts and mapping ratios are listed in Additional file 2) were annotated and classified by alignment against ncRNAs. Among the diverse sequences of ovine pineal ncRNAs (including miRNAs, rRNAs, sRNAs, tRNAs and other RfamRNAs), the proportion of miRNAs was always the highest in each stage (Fig. 2b), and their values were also similar among different stages. However, the proportions of rRNAs, sRNA and other RfamRNA were relatively higher during anestrus than during the breeding season. In total, 427 miRNAs were identified in the ovine pineal gland, including 91 known miRNAs, 311 conserved miRNAs and 25 predicted novel miRNAs (Fig. 2c). Compared with the two stages (luteal and follicular phases) in the breeding season, expressed miRNAs (including known, conserved and novel miRNAs) were the least abundant in ovine anestrus (Fig. 2d).
Next, the functions of miRNAs that were specifically expressed in anestrus or the breeding season were predicted. KEGG pathway analysis (Additional file 3) showed that the target genes of miRNAs that are expressed specifically in anestrus were predominantly enriched in endocytosis, mTOR and MAPK signaling pathways. These pathways are mainly associated with hormone uptake, protein synthesis, and cell proliferation and differentiation. On the other hand, the target genes of miRNAs that were expressed specifically in the breeding season were predominantly involved in pathways such as the mTOR signaling pathway, apoptosis and axon guidance (Additional file 3). These pathways are mainly associated with protein synthesis, cell growth and death, and the formation of neuronal networks.
Meanwhile, the expression of miRNAs was also ranked in each reproductive stage, and the 20 most highly expressed miRNAs are displayed in Table 1. The results indicated that the top 20 miRNAs were similar between the two stages (luteal and follicular phases) in the breeding season; however, they were significantly different between the breeding season and anestrus. In anestrus, miR-142 (homology ID: aca-miR-5441) was the most abundant miRNA, accounting for 86% of the total expressed miRNA. KEGG analysis showed that the target genes of miR-142 were predominantly enriched in oxidative phosphorylation, glycerolipid metabolism and phosphatidylinositol signaling pathways. In addition to miR-142, high expression of miR-202 (homology ID: tae-miR-5086) and miR-2 (homology ID: cel-miR-51-5p) was also observed during anestrus. Oar-miR-181a, Oar-miR-26a and Oar-miR-143 showed the highest levels of expression in the breeding season. Additionally, Oar-let-7a was highly expressed in all reproductive stages.
Differentially expressed (DE) miRNAs among different ovine reproductive stages and their probable functions in the pineal gland
We determined the DE miRNAs among three reproductive stages (anestrus, follicular phase and luteal phase). The largest number of DE miRNAs was detected between anestrus and the follicular phase (Fig. 3a). Hierarchical clustering of miRNAs (Fig. 3b) also indicated that the differences in miRNA expression between anestrus and the follicular phase were greatest among the three stages analyzed. Furthermore, the majority of the DE miRNAs between anestrus and the two stages of the breeding season overlapped (Fig. 3c). Therefore, these overlapping miRNAs could be considered DE miRNAs between anestrus and the breeding season. To determine the probable biological functions of the overlapping miRNAs, we performed a pathway analysis of the target genes of these miRNAs. The analysis revealed that these miRNAs were mainly enriched in pathways related to protein biosynthesis, secretion and uptake (such as biosynthesis of amino acids, ribosome, cAMP signaling pathway, vascular smooth muscle contraction, axon guidance, dopaminergic synapses, and endocytosis pathway) and the phototransduction pathway (P < 0.01) (Fig. 3d, Additional file 4). Moreover, the results of the transcriptome analysis (Additional file 5) showed that the majority of the target genes in these pathways exhibited differential expression between the seasons. For example, RPLP1, RPLP2, RPL18A, RPL35, RPS5, RPS13 and RPSA in the ribosome pathway showed significantly lower expression levels in the breeding season than in anestrus (Additional file 6).
In addition, the overlapping differentially expressed genes between anestrus and the luteal phase and between anestrus and the follicular phase were also screened out to represent the expression differences in genes between nonbreeding and breeding seasons. The highly expressed genes during anestrus and the breeding season in the pineal gland of sheep are shown in Additional file 7. Some of the highly expressed genes in anestrus were related to protein synthesis and hormone secretion. Highly expressed genes in the breeding season were involved in the ribosome pathway, cAMP signaling pathway and other pathways.
Prediction of important miRNA–target gene pairs
The joint analysis of negatively correlated miRNA–mRNA pairs and miRNA target binding prediction has been demonstrated to be a feasible approach for predicting miRNA-target gene pairs [31, 32]. We therefore measured pineal mRNA profiles at different reproductive stages to examine miRNA–mRNA correlations at the whole-genome scale. Among the negatively correlated pairs, many miRNA-target gene pairs with binding sites were predicted. We first investigated the transcriptional regulatory role of miRNAs on key rate-limiting enzyme genes in melatonin synthesis. The expression of AANAT mRNA showed significant variation at different reproductive stages (Fig. 4). Therefore, the miRNAs that were significantly and negatively correlated with the AANAT expression pattern were predicted and summarized in Table 2. To validate the results from high-throughput sequencing, real-time quantitative PCR was performed for the five miRNAs in Table 2 and the AANAT gene. The results of quantitative PCR (Fig. 4) were consistent with the sequencing data, and these miRNAs exhibited an inverted expression to that of the AANAT gene. Additionally, for the expression of HIOMT mRNA, there was no significant difference among the reproductive stages. Therefore, miRNAs targeting the gene were not predicted in this study. In addition to AANAT, miRNAs potentially targeting several differentially expressed genes in the ribosome pathway were also predicted (Additional file 8).
Validation of the target relationships between AANAT and candidate miRNAs
To further verify the target relationships between AANAT and candidate miRNAs, one miRNA (miR-89) that was predicted to show a target relationship (completely binding) and one (miR-201) that was predicted to show just partial binding but exhibited a high negative correlation coefficient were selected for validation in vitro. Among the miRNAs with a complete binding relationship, miR-89 was selected for validation because it was previously detected in sheep (defined as Oar-miR-214-3p) and suggested to be involved in the regulation of breeding activities . Next, the 3′ UTR sequence of the AANAT gene in Tan sheep was obtained, and the sequence was consistent with the corresponding region of the reported AANAT mRNA sequence (NM_001009461.1) in a cross of the Dorset × Rambouillet sheep [34, 35]. The predicted target site of miR-89 in the 3’UTR of the AANAT gene is shown in Fig. 5a. The possible binding site of miR-201 with the 3’UTR of AANAT was also predicted by RNAhybrid software. The 2nd and 8th bases at the 5′ end of the miRNA showed semibinding states with the 3’UTR of AANAT (Fig. 5b). Then, the dual luciferase reporter system was used to further verify the targeting role of the miRNAs associated with the 3’UTR of AANAT in vitro. The miRNAs were successfully transfected into 293FT cells, and the efficiency of transfection is indicated in Fig. 6. As shown in Fig. 5c, the relative luciferase activity of luc2/hRluc in the group cotransfected with miR-89 and the wild-type 3′ UTR of AANAT was significantly lower (P < 0.01) than in the group transfected with the wild-type 3′ UTR of AANAT alone or the group cotransfected with miR-89 and the mutant-type 3′ UTR of AANAT. This result verified that miR-89 can target the 3′ UTR of AANAT and that the binding site is unique. Taken together, the results of validation of the target relationship and the observed expression of miRNA and AANAT mRNA during different stages implied that miR-89 participates in the negative regulation of AANAT mRNA expression by targeting its 3′ UTR. However, the data demonstrated that miR-201 had no apparent target effect on the 3′ UTR of AANAT (Fig. 5d). Therefore, the results of in vitro analysis demonstrated that the predicted miRNA–target mRNA pairs were accurate.
Many recent studies have shown that miRNAs play important roles in the regulation of reproduction [3,4,5,6,7,8,9], including seasonal reproduction [10,11,12,13]. The pineal gland is a crucial organ for seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons. In this study, we investigated the miRNA expression profiles of the ovine pineal gland for the first time.
First, a higher number of expressed miRNAs were detected in the pineal gland of sheep compared with those identified miRNAs in zebrafish  and rats [36, 37]. Second, the results suggested that there are some differences in the expression characteristics of miRNAs between reproductive stages. Our previous and current results have demonstrated that the number of expressed miRNAs is lowest during anestrus in both the ovary  and pineal gland of sheep. Additionally, the most highly expressed miRNAs were also distinct among reproductive stages. Specifically, miR-142 was detected as the most highly expressed miRNA during anestrus. Similarly, high expression of miR-142 was also observed in the gonad during ovine anestrus . Previous studies revealed the roles of miR-142 in suppressing proliferation, inducing apoptosis [38, 39] and regulating neurotransmission . In this study, the synaptic vesicle cycle was one of the enrichment pathways for miR-142 target genes, which is related to neurotransmission. In addition, the enriched pathways (oxidative phosphorylation, glycerolipid metabolism and phosphatidylinositol signaling system) for the target genes of miR-142 were mainly associated with energy and lipid metabolism. Previous studies observed that the volume and number of mitochondria and lipid droplets in pinealocytes of seasonal breeders varied between short and long photoperiods [41, 42], implying the existence of differences in energy and lipid metabolism in pinealocytes between different seasons for seasonal breeders. In addition to miR-142, miR-202 is the second ranked miRNA that is expressed in anestrus. Several studies have shown, such as miR-142, that miR-202 is involved in cell proliferation [43,44,45]. For miR-26, which was highly expressed in the pineal gland of the breeding season, previous studies have suggested that it may play a role in the transition of the reproductive stages in goats , and its high expression was also detected in brain tissues of beef cattle . For Oar-let-7a, which was highly expressed in each reproductive stage, several studies have shown that it is crucial for reproduction in the ovaries of sheep during different reproductive stages [10, 48]. Specifically, let-7 is involved in regulating important aspects of reproduction, including steroidogenesis  and follicular development . In addition, miR-7 was listed in the top 20 expressed genes of the three stages simultaneously. miR-7 has been shown to regulate melatonin synthesis by acting as a linking molecule between leptin and norepinephrine signaling pathways in pig pineal .
A clear difference in miRNA expression existed between anestrus and the breeding season, as revealed in Fig. 3, implying distinct roles for pineal miRNAs between the seasons. KEGG analysis indicated that the target genes of both miRNAs that were expressed specifically during anestrus or the breeding season and the DE miRNAs between anestrus and the breeding season were highly enriched in pathways that were related to protein synthesis, secretion and uptake. Furthermore, based on the results of transcriptome analysis, many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season, which might lead to a change in protein synthesis. Previous studies indicated that the morphology of Golgi bodies and endoplasmic reticulum, which represent the ultrastructural characteristics of protein synthesis, showed seasonal-specific changes in the pinealocytes of seasonal breeders. For instance, when exposed to an artificially long photoperiod, rams exhibited some morphological changes in the endoplasmic reticulum and microtubules, and the volumes of mitochondria and lipid droplets in pinealocytes increased compared to rams under naturally short photoperiods . Our results also suggest that protein synthesis in pinealocytes is more active under long photoperiods. However, for long-day breeders, the cisternae of the endoplasmic reticulum and Golgi bodies tend to dilate under short photoperiods in ground squirrels , and the Golgi body-endoplasmic reticulum-lysosome complex will be clearer under short photoperiods in Parakeet . These studies, together with ours, suggest that protein biosynthesis in pinealocytes has seasonal features in seasonal breeders, and miRNAs might participate in the regulation of this variation.
Additionally, our study demonstrated that AANAT mRNA expression varied and may be regulated by miRNAs during different stages of reproduction. To date, few reports have addressed the differences in AANAT gene expression in the pineal glands of seasonal breeders among different reproductive stages. In the estrus cycle, it was found that the expression levels of AANAT in oocytes were higher in the follicular phase than in the luteal phase in rat ovaries , which is similar to the trend of AANAT variation in the pineal gland observed in this study. Analyses that have combined miRNA-gene sequencing data with target relationship validation in vitro have become popular in recent years [54, 55]. Using such analysis, it was found that miR-483 could negatively regulate AANAT expression in the rat pineal gland during development  and that the upregulation of miR-325-3p suppressed pineal AANAT expression after neonatal hypoxia–ischemia brain injury in rats . In this study, the joint analysis of the expression of AANAT and miRNA, along with validation testing in vitro, suggested that miR-89 may negatively regulate AANAT mRNA expression at different reproductive stages by targeting its 3’UTR. Meanwhile, the miRNA without a predicted complete target relationship with the AANAT gene was proven to be unable to target the 3’UTR of the gene.
The expression characteristics of pineal miRNAs in different stages of reproduction were revealed for the first time using domestic sheep. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season. KEGG analysis of differentially expressed miRNAs between anestrus and the breeding season indicated that these miRNAs are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, the transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway were expressed at relatively low levels in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of AANAT mRNA expression by targeting its 3’UTR at a unique binding site.
Sample collection and deep sequencing
The Tan sheep in northern China are a typical seasonal breeder with a long anestrus period, which usually extends from April to July . Therefore, the Tan sheep used in this study were selected from the Yanchi Tan sheep conservation farm in the Ningxia Autonomous Region, China. First, nine healthy and nonpregnant Tan ewes (three years old, body weight: 45–50 kg) were examined daily for estrous activity with a teaser ram during anestrus and the breeding season of 2016. Anestrus was judged to be a period without obvious signs of estrus for more than 36 days (equivalent to the period of two estrous cycles). On the 50th day during anestrus, three ewes (anestrus, n = 3) were simultaneously euthanized in a separate room. Before euthanasia, the overdose of pentobarbital (a barbiturate which depresses the central nervous system causing respiratory and cardiac arrest) and phenytoin (an antiepileptic drug which further decreases the central nervous system via sodium and calcium channel blockade) were administered at 0.22 mL/kg (85.8 mg/kg of pentobarbital and 11 mg/kg of phenytoin) over 10 s referring to the previous report . Within 15 min after death, the entire pineal gland of ewes was removed from the brain and immediately collected in liquid nitrogen. Meanwhile, their ovaries were placed in 4% paraformaldehyde fixative to prepare paraffin sections in order to verify the accuracy of the sampling stages. Ovarian paraffin sections were stained with hematoxylin-eosin. In the breeding season, an estrous cycle includes the luteal and follicular phases. The date of estrus was recorded as the 1st day in each estrous cycle. Two continuous estrous cycles were recorded. The mean estrus cycle among six Tan ewes was 17.33 ± 0.26 days. Therefore, in the third estrous cycle, ewes were euthanized using the same method as above on the 9th day (luteal phase, n = 3) and during estrus (follicular phase, n = 3), and their pineal glands and ovaries were collected immediately as described above. Total RNA extraction and the construction of small RNA and mRNA libraries were performed using the same method as that in the literature . The Illumina HiSeq 2000 system was used to sequence the library products. The sequence datasets have been deposited in the BioProject (Biological Project Library) database of Genome Sequence Archive  in National Genomics Data Center , Beijing Institute of Genomics (China National Center for Bioinformation), Chinese Academy of Sciences, with the accession number PRJCA000331.
Bioinformatic analysis for small RNAs and mRNA
For miRNA analysis, clean reads of 17 ~ 36 bp were obtained using a method we described in a previous study . Then, all clean reads were collapsed into tags based on sequence identity using an in-house built pipeline . The unique tags were aligned to the Ovine reference genome (v3.1) by the mapper.pl of miRDeep2 (version 188.8.131.52) (no nucleotide changes; genome hits were less than 5), and then the miRNAs were predicted with the mirDeep2.pl program . By comparison with the known ncRNAs (rRNAs, tRNAs, snRNAs, and snoRNA) from the NCBI GenBank and Rfam10.1 database, all the mapped tags were classed and annotated. In addition, using miRBase v21 [62, 63], all predicted miRNAs were classified as novel, known and conserved miRNAs using a method we previously described . Finally, the expression levels of miRNAs and their difference between different reproductive stages were analyzed with the same method as in our previous study . For the mRNA, quality control of the raw data, mapping of high-quality reads and expression analysis of genes were performed using a method we previously described . Genes with an FDR < 0.05 were defined as differentially expressed genes .
Based on the 3′-UTR sequences that were extracted according to gene annotation information, miRNAs potentially targeting candidate genes were simultaneously predicted using the miRanda and RNAhybrid algorithms. Pairs of DE miRNAs and genes showing negative correlations were identified using R software. The function of target genes was predicted by KEGG pathway analysis. The enrichment p-values were calculated using a hypergeometric distribution.
Validation of the expression of AANAT and its associated miRNAs
Real-time quantitative PCR (Q-PCR) was used to quantify the expression levels of AANAT mRNA and miRNAs that were randomly selected among those that were negatively correlated with AANAT in sequencing data. The stem-loop primers employed in RT-PCR of miRNAs and the primers employed in Q-PCR of AANAT and miRNAs are summarized in Additional file 9. Reverse transcription and Q-PCR were performed with the same method as in our previous study . The data were analyzed with the 2–ΔΔCt method .
Validation of the target relationship between candidate miRNAs and AANAT
The 3′ UTR sequence of the AANAT gene of Tan sheep was obtained using the 3′-Full RACE Core Set with PrimeScript™ RTase kit (TaKaRa) and the specific primers GSP1 (5′-AGCGTCCACTGCCTGAAACC-3′) and GSP2 (5′-AGGTTTGGCTTCCATCCCG-3′). Following sequencing, the potential target site bound by candidate miRNAs was predicted using miRanda and RNAhybrid algorithms. For candidate miRNA-gene pairs, both a miRNA expression vector (pcDNA6.2-GW/EmGFP-miR) and a reporter vector (pmirGLO) in which the AANAT 3’UTR was fused to a firefly luciferase reporter were obtained to perform the dual luciferase reporter assay. The sequence that was inserted in the pcDNA6.2-GW/EmGFP-miR vector was designed based on the mature sequence of the miRNA (Additional file 10) and was synthesized by Invitrogen Trading (Shanghai) Company. The AANAT 3’UTR sequence that was inserted into the pmirGLO vector was synthesized by Shanghai Generay Biotech Company. The two vectors were cotransfected into 293FT cells using Lipofectamine 3000. The luciferase activity of luc2/hRluc in transfected cells will be significantly reduced when the miRNA targets the 3’UTR of the candidate gene. To determine whether there were other target sites besides the predicted one, a mutant-type 3’UTR of AANAT mutated at the predicted target site was also designed.
Availability of data and materials
The raw datasets generated during the current study are available in the BioProject (Biological Project Library) database at National Genomics Data Center, Beijing Institute of Genomics (China National Center for Bioinformation) with the accession number PRJCA000331.
Kim D, Sung YM, Park J, Kim S, Kim J, Park J, Ha H, Bae JY, Kim S, Baek D. General rules for functional microRNA targeting. Nat Genet. 2016;48(12):1517–26. https://doi.org/10.1038/ng.3694.
Hansen TB, Veno MT, Jensen TI, Schaefer A, Damgaard CK, Kjems J. Argonaute-associated short introns are a novel class of gene regulators. Nat Commun. 2016;7(1):11538. https://doi.org/10.1038/ncomms11538.
Zhang Y, Zhao B, Roy S, Saha TT, Kokoza VA, Li M, Raikhel AS. microRNA-309 targets the Homeobox gene SIX4 and controls ovarian development in the mosquito Aedes aegypti. Proc Natl Acad Sci U S A. 2016;113(33):E4828–36. https://doi.org/10.1073/pnas.1609792113.
Inyawilert W, Liao YJ, Tang PC. Superovulation at a specific stage of the estrous cycle determines the reproductive performance in mice. Reprod Biol. 2016;16(4):279–86. https://doi.org/10.1016/j.repbio.2016.10.004.
Saul MC, Zhao C, Driessen TM, Eisinger BE, Gammie SC. MicroRNA expression is altered in lateral septum across reproductive stages. Neuroscience. 2016;312:130–40. https://doi.org/10.1016/j.neuroscience.2015.11.019.
Wu N, Gaur U, Zhu Q, Chen B, Xu Z, Zhao X, Yang M, Li D. Expressed microRNA associated with high rate of egg production in chicken ovarian follicles. Anim Genet. 2017;48(2):205–16. https://doi.org/10.1111/age.12516.
Wang Y, Yang HM, Cao W, Li YB, Wang ZY. Deep sequencing identification of miRNAs in pigeon ovaries illuminated with monochromatic light. BMC Genomics. 2018;19(1):446. https://doi.org/10.1186/s12864-018-4831-6.
Wu N, Zhu Q, Chen BL, Gao J, Xu ZX, Li DY. High-throughput sequencing of pituitary and hypothalamic microRNA transcriptome associated with high rate of egg production. BMC Genomics. 2017;18(1):255. https://doi.org/10.1186/s12864-017-3644-3.
Zou X, Lu TT, Zhao ZF, Liu GB, Lian ZQ, Guo YQ, Sun BL, Liu DW, Li YK. Comprehensive analysis of mRNAs and miRNAs in the ovarian follicles of uniparous and multiple goats at estrus phase. BMC Genomics. 2020;21(1):267. https://doi.org/10.1186/s12864-020-6671-4.
Di R, He J, Song S, Tian D, Liu Q, Liang X, Ma Q, Sun M, Wang J, Zhao W, et al. Characterization and comparative profiling of ovarian microRNAs during ovine anestrus and the breeding season. BMC Genomics. 2014;15(1):899. https://doi.org/10.1186/1471-2164-15-899.
Surbhi, Schatz KC, Kyne RF, Nelson RJ, Paul MJ. Photoperiod regulates hypothalamic miR-155 gene expression in female, but not male, Siberian hamsters (Phodopus sungorus). Behav Neurosci. 2019;133(2):240–6. https://doi.org/10.1037/bne0000296.
Yang H, Lin S, Lei XP, Yuan C, Yu YS, Zhao ZS, Chen JB. Nutritional status affects the microRNA profile of the hypothalamus of female sheep. Reprod Fertil Dev. 2018;30(7):946–57. https://doi.org/10.1071/RD17179.
Yang H, Liu XX, Hu GD, Xie YF, Lin S, Zhao ZS, Chen JB. Identification and analysis of microRNAs-mRNAs pairs associated with nutritional status in seasonal sheep. Biochem Biophys Res Commun. 2018;499(2):321–7. https://doi.org/10.1016/j.bbrc.2018.03.155.
Malpaux B, Thiéry JC, Chemineau P. Melatonin and the seasonal control of reproduction. Reprod Nutr Dev. 1999;39(3):355–66. https://doi.org/10.1051/rnd:19990308.
Robinson JE, Karsch FJ. Timing the breeding season of the ewe: what is the role of daylength? Reprod Nutr Dev. 1988;28(2b):365–74. https://doi.org/10.1051/rnd:19880303.
Dardente H, Hazlerigg DG, Ebling FJP. Thyroid hormone and seasonal rhythmicity. Front Endocrinol. 2014;5(26):19.
Dardente H, Wyse CA, Birnie MJ, Dupre SM, Loudon ASI, Lincoln GA, Hazlerigg DG. A molecular switch for photoperiod responsiveness in mammals. Curr Biol. 2010;20(24):2193–8. https://doi.org/10.1016/j.cub.2010.10.048.
Lomet D, Cognie J, Chesneau D, Dubois E, Hazlerigg D, Dardente H. The impact of thyroid hormone in seasonal breeding has a restricted transcriptional signature. Cell Mol Life Sci. 2018;75(5):905–19. https://doi.org/10.1007/s00018-017-2667-x.
Lomet D, Druart X, Hazlerigg D, Beltramo M, Dardente H. Circuit-level analysis identi fies target genes of sex steroids in ewe seasonal breeding. Mol Cell Endocrinol. 2020;512:110825. https://doi.org/10.1016/j.mce.2020.110825.
Wood SH, Christian HC, Miedzinska K, Saer BRC, Johnson M, Paton B, Yu L, McNeilly J, Davis JRE, McNeilly AS, et al. Binary switching of calendar cells in the pituitary defines the phase of the circannual cycle in mammals. Curr Biol. 2015;25(20):2651–62. https://doi.org/10.1016/j.cub.2015.09.014.
Lee SE, Kim SJ, Youn JP, Hwang SY, Park CS, Park YS. MicroRNA and gene expression analysis of melatonin-exposed human breast cancer cell lines indicating involvement of the anticancer effect. J Pineal Res. 2011;51(3):345–52. https://doi.org/10.1111/j.1600-079X.2011.00896.x.
Agarwal S, Nagpure NS, Srivastava P, Kushwaha B, Kumar R, Pandey M, Srivastava S. In silico genome wide mining of conserved and novel miRNAs in the brain and pineal gland of Danio rerio using small RNA sequencing data. Genom Data. 2016;7:46–53. https://doi.org/10.1016/j.gdata.2015.11.013.
Wang L, Huang W, Ren C, Zhao M, Jiang X, Fang X, Xia X. Analysis of serum microRNA profile by Solexa sequencing in women with endometriosis. Reprod Sci. 2016;23(10):1359–70. https://doi.org/10.1177/1933719116641761.
Xie J, Yang X, Song Y, Du Q, Li Y, Chen J, Zhang D. Adaptive evolution and functional innovation of Populus-specific recently evolved microRNAs. New Phytol. 2017;213(1):206–19. https://doi.org/10.1111/nph.14046.
Lennox KA, Vakulskas CA, Behlke MA. Non-nucleotide modification of anti-miRNA oligonucleotides. Methods Mol Biol. 2017;1517:51–69. https://doi.org/10.1007/978-1-4939-6563-2_3.
Li X, Chen W, Zeng W, Wan C, Duan S, Jiang S. microRNA-137 promotes apoptosis in ovarian cancer cells via the regulation of XIAP. Br J Cancer. 2017;116(1):66–76. https://doi.org/10.1038/bjc.2016.379.
Vriend J, Reiter RJ. Melatonin feedback on clock genes: a theory involving the proteasome. J Pineal Res. 2015;58(1):1–11. https://doi.org/10.1111/jpi.12189.
Chesneau D, Guillaume D, Chemineau P, Malpaux B. Continuous light after 2 months of long days stimulates ram testis volume and increases fertility in spring. Animal. 2017;11(7):1189–95. https://doi.org/10.1017/S1751731116002299.
Nishiwaki-Ohkawa T, Yoshimura T. Molecular basis for regulating seasonal reproduction in vertebrates. J Endocrinol. 2016;229(3):R117–27. https://doi.org/10.1530/JOE-16-0066.
Barrett P, Bolborea M. Molecular pathways involved in seasonal body weight and reproductive responses governed by melatonin. J Pineal Res. 2012;52(4):376–88. https://doi.org/10.1111/j.1600-079X.2011.00963.x.
Dini P, Daels P, Loux SC, Esteller-Vico A, Carossino M, Scoggin KE, Ball BA. Kinetics of the chromosome 14 microRNA cluster ortholog and its potential role during placental development in the pregnant mare. BMC Genomics. 2018;19(1):954. https://doi.org/10.1186/s12864-018-5341-2.
Fernandes JG, Borrego A, Jensen JR, Cabrera WHK, Correa MA, Starobinas N, Ribeiro OG, Ibañez OM, De Franco M. miRNA expression and interaction with genes involved in susceptibility to Pristane-induced arthritis. J Immunol Res. 2018;2018:1928405.
Chang W, Wang J, Tao D, Zhang Y, He J, Shi C. Identification of a novel miRNA from the ovine ovary by a combinatorial approach of bioinformatics and experiments. J Vet Med Sci. 2016;77(12):1617–24. https://doi.org/10.1292/jvms.15-0289.
Coon SL, Roseboom PH, Baler R, Weller JL, Namboodiri MA, Koonin EV, Klein DC. Pineal serotonin N-acetyltransferase: expression cloning and molecular analysis. Science. 1995;270(5242):1681–3. https://doi.org/10.1126/science.270.5242.1681.
Herman AP, Bochenek J, Skipor J, Krol K, Krawczynska A, Antushevich H, Pawlina B, Marciniak E, Tomaszewska-Zaremba D. Interleukin-1 beta modulates melatonin secretion in ovine pineal gland: ex vivo study. Biomed Res Int. 2015;2015:526464.
Ding X, Sun B, Huang J, Xu L, Pan J, Fang C, Tao Y, Hu S, Li R, Han X, Miao P, Wang Y, Yu J, Feng X. The role of miR-182 in regulating pineal CLOCK expression after hypoxia-ischemia brain injury in neonatal rats. Neurosci Lett. 2015;591:75–80. https://doi.org/10.1016/j.neulet.2015.02.026.
Clokie SJ, Lau P, Kim HH, Coon SL, Klein DC. MicroRNAs in the pineal gland: miR-483 regulates melatonin synthesis by targeting arylalkylamine N-acetyltransferase. J Biol Chem. 2012;287(30):25312–24. https://doi.org/10.1074/jbc.M112.356733.
Gao YF, Zhang QJ, Yu Z, Liu SH, Liang J. miR-142 suppresses proliferation and induces apoptosis of osteosarcoma cells by upregulating Rb. Oncol Lett. 2018;16(1):733–40. https://doi.org/10.3892/ol.2018.8761.
Pahlavan Y, Mohammadi Nasr M, Dalir Abdolahinia E, Pirdel Z, Razi Soofiyani S, Siahpoush S, Nejati K. Prominent roles of microRNA-142 in cancer. Pathol Res Pract. 2020;216(11):153220.
Chaudhuri AD, Yelamanchili SV, Fox HS. MicroRNA-142 reduces monoamine oxidase a expression and activity in neuronal cells by downregulating SIRT1. PLoS One. 2013;8(11):e79579. https://doi.org/10.1371/journal.pone.0079579.
Lewczuk B, Przybylska B, Udala J, Koncicka A, Wyrzykowski Z. Ultrastructure of the ram pinealocytes during breeding season in natural short and in artificial long photoperiods. Folia Morphol (Warsz). 1993;52(3):133–42.
Prasadan TN, Kotak VC. Fine structure of the free-living parakeet pineal in relation to the breeding cycle. J Pineal Res. 1993;15(3):122–31. https://doi.org/10.1111/j.1600-079X.1993.tb00518.x.
Jiang J, Huang J, Wang XR, Quan YH. MicroRNA-202 induces cell cycle arrest and apoptosis in lung cancer cells through targeting cyclin D1. Eur Rev Med Pharmacol Sci. 2016;20(11):2278–84.
Li YG, He JH, Yu L, Hang ZP, Li W, Shun WH, Huang GX. microRNA-202 suppresses MYCN expression under the control of E2F1 in the neuroblastoma cell line LAN-5. Mol Med Rep. 2014;9(2):541–6. https://doi.org/10.3892/mmr.2013.1845.
Yang J, Fan B, Zhao Y, Fang J. MicroRNA-202 inhibits cell proliferation, migration and invasion of glioma by directly targeting metadherin. Oncol Rep. 2017;38(3):1670–8. https://doi.org/10.3892/or.2017.5815.
Ling YH, Guo XF, Chen T, Ding JP, Ma YH, Chu MX, Di R, Zhang YH, Zhang XR. Characterization and analysis of differentially expressed microRNAs in hircine ovaries during the follicular and luteal phases. Anim Reprod Sci. 2016;166:47–57. https://doi.org/10.1016/j.anireprosci.2016.01.003.
Jin W, Grant JR, Stothard P, Moore SS, Guan LL. Characterization of bovine miRNAs by sequencing and bioinformatics analysis. BMC Mol Biol. 2009;10(1):90. https://doi.org/10.1186/1471-2199-10-90.
McBride D, Carre W, Sontakke SD, Hogg CO, Law A, Donadeu FX, Clinton M. Identification of miRNAs associated with the follicular-luteal transition in the ruminant ovary. Reproduction. 2012;144(2):221–33. https://doi.org/10.1530/REP-12-0025.
Men Y, Fan Y, Shen Y, Lu L, Kallen AN. The Steroidogenic acute regulatory protein (StAR) is regulated by the H19/let-7 Axis. Endocrinology. 2017;158(2):402–9. https://doi.org/10.1210/en.2016-1340.
Cao R, Wu WJ, Zhou XL, Xiao P, Wang Y, Liu HL. Expression and preliminary functional profiling of the let-7 family during porcine ovary follicle atresia. Mol Cells. 2015;38(4):304–11. https://doi.org/10.14348/molcells.2015.2122.
Qiu J, Zhang J, Zhou Y, Li X, Li H, Liu J, Gou K, Zhao J, Cui S. MicroRNA-7 inhibits melatonin synthesis by acting as a linking molecule between leptin and norepinephrine signaling pathways in pig pineal gland. J Pineal Res. 2019;66(3):e12552. https://doi.org/10.1111/jpi.12552.
McNulty JA, Dombrowski TA. Ultrastructural evidence for seasonal changes in pinealocytes of the 13-lined ground squirrel, Spermophilus tridecemlineatus: a qualitative and quantitative study. Anat Rec. 1980;196(4):387–400. https://doi.org/10.1002/ar.1091960404.
Sakaguchi K, Itoh MT, Takahashi N, Tarumi W, Ishizuka B. The rat oocyte synthesises melatonin. Reprod Fertil Dev. 2013;25(4):674–82. https://doi.org/10.1071/RD12091.
Fang YN, Zheng BB, Wang L, Yang W, Wu XM, Xu Q, Guo WW. High-throughput sequencing and degradome analysis reveal altered expression of miRNAs and their targets in a male-sterile cybrid pummelo (Citrus grandis). BMC Genomics. 2016;17(1):591. https://doi.org/10.1186/s12864-016-2882-0.
Lin R, He L, He J, Qin P, Wang Y, Deng Q, Yang X, Li S, Wang S, Wang W, Liu H, Li P, Zheng A. Comprehensive analysis of microRNA-Seq and target mRNAs of rice sheath blight pathogen provides new insights into pathogenic regulatory mechanisms. DNA Res. 2016;23(5):415–25. https://doi.org/10.1093/dnares/dsw024.
Yang Y, Sun B, Huang J, Xu L, Pan J, Fang C, Li M, Li G, Tao Y, Yang X, et al. Up-regulation of miR-325-3p suppresses pineal aralkylamine N-acetyltransferase (Aanat) after neonatal hypoxia-ischemia brain injury in rats. Brain Res. 2017;1668:28–35.
Barletta M, Hofmeister EH, Peroni JF, Thoresen M, Scharf AM, Quandt JE. Influence of sedation on onset and quality of euthanasia in sheep. Res Vet Sci. 2018;117:57–9. https://doi.org/10.1016/j.rvsc.2017.11.012.
Wang YQ, Song FH, Zhu JW, Zhang SS, Yang YD, Chen TT, Tang BX, Dong LL, Ding N, Zhang Q, Bai Z, Dong X, Chen H, Sun M, Zhai S, Sun Y, Yu L, Lan L, Xiao J, Fang X, Lei H, Zhang Z, Zhao W. GSA: genome sequence archive. Genomics Proteomics Bioinformatics. 2017;15(1):14–8. https://doi.org/10.1016/j.gpb.2017.01.001.
Zhang Z, Zhao WM, Xiao JF, Bao YM, He SM, Zhang GQ, Li YX, Zhao GP, Chen RS, Gao Y, et al. Database resources of the National Genomics Data Center in 2020. Nucleic Acids Res. 2020;48(D1):D24–33.
Zhao WM, Liu WF, Tian DM, Tang BX, Wang YQ, Yu CX, Li RJ, Ling YC, Wu JY, Song SH, et al. wapRNA: a web-based application for the processing of RNA sequences. Bioinformatics. 2011;27(21):3076–7. https://doi.org/10.1093/bioinformatics/btr504.
Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52. https://doi.org/10.1093/nar/gkr688.
Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019;47(D1):D155–62. https://doi.org/10.1093/nar/gky1141.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(D1):D68–73. https://doi.org/10.1093/nar/gkt1181.
Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Ser B Stat Methodol. 1995;57(1):289–300.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
We thank JX Wang, XY Pang and SJ Wu for sample collection. We thank M Sun for construction of the RNA libraries.
This work was supported by the National Natural Science Foundation of China (Grant No. 31472078, 31572371, 31402041), the Earmarked Fund for China Agriculture Research System (Grant No. CARS-38), Agricultural Science and Technology Innovation Program of China (Grant No. ASTIP-IAS13), Central Public-interest Scientific Institution Basal Research Fund (Grant No. Y2017JC24), Science and Technology Innovation Guide Program of Ningxia Academy of Agriculture and Forestry Sciences (Grant No. DWJLC-2016001), China Agricultural Scientific Research Outstanding Talents and Their Innovative Teams Program, China High-level Talents Special Support Plan Scientific and Technological Innovation Leading Talents Program (Grant No. W02020274), Tianjin Agricultural Science and Technology Achievements Transformation and Popularization Program (Grant No. 201704020). The funding bodies didn’t play any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All procedures in this study involving animals were carried out in accordance with the guidelines for the care and use of experimental animals established by the Ministry of Science and Technology of China (Approval Number 2006–398), and the work was approved by the Animal Care and Use Committee of Chinese Academy of Agricultural Sciences and Ningxia Academy of Agriculture and Forestry Sciences in China. The written informed consent to participate was obtained from Yanchi Tan sheep conservation farm in the Ningxia Autonomous Region, China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Percentage of reads number with different lengths at three reproductive stages of sheep.
Expression information of small RNA reads in ovine pineal gland.
The signaling pathways which the stage-specific expressed miRNAs were predominantly enriched in.
The genes number and P value of pathways which target genes for differentially expressed miRNAs between anestrus and breeding season were enriched in.
The expression level of target genes for DE miRNAs between anestrus and breeding season in different reproductive stages.
Expression pattern of several genes in the ribosome pathway during different reproductive stages. (EPS 16678 kb)
Highly expressed genes during anestrus or breeding season in pineal gland of sheep.
Predicted miRNAs potentially targeting genes in ribosome pathway.
Stem-loop primers employed for RT-PCR of miRNAs and primers employed for Q-PCR of AANAT and miRNAs.
Mature sequences of miRNAs and inserted sequences in the pcDNA6.2-GW/EmGFP-miR vector.
About this article
Cite this article
Di, R., Liu, QY., Song, SH. et al. Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene. BMC Genomics 22, 217 (2021). https://doi.org/10.1186/s12864-021-07536-y
- Pineal gland
- Breeding season