Identification and characterization of microRNAs in the ovaries of multiple and uniparous goats (Capra hircus) during follicular phase

Identification and analysis of differentially expressed microRNAs in endometrium to explore the regulation of sheep fecundity

BACKGROUND

MicroRNAs (miRNAs) play an important regulatory role in mammalian reproduction. Currently, most studies are primarily concentrated on ovarian miRNAs, ignoring the influence of endometrial miRNAs on the fecundity of female sheep. To uncover potential regulators of sheep fecundity, RNA-seq was used to comparatively analyze miRNA expression profiles of endometrium between high prolificacy sheep (HP, litter size = 3) and low prolificacy sheep (LP, litter size = 1) with FecB genotype.

RESULTS

Firstly, genomic features of miRNAs from endometrium were analyzed. Furthermore, 58 differentially expressed (DE) miRNAs were found in the endometrium of Hu sheep with different litter size. A co-expression network of DE miRNAs and target genes has been constructed, and hub genes related litter size are included, such as DE miRNA unconservative_NC_019472.2_1229533 and unconservative_NC_019481.2_1637827 target to estrogen receptor α (ESR1) and unconservative_NC_019481.2_1637827 targets to transcription factor 7 (TCF7). Moreover, functional annotation analysis showed that the target genes (NRCAM and NEGR1) of the DE miRNAs were significantly enriched in cell adhesion molecules (CAMs) signaling pathway, which was related to uterine receptivity.

CONCLUSION

Taken together, this study provides a new valuable resource for understanding the molecular mechanisms underlying Hu sheep prolificacy.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

Identification and profiling of microRNAs during yak's testicular development

BACKGROUND

Normal testicular development is highly crucial for male reproduction and is a precondition for spermatogenesis that is the production of spermatozoa in the testes. MiRNAs have been implicated in several testicular biological processes, including cell proliferation, spermatogenesis, hormone secretion, metabolism and reproductive regulation. In the present study, we used deep sequencing data to study the functions of miRNAs during testicular development and spermatogenesis, by analyzing the expression patterns of small RNAs in 6-, 18- and 30-month-old yak testis tissues.

RESULTS

A total of 737 known and 359 novel miRNAs were obtained from 6-, 18- and 30-month-old yak testes. In all, we obtained 12, 142 and 139 differentially expressed (DE) miRNAs in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of all DE miRNA target genes revealed BMP2, TGFB2, GDF6, SMAD6, TGFBR2 and other target genes as participants in different biological processes, including TGF-β, GnRH, Wnt, PI3K-Akt, MAPK signaling pathways and several other reproductive pathways. In addition, quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) was used to detect the expression of seven randomly selected miRNAs in 6-, 18- and 30-month-old testes, and the results were consistent with the sequencing data.

CONCLUSIONS

The differential expression of miRNAs in yak testes at different development stages was characterized and investigated using deep sequencing technology. We believe that the results will contribute to further understanding the functions of miRNAs in regulating the development of yak testes and improving the reproductive performance of male yaks.

Identification of differentially expressed microRNAs in ovulatory and subordinate follicles in Dazu black goats

How multiple ovulations happen in prolific goats is still unknown even though studies on ovarian physiology and folliculogenesis have made extensive progress. MicroRNAs (miRNAs) are a class of endogenous regulatory factors and regulate structural gene expression mainly at the post-transcriptional level. In this study, in the follicular phase, six large follicles were collected from prolific Dazu black goat and used to generate RNA libraries for RNA sequencing. Based on the litter size and average number of ovulatory follicles in Dazu black goats, the largest three follicles were sorted as ovulatory follicles, and the remaining as subordinate ones. In total, 418 known miRNAs and 110 novel miRNAs were found, and the expression of six randomly selected miRNAs was validated by quantitative PCR analysis. Nine miRNAs were differently expressed between the ovulatory and subordinate follicles (p < 0.01). Chi-miR-582-5p, novel-130, chi-miR-214-3p, and chi-miR-500-5p were upregulated in the ovulatory group, and chi-miR-383, chi-miR-130b-5p, chi-miR-92a-3p, chi-miR-125b-5p, and novel-9 were downregulated in the same group. Chi-miR-130b-5p and chi-miR-214-3p were predicted to target at LHR (XM_013967581.1), GDF9 (NM_001285708.1), BMP15 (NM_001285588.1), and CYP19A1 (XM_013967046.1). In conclusion, nine miRNAs were differently expressed between ovulatory and subordinate follicles, and chi-miR-130b-5p and chi-miR-214-3p were predicted to regulate the expression of genes involved in gonadotropin hormone signaling and oocyte-derived growth factors. Additional studies are needed to confirm these findings.

Expression characteristics of piRNAs in ovine luteal phase and follicular phase ovaries

PIWI-interacting RNAs (piRNAs), as a novel class of small non-coding RNAs that have been shown to be indispensable in germline integrity and stem cell development. However, the expressed characteristics and regulatory roles of piRNAs during different reproductive phases of animals remain unknown. In this study, we investigated the piRNAs expression profiles in ovaries of sheep during the luteal phase (LP) and follicular phase (FP) using the Solexa sequencing technique. A total of 85,219 and 1,27,156 piRNAs tags were identified in ovine ovaries across the two phases. Most expressed piRNAs start with uracil. piRNAs with a length of 24 nt or 27–29 nts accounted for the largest proportion. The obvious ping-pong signature appeared in the FP ovary. The piRNA clusters in the sheep ovary were unevenly distributed on the chromosomes, with high density on Chr 3 and 1. For genome distribution, piRNAs in sheep ovary were mainly derived from intron, CDS, and repeat sequence regions. Compared to the LP ovary, a greater number of expressed piRNA clusters were detected in the FP ovary. Simultaneously, we identified 271 differentially expressed (DE) piRNAs between LP and FP ovaries, with 96 piRNAs upregulated and 175 piRNAs downregulated, respectively. Functional enrichment analysis (GO and KEGG) indicated that their target genes were enriched in reproduction-related pathways including oocyte meiosis, PI3K-Akt, Wnt, and TGF-β signaling pathways. Together, our results highlighted the sequence and expression characteristics of the piRNAs in the sheep ovary, which will help us understand the roles of piRNAs in the ovine estrus cycle.

Integration of small RNA, degradome, and transcriptome sequencing data illustrates the mechanism of low phosphorus adaptation in Camellia oleifera

Phosphorus (P) is an indispensable macronutrient for plant growth and development, and it is involved in various cellular biological activities in plants. Camellia oleifera is a unique high-quality woody oil plant that grows in the hills and mountains of southern China. However, the available P content is deficient in southern woodland soil. Until now, few studies focused on the regulatory functions of microRNAs (miRNAs) and their target genes under low inorganic phosphate (Pi) stress. In this study, we integrated small RNA, degradome, and transcriptome sequencing data to investigate the mechanism of low Pi adaptation in C. oleifera. We identified 40,689 unigenes and 386 miRNAs by the deep sequencing technology and divided the miRNAs into four different groups. We found 32 miRNAs which were differentially expressed under low Pi treatment. A total of 414 target genes of 108 miRNAs were verified by degradome sequencing. Gene ontology (GO) functional analysis of target genes found that they were related to the signal response to the stimulus and transporter activity, indicating that they may respond to low Pi stress. The integrated analysis revealed that 31 miRNA-target pairs had negatively correlated expression patterns. A co-expression regulatory network was established based on the profiles of differentially expressed genes. In total, three hub genes (ARF22, WRKY53, and SCL6), which were the targets of differentially expressed miRNAs, were discovered. Our results showed that integrated analyses of the small RNA, degradome, and transcriptome sequencing data provided a valuable basis for investigating low Pi in C. oleifera and offer new perspectives on the mechanism of low Pi tolerance in woody oil plants.

Integrated analyses of miRNA-mRNA expression profiles of ovaries reveal the crucial interaction networks that regulate the prolificacy of goats in the follicular phase

BACKGROUND

Litter size is an important index of mammalian prolificacy and is determined by the ovulation rate. The ovary is a crucial organ for mammalian reproduction and is associated with follicular development, maturation and ovulation. However, prolificacy is influenced by multiple factors, and its molecular regulation in the follicular phase remains unclear.

METHODS

Ten female goats with no significant differences in age and weight were randomly selected and divided into either the high-yielding group (n = 5, HF) or the low-yielding group (n = 5, LF). Ovarian tissues were collected from goats in the follicular phase and used to construct mRNA and miRNA sequencing libraries to analyze transcriptomic variation between high- and low-yield Yunshang black goats. Furthermore, integrated analysis of the differentially expressed (DE) miRNA-mRNA pairs was performed based on their correlation. The STRING database was used to construct a PPI network of the DEGs. RT-qPCR was used to validate the results of the predicted miRNA-mRNA pairs. Luciferase analysis and CCK-8 assay were used to detect the function of the miRNA-mRNA pairs and the proliferation of goat granulosa cells (GCs).

RESULTS

A total of 43,779 known transcripts, 23,067 novel transcripts, 424 known miRNAs and 656 novel miRNAs were identified by RNA-seq in the ovaries from both groups. Through correlation analysis of the miRNA and mRNA expression profiles, 263 negatively correlated miRNA-mRNA pairs were identified in the LF vs. HF comparison. Annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes such as "estrogen receptor binding (GO:0030331)", "oogenesis (GO:0048477)", "ovulation cycle process (GO:0022602)" and "ovarian follicle development (GO:0001541)". Subsequently, five KEGG pathways (oocyte meiosis, progesterone-mediated oocyte maturation, GnRH signaling pathway, Notch signaling pathway and TGF-β signaling pathway) were identified in the interaction network related to follicular development, and a PPI network was also constructed. In the network, we found that CDK12, FAM91A1, PGS1, SERTM1, SPAG5, SYNE1, TMEM14A, WNT4, and CAMK2G were the key nodes, all of which were targets of the DE miRNAs. The PPI analysis showed that there was a clear interaction among the CAMK2G, SERTM1, TMEM14A, CDK12, SYNE1 and WNT4 genes. In addition, dual luciferase reporter and CCK-8 assays confirmed that miR-1271-3p suppressed the proliferation of GCs by inhibiting the expression of TXLNA.

CONCLUSIONS

These results increase the understanding of the molecular mechanisms underlying goat prolificacy. These results also provide a basis for studying interactions between genes and miRNAs, as well as the functions of the pathways in ovarian tissues involved in goat prolificacy in the follicular phase.

Epigenetic transgenerational inheritance, gametogenesis and germline development†

One of the most important developing cell types in any biological system is the gamete (sperm and egg). The transmission of phenotypes and optimally adapted physiology to subsequent generations is in large part controlled by gametogenesis. In contrast to genetics, the environment actively regulates epigenetics to impact the physiology and phenotype of cellular and biological systems. The integration of epigenetics and genetics is critical for all developmental biology systems at the cellular and organism level. The current review is focused on the role of epigenetics during gametogenesis for both the spermatogenesis system in the male and oogenesis system in the female. The developmental stages from the initial primordial germ cell through gametogenesis to the mature sperm and egg are presented. How environmental factors can influence the epigenetics of gametogenesis to impact the epigenetic transgenerational inheritance of phenotypic and physiological change in subsequent generations is reviewed.

Preliminary screening of fluorine-stained osteoblastic apoptosis-related microRNA

Endemic fluorosis is a chronic systemic disease that seriously endangers human health. In high fluoride areas, people consume excessive fluoride for a long time through drinking water or food, which leads to chronic cumulative fluorosis in the body. Fluorosis can cause changes in the expression of some miRNA in cells, and the miRNA can participate in fluoride-induced osteoblast activation through various signal pathways. To observe the differential expression of apoptosis-related microRNA (miRNA) in mouse osteoblasts under the action of excessive fluoride. Primary cultured mouse osteoblasts, identified by osteocalcin (OC) and alkaline phosphatase (ALP) staining, were treated with 20 mg/L sodium fluoride and 40 mg/L sodium fluoride for 12/24 hr, respectively, to establish the fluoride staining model for comparing and analyzing the sequence of miRNA among groups by bioinformatics methods; four miRNA chains were verified by fluorescence quantitative PCR. After treatment with 20 mg/L sodium fluoride for 12 hr and 24 hr, 128 miRNA expressions were up-regulated while 36 miRNA expressions were down-regulated. In Group 40 mg/L, 130 miRNA expressions were up-regulated while 29 miRNA expressions were down-regulated after 12 hr and 24 hr; 72 miRNA were up-regulated and 2 miRNA were down-regulated at the two time points. 10 up-regulated miRNA and 2 down-regulated miRNA with higher scores in Bioinformatics software were analyzed the target genes. Fluorescence quantitative PCR verified that the expressions of four miRNA were up-regulated. Target gene analysis of the 10 selected mouse osteoblastic apoptosis-related miRNA reveals their involvement of the functions of inhibiting or promoting apoptosis, which has certain theoretical significance for early identification of skeletal fluorosis. The involved signaling pathways include the Wnt signaling pathway, ubiquitin-regulated proteolysis, Toll signaling pathway, TNF signaling pathway, pluripotent stem cell signaling pathway, MAPK signaling pathway, phosphatidylinositide metabolism, FoxO signaling pathway, ErbB signaling pathway, autophagy, and so forth.

Differential Expression and Functional Analysis of CircRNA in the Ovaries of Low and High Fecundity Hanper Sheep

Simple Summary

Litter size is an important trait affecting reproductive capacity and breeding economics in meat sheep. Consequently, revealing its molecular mechanism helps us understand multiple lambs from the genetic perspective. In this study, we provide a genome-wide expression profile of circular RNAs (circRNAs) expression in Hanper sheep, which is a new breed of meat sheep raised by cross and self-group breeding for 15 years. In this study, ovarian circular RNAs and miRNAs associated with high and low fertility Hanper sheep are identified during the follicular and luteal phases of the estrous cycle, and their potential biological functions are predicted through Gene Ontology (GO), KEGG, GSEA, STEM, WGCNA analysis.

Abstract

Litter size is a considerable quality that determines the production efficiency of mutton sheep. Therefore, revealing the molecular regulation of high and low fertility may aid the breeding process to develop new varieties of mutton sheep. CircRNAs are the important factors regulating follicular development, but their mechanism role in the regulation of litter size in Hanper sheep is not clear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either consecutive monotocous (M) or polytocous were collected. Then, we performed transcriptome sequencing to screen for differentially expressed circRNAs (DE-circRNAs) and elucidate their function. In total, 4256 circRNA derived from 2184 host genes were identified in which 183 (146 were upregulated, while 37 were downregulated) were differentially expressed in monotocous sheep in the follicular phase versus polytocous sheep in the follicular phase (MF vs. PF). Moreover, 34 circRNAs (14 were upregulated, while 20 were downregulated) were differentially expressed in monotocous sheep in the luteal phase versus polytocous sheep in the luteal sheep (ML vs. PL). This was achieved through DE-circRNAs function enrichment annotation analysis by GESA, GO, and KEGG, which function through the EGF-EGFR-RAS-JNK, TGF-β and thyroid hormone signaling pathway to affect the litter size of Hanper sheep in MF vs. PF and ML vs. PL. STEM results showed that MAPK signaling pathways play a key role in MF vs. PF and ML vs. PL. Through WGCNA analysis, AKT3 was a core gene in MF vs. PF and ML vs. PL. Moreover, competitive endogenous RNA (ceRNA) network analysis revealed the target binding sites for miRNA such as oar-miR-27a, oar-miR-16b, oar-miR-200a/b/c, oar-miR-181a, oar-miR-10a/b, and oar-miR-432 in the identified DE-cirRNAs.

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

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.

SMARCA2 is regulated by NORFA-miR-29c, a novel pathway that controls granulosa cell apoptosis and is related to female fertility

SMARCA2, an evolutionarily conserved catalytic ATPase subunit of SWI/SNF complexes, has been implicated in development and diseases; however, its role in mammalian ovarian function and female fertility is unknown. Here, we identified and characterized the 3'-UTR of the porcine SMARCA2 gene and identified a novel adenylate number variation. Notably, this mutation was significantly associated with sow litter size traits and SMARCA2 levels, due to its influence on the stability of SMARCA2 mRNA in ovarian granulosa cells (GCs). Immunohistochemistry and functional analysis showed that SMARCA2 is involved in the regulation of follicular atresia by inhibiting GC apoptosis. In addition, miR-29c, a pro-apoptotic factor, was identified as a functional miRNA that targets SMARCA2 in GCs and mediates regulation of SMARCA2 expression via the NORFA-SMAD4 axis. Although a potential miR-29c-responsive element was identified within NORFA, negative regulation of miR-29c expression by NORFA was not due to activity as a competing endogenous RNA. In conclusion, our findings demonstrate that SMARCA2 is a candidate gene for sow litter size traits, because it regulates follicular atresia and GC apoptosis. Additionally, we have defined a novel candidate pathway for sow fertility, the NORFA-TGFBR2-SMAD4-miR-29c-SMARCA2 pathway.This article has an associated First Person interview with the first author of the paper.

A Preliminary Study on the Characteristics of microRNAs in Ovarian Stroma and Follicles of Chuanzhong Black Goat during Estrus

microRNAs (miRNAs) play a significant role in ovarian follicular maturity, but miRNA expression patterns in ovarian stroma (OS), large follicles (LF), and small follicles (SF) have been rarely explored. We herein aimed to identify miRNAs, their target genes and signaling pathways, as well as their interaction networks in OS, LF, and SF of Chuanzhong black goats at the estrus phase using small RNA-sequencing. We found that the miRNA expression profiles of LF and SF were more similar than those of OS—32, 16, and 29 differentially expressed miRNAs were identified in OS vs. LF, OS vs. SF, and LF vs. SF, respectively. Analyses of functional enrichment and the miRNA-targeted gene interaction network suggested that miR-182 (SMC3), miR-122 (SGO1), and miR-206 (AURKA) were involved in ovarian organogenesis and hormone secretion by oocyte meiosis. Furthermore, miR-202-5p (EREG) and miR-485-3p (FLT3) were involved in follicular maturation through the MAPK signaling pathway, and miR-2404 (BMP7 and CDKN1C) played a key role in follicular development through the TGF-β signaling pathway and cell cycle; nevertheless, further research is warranted. To our knowledge, this is the first study to investigate miRNA expression patterns in OS, LF, and SF of Chuanzhong black goats during estrus. Our findings provide a theoretical basis to elucidate the role of miRNAs in follicular maturation. These key miRNAs might provide candidate biomarkers for the diagnosis of follicular maturation and will assist in developing new therapeutic targets for female goat infertility.

Unconservative_15_2570409 suppresses progesterone receptor expression in the granulosa cells of Hu sheep

Female fertility potential depends on the number of mature follicles; however, the underlying molecular mechanisms remain unclear. Based on previously generated miRNA and mRNA sequencing data of healthy ovarian follicles (>5 mm in diameter) isolated from Hu sheep with different prolificacy, we investigated the roles of a novel miRNA (unconservative_15_2570409) and the progesterone receptor (PGR) gene in follicular development. During the periovulatory phase, the expression of unconservative_15_2570409 and PGR was lower and higher, respectively, in the >5 mm follicles of high prolificacy (HP) ewes than in those of low prolificacy (LP) ewes and in the >3 mm follicles than in the smaller follicles of the HP ewes. Subsequently, the granulosa cells (GCs) of Hu sheep were used as an in vitro model. PGR overexpression significantly increased the mRNA expression of steroidogenic acute regulatory protein (StAR), 3-beta-hydroxysteroid dehydrogenase/isomerase (3β-HSD), and cytochrome P450 family 19 subfamily A member 1 (CYP19A1), which promoted the secretion of progesterone (P₄) and estradiol (E₂). PGR knockdown significantly downregulated the levels of StAR and 3β-HSD mRNA and decreased the production of P₄, whereas no effects on CYP19A1 mRNA expression and E₂ levels were observed. We also found the negative regulatory effect of unconservative_15_2570409 on the mRNA and protein expression of PGR by targeting the 3'-untranslated region. The regulation of PGR levels resulted in a corresponding change in the ADAMTS1, PPAR-γ, and CTSL gene transcripts, which are important for follicular maturation and ovulation. Additionally, PGR, ADAMTS1, and PPAR-γ were predominantly localized in the GCs. Collectively, our results suggest that by regulating PGR expression and consequently affecting the expression of target genes and steroidogenesis, unconservative_15_2570409 plays a role in follicular development during the periovulatory stage, which provides novel insights into the molecular mechanisms affecting Hu sheep prolificacy.

The diagnostic, prognostic and therapeutic potential of circulating microRNAs in ovarian cancer

Ovarian cancer (OC) is often diagnosed at an advanced stage because of the late onset of symptoms, and this together with the lack of effective treatments, has meant it is associated with a very high mortality. The aberrant expression of MicroRNA (miRNA) contributes to the initiation and development of human tumors including OC. Several miRNAs are secreted by tumor cells and can be identified in body fluids. Serum miRNAs levels are associated with several clinical conditions, and may be used to predict prognosis and response to treatments in some cancers including OC. This review summarizes the current progresses regarding the potential applications of circulating miRNA as innovative biomarkers in OC.

Comprehensive analysis of mRNAs and miRNAs in the ovarian follicles of uniparous and multiple goats at estrus phase

Background

Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique.

Results

The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells.

Conclusions

This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.

The expression of miR-129-5p and its target genes in the skin of goats

Coat color is one of the major quality traits of animals, and miR-129-5p acts as an important regulator for melanin biosynthesis in mammals. In this study, real-time PCR and western blotting were used to examine the expression of miR-129-5p and its targets genes in the skin of different coat color goats. The results showed that the expression of miR-129-5p in the skin samples of Inner Mongolia cashmere goats (IMCG) was higher than that of Dazu black goat (DBG). Also, the target genes (tyrosinase (TYR), frizzled 6 (FZD6) and glycogen synthase kinase 3β (GSK3β)) of miR-129-5p was highly expressed in the skin samples of DBG. The expression of miR-129-5p firstly increased and then decreased with age in F1 hybrid generation of DBG and IMCG. In addition, the expression of TYR decreased with age, while the expression of MITF increased with age but then decreased. The expression of FZD6 and GSK3β in the skin samples of F1 of different ages were irregular. Our results indicated that miR-129-5p mainly affects the formation of coat color of goats by decreasing the expression of TYR. This study suggests that miR-129-5p can act as a suppressor in the formation of coat color to lay the foundation for studying the effect of miR-129-5p on melanin synthesis.

Effects of Decreasing Dietary Crude Protein Level on Growth Performance, Nutrient Digestion, Serum Metabolites, and Nitrogen Utilization in Growing Goat Kids (Capra hircus)

Simple Summary

Reducing the dietary protein content could potentially reduce losses of nitrogen from ruminant farms and mitigate pressure on the protein ingredient supply. However, there is little information in the literature on the effect of low-protein diets in growing Anhui white goat kids. We demonstrated that decreasing the dietary crude protein level in Anhui white goat kids affected growth performance, improved nitrogen utilization, and reduced environmental nitrogen pollution. The key finding of this study was that a diet containing 13.4% crude protein supplied adequate protein to improve nitrogen utilization in white goat kids without any adverse effect on growth performance.

Abstract

The effects of decreasing dietary crude protein (CP) level on growth performance, nutrient digestion, serum metabolites, and nitrogen utilization in growing goat kids were investigated in the current study. Thirty-six male Anhui white goat kids were randomly assigned to one of three CP content diets: 14.8% (control), 13.4%, and 12.0% of dry matter, respectively. Diets were isoenergetic. The experiment lasted for 14 weeks, with the first two weeks being for adaptation. Results showed that the low-CP diet decreased average daily gain, feed efficiency, digestibility of dry matter, organic matter, crude protein, and fiber. No significant changes were observed in dry-matter intake. With a decrease in dietary CP level, fecal nitrogen excretion (% of nitrogen intake) increased linearly, whereas CP intake, blood urea nitrogen, urinary nitrogen excretion (% of nitrogen intake), and total nitrogen excretion (% of nitrogen intake) decreased. Serum glucose concentration decreased, while concentrations of low-density lipoproteins and non-esterified fatty acids increased with the low-CP diet. In conclusion, decreasing the dietary CP level decreased goats’ nitrogen excretion, but with restrictive effects on growth performance. A diet containing 13.4% CP is optimal for reducing nitrogen excretion without any adverse effect on growth performance of Anhui white goat kids. This concentration is 1.4% points lower than the NRC recommendations and thus is also environmentally beneficial on the input side because it decreases the use of feed (soy) protein.

Identification and characterization of a specific 13-miRNA expression signature during follicle activation in the zebrafish ovary

Ovarian folliculogenesis is always of great interest in reproductive biology. However, the molecular mechanisms that control follicle development, particularly the early phase of follicle activation or recruitment, still remain poorly understood. In an attempt to decipher the gene networks and signaling pathways involved in such transition, we conducted a transcriptomic analysis (RNA-seq) on zebrafish primary growth (PG, stage I; inactive) and previtellogenic (PV, stage II; activated) follicles. A total of 118 unique microRNAs (miRNAs) (11 downregulated and 83 upregulated during PG/PV transition) and 56711 unique messenger RNAs (mRNAs) (1839 downregulated and 7243 upregulated during PG/PV transition) were identified. Real-time quantitative polymerase chain reaction analysis confirmed differential expression of 46 miRNAs from 66 candidates (66.67%). Among which, we chose to focus on 13 miRNAs (let-7a, -7b, -7c-5p, -7d-5p, -7h, -7i; miR-21, -23a-3p, -27c-3p, -107a-3p, -125b-5p, -145-3p, and -202-5p) that exhibited significant differential expression between PG and PV follicles (P ≤ 0.045*). With this 13-miRNA expression signature alone, PG follicles can be well differentiated from PV follicles by hierarchical clustering, suggesting their functional relevance during PG-to-PV transition. By overlaying predicted target genes and the differentially expressed mRNAs revealed by the RNA-seq analysis, especially those showing reciprocal miRNA-mRNA expression patterns, we shortlisted a panel of miRNA downstream targets for luciferase reporter validation. The reporter assay confirmed the interactions of let-7i:: atg4a (P = 0.01*), miR-202-5p::c23h20orf24 (P = 0.0004***), and miR-144-5p::ybx1 (P = 0.003**), implicating these potential miRNA-mRNA gene pairs in follicle activation during folliculogenesis. Our transcriptomic data analyses suggest that miRNA-mediated post-transcriptional control may represent an important mechanism underlying follicle activation.

Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq

Background

A multitude of genes tightly regulate ovarian follicular development and hormone secretion. These complex and coordinated biological processes are altered during pregnancy. In order to further understand the regulatory role of these genes during pregnancy, it is important to screen the differentially expressed genes (DEGs) in the ovaries of pregnant and nonpregnant mammals. To detect the genes associated with the development of pregnancy in goats, DEGs from the ovaries from pregnant and nonpregnant Anhui white goats (pAWGs and nAWGs, respectively) were analyzed using RNA sequencing technology (RNA-Seq).

Results

In this study, 13,676,394 and 13,549,560 clean reads were generated from pAWGs and nAWGs, respectively, and 1724 DEGs were identified between the two libraries. Compared with nAWGs, 1033 genes were upregulated and 691 genes were downregulated in pAWGs, including PGR, PRLR, STAR and CYP19A1, which play important roles in goat reproduction. Gene Ontology analysis showed that the DEGs were enriched for 49 functional GO terms. Kyoto Encyclopedia of Genes and Genomes analysis revealed that 397 DEGs were significantly enriched in 13 pathways, including “cell cycle”, “cytokine–cytokine receptor interaction” and “steroid biosynthesis”, suggesting that the genes may be associated with cell cycle regulation, follicular development and hormone secretion to regulate the reproduction process. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data.

Conclusions

The data obtained in this work enrich the genetic resources of goat and provide a further understanding of the complex molecular regulatory mechanisms occurring during the development of pregnancy and reproduction in goats. The DEGs screened in this study may play an important role in follicular development and hormone secretion and they would provide scientific basis for related research in the future.

Electronic supplementary material

The online version of this article (10.1186/s40709-019-0095-9) contains supplementary material, which is available to authorized users.

Non-coding RNA regulation in reproduction: Their potential use as biomarkers

Non-coding RNAs (ncRNAs) are crucial regulatory elements in most biological processes and reproduction is also controlled by them. The different types of ncRNAs, as well as the high complexity of these regulatory pathways, present a complex scenario; however, recent studies have shed some light on these questions, discovering the regulatory function of specific ncRNAs on concrete reproductive biology processes. This mini review will focus on the role of ncRNAs in spermatogenesis and oogenesis, and their potential use as biomarkers for reproductive diseases or for reproduction success.

Analyses of circRNA profiling during the development from pre-receptive to receptive phases in the goat endometrium

BACKGROUND

Recent studies have revealed that noncoding RNAs play important regulatory roles in the formation of endometrial receptivity. Circular RNAs (circRNAs) are a universally expressed noncoding RNA species that have been recently proposed to act as miRNA sponges that directly regulate expression of target genes or parental genes.

RESULTS

We used Illumina Solexa technology to analyze the expression profiles of circRNAs in the endometrium from three goats at gestational day 5 (pre-receptive endometrium, PE) and three goats at gestational day 15 (receptive endometrium, RE). Overall, 21,813 circRNAs were identified, of which 5,925 circRNAs were specific to the RE and 9,078 were specific to the PE, which suggested high stage-specificity. Further analysis found 334 differentially expressed circRNAs in the RE compared with PE (P < 0.05). The analysis of the circRNA-miRNA interaction network further supported the idea that circRNAs act as miRNA sponges to regulate gene expression. Moreover, some circRNAs were regulated by estrogen (E2)/progesterone (P4) in endometrial epithelium cell lines (EECs) and endometrial stromal cell line (ESCs), and each circRNA molecule exhibited unique regulation characteristics with respect to E2 and P4.

CONCLUSIONS

These data provide an endometrium circRNA expression atlas corresponding to the biology of the goat receptive endometrium during embryo implantation.

Identification of key pathways and candidate genes in pancreatic ductal adenocarcinoma using bioinformatics analysis

Pancreatic ductal adenocarcinoma (PDAC) is a malignant tumor with a high degree of malignancy that is difficult to diagnose and treat. The present study integrated PDAC cohort profile datasets to identify key candidate genes and pathways involved in the pathogenesis of the disease. The expression profiles of GSE28735 included 45 PDCA and matching pairs of adjacent non-tumor tissue. Differentially expressed genes (DEGs) were sorted and candidate genes and pathway enrichment were analyzed. A DEG-associated protein-protein interaction (PPI) network was constructed. A total of 424 DEGs were identified in PDAC, including 159 upregulated genes and 265 downregulated genes. Gene Ontology analysis results indicated that upregulated DEGs were significantly enriched in biological process, molecular function and cellular component categories. Kyoto Encyclopedia of Genes and Genomes pathway analysis demonstrated that the upregulated DEGs were enriched in ‘pancreatic secretion’, ‘protein digestion’ and ‘absorption’. Downregulated DEGs were enriched in ‘ECM-receptor interaction’, ‘focal adhesion’ and ‘PI3K/AKT’ signaling pathways. The PPI network revealed that these genes were involved in significant pathways, including ‘ECM organization’ signaling pathways (Hippo signaling pathway, TGF-β signaling pathway, Hedgehog signaling pathway and Wnt signaling pathway), ‘serine-type peptidase activity’ signaling pathway (PI3K-Akt signaling pathway, TNF-α signaling pathway and Wnt signaling pathway) and ‘extracellular region’ signaling pathways (RTP signaling pathway, G protein-coupled receptor signaling pathway and RAS-RAF-MAPK signaling pathway). The identification of these candidate genes and pathways sheds light on the etiology and molecular mechanisms of PDAC and may guide the development of novel therapies for pancreatic cancer.

MicroRNAs in ovarian follicular atresia and granulosa cell apoptosis

MicroRNAs (miRNAs) are short, noncoding RNAs that posttranscriptionally regulate gene expression. In the past decade, studies on miRNAs in ovaries have revealed the key roles of miRNAs in ovarian development and function. In this review, we first introduce the development of follicular atresia research and then summarize genome-wide studies on the ovarian miRNA profiles of different mammalian species. Differentially expressed miRNA profiles during atresia and other biological processes are herein compared. In addition, current knowledge on confirmed functional miRNAs during the follicular atresia process, which is mostly indicated by granulosa cell (GC) apoptosis, is presented. The main miRNA families and clusters, including the let-7 family, miR-23-27-24 cluster, miR-183-96-182 cluster and miR-17-92 cluster, and related pathways that are involved in follicular atresia are thoroughly summarized. A deep understanding of the roles of miRNA networks will not only help elucidate the mechanisms of GC apoptosis, follicular development, atresia and their disorders but also offer new diagnostic and treatment strategies for infertility and other ovarian dysfunctions.

Identification of microRNAs in granulosa cells from patients with different levels of ovarian reserve function and the potential regulatory function of miR-23a in granulosa cell apoptosis

This study aimed to determine the microRNA (miRNA) profiles in granulosa cells (GCs) from the follicular fluid (FF) of patients with varying levels of ovarian reserve function. We included 45 women undergoing in vitro fertilization (IVF) and intracytoplasmic sperm injection (ICSI) treatment. After collecting GCs from each patient, total RNA was extracted from 12 samples. Using Illumina/deep-sequencing technology, we analyzed the small RNAs in each group. Using the R package, we identified the differentially expressed (DE) miRNAs among patients with varying levels of ovarian reserve function. We identified 20 conserved and 3 novel miRNAs that were upregulated in the poor ovarian response (POR) group and 30 conserved miRNAs and 1 novel miRNA that were upregulated in the polycystic ovary syndrome (PCOS) group. Bioinformatics analysis revealed complementary pairing between miR-23a and the 3'-untranslated region (UTR) of the Sirt1 mRNA. miR-23a can regulate SIRT1 protein expression at the posttranscriptional level in GCs. Overexpressing miR-23a can inhibit the expression of SIRT1, decrease the stimulatory effect of SIRT1 on the ERK1/2 pathway, inhibit the expression of p-ERK1/2, and increase apoptosis in GCs. Previous studies confirmed that miR-23a targets SIRT1 and promotes apoptosis in GCs by inhibiting the ERK1/2 signaling pathway. This study provides a novel perspective regarding the role of miRNAs in the regulation of human GC apoptosis in vitro.

Comparative profiling of differentially expressed microRNAs in estrous ovaries of Kazakh sheep in different seasons

Seasonal estrus is a critical limiting factor for animal fecundity. However, estrus occurs in some seasonally estrous sheep in the non-breeding season, and this phenomenon involves changes in ovarian biology. Previous studies indicated that small RNAs, such as microRNAs (miRNAs), play important regulatory roles in ovarian biology. Differentially expressed miRNAs in the ovaries of estrous sheep were identified using Solexa sequencing technology. A total of 423 known miRNAs were identified in ovaries of estrous sheep in the breeding season and non-breeding season. In the comparison of these two groups, 48 miRNAs were identified that were differentially expressed between the two groups (including 5 up-regulated and 43 down-regulated miRNAs). KEGG pathway analysis revealed that the target genes of some differentially expressed miRNAs were involved in pathways related to reproductive hormone signaling and follicular development. Furthermore, the levels of estradiol (E2), progesterone (P4), luteinizing hormone (LH) and follicle-stimulating hormone (FSH) were lower in anestrus sheep than in sheep during the breeding season. Upon combining the pathway enrichment analysis, target gene expression and hormone measurement results, we suggest that these differentially expressed miRNAs might influence ovarian activity in the non-breeding season by regulating the above pathways. The identification of miRNAs that are differentially expressed between ovines in the breeding season and non-breeding season will contribute to our understanding of the role of miRNAs in estrus regulation, and these data may provide a basis for regulating estrus in sheep during the non-breeding season.

MicroRNAs: tiny molecules with a significant role in mammalian follicular and oocyte development

The genetic regulation of female fertility (follicular development, oocyte maturation and early preimplantation embryo development) involves the spatio-temporal regulation of those genes that play key roles in various stages of the female reproductive axis. MicroRNAs (miRNAs), a class of small non-coding RNAs, are known to regulate the expression of a large proportion of such genes. In recent decades, multiple studies have aimed to determine the roles of these non-coding RNAs in mammalian follicular development, oocyte growth and embryo development. These studies have applied a variety of approaches, including conditional knockout of miRNA biogenesis genes, high-throughput sequencing technologies for pattern recognition in miRNA expression and loss- and gain-of-function of miRNAs in various animal models. In addition to the cellular miRNAs, a large variety of RNAs are found in circulation, being coupled with extracellular vesicles, proteins and lipids. Because of their potential as diagnostic markers for abnormal physiologies, there is increasing interest in the identification of extracellular miRNAs in various biological fluids and spent in vitro culture media. This review focuses on studies addressing the expression and potential role of cellular and extracellular miRNAs in mammalian follicular cell physiology and subsequent ovarian functionality and oocyte maturation.

Integrated analysis of mRNA and viral miRNAs in the kidney of Carassius auratus gibelio response to cyprinid herpesvirus 2

MicroRNAs (miRNAs) are small, non-coding single stranded RNAs that play crucial roles in numerous biological processes. Vertebrate herpesviruses encode multiple viral miRNAs that modulate host and viral genes. However, the roles of viral miRNAs in lower vertebrates have not been fully determined. Here, we used high-throughput sequencing to analyse the miRNA and mRNA expression profiles of Carassius auratus gibelio in response to infection by cyprinid herpesvirus 2 (CyHV-2). RNA sequencing obtained 26,664 assembled transcripts, including 2,912 differentially expressed genes. Based on small RNA sequencing and secondary structure predictions, we identified 17 CyHV-2 encoded miRNAs, among which 14 were validated by stem-loop quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) and eight were validated by northern blotting. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of miRNAs-mRNA pairs revealed diverse affected immune signalling pathways, including the RIG-I-like receptor and JAK-STAT pathways. Finally, we presented four genes involved in RIG-I-like pathways, including host gene IRF3, RBMX, PIN1, viral gene ORF4, which are negatively regulated by CyHV-2 encoded miRNA miR-C4. The present study is the first to provide a comprehensive overview of viral miRNA-mRNA co-regulation, which might have a key role in controlling post-transcriptomic regulation during CyHV-2 infection.

Genome-wide identification of microRNAs responsive to Ectropis oblique feeding in tea plant (Camellia sinensis L.)

The tea plant (Camellia sinensis L.) is vulnerable to the geometrid Ectropis oblique; although microRNAs (miRNAs) are important for plant growth, development and stress response, the function of miRNAs in the response of C. sinensis to stress from E. oblique is unclear. To identify E. oblique stress-responsive miRNAs and their target genes in tea plant, three small RNA libraries were constructed from leaves subjected to mechanical wounding (MW), geometrid attack (GA) and from healthy control (CK) leaves. Using high-throughput sequencing, 130 known miRNAs and 512 novel miRNAs were identified; of these, differential expression under GA stress was observed for 36 known and 139 novel miRNAs. Furthermore, 169 GA-responsive and 173 MW-responsive miRNAs were detected by miRNA microarray. The expression patterns of six GA-responsive miRNAs were validated by qRT-PCR. Several target genes for these miRNAs encode various transcription factors, including ethylene-responsive transcription factors and squamosa promoter-binding-like proteins, which suggests that these miRNAs may regulate stress-responsive transcriptional processes in tea plant. The present findings provide novel insights into miRNA-mediated regulatory mechanisms underlying the response to GA stress, and also offer valuable information for development of pest resistance using RNA interference-based strategies in tea plants.

Identification and analysis of differentially expressed long non-coding RNAs between multiparous and uniparous goat (Capra hircus) ovaries

Long non-coding RNAs (lncRNAs) play important roles in almost all biological processes. However, there is little information on the effects of lncRNAs on ovulation and lambing rates. In the present study, we used high-throughput RNA sequencing to identify differentially expressed lncRNAs between the ovaries of multiparous (Mul) and uniparous (Uni) Anhui White goats. Among the 107,255,422 clean reads, 183,754 lncRNAs were significantly differentially expressed between the Uni and Mul. Among them, 455 lncRNAs were co-expressed between the two samples, whereas, 157,523 lncRNAs were uniquely expressed in the Uni, and 25,776 uniquely lncRNAs were expressed in the Mul. Through Cis role analysis, 24 lncRNAs were predicted to overlap with cis-regulatory elements, which involved in Progesterone-mediated oocyte maturation, Steroid biosynthesis, Oocyte meiosis, and gonadotropin-releasing hormone (GnRH) signaling pathway. These 4 pathways were related to ovulation, and the KEGG pathway analysis on target genes of the differentially expressed lncRNAs confirmed this results. In addition, 10 lncRNAs harbored precursors of 40 miRNAs, such as TCONS_00320849 related to a mature miRNA sequence, miR-365a, which was reported to be related to proliferation, were annotated in the precursor analysis of miRNAs. The present expand the understanding of lncRNA biology and contribute to the annotation of the goat genome. The study will provide a resource for lncRNA studies of ovulation and lambing.

The developmental transcriptome landscape of receptive endometrium during embryo implantation in dairy goats

Under natural conditions, some embryos cannot implant successfully because of the dysfunction of receptive endometrium (RE). Thus, it is imperative for us to study the molecular mechanisms involved in the formation of the RE from pre-receptive endometrium (PE). In this study, the endometrium from gestational day 5 (D5, PE) and gestational day 15 (D15, RE) dairy goats were selected to systematically analyze the transcriptome using strand-specific Ribo-Zero RNA-Seq, >120 million high-quality paired-end reads were generated and 47,616 transcripts were identified in the endometrium of dairy goats. A total of 810 mRNAs were differentially expressed genes (DEGs) between the RE and PE meeting the criteria of P-values<0.05. Bioinformatics analysis of the DEGs revealed that a number of biological processes and pathways were potentially involved in the establishment of the RE, notably energy metabolism and amino acid metabolism. Furthermore, we speculated that CXCL14, IGFBP3, and LGALS15 potentially participated in the development of endometrium. What's more, putative SNPs, InDels and AS events were identified and analyzed in the endometrium. In a word, this resulting view of the transcriptome greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during the formation of receptive endometrium in dairy goats.

Identification and comparative analysis of the ovarian microRNAs of prolific and non-prolific goats during the follicular phase using high-throughput sequencing

The kidding rate is one of the most important economic traits for goat production, but the genetic mechanism that is associated with ovulation rate is poorly understood. Recently, increasing evidence has suggested that microRNAs (miRNAs) influence ovarian biological processes. The present study provides the first comparison of the ovarian miRNAs of prolific Jintang black goats (JTGs) and non-prolific Tibetan goats (TBGs) during the follicular phase using RNA-Seq technology. We generated 11.19 million (M) and 11.34 M clean reads from the TBG and JTG libraries, respectively, from which a total of 389 known miRNAs were identified and 142 novel miRNAs were predicted. A total of 191 miRNAs were differentially expressed between the two breeds. Among the 10 most abundant miRNAs, miR-21-5p was defined as differentially expressed miRNA with a higher level in the JTG library than in the TBG library, but the other miRNAs were not different between the breeds. The predicted miRNA-targeted genes were further analyzed by Gene Ontology and KEGG pathway analyses. The results revealed that miR-21, miR-99a, miRNA-143, let-7f, miR-493 and miR-200b may affect follicular development. These findings will increase the current understanding of the role of ovarian miRNAs in the regulation of ovulation rate in goats.

Identification of differentially expressed miRNAs through high-throughput sequencing in the chicken lung in response to Mycoplasma gallisepticum HS

Mycoplasma gallisepticum (MG) infects chickens, causes chronic respiratory diseases (CRD) and severely damages the poultry industry. It has been suggested that micro-ribonucleic acids (miRNAs) are involved in microbial pathogenesis. Here, we identified miRNAs that are associated with MG infection in chicken lungs at 3 and 10days post-infection by deep sequencing. Thirty-six down-regulated and 9 up-regulated miRNAs belonging to 31 miRNA families were detected at 3days post-infection, whereas 50 down-regulated and 18 up-regulated miRNAs belonging to 41 miRNA families were found at 10days post-infection. The 45 and 68 differentially expressed miRNAs at 3 and 10days target 6280 and 7181 genes, respectively. In this study, 8 candidate novel chicken miRNAs were identified. Analyses via GO, KEGG, miRNA-GO-network, path-net and gene-net showed that these altered miRNAs might be involved in regulating the host response to MG infection by targeting genes in many pathways, such as the MAPK pathway, focal adhesion, Wnt pathway, endocytosis, Jak/STAT pathway, phosphatidylinositol pathway, adherens junctions, regulation of actin cytoskeleton among others. These analyses indicate that the MAPK pathway may be a key regulatory route. Also, the miR-8 family, miR-499 family, miR-17 family, and PIK3 family genes, as well as the MAP2K1 and RAC1 genes, might be important in MG infection. miR-20 of the miR-17 family was further confirmed by RT-qPCR. The important miRNAs, mRNAs and pathways associated with MG infection in chicken are valuable for further research. Our data provide new insights into the mechanism of these miRNAs on the regulation of host-MG interactions.

MicroRNA Mediating Networks in Granulosa Cells Associated with Ovarian Follicular Development

Ovaries, which provide a place for follicular development and oocyte maturation, are important organs in female mammals. Follicular development is complicated physiological progress mediated by various regulatory factors including microRNAs (miRNAs). To demonstrate the role of miRNAs in follicular development, this study analyzed the expression patterns of miRNAs in granulosa cells through investigating three previous datasets generated by Illumina miRNA deep sequencing. Furthermore, via bioinformatic analyses, we dissected the associated functional networks of the observed significant miRNAs, in terms of interacting with signal pathways and transcription factors. During the growth and selection of dominant follicles, 15 dysregulated miRNAs and 139 associated pathways were screened out. In comparison of different styles of follicles, 7 commonly abundant miRNAs and 195 pathways, as well as 10 differentially expressed miRNAs and 117 pathways in dominant follicles in comparison with subordinate follicles, were collected. Furthermore, SMAD2 was identified as a hub factor in regulating follicular development. The regulation of miR-26a/b on smad2 messenger RNA has been further testified by real time PCR. In conclusion, we established functional networks which play critical roles in follicular development including pivotal miRNAs, pathways, and transcription factors, which contributed to the further investigation about miRNAs associated with mammalian follicular development.

MicroRNAs: New Insight in Modulating Follicular Atresia: A Review

Our understanding of the post-transcriptional mechanisms involved in follicular atresia is limited; however, an important development has been made in understanding the biological regulatory networks responsible for mediating follicular atresia. MicroRNAs have come to be seen as a key regulatory actor in determining cell fate in a wide range of tissues in normal and pathological processes. Profiling studies of miRNAs during follicular atresia and development have identified several putative miRNAs enriched in apoptosis signaling pathways. Subsequent in vitro and/or in vivo studies of granulosa cells have elucidated the functional role of some miRNAs along with their molecular pathways. In particular, the regulatory roles of some miRNAs have been consistently observed during studies of follicular cellular apoptosis. Continued work should gradually lead to better understanding of the role of miRNAs in this field. Ultimately, we expect this understanding will have substantial benefits for fertility management at both the in vivo or/and in vitro levels. The stable nature of miRNA holds remarkable promise in clinical use as a diagnostic tool and in reproductive medicine to solve the ever-increasing fertility problem. In this review, we summarize current knowledge of the involvement of miRNAs in follicular atresia, discuss the challenges for further work and pinpoint areas for future research.

Effect of dietary Astragalus Polysaccharide supplements on testicular miRNA expression profiles and enzymatic changes of breeder cocks

Astragalus Polysaccharide (APS) is an important feed additive due to its immunomodulatory functions. Previous studies have proven that miRNAs play important roles in posttranscriptional gene regulation. Our goals were to identify differentially expressed miRNAs in testes in responses to APS dietary supplements and to find the effects of APS on breeder cock testes. We measured several enzymatic activities in testes and sperm samples and further generated miRNA expression profiles of testes from breeder cocks fed with control diets and extra APS. As a result, we found APS could increase testicular functional activities of marker enzymes. Meanwhile, there were 16 up-regulated and 17 down-regulated miRNAs in APS group, compared with the control group meeting the criteria of P-values < 0.05. Meanwhile, twelve differentially expressed miRNAs were validated by Mir-XTM miRNA RT-qPCR. Further GO and KEGG analyses of target genes for differentially expressed miRNAs revealed that some miRNAs may be involved in testicular nutrient metabolisms and NK cell mediated cytotoxicity pathway. Moreover, the effect of dietary APS supplements on NK cell mediated cytotoxicity pathway was also validated by RT-qPCR. Our results provided a novel insight into the effect of dietary APS supplements on testicular miRNA expression profiles and enzymatic changes of breeder cocks.

Comparative profiling of differentially expressed microRNAs between the follicular and luteal phases ovaries of goats

To explore if the regulation at post-transcriptional level of follicular phase (Fols) to luteal phase (Luts) transition occurs in the ovaries of Anhuai goats, the differentially expressed microRNAs (miRNAs) of ovaries in the Fols and Luts were analyzed using Solexa sequencing in the study. In total, 320 known miRNAs were co-expressed in the two phases, 339 and 353 known miRNAs were expressed in the ovary in the Fols and Luts, respectively. In addition, 45 novel miRNAs were co-expressed in the two phases, 70 and 94 novel miRNAs were expressed in the ovary in the Fols and Luts, respectively. Let-7f was the highest expressed significantly different known miRNA in the two phases, and mir-159 was the highest expressed significantly different novel miRNA in the two phases, which may participate in the follicular-luteal transition of Anhuai goats. GO annotation and KEGG pathway analysis were applied to analyze the target genes of differentially expressed miRNAs detected in the two phases. The results will help to further understand the role of miRNAs in the regulation of follicular to luteal transition in goat ovaries.

High-Throughput Sequencing Reveals Circulating miRNAs as Potential Biomarkers for Measuring Puberty Onset in Chicken (Gallus gallus)

There are still no highly sensitive and unique biomarkers for measurement of puberty onset. Circulating miRNAs have been shown to be promising biomarkers for diagnosis of various diseases. To identify circulating miRNAs that could be served as biomarkers for measuring chicken (Gallus gallus) puberty onset, the Solexa deep sequencing was performed to analyze the miRNA expression profiles in serum and plasma of hens from two different pubertal stages, before puberty onset (BO) and after puberty onset (AO). 197 conserved and 19 novel miRNAs (reads > 10) were identified as serum/plasma-expressed miRNAs in the chicken. The common miRNA amounts and their expression changes from BO to AO between serum and plasma were very similar, indicating the different treatments to generate serum and plasma had quite small influence on the miRNAs. 130 conserved serum-miRNAs were showed to be differentially expressed (reads > 10, P < 0.05) from BO to AO, with 68 up-regulated and 62 down-regulated. 4829 putative genes were predicted as the targets of the 40 most differentially expressed miRNAs (|log2(fold-change)|>1.0, P < 0.01). Functional analysis revealed several pathways that were associated with puberty onset. Further quantitative real-time PCR (RT-qPCR) test found that a seven-miRNA panel, including miR-29c, miR-375, miR-215, miR-217, miR-19b, miR-133a and let-7a, had great potentials to serve as novel biomarkers for measuring puberty onset in chicken. Due to highly conserved nature of miRNAs, the findings could provide cues for measurement of puberty onset in other animals as well as humans.

MicroRNA in ovarian function

The mammalian ovary is a dynamic organ. The coordination of follicle recruitment, selection, and ovulation and the timely development and regression of the corpus luteum are essential for a functional ovary and fertility. Deregulation of any of these processes results in ovarian dysfunction and potential infertility. MicroRNA (miRNA) are short noncoding RNA that regulate developmental processes and time-sensitive functions. The expression of miRNA in the ovary varies with cell type, function, and stage of the estrous cycle. miRNA are involved in the formation of primordial follicles, follicular recruitment and selection, follicular atresia, oocyte-cumulus cell interaction, granulosal cell function, and luteinization. miRNA are differentially expressed in luteal cells at the various stages of the estrous cycle and during maternal recognition of pregnancy, suggesting a role in luteal development, maintenance, and regression. An understanding of the patterns of expression and functions of miRNA in the ovary will lead to novel therapeutics to treat ovarian dysfunction and improve fertility and, potentially, to the development of better contraceptives.

Identification and profiling of microRNAs in goat endometrium during embryo implantation

Background

MicroRNAs (miRNAs) are short, highly conserved small noncoding RNAs that had fundamental roles in post-transcriptional gene expression, and they are crucial for proper control of biological processes and known to participate in embryo implantation. However, miRNA expression profiles in the pre-receptive and receptive phases of the goat endometrium during embryo implantation are unknown.

Results

A total of 1,069 and 847 miRNAs were expressed in receptive (R) and pre-receptive (P) goat endometrium, and 632 miRNAs were co-expressed in both phases. We identified 545 (50.98%) known miRNAs in the R library and 522 (61.63%) in the P library. There were 110 up-expressed miRNAs and 33 down-expressed miRNAs in receptive endometrium compared with the pre-receptive endometrium meeting the criteria of P-values< 0.05. Moreover, GO and KEGG analysis of the target genes of the differentially expressed miRNAs revealed some candidate miRNAs, genes and pathways that may involve in the formation of the receptive endometrium. Based on stem-loop RT-qPCR, 15 miRNAs were detected and the results suggested that the majority of the miRNA expression data measured by Solexa deep sequencing could represent actual miRNA expression levels.

Conclusions

Our data revealed the first miRNA profile related to the biology of the goat receptive endometrium during embryo implantation, and the results suggested that a subset of miRNAs might play important roles in the formation of endometrial receptivity. Thus, elucidating the physiological roles of endometrial miRNAs will help us better understand the genetic control of embryo implantation in goats.

Annotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches

Background

MicroRNAs (miRNA) are small endogenous non-coding RNA involved in the post-transcriptional regulation of specific mRNA targets. The first whole goat genome sequence became available in 2013, with few annotations. Our goal was to establish a list of the miRNA expressed in the mammary gland of lactating goats, thus enabling implementation of the goat miRNA repertoire and considerably enriching annotation of the goat genome.

Results

Here, we performed high throughput RNA sequencing on 10 lactating goat mammary glands. The bioinformatic detection of miRNA was carried out using miRDeep2 software. Three different methods were used to predict, quantify and annotate the sequenced reads. The first was a de novo approach based on the prediction of miRNA from the goat genome only. The second approach used bovine miRNA as an external reference whereas the last one used recently available goat miRNA. The three methods enabled the prediction and annotation of hundreds of miRNA, more than 95% were commonly identified. Using bovine miRNA, 1,178 distinct miRNA were detected, together with the annotation of 88 miRNA for which corresponding precursors could not be retrieved in the goat genome, and which were not detected using the de novo approach or with the use of goat miRNA. Each chromosomal coordinate of the precursors determined here were generated and depicted on a reference localisation map. Forty six goat miRNA clusters were also reported. The study revealed 263 precursors located in goat protein-coding genes, amongst which the location of 43 precursors was conserved between human, mouse and bovine, revealing potential new gene regulations in the goat mammary gland. Using the publicly available cattle QTL database, and cow precursors conserved in the goat and expressed in lactating mammary gland, 114 precursors were located within known QTL regions for milk production and composition.

Conclusions

The results reported here represent the first major identification study on miRNA expressed in the goat mammary gland at peak lactation. The elements generated by this study will now be used as references to decipher the regulation of miRNA expression in the goat mammary gland and to clarify their involvement in the lactation process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1471-y) contains supplementary material, which is available to authorized users.

Characterization and comparative profiling of ovarian microRNAs during ovine anestrus and the breeding season

Background

Seasonal estrus is a critical limiting factor of animal fecundity, and it involves changes in both ovarian biology and hormone secretion in different seasons. Previous studies indicate that two classes of small RNAs (miRNAs and piRNAs) play important regulatory roles in ovarian biology. To understand the roles of small RNA-mediated post-transcriptional regulation in ovine seasonal estrus, the variation in expression patterns of ovarian small RNAs during anestrus and the breeding season were analyzed using Solexa sequencing technology. In addition, reproductive hormone levels were determined during ovine anestrus and the breeding season.

Results

A total of 483 miRNAs (including 97 known, 369 conserved and 17 predicated novel miRNAs), which belong to 183 different miRNA families, were identified in ovaries of Tan sheep and Small Tail Han (STH) sheep. Compared with the three stages of the breeding season, 25 shared significantly differentially expressed (including 19 up- and six down-regulated) miRNAs were identified in ovine anestrus. KEGG Pathway analysis revealed that the target genes for some of the differentially expressed miRNAs were involved in reproductive hormone related pathways (e.g. steroid biosynthesis, androgen and estrogen metabolism and GnRH signaling pathway) as well as follicular/luteal development related pathways. Moreover, the expression of the differentially expressed miRNAs and most of their target genes were negatively correlated in the above pathways. Furthermore, the levels of estrogen, progesterone and LH in ovine anestrus were significantly lower than those in the breeding season. Combining the results of pathway enrichment analysis, expression of target genes and hormone measurement, we suggest that these differentially expressed miRNAs in anestrus might participate in attenuation of ovarian activity by regulating the above pathways. Besides miRNAs, a large and unexpectedly diverse set of piRNAs were also identified.

Conclusions

The miRNA profiles of ovine ovaries in anestrus were presented for the first time. The identification and characterization of miRNAs that are differentially expressed between ovine anestrus and the breeding season will help understanding of the role of miRNAs in the regulation of seasonal estrus, and provides candidates for determining miRNAs which could be potentially used to regulate ovine seasonal estrus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-899) contains supplementary material, which is available to authorized users.