Characterization and differential expression of microRNAs in the ovaries of pregnant and non-pregnant goats (Capra hircus)
- Xiao-Dong Zhang†1, 2,
- Yun-Hai Zhang†1, 2,
- Ying-Hui Ling1, 2,
- Ya Liu1, 2,
- Hong-Guo Cao1, 2,
- Zong-Jun Yin1, 2,
- Jian-Ping Ding1, 2Email author and
- Xiao-Rong Zhang1, 2Email author
© Zhang et al; licensee BioMed Central Ltd. 2013
Received: 4 August 2012
Accepted: 22 February 2013
Published: 7 March 2013
Ovarian follicular development and hormone secretion are complex and coordinated biological processes which will usually be altered during pregnancy. Ovarian function is tightly regulated by a multitude of genes, and also by some specific miRNAs. It is necessary to identify the differentially expressed miRNAs in the ovaries of pregnant and non-pregnant mammals, in order to further understand the role of miRNA-mediated post-transcriptional regulation in mammalian reproduction. Here, we performed a comprehensive search for hircine miRNAs using two small RNA sequencing libraries prepared from the ovaries of pregnant and non-pregnant goats.
617 conserved and 7 putative novel miRNAs were identified in the hircine ovaries. A total of 471 conserved miRNAs (76.34%) were co-expressed in both pregnant and non-pregnant libraries, and 90 pregnancy-specific and 56 non-pregnancy-specific conserved miRNAs were identified. Additionally, 407 unique miRNAs (65.96%) were significantly differentially expressed in the pregnant and non-pregnant libraries, of which 294 were upregulated and 113 were downregulated in the pregnant library compared to the non-pregnant library. Further analysis showed that miR-143 was predicted to bind to the target sequences of Frizzled-6 and -3 receptor genes in the Wnt/beta-catenin signaling pathway, and let-7b may target the Activin receptor I and Smad 2/3 genes in the TGF-beta signaling pathway. The expression level of 5 randomly selected miRNAs were analyzed by quantitative real-time PCR (q-PCR), and the results demonstrated that the expression patterns were consistent with the Solexa sequencing results.
The identification and characterization of differentially expressed miRNAs in the ovaries of pregnant and non-pregnant goats provides important information on the role of miRNA in the regulation of the ovarian development and function. This data will be helpful to facilitate studies on the regulation of miRNAs during mammalian reproduction.
KeywordsMicroRNA Solexa sequencing Ovary Anhui white goat
MicroRNAs (miRNAs) are a group of single-stranded noncoding 21–24 nt RNAs which are involved in diverse aspects of eukaryotic biology including reproduction, development, pathogenesis, cell proliferation, apoptosis and lipometabolism by pairing to mRNAs which mainly results in target-specific post-transcriptional repression [1–4]. Recent research has demonstrated that miRNAs are involved in the regulation of mammalian reproduction, especially the regulation of ovarian function [4, 5]. The Dicer1-deficient mice which lost ~75% of Dicer1 mRNA expression demonstrated female infertility, and were predicted to reduce angiogenesis in the corpus luteum because of the loss of miR-17-5p and let-7b which regulate the expression of tissue inhibitor of metalloproteinase 1 [6, 7]. When the mouse granulosa cells were treated with the human chorionic gonadotropin (hCG), the miR-132 and miR-212 expression were upregulated, and knockdown of these miRNAs increased the expression of C-terminal binding protein 1 (CtBP1) in granulosa cells . CtBP1 was recently shown to regulate adrenal steroidogenesis, in conjunction with steroidogenic factor-1 . Additionally, massive parallel sequencing in the newborn mouse ovaries showed the expression of 398 known miRNAs, among which the X-linked mir-503, mir-672 and mir-465 family were found to be preferentially expressed in the testes and ovaries which play important roles in folliculogenesis and spermatogenesis, respectively .
Pregnancy is a complex reproductive process, which is tightly regulated by various endocrine factors and a large number of genes. The ovaries play an important role during pregnancy. There are significant differences in the activity and endocrine characteristics of the ovary during pregnancy and non-pregnancy [11–13]. In the non-pregnant phase, ovulation is normal and estrogen secretion dominates; whereas ovulation is temporarily suspended during pregnancy and progesterone secretion gradually increases to maintain pregnancy. In the present study, we characterized and investigated the differential expression of miRNAs in the ovaries of pregnant and non-pregnant goats using deep sequencing technology. The result will help to further understand the role of miRNAs in reproductive biological processes, including follicular development, hormone secretion, luteinization and the maintenance of pregnancy, and also may help to identify miRNAs which could be potentially used to regulate hircine reproduction and breeding practice in the future.
Overview of sequencing data
Conserved microRNAs and microRNA*s
To identify conserved miRNAs in hircine ovaries, the dataset was compared to the known mammalian miRNAs (miRNA precursors and mature miRNAs) in miRBase 18.0 (http://www.mirbase.org). A total of 36501 and 47316 unique sequences in the non-pregnant and pregnant libraries were mapped to known mammalian miRNAs in miRBase 18.0, respectively. Considering one or two mismatches between sequences, 508 and 535 conserved miRNAs were identified in the non-pregnant and pregnant libraries, respectively.
Comparison of the read counts for miRNAs and the corresponding miRNA*s
Identification of putative miRNAs
Though the substantial goat genome sequence data is currently short, it is feasible to identify putative miRNAs by means of alignment with goat expressed sequence tags (ESTs). Two goat putative miRNAs with a total of 24 read counts, and five goat putative miRNAs with a total of 43 read counts were predicted in pregnant and non-pregnant libraries. The RNA stem loop structures and the positions of the reads were shown in Additional file 5 and Additional file 4: Figure S2. These seven putative miRNAs were not analyzed further, as their expression levels were too low in the two small RNA libraries.
Differential expression of miRNAs in the ovaries of pregnant and non-pregnant goats
MiRNA target prediction
Validation of hircine miRNAs
Goats are one of the important economical animals, as they provide high-quality wool, meat and other products. However, hircine fecundity is relatively low, and this is a major constraint which prevents the development of the goat industry. Fecundity is a low heritability trait and conventional selection only supports slow improvement. Therefore, new methods to improve hircine fertility are highly desirable [21–23]. Molecular assisted breeding technology has been successfully used in the goat breeding industry, and the resource population of high fecundity of Anhui White goats was created successfully by our research group. In this study, we sequenced the small RNAs in the ovarian tissues of pregnant and non-pregnant Anhui White goats using Illumina Solexa technology to identify ovary-specific and differentially expressed miRNAs. Over eighty percent of known conserved mammalian miRNAs were found to be expressed in the ovaries of pregnant and non-pregnant goats (Figure 2A and C).
As the complete goat genome sequence has not yet been published, and goat-specific miRNAs are not recorded in the miRBase 18.0 database , we can currently only research hircine miRNAs using the miRNA sequences of closely related species (including mammals such as Ovis aries, Bos taurus, Sus scrofa, Canis familiaris and Equus caballus), or using goat EST sequences and Unigene information. Although 617 conserved miRNAs were found to be expressed in hircine ovaries using BLASTN in the study, the identification was highly dependent on the quantity and quality of the data in miRBase 18.0. Q-PCR was preformed to analyze the expression of 5 selected differentially expressed miRNAs in pregnant and non-pregnant hircine ovaries, and the results were consistent with the Solexa sequencing data. In theory, the expression levels of each candidate miRNA need to be validated using q-PCR; therefore, the 617 miRNAs identified in this study can only be regarded as a hircine ovary-specific miRNA reference dataset, and further research should be performed to validate the expression of miRNAs of interest in a larger number of samples. On the other hand, only one sample per group (pregnant or non-pregnant) was sequenced in this study, so no expression variances within the groups and intro-species could be estimated. To overcome this deficiency, we will add 1-2 samples in each physio-stage, or add samples with these conditions (before pregnancy, in pregnancy, after parturition) in the future, to make the library a complete process of pregnancy. Furthermore, we predicted 7 putative miRNAs using the goat EST database; however, the expression levels of these miRNAs in the libraries were very low. This discrepancy is due to the shortage of goat EST sequence information in the Genbank database, and is consistent with the previous studies [25, 26]. With continual improvements in goat genome sequence information, we will be able to obtain richer, more accurate data on hircine miRNAs.
In the view of the different ovarian activity and endocrine function in pregnant and non-pregnant goats, differentially expressed miRNAs were identified in the ovaries of pregnant and non-pregnant goats which enabled to analyze the relationship between hircine miRNAs and reproductive traits such as follicular development, hormone secretion, luteinization and pregnancy maintenance. We identified 90 and 56 miRNAs which were specifically expressed in the ovaries of pregnant and non-pregnant goats, respectively. Some miRNAs which were not expressed in the ovaries of non-pregnant goats were expressed during pregnancy (Figure 3 and Additional file 6). MiRNAs usually regulate protein expression by binding to and repressing translation or promoting the degradation of their target mRNAs [27, 28]. The present study indicated that upregulation of miRNAs in the ovaries of pregnant goats may inhibit the expression of target genes associated with follicular development, ovulation and estrogen secretion; thereby inhibiting ovulation and estrogen secretion. On the contrary, downregulation of other miRNAs in the ovaries of pregnant goats may relieve the repression of target genes associated with luteinization, progesterone secretion and pregnancy maintenance; thereby promoting luteinization and progesterone secretion.
Although much research has been performed on miRNAs in the gonads of pigs, cattle and other mammals, research of goat ovarian miRNAs has seldom been reported. In previous, miR-143 was found to be the most highly expressed miRNA in the testis and ovaries of Holstein cows, and 10 putative miRNA target genes involved in the GnRH-signaling and insulin-signaling pathways, which are associated with endocrine system function, were also identified . Interestingly, in the study, miR-143 was also found to be the most highly expressed miRNA in the ovaries of pregnant goats (normalized expression level of 276,985). Target prediction indicated that miR-143 may bind to the Frizzled-6 and -3 receptor genes in the Wnt signaling pathway, thereby affecting the binding of Wnt-4 to its receptor (Figure 6A and B). Moreover, previous studies have suggested that the Wnt-4 gene, one of the most important members of the Wnt family, may regulate the function of ovarian granulosa cells and luteal cells by binding to specific members of the Frizzled receptor family [30, 31]. Therefore, miR-143 may play an important role in reproduction, such as mammalian gonadal endocrine function and pregnancy maintenance. Other studies have also shown that miR-143 promotes normal adipocyte differentiation and fat deposition [2, 32–34], indicating that miR-143 may play an important role in normal mammalian physical development, as well as mammalian reproduction.
It is a challenge to verify the complex functions of miRNAs. In the previous studies, miR-21 exerted an anti-apoptotic effect during the transformation of ovarian granulosa cells into luteal cells, and repression of miR-21 expression induced granulosa cell apoptosis and significantly reduced the rate of ovulation, via a mechanism dependent on luteinizing hormone secretion . Also, miR-21 was demonstrated to express at significantly higher levels in the ovaries of Holstein cows compared to the testis (1.97-fold), indicating that miR-21 may play an important role in ovarian function . Interestingly, the present study also showed that miR-21 was significantly upregulated in pregnant goat ovaries, compared to non-pregnant goat ovaries (1.92-fold). Increased LH secretion after becoming pregnant, followed by upregulation of miR-21 expression, may contribute to the transformation of ovarian granulosa cells into luteal cells. And the mechanisms regulating the anti-apoptotic effects of miR-21 in ovarian granulosa cells still need to be investigated further .
Recent research has indicated that specific members of the let-7 family can affect mammalian reproduction, development, cell proliferation and apoptosis [37–40]. In the present study, eight members of the let-7 family (let-7a, let-7b, let-7c, let-7d, let-7e, let-7f, let-7 g and let-7i) were expressed at high levels in both libraries, and let-7b was the most highly expressed miRNA in the ovaries of non-pregnant goats (normalized expression level of 246,928). Bioinformatics analysis indicated that let-7b may bind to the Activin receptor I and Smad2/3 genes in the TGF-beta signaling pathway, and may affect follicular development and estrogen secretion (Figure 6C and D).
MiRNA expression patterns vary in the ovaries of pregnant and non-pregnant goats. In total, 617 conserved and 7 putative miRNAs were detected, and 407 differentially expressed miRNAs were identified in pregnant and non-pregnant ovaries, suggesting that miRNAs may play an important role in the regulation of goat ovarian function. For example, miRNA-143 was strongly associated with reproduction, and may potentially regulate the Wnt/beta-catenin signaling pathway by targeting Frizzled-6 and -3. Let-7b may regulate the TGF-beta signaling pathway by targeting the Activin receptor І and Smad2/3. Future work to characterize the expression of ovarian miRNAs at different stages of reproduction and in different breeds of goat, or in specific cell lines derived from ovarian tissues, is necessary to fully elucidate the functions of miRNAs in goat follicular development and hormone secretion, which will help to understand the relationships between miRNAs and mammalian reproduction, while enhancing the development of artificial reproduction and marker assisted selection (MAS) techniques in goats.
Anhui White goats (a Chinese indigenous breed) were obtained from the College of Animal Science and Technology, Anhui Agricultural University, Hefei, China. Experiments were performed according to the Regulations for the Administration of Affairs Concerning Experimental Animals (Ministry of Science and Technology, China; revised in June 2004) and approved by the ethics committee of Anhui Agricultural University, Anhui, China, under permit No. AHAU20101025. The animals were allowed access to feed and water ad libitum under normal condition and were sacrificed humanely to minimize suffering.
Ovary collection and total RNA isolation
The ovaries of Anhui White goats were collected and used to generate small RNA libraries. The collected ovaries were divided into two groups: six ovaries were from three 24-month old non-pregnant goats and six ovaries were from three 24 -month old pregnant goats. The ovaries were immersed in liquid nitrogen immediately after collection and stored at -80°C. Total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions. The quality of the total RNA was checked using the Agilent 2100 Bioanalyzer system (Santa Clara, CA, USA) and the samples were stored at -80°C until analysis.
Small RNA library construction and sequencing
Two groups of total RNA were used for library preparation and sequencing by pooling equal quantity (10 μg) of total RNA isolated from six individual pregnant or non-pregnant goat ovaries. Briefly, total RNA were purified by polyacrylamide gel electrophoresis (PAGE) to enrich 15–35 nt molecules, then proprietary adapters were ligated to the 5′ and 3′termini of the RNAs and the samples were used as templates for cDNA synthesis. The cDNA was amplified using the appropriate number of PCR cycles to produce sequencing libraries, which were subsequently subjected to the proprietary Solexa sequencing-by-synthesis method using the Illumina Genome Analyzer (SanDiego, CA, USA). Sequencing was carried out at the Beijing Genomics Institute (BGI, Beijing, China).
According to the principle of bioinformatics analysis, low-quality reads were removed from the raw reads. After trimming the 3′adaptor sequence, removing 5′ adaptor contaminants and counting the total, unique and length of reads, all valid sequences were obtained for further analysis. The overall flow of the sequencing data analysis is represented schematically in Additional file 4: Figure S3. All unique sequences were used to search the ncRNA data (Genbank, Repeat sequence and Rfam) with BLASTN to remove non-miRNA sequences (rRNA, tRNA, snoRNA, snRNA, etc.). Subsequently, the remaining sequences were analyzed using a BLAST search against miRBase 18.0. Sequences in the libraries with identical or related sequences (1~2 nucleotide substitutions permitted) to Ovis aries or other mammals (Bos taurus, Sus scrofa, Canis familiaris, Equus caballus) were identified as conserved miRNAs. Although the full goat genome sequence has not yet been published, we integrated data from the small RNA libraries with the goat EST sequences (http://www.ncbi.nlm.nih.gov/nucest) to identify goat putative miRNAs using the Mireap software (http://sourceforge.net/projects/mireap). Sequences with a perfect match or one mismatch were retained for further analysis. Subsequently, 60-80 nt of the EST sequences were extracted, and secondary structure was predicted and analyzed with Mireap using specific parameter settings .
The x and y represent normalized expression level, and the N1 and N2 represent total count of clean reads of a given miRNA in small RNA library of ovaries of pregnant and non-pregnant goats, respectively . Additionally, Hircine Unigene sequences (http://www.ncbi.nlm.nih.gov/unigene) were chosen to predict miRNA targets with RNAhybrid, using the parameter settings described by Rehmsmeier et al. (2004) .
MiRNA validation via q-PCR
Quantitative PCR (q-PCR) was used to validate 5 randomly selected miRNAs that were differentially expressed by Solexa sequencing. One microgram of total RNA from each sample were reverse-transcript into cDNA using the miScript Reverse Transcription Kit (Qiagen, Dusseldorf, Germany) according to the manufacturer’s instructions. After incubation at 37°C for 1 h and deactivation at 95°C for 10 min, the mix was used as the template for q-PCR. Q-PCR was performed using standard protocols on the Roche LightCycler 480 II Real-Time PCR Detection System (Roche; LC480 II, Basel, Switzerland). In each assay, 1 μl of cDNA was added to 19 μl of mix containing 10 μl 2×SYBR green SuperReal PreMix (TIANGEN, Beijing, China), 0.4 μl of each primer and 8.2 μl RNase free H2O. The reaction was amplified at 95°C for 15 min, followed by 40 cycles of 95°C 10 s and 60°C 30 s. For each miRNA, three biological replicates were performed, and all of the reactions were run in triplicate. The cycle threshold was collected from each reaction, and the relative expression level of each miRNA to U6 snRNA was evaluated using the equation 2-(CTmiRNA-CTU6RNA), and the fold-change (log2-ratio) and P value were used to show the differential expression of miRNA in the two samples [43–45]. The miRNA-specific primers were presented in Additional file 7.
This work was supported by grants from the National High Science and Technology Foundation of China “863” (No. 2011AA100307-4), the Technology Innovation Project Special Program and the Science and Technology Program of Anhui Province (No. 11Z0101095 and No. 11010302108). We are grateful to Yunsheng Li and Chunhuan Ren for construction of small RNA libraries, and Yu Zhang, Yong Pu and Meiling Zhang for assistance with data analysis.
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871.View ArticlePubMedGoogle Scholar
- Esau C, Kang X, Peralta E, Hanson E, Marcusson EG, Ravichandran LV, Sun Y, Koo S, Perera RJ, Jain R: MicroRNA-143 regulates adipocyte differentiation. J Biol Chem. 2004, 279 (50): 52361-52365. 10.1074/jbc.C400438200.View ArticlePubMedGoogle Scholar
- Inoue K: MicroRNA function in animal development. Tanpakushitsu kakusan koso Protein, nucleic acid, enzyme. 2007, 52 (3): 197-204.PubMedGoogle Scholar
- Baley J, Li J: MicroRNAs and ovarian function. J Ovarian Res. 2012, 5: 8-10.1186/1757-2215-5-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Cutting AD, Bannister SC, Doran TJ, Sinclair AH, Tizard MV, Smith CA: The potential role of microRNAs in regulating gonadal sex differentiation in the chicken embryo. Chromosome research: an international journal on the molecular, supramolecular and evolutionary aspects of chromosome biology. 2012, 20 (1): 201-213. 10.1007/s10577-011-9263-y.View ArticleGoogle Scholar
- Otsuka M, Jing Q, Georgel P, New L, Chen J, Mols J, Kang YJ, Jiang Z, Du X, Cook R, Das SC, Pattnaik AK, Beutler B, Han J: Hypersusceptibility to vesicular stomatitis virus infection in Dicer1-deficient mice is due to impaired miR24 and miR93 expression. Immunity. 2007, 27 (1): 123-134. 10.1016/j.immuni.2007.05.014.View ArticlePubMedGoogle Scholar
- Otsuka M, Zheng M, Hayashi M, Lee JD, Yoshino O, Lin S, Han J: Impaired microRNA processing causes corpus luteum insufficiency and infertility in mice. J Clin Invest. 2008, 118 (5): 1944-1954. 10.1172/JCI33680.PubMed CentralView ArticlePubMedGoogle Scholar
- Fiedler SD, Carletti MZ, Hong XM, Christenson LK: Hormonal regulation of microRNA expression in periovulatory mouse mural granulosa cells. Biol Reprod. 2008, 79 (6): 1030-1037. 10.1095/biolreprod.108.069690.PubMed CentralView ArticlePubMedGoogle Scholar
- Dammer EB, Sewer MB: Phosphorylation of CtBP1 by cAMP-dependent protein kinase modulates induction of CYP17 by stimulating partnering of CtBP1 and 2. J Biol Chem. 2008, 283 (11): 6925-6934. 10.1074/jbc.M708432200.PubMed CentralView ArticlePubMedGoogle Scholar
- Ahn HW, Morin RD, Zhao H, Harris RA, Coarfa C, Chen ZJ, Milosavljevic A, Marra MA, Rajkovic A: MicroRNA transcriptome in the newborn mouse ovaries determined by massive parallel sequencing. Mol Hum Reprod. 2010, 16 (7): 463-471. 10.1093/molehr/gaq017.PubMed CentralView ArticlePubMedGoogle Scholar
- Judd HL, Judd GE, Lucas WE, Yen SS: Endocrine function of the postmenopausal ovary: concentration of androgens and estrogens in ovarian and peripheral vein blood. J Clin Endocrinol Metab. 1974, 39 (6): 1020-1024. 10.1210/jcem-39-6-1020.View ArticlePubMedGoogle Scholar
- Richards JS: Hormonal control of gene expression in the ovary. Endocr Rev. 1994, 15 (6): 725-751.View ArticlePubMedGoogle Scholar
- Stephens SM, Moley KH: Follicular origins of modern reproductive endocrinology. Am J Physiol Endocrinol Metab. 2009, 297 (6): E1235-E1236. 10.1152/ajpendo.00575.2009.View ArticlePubMedGoogle Scholar
- Xie SS, Li XY, Liu T, Cao JH, Zhong Q, Zhao SH: Discovery of porcine microRNAs in multiple tissues by a Solexa deep sequencing approach. PLoS One. 2011, 6 (1): e16235-10.1371/journal.pone.0016235.PubMed CentralView ArticlePubMedGoogle Scholar
- Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18 (6): 957-964. 10.1101/gr.074740.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Li T, Wu R, Zhang Y, Zhu D: A systematic analysis of the skeletal muscle miRNA transcriptome of chicken varieties with divergent skeletal muscle growth identifies novel miRNAs and differentially expressed miRNAs. BMC Genomics. 2011, 12: 186-10.1186/1471-2164-12-186.PubMed CentralView ArticlePubMedGoogle Scholar
- Jadhav V, Hackl M, Bort JAH, Wieser M, Harreither E, Kunert R, Borth N, Grillari J: A screening method to assess biological effects of microRNA overexpression in Chinese hamster ovary cells. Biotechnol Bioeng. 2012, 109 (6): 1376-1385. 10.1002/bit.24490.View ArticlePubMedGoogle Scholar
- Xu SY, Linher-Melville K, Yang BB, Wu D, Li JL: Micro-RNA378 (miR-378) regulates ovarian estradiol production by targeting aromatase. Endocrinology. 2011, 152 (10): 3941-3951. 10.1210/en.2011-1147.PubMed CentralView ArticlePubMedGoogle Scholar
- Reid JG, Nagaraja AK, Lynn FC, Drabek RB, Muzny DM, Shaw CA, Weiss MK, Naghavi AO, Khan M, Zhu HF, Tennakoon J, Gunaratne GH, Corry DB, Miller J, McManus MT, German MS, Gibbs RA, Matzuk MM, Gunaratne PH: Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5′-seed/cleavage/anchor regions and stabilize predicted mmu-let-7a: mRNA duplexes. Genome Res. 2008, 18 (10): 1571-1581. 10.1101/gr.078246.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Pandey DP, Picard D: miR-22 inhibits estrogen signaling by directly targeting the estrogen receptor alpha mRNA (vol 29, pg 3783, 2009). Mol Cell Biol. 2009, 29 (17): 4873-4873. 10.1128/MCB.00830-09.PubMed CentralView ArticleGoogle Scholar
- Amoah EA, Gelaye S: Biotechnological advances in goat reproduction. J Anim Sci. 1997, 75 (2): 578-585.PubMedGoogle Scholar
- Zhang CY, Chen SL, Li X, Xu DQ, Zhang Y, Yang LG: Genetic and phenotypic parameter estimates for reproduction traits in the Boer dam. Livest Sci. 2009, 125 (1): 60-65. 10.1016/j.livsci.2009.03.002.View ArticleGoogle Scholar
- Leboeuf B, Delgadillo JA, Manfredi E, Piacere A, Clement V, Martin P, Pellicer M, Boue P, de Cremoux R: Management of goat reproduction and insemination for genetic improvement in France. Reproduction in domestic animals = Zuchthygiene. 2008, 43 (Suppl 2): 379-385.View ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-D157.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei Y, Chen S, Yang P, Ma Z, Kang L: Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Biol. 2009, 10 (1): R6-10.1186/gb-2009-10-1-r6.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Y, Shi Z, Wu M, Zhang J, Jia L, Chen X: Identification and differential expression of microRNAs during metamorphosis of the Japanese flounder (Paralichthys olivaceus). PLoS One. 2011, 6 (7): e22957-10.1371/journal.pone.0022957.PubMed CentralView ArticlePubMedGoogle Scholar
- Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett. 2005, 579 (26): 5911-5922. 10.1016/j.febslet.2005.07.070.View ArticlePubMedGoogle Scholar
- Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301 (5631): 336-338. 10.1126/science.1085242.View ArticlePubMedGoogle Scholar
- Huang J, Ju Z, Li Q, Hou Q, Wang C, Li J, Li R, Wang L, Sun T, Hang S, Gao Y, Hou M, Zhong J: Solexa sequencing of novel and differentially expressed microRNAs in testicular and ovarian tissues in Holstein cattle. Int J Biol Sci. 2011, 7 (7): 1016-1026.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh M, Johnson MA, Greenberg NM, Richards JS: Regulated expression of Wnts and frizzleds at specific stages of follicular development in the rodent ovary. Endocrinology. 2002, 143 (3): 898-908. 10.1210/en.143.3.898.View ArticlePubMedGoogle Scholar
- Ricken A, Lochhead P, Kontogiannea M, Farookhi R: Wnt signaling in the ovary: identification and compartmentalized expression of wnt-2, wnt-2b, and frizzled-4 mRNAs. Endocrinology. 2002, 143 (7): 2741-2749. 10.1210/en.143.7.2741.View ArticlePubMedGoogle Scholar
- Takanabe R, Ono K, Abe Y, Takaya T, Horie T, Wada H, Kita T, Satoh N, Shimatsu A, Hasegawa K: Up-regulated expression of microRNA-143 in association with obesity in adipose tissue of mice fed high-fat diet. Biochem Biophys Res Commun. 2008, 376 (4): 728-732. 10.1016/j.bbrc.2008.09.050.View ArticlePubMedGoogle Scholar
- Wang T, Li M, Guan J, Li P, Wang H, Guo Y, Shuai S, Li X: MicroRNAs miR-27a and miR-143 regulate porcine adipocyte lipid metabolism. Int J Mol Sci. 2011, 12 (11): 7950-7959.PubMed CentralView ArticlePubMedGoogle Scholar
- Yi C, Xie WD, Li F, Lv Q, He J, Wu J, Gu D, Xu N, Zhang Y: MiR-143 enhances adipogenic differentiation of 3 T3-L1 cells through targeting the coding region of mouse pleiotrophin. FEBS Lett. 2011, 585 (20): 3303-3309. 10.1016/j.febslet.2011.09.015.View ArticlePubMedGoogle Scholar
- Carletti MZ, Fiedler SD, Christenson LK: MicroRNA 21 blocks apoptosis in mouse periovulatory granulosa cells. Biol Reprod. 2010, 83 (2): 286-295. 10.1095/biolreprod.109.081448.PubMed CentralView ArticlePubMedGoogle Scholar
- Christenson LK: MicroRNA control of ovarian function. Animal reproduction / Colegio Brasileiro de Reproducao Animal. 2010, 7 (3): 129-133.PubMed CentralPubMedGoogle Scholar
- Roush S, Slack FJ: The let-7 family of microRNAs. Trends Cell Biol. 2008, 18 (10): 505-516. 10.1016/j.tcb.2008.07.007.View ArticlePubMedGoogle Scholar
- Fu TY, Lin CT, Tang PC: Steroid hormone-regulated let-7b mediates cell proliferation and basigin expression in the mouse endometrium. J Reprod Dev. 2011, 57 (5): 627-635. 10.1262/jrd.11-018E.View ArticlePubMedGoogle Scholar
- Saleh AD, Savage JE, Cao L, Soule BP, Ly D, DeGraff W, Harris CC, Mitchell JB, Simone NL: Cellular stress induced alterations in microRNA let-7a and let-7b expression are dependent on p53. PLoS One. 2011, 6 (10): e24429-10.1371/journal.pone.0024429.PubMed CentralView ArticlePubMedGoogle Scholar
- Miles JR, McDaneld TG, Wiedmann RT, Cushman RA, Echternkamp SE, Vallet JL, Smith TP: MicroRNA expression profile in bovine cumulus-oocyte complexes: possible role of let-7 and miR-106a in the development of bovine oocytes. Anim Reprod Sci. 2012, 130 (1–2): 16-26.View ArticlePubMedGoogle Scholar
- Bentwich I: Prediction and validation of microRNAs and their targets. FEBS Lett. 2005, 579 (26): 5904-5910. 10.1016/j.febslet.2005.09.040.View ArticlePubMedGoogle Scholar
- Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10 (10): 1507-1517. 10.1261/rna.5248604.PubMed CentralView ArticlePubMedGoogle Scholar
- Deo A, Carlsson J, Lindlof A: How to choose a normalization strategy for miRNA quantitative real-time (qPCR) arrays. J Bioinform Comput Biol. 2011, 9 (6): 795-812. 10.1142/S0219720011005793.View ArticlePubMedGoogle Scholar
- Mestdagh P, Feys T, Bernard N, Guenther S, Chen C, Speleman F, Vandesompele J: High-throughput stem-loop RT-qPCR miRNA expression profiling using minute amounts of input RNA. Nucleic Acids Res. 2008, 36 (21): e143-10.1093/nar/gkn725.PubMed CentralView ArticlePubMedGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8 (2): R19-10.1186/gb-2007-8-2-r19.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.