Identification and differential expression of microRNAs in the ovaries of pigs (Sus scrofa) with high and low litter sizes

Comparative miRNA expression profile analysis of porcine ovarian follicles: new insights into the initiation mechanism of follicular atresia

Follicular atresia occurs in every stage of ovarian development, which is relevant to female fertility. In the past decade, increasing studies have confirmed that miRNAs, a class of short non-coding RNAs, play an important role in follicular atresia by post-transcription regulation of their target genes. However, the function of miRNAs on follicular atresia initiation is unknown. In the present study, high-throughput small RNA sequencing was performed to analyze differential miRNA expression profiles between healthy (HF) follicles and early atretic (EAF) follicles. A total of 237 conserved miRNA were detected, and the miR-143 is the highest expressed in follicles. Meanwhile, we also found wide sequence variations (isomiRs) in porcine ovarian miRNA, including in 5'un-translation region, core seed sequences and 3'untranslation region. Furthermore, we identified 22 differentially expressed miRNAs in EAF groups compared to HF group, of which 3 miRNAs were upregulated, as well as 19 miRNAs were downregulated, and then the RT-PCR was performed to validate these profiles. The target genes of these differentially expressed miRNAs were predicted by using miRwalk, miRDB, and Targetscan database, respectively. Moreover, the gene ontology and KEGG pathway enrichment established that the regulating functions and signaling pathways of these miRNAs contribute to follicular atresia initiation and cell fate. In conclusion, this study provides new insights into the changes of miRNAs in early atretic follicles to demonstrate their molecular regulation in ovarian follicular atretic initiation.

CeRNA Network Reveals the Circular RNA Characterization in Goat Ear Fibroblasts Reprogramming into Mammary Epithelial Cells

Circular RNAs (circRNAs) are a type of non-coding RNA that play a crucial role in the development and lactation of mammary glands in mammals. A total of 107 differentially expressed circRNAs (DE circRNAs) were found, of which 52 were up-regulated and 55 were down-regulated. We also found that DE circRNA host genes were mainly involved in GO terms related to the development process of mammary epithelial cells and KEGG pathways were mostly related to mammary epithelial cells, lactation, and gland development. Protein network analysis found that DE circRNAs can competitively bind to miRNAs as key circRNAs by constructing a circRNA-miRNA-mRNA network. CircRNAs competitively bind to miRNAs (miR-10b-3p, miR-671-5p, chi-miR-200c, chi-miR-378-3p, and chi-miR-30e-5p) involved in goat mammary gland development, mammary epithelial cells, and lactation, affecting the expression of core genes (CDH2, MAPK1, ITGB1, CAMSAP2, and MAPKAPK5). Here, we generated CiMECs and systematically explored the differences in the transcription profile for the first time using whole-transcriptome sequencing. We also analyzed the interaction among mRNA, miRNA, and cirRNA and predicted that circRNA plays an important role in the maintenance of mammary epithelial cells.

Identification of hub genes associated with follicle development in multiple births sheep by WGCNA

Sheep exhibit a distinct estrous cycle that includes four different phases: proestrus, estrus, late estrus, and luteal phase. As the estrous cycle repeats, follicular development regularly alternates. We thus investigated ovarian transcriptome data from each of the four phases using weighted gene co-expression network analysis (WGCNA) to identify modules, pathways, and genes essential to follicle growth and development. We clustered mRNA and long non-coding RNA (lncRNA) into different modules by WGCNA, and calculated correlation coefficients between genes and Stages of the estrous cycle. Co-expression of the black module (cor = 0.81, P<0.001) and the yellow module (cor = 0.61, P<0.04) was found to be critical for follicle growth and development. A total of 2066 genes comprising the black and yellow modules was used for functional enrichment. The results reveal that these genes are mainly enriched in Cell cycle, PI3K-Akt signaling pathway, Oocyte meiosis, Apoptosis, and other important signaling pathways. We also identified seven hub genes (BUB1B, MAD2L1, ASPM, HSD3B1, WDHD1, CENPA, and MXI1) that may play a role in follicle development. Our study may provide several important new markers allowing in depth exploration of the genetic basis for multiparous reproduction in sheep.

Totipotency of miR-184 in porcine granulosa cells

Estradiol (E₂) synthesis, cell proliferation and the apoptosis of porcine granulosa cells (GCs) affect follicular growth and development. The miR-184 level in ovary tissues of Yorkshire × Landrace sows was significantly higher in high-yielding sows than that in low-yielding sows, which was the same as in Yorkshire sows. However, the roles of miR-184 on E₂ granulosa cells (GCs) are still unclear. We found that miR-184 promoted E₂ synthesis and proliferation but inhibited apoptosis in GCs by targeting nuclear receptor subfamily 1 group D member 1 (NR1D1), cyclin dependent kinase inhibitor 1A (P21,CDKN1A) and homeodomain interacting protein kinase 2 (HIPK2) respectively. These findings indicated that miR-184 is a novel key factor that regulates the physiological functions of GCs.

Investigation of Sperm and Seminal Plasma Candidate MicroRNAs of Bulls with Differing Fertility and In Silico Prediction of miRNA-mRNA Interaction Network of Reproductive Function

Simple Summary

The objective of this study was to identify differentially expressed (DE) sperm and seminal plasma microRNAs (miRNAs) in high- and low-fertile Holstein bulls (four bulls per group), integrate miRNAs to their target genes, and categorize target genes based on predicted biological processes. Out of 84 bovine-specific, prioritized miRNAs analyzed by RT-PCR, 30 were differentially expressed in high-fertile sperm and seminal plasma compared to low-fertile sperm and seminal plasma, respectively (p ≤ 0.05, fold regulation ≥5 magnitudes). Interestingly, expression levels of DE-miRNAs in sperm and seminal plasma followed a similar pattern. Highly scored integrated genes of DE-miRNAs predicted various biological and molecular functions, cellular process, and pathways. Further in silico analysis revealed categorized genes may have a plausible association with pathways regulating sperm structure and function, fertilization, and embryo and placental development. In conclusion, highly DE-miRNAs in bovine sperm and seminal plasma could be used as a tool for predicting reproductive functions. Since the identified miRNA-mRNA interactions were mostly based on predictions from public databases, the causal regulations of miRNA-mRNA and the underlying mechanisms require further functional characterization in future studies.

Abstract

Recent advances in high-throughput in silico techniques portray experimental data as exemplified biological networks and help us understand the role of individual proteins, interactions, and their biological functions. The objective of this study was to identify differentially expressed (DE) sperm and seminal plasma microRNAs (miRNAs) in high- and low-fertile Holstein bulls (four bulls per group), integrate miRNAs to their target genes, and categorize the target genes based on biological process predictions. Out of 84 bovine-specific, prioritized miRNAs analyzed by RT-PCR, 30 were differentially expressed in high-fertile sperm and seminal plasma compared to low-fertile sperm and seminal plasma, respectively (p ≤ 0.05, fold regulation ≥ 5 magnitudes). The expression levels of DE-miRNAs in sperm and seminal plasma followed a similar pattern. Highly scored integrated genes of DE-miRNAs predicted various biological and molecular functions, cellular process, and pathways. Further, analysis of the categorized genes showed association with pathways regulating sperm structure and function, fertilization, and embryo and placental development. In conclusion, highly DE-miRNAs in bovine sperm and seminal plasma could be used as a tool for predicting reproductive functions. Since the identified miRNA-mRNA interactions were mostly based on predictions from public databases, the causal regulations of miRNA-mRNA and the underlying mechanisms require further functional characterization in future studies.

Comprehensive Transcriptome Analysis of Follicles from Two Stages of the Estrus Cycle of Two Breeds Reveals the Roles of Long Intergenic Non-Coding RNAs in Gilts

Simple Summary

This study provides new perspectives about the roles of lincRNAs in the estrus expression of gilts, which is correlated with ovarian steroid hormone and follicular development. Follicular tissues from two stages of the estrus cycle of Large White and Mi gilts were used for RNA-seq. Some genes and lincRNAs related to estrus expression in pigs were discovered. PPI and ceRNA networks related to the estrus expression were constructed. These results suggest that the estrus expression may be affected by lincRNAs and their target genes.

Abstract

Visible and long-lasting estrus expression of gilts and sows effectively sends a mating signal. To reveal the roles of Long Intergenic Non-coding RNAs (lincRNAs) in estrus expression, RNA-seq was used to investigate the lincRNAs expression of follicular tissues from Large White gilts at diestrus (LD) and estrus (LE), and Chinese Mi gilts at diestrus (MD) and estrus (ME). Seventy-three differentially expressed lincRNAs (DELs) were found in all comparisons (LE vs. ME, LD vs. LE, and MD vs. ME comparisons). Eleven lincRNAs were differentially expressed in both LD vs. LE and MD vs. ME comparisons. Fifteen DELs were mapped onto the pig corpus luteum number Quantitative Trait Loci (QTL) fragments. A protein–protein interaction (PPI) network that involved estrus expression using 20 DEGs was then constructed. Interestingly, three predicted target DEGs (PTGs) (CYP19A1 of MSTRG.10910, CDK1 of MSTRG.10910 and MSTRG.23984, SCARB1 of MSTRG.1559) were observed in the PPI network. A competitive endogenous RNA (ceRNA) network including three lincRNAs, five miRNAs, and five genes was constructed. Our study provides new insight into the lincRNAs associated with estrus expression and follicular development in gilts.

Mechanisms of Oogenesis-Related Long Non-coding RNAs in Porcine Ovaries Treated With Recombinant Pig Follicle-Stimulating Hormone

Reproductive efficiency is of significant importance in pork production for it has a great impact on economic success. Ovulation rate is an early component of reproduction efficiency of pigs, and it contributes to the upper limit of litter size. In this study, we used the newly developed recombinant pig follicle stimulating hormone (rpFSH) instead of traditional PMSG to increase ovulation rate of pigs in order to achieve higher litter size, for it was better at stimulating ovulation, and showed more cheaper and greener. However, relatively little is known about the underlying genetic bases and molecular mechanisms. Consequently, an experiment was carried out in ovaries of replacement gilts to screen the key genes and lncRNAs that affect the fecundity of pigs by RNA-seq technology. Twenty gilts were divided into two groups, including 10 rpFSH treatment pigs and 10 control animals. After slaughtering and collecting the phenotypic data, ovaries of five pigs in each group were selected for RNA-seq. Total RNA was extracted to construct the library and then sequence on an Illumina Hiseq 4000 system. A comprehensive analysis of mRNAs and long non-coding RNAs (lncRNAs) from 10 samples was performed with bioinformatics. The phenotypic data showed that rpFSH treatment groups had the higher (P < 0.01) ovarian weight and more mature follicles. The RNA-seq results showed that a total of 43,499 mRNAs and 21,703 lncRNAs were identified, including 21,300 novel lncRNAs and 403 known lncRNAs, of which 585 mRNAs and 398 lncRNAs (P < 0.05) were significantly differentially expressed (DE) between the two groups of rpFSH treatment group and controlled group. GO and KEGG annotation analysis indicated that the target genes of DE lncRNAs and DE mRNAs were related to prolactin receptor activity, mitophagy by induced vacuole formation, and meiotic spindle. Moreover, we found that NR5A2 (nuclear receptor subfamily 5, group A, member 2), a target gene of lncRNA MSTRG.3902.1, was involved in regulating follicular development, ovulation, and estrogen production. Our study provided a catalog of lncRNAs and mRNAs associated with ovulation of rpFSH treatment, and they deserve further study to deepen the understanding of biological processes in the regulation of ovaries of rpFSH treatment pigs.

Discovery of Differentially Expressed MicroRNAs in Porcine Ovaries With Smaller and Larger Litter Size

The number of live births in a litter is an important reproductive trait, and is one of the main indicators which reflect the production level and economic benefit of a pig farm. The ovary is an important reproductive organ of the sow, and it undergoes a series of biological processes during each estrous cycle. A complex transcriptional network containing coding and non-coding RNAs in the ovary closely regulates the reproductive capability of sows. However, the molecular regulation mechanisms affecting sow litter size are still unclear. We investigated the expression profiles of microRNAs (miRNAs) in porcine ovaries from sows with smaller than average litter sizes (SLS) and those with larger litter sizes (LLS). In total, 411 miRNAs were identified, and of these 17 were significantly down-regulated and 16 miRNAs were up-regulated when comparing sows with LLS and SLS, respectively. We further characterized the role of miR-183 which was one of the most up-regulated miRNAs. CCK-8, EdU incorporation and western blotting assays demonstrated that miR-183 promoted the proliferation of granulosa cells (GCs) in pig ovaries. Moreover, miR-183 inhibited the synthesis of estradiol in GCs and promoted the synthesis of progesterone. These results will help in gaining understanding of the role of miRNAs in regulating porcine litter size.

Insights Into Extracellular Vesicle/Exosome and miRNA Mediated Bi-Directional Communication During Porcine Pregnancy

Spontaneous fetal loss is one of the most important challenges that commercial pig industry is still facing in North America. Research over the decade provided significant insights into some of the associated mechanisms including uterine capacity, placental efficiency, deficits in vasculature, and immune-inflammatory alterations at the maternal-fetal interface. Pigs have unique epitheliochorial placentation where maternal and fetal layers lay in opposition without any invasion. This has provided researchers opportunities to accurately tease out some of the mechanisms associated with maternal-fetal interface adaptations to the constantly evolving needs of a developing conceptus. Another unique feature of porcine pregnancy is the conceptus derived recruitment of immune cells during the window of conceptus attachment. These immune cells in turn participate in pregnancy associated vascular changes and contribute toward tolerance to the semi-allogeneic fetus. However, the precise mechanism of how maternal-fetal cells communicate during the critical times in gestation is not fully understood. Recently, it has been established that bi-directional communication between fetal trophoblasts and maternal cells/tissues is mediated by extracellular vesicles (EVs) including exosomes. These EVs are detected in a variety of tissues and body fluids and their role has been described in modulating several physiological and pathological processes including vascularization, immune-modulation, and homeostasis. Recent literature also suggests that these EVs (exosomes) carry cargo (nucleic acids, protein, and lipids) as unique signatures associated with some of the pregnancy associated pathologies. In this review, we provide overview of important mechanisms in porcine pregnancy success and failure and summarize current knowledge about the unique cargo containing biomolecules in EVs. We also discuss how EVs (including exosomes) transfer their contents into other cells and regulate important biological pathways critical for pregnancy success.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Comparative profiling of the resistance of different genotypes of mannose-binding lectin to Mycoplasma pneumoniae infection in Chinese Merino sheep based on high-throughput sequencing technology

Mannose-binding lectin (MBL) glycoproteins in blood can selectively recognise lectins on the surface of bacteria, and play an important role in natural immunity. Micro RNAs (miRNAs) are key molecules that regulate gene expression at the post-transcriptional level in vivo, and their pathways are specific and effective. Previous studies indicate that small RNAs such as miRNAs perform regulatory roles in immunology. Herein, we investigated differential expression of miRNAs during MBL protein immunotherapy in sheep following treatment with different MBL genotypes (resistant and susceptible), and identified miRNAs linked to different target genes and pathways. RNA was extracted from liver tissue of resistant and susceptible sheep, miRNAs were identified by high-throughput sequencing, and differentially expressed miRNAs were analysed by SOAP to predict target genes and biological pathways. Results: Some miRNAs (oar-mir-143, oar-mir-10b, oar-mir-382, oar-mir-432 and oar-mir-379) were up-regulated, while others were down-regulated. GPATCH3 and DNAJC5 were predicted target genes of oar-mir-379, DMRT1 and GATA4 were linked to oar-mir-382, and oar-mir-432 was associated with STAT2, DMRT1 and ATG16L1. Identification of miRNAs differentially expressed in resistant and susceptible sheep may expand our understanding of miRNAs in immune regulation, and the role of MBL in innate immunity.

MiR-214-3p promotes proliferation and inhibits estradiol synthesis in porcine granulosa cells

Background

Granulosa cells (GCs) proliferation and estradiol synthesis significantly affect follicular development. The miR-214-3p expression in the ovarian tissues of high-yielding sows is higher than that in low-yielding sows, indicating that miR-214-3p may be involved in sow fertility. However, the functions and mechanisms of miR-214-3p on GCs are unclear. This study focuses on miR-214-3p in terms of the effects on GCs proliferation and estradiol synthesis.

Results

Our findings revealed that miR-214-3p promotes proliferation and inhibits estradiol synthesis in porcine GCs. MiR-214-3p can increase the percentage of S-phase cells, the number of EdU labeled positive cells, and cell viability. However, E₂ concentration was reduced after miR-214-3p agomir treatment. We also found that miR-214-3p up-regulates the expression of cell cycle genes including cell cycle protein B (Cyclin B), cell cycle protein D (Cyclin D), cell cycle protein E (Cyclin E), and cyclin-dependent kinase 4 (CDK4) at the transcription and translation levels, but down-regulates the mRNA and protein levels of cytochrome P450 family 11 subfamily A member 1 (CYP11A1), cytochrome P450 family 19 subfamily A member 1 (CYP19A1), and steroidogenic acute regulatory protein (StAR) (i.e., the key enzymes in estradiol synthesis). On-line prediction, bioinformatics analysis, a luciferase reporter assay, RT-qPCR, and Western blot results showed that the target genes of miR-214-3p in proliferation and estradiol synthesis are Mfn2 and NR5A1, respectively.

Conclusions

Our findings suggest that miR-214-3p plays an important role in the functional regulation of porcine GCs and therefore may be a target gene for regulating follicular development.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.

Identification of Ovarian Circular RNAs and Differential Expression Analysis between MeiShan and Large White Pigs

MeiShan and Large White pigs differ in their female fecundity. However, the mechanisms behind the gene expression and regulation that cause these differences remain unclear. In this study, we profiled circRNAs and identified 5,879 circRNAs from the ovaries of MeiShan and Large White pigs. Eighty-five circRNAs were differentially expressed between the two pig breeds. Of these, 37 were up-regulated and 48 were down-regulated in MeiShan pigs. Gene ontology enrichment analysis suggested that the differentially expressed circRNA were involved in the hormone-mediated signaling pathway. We verified that circSCIN and its parent gene, scinderin (SCIN), were differentially expressed by reverse transcription and quantitative PCR (RT-qPCR). Luciferase assays demonstrated that circSCIN can target and sponge miR-133 and miR-148a/b. The identification of differentially expressed circRNAs (DECs) and their regulatory functions increased our understanding of the differences in reproductive efficiency between MeiShan and Large White pigs.

Ovarian Circular RNAs Associated with High and Low Fertility in Large White Sows during the Follicular and Luteal Phases of the Estrous Cycle

In this study, the ovarian tissues of Large White pigs were mined for novel circular RNAs (circRNAs), following which, their molecular characteristics and potential mechanisms for fertility regulation were examined. RNA sequencing was used for transcriptome analysis of ovarian follicles and corpora lutea in Large White sows with high (H) and low (L) fertility during the follicular (F) and luteal (L) phases of the estrous cycle. In total, 21,386 circRNA derived from 4535 host genes were identified. Differentially expressed circRNAs were detected in the LH vs. LL (1079) and in the FH vs. FL (1077) comparisons, and their host genes were enriched in steroid biosynthesis and forkhead box O (FOXO), thyroid hormone, cell cycle, and tumor growth factor (TGF)-beta signaling pathways. Protein–protein interaction networks were constructed on the basis of the host genes that were significantly enriched in pathways related to reproductive processes, with AKT3 and PP2CB serving as the hub genes in the networks of the LH vs. LL and FH vs. FL comparisons, respectively. The microRNA (miRNA) binding sites of the differentially expressed circRNAs were predicted, and 128 (LH vs. LL) and 113 (FH vs. FL) circRNA–miRNA pairs were identified. Finally, circRNA–miRNA negative regulatory networks were established on the basis of the gene expression profiles and bioinformatic analyses. In the current study, differentially expressed circRNAs were observed in ovarian tissues between the H and L fertility groups in both F and L phases of the estrous cycle, which suggested roles in pig fertility regulation. These findings provide new clues for elucidating fertility differences in pigs.

Genome-Wide Differential Expression Profiling of Ovarian circRNAs Associated With Litter Size in Pigs

Circular RNAs (circRNAs) have been emerging as an important regulator in mammalian reproduction via acting as miRNA sponges. However, the circRNAs in porcine ovaries related with litter size remains largely unknown. In this study, porcine ovaries with smaller or larger litter size (LLS) were subjected to high-throughput RNA sequencing. In total, 38,722 circRNAs were identified, of which 1,291 circRNAs were commonly expressed in all samples. There were 56 circRNAs significantly down-regulated and 54 circRNAs up-regulated in LLS pig (|log2 (fold change) | > 1, FDR < 0.05). Bioinformatics predicted that most of circRNAs harbored miRNA binding sites, and the expression patterns of circRNAs and their putative binding miRNAs were validated by qPCR. Moreover, the expression of circ-TCP11/miR-183 was significantly reversely correlated and their direct interaction was confirmed by dual-luciferase assay. Our study indicates that circRNAs may play potential effects on modulating porcine litter size.

Transcriptome analysis of the effect of pyrroloquinoline quinone disodium (PQQ·Na2) on reproductive performance in sows during gestation and lactation

Background

Pyrroloquinoline quinone (PQQ), which is a water soluble, thermo-stable triglyceride-quinone, is widely distributed in nature and characterized as a mammalian vitamin-like redox cofactor. The objective of this study was to investigate the effects of pyrroloquinoline quinone disodium (PQQ·Na₂) on reproductive performance in sows.

Results

Dietary supplementation with PQQ·Na₂ significantly increased the total number of piglets born, the number of piglets born alive and the born alive litter weight. It also increased the antioxidant status in the placenta, plasma and milk. The concentration of NO was significantly increased in the plasma and placenta. RNA-seq analysis showed that 462 unigenes were differentially expressed between the control (Con) treatment and PQQ treatment groups. Among these unigenes, 199 were upregulated, while 263 unigenes were downregulated. The assigned functions of the unigenes covered a broad range of GO categories. Reproduction (27, 7.03%) and the reproduction process (27, 7.03%) were assigned to the biological process category. By matching DEGs to the KEGG database, we identified 29 pathways.

Conclusions

In conclusion, dietary supplementation with PQQ·Na₂ in gestating and lactating sows had positive effects on their reproductive performance.

Electronic supplementary material

The online version of this article (10.1186/s40104-019-0369-y) contains supplementary material, which is available to authorized users.

MiR-26a promotes apoptosis of porcine granulosa cells by targeting the 3β-hydroxysteroid-Δ24-reductase gene

Objective

Apoptosis of ovarian granulosa cells (GCs) affects mammalian follicular development and fecundity. This study aimed to explore the regulatory relationship between microRNA-26a (miR-26a) and the 3β-hydroxysteroid-Δ24-reductase gene (DHCR24) gene in porcine follicular granular cells (pGCs), and to provide empirical data for the development of methods to improve the reproductive capacity of pigs.

Methods

The pGCs were transfected with miR-26a mimic, miR-26a inhibitor and DHCR24-siRNA in vitro. The cell apoptosis rate of pGCs was detected by the flow cytometry. The secretion levels of estradiol (E2) and progesterone (P) in pGCs were detected by enzyme-linked immunosorbent assay. Double luciferase validation system was used to detect the binding sites between miR-26a and DHCR24 3′-UTR region. Qualitative real-time polymerase chain reaction and Western blotting were used to verify the DHCR24 mRNA and protein expression in pGCs, respectively, after transfecting with miR-26a mimic and miR-26a inhibitor.

Results

Results showed that enhancement of miR-26a promoted apoptosis, and inhibited E2 and P secretion in pGCs. Meanwhile, inhibition of DHCR24 also upregulated the Caspase-3 expression, reduced the BCL-2 expression, promoted pGCs apoptosis, and inhibited E2 and P secretion in pGCs. There were the binding sites of miR-26a located within DHCR24 3′-UTR. Up-regulation of miR-26a inhibited DHCR24 mRNA and protein expression in pGCs.

Conclusion

This study demonstrates that miR-26a can promote cell apoptosis and inhibit E2 and P secretion by inhibiting the expression of DHCR24 in pGCs.

MicroRNAs in porcine uterus and serum are affected by zearalenone and represent a new target for mycotoxin biomarker discovery

The mycotoxin zearalenone (ZEN) poses a risk to animal health because of its estrogenic effects. Diagnosis of ZEN-induced disorders remains challenging due to the lack of appropriate biomarkers. In this regard, circulating microRNAs (small non-coding RNAs) have remarkable potential, as they can serve as indicators for pathological processes in tissue. Thus, we combined untargeted and targeted transcriptomics approaches to investigate the effects of ZEN on the microRNA expression in porcine uterus, jejunum and serum, respectively. To this end, twenty-four piglets received uncontaminated feed (Control) or feed containing 0.17 mg/kg ZEN (ZEN low), 1.46 mg/kg ZEN (ZEN medium) and 4.58 mg/kg ZEN (ZEN high). After 28 days, the microRNA expression in the jejunum remained unaffected, while significant changes in the uterine microRNA profile were observed. Importantly, 14 microRNAs were commonly and dose-dependently affected in both the ZEN medium and ZEN high group, including microRNAs from the miR-503 cluster (i.e. ssc-miR-424-5p, ssc-miR-450a, ssc-miR-450b-5p, ssc-miR-450c-5p, ssc-miR-503 and ssc-miR-542-3p). Predicted target genes for those microRNAs are associated with regulation of gene expression and signal transduction (e.g. cell cycle). Although the effects in serum were less pronounced, receiver operating characteristic analysis revealed that several microRNA ratios were able to discriminate properly between non-exposed and ZEN-exposed pigs (e.g. ssc-miR-135a-5p/ssc-miR-432-5p, ssc-miR-542-3p/ssc-miR-493-3p). This work sheds new light on the molecular mechanisms of ZEN, and fosters biomarker discovery.

Combined microRNAome and transcriptome analysis of follicular phase and luteal phase in porcine ovaries

The aim of this study was to explore the expression difference of miRNAs and mRNAs between the follicular phase (FP) and luteal phase (LP) in porcine ovaries and provide a theoretical basis for the research on mammalian reproductive regulation. RNA-Seq and miRNA-Seq were used to identify differentially expressed genes (DEGs) and miRNAs (DEMs) between the FP and LP in ovaries of six sows (3-year-old Yorkshire pigs with similar weights and same parities). Bioinformatic analysis was used to screen potential genes and miRNAs related to porcine ovarian function. Real-time qualitative PCR was used to validate the sequencing results. RNA-Seq results showed that 3,078 genes were up-regulated, and 1,444 genes were down-regulated in the LP compared with the FP, and DEGs were significantly enriched in 242 Gene Ontology (GO) terms and 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. miRNA-Seq identified 112 DEMs, of which 25 were up-regulated and 87 were down-regulated in the LP compared with the FP. We obtained 186 intersection genes (IGs) between the 4,522 DEGs and 2,444 target genes predicted from the 112 DEMs. After constructing a miRNA-gene-pathway network, we identified key miRNAs and genes including miR-17-3p, miR-214, miR-221-5p, miR-125b, FGF1, YWHAG, YWHAZ, FDFT1 and DHCR24, which are enriched in Hippo and PI3K-Akt signalling pathways, and various metabolic pathways. These results indicate that these key genes and miRNAs may play important roles in the developmental transition from FP to LP in porcine ovaries and represent candidate targets for further study.

MiR-126-3p promotes the cell proliferation and inhibits the cell apoptosis by targeting TSC1 in the porcine granulosa cells

In mammalian ovaries, many studies demonstrated that the proliferation and apoptosis of granulosa cells are involved in folliculogenesis. Previous evidence suggests that miR-126-3p might get involved in the proliferation and apoptosis of granulosa cells, and tuberous sclerosis complex 1 (TSC1) gene was predicted as one target of miR-126-3p, and moreover, granulosa cell-specific TSC1 knockout stimulated folliculogenesis in mice. However, the molecular regulation of miR-126-3p on TSC1 and its effects on cell proliferation and apoptosis remain virtually unexplored in granulosa cells. Using porcine granulosa cells as a model, the luciferase report assay, mutation, deletion, Annexin-V/PI staining, and EdU assays were applied to investigate the molecular mechanism for miR-126-3p regulating the expression of TSC1 and their effects on the cell proliferation and apoptosis. We found that miR-126-3p showed a positive effect on cell proliferation and a negative effect on cell apoptosis in porcine granulosa cells, and knockdown of TSC1 significantly promoted cell proliferation and significantly inhibited cell apoptosis in porcine granulosa cells. Furthermore, miR-126-3p might target and repress the expressions of TSC1 at the post-transcriptional level, thereby promoting cell proliferation and inhibiting cell apoptosis of granulosa cells. These findings would provide of great insight in further exploring the molecular regulation of miR-126-3p and TSC1 on the functions of granulosa cells during the folliculogenesis in mammals.

Identification and Characterization of MicroRNAs in Skin of Chinese Giant Salamander (Andrias davidianus) by the Deep Sequencing Approach

MicroRNAs (miRNA) play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. However, the miRNAs from skin of Andrias davidianus have not been reported. In this study, a small-RNA cDNA library was constructed and sequenced from skin of A. davidianus. A total of 513 conserved miRNAs belonging to 174 families were identified. The remaining 108 miRNAs we identified were novel and likely to be skin tissue-specific but were expressed at low levels. The presence of randomly selected 15 miRNAs identified and their expression in eight different tissues from A. davidianus were validated by stem-loop qRT-PCR. For better understanding the functions of miRNAs, 129,791 predicated target genes were analyzed by GO and their pathways illustrated by KEGG pathway analyses. The results show that these identified miRNAs from A. davidianus skin are involved in a broad range of physiological functions including metabolism, growth, development, and immune responses. This study exhaustively identifies miRNAs and their target genes, which will ultimately pave the way for understanding their role in skin of A. davidianus and other amphibians. Further studies are necessary to better understand miRNA-mediated gene regulation.

MicroRNA-21 and microRNA-214 play important role in reproduction regulation during porcine estrous

Normal estrous cycle is crucial for porcine reproduction, and microRNA is closely related to regulation of estrous cycle in porcine ovaries. In this study, we found that the expression of miR-214 in porcine ovaries was higher than in many other tissues, and miR-21 expression in ovaries was significantly higher than in the uterus and pituitary. Meanwhile, miR-21 was upregulated and miR-214 was downregulated in the ovaries of high litter size (YH) pigs compared with low litter size (YL) pigs. Moreover, the lowest expression of miR-21 and miR-214 occurred on Days 14 and 7 of the estrous cycle and was expressed at greater levels in the granulosa cells of subordinate follicles than in dominant follicles on Day 3 of the estrous cycle. Bioinformatics analysis showed that miR-21 and miR-214 might target several genes that involved in the mTOR signaling, apoptosis, and steroid biosynthesis pathways, and they play important roles in maintaining the porcine estrous cycle. The qPCR and western blot analysis indicated that miR-214 inhibited the expression of SCARB1 gene in the transcriptional level, but not affected the SCARB1 gene's protein level. Our research findings indicated that miR-21 and miR-214 played important roles in reproduction regulation during porcine estrous.

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.

Pregnancy and spontaneous fetal loss: A pig perspective

Pigs have a unique, non-invasive epitheliochorial placenta where maternal and fetal layers lay in apposition. Indentation of fetal capillaries into the trophoblasts and maternal capillaries into the uterine epithelium reduce the distance between the fetal and maternal blood, ensuring nutrient transfer for proper conceptus development. Another unique feature of pig pregnancy is conceptus-mediated immune cell enrichment during the early stages of conceptus attachment (around gestation Day 15). This period coincides with the development of vasculature networks at the maternal-fetal interface, which is critical for successful conceptus growth. Specific chemokines, their receptors, and chemokine decoy receptor networks coordinate this immune cell enrichment and the positioning at the maternal-fetal interface. The recruited immune cells, in turn, adopt a specialized phenotype to support key processes of maternal-fetal adaptations, including tolerance to the semi-allogeneic fetus and supporting vascularization. Disturbance in coordinated cross talk between the conceptus and maternal endometrium is an important mechanism associated with spontaneous fetal loss. The exact mechanism of fetal loss is still not yet identified, although research in the last two decades point to various factors including genetics, nutrition, uterine capacity, placental efficiency, and imbalanced immune factors at the maternal-fetal interface. In this review, we summarize some of the recent advances in endometrial immune cell functions and their regulation. We also provide insights into endometrial/placental transcriptome, microRNA biology, and extravesicular transport across the maternal-fetal interface, as well as their potential implications in porcine pregnancy success or failure.

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.