Characterization of the Transcriptional Complexity of the Receptive and Pre-receptive Endometria of Dairy Goats

SPOCK: Master regulator of malignant tumors (Review)

SPARC/osteonectin, CWCV and Kazal‑like domain proteoglycan (SPOCK) is a family of highly conserved multidomain proteins. In total, three such family members, SPOCK1, SPOCK2 and SPOCK3, constitute the majority of extracellular matrix glycoproteins. The SPOCK gene family has been demonstrated to serve key roles in tumor regulation by affecting MMPs, which accelerates the progression of cancer epithelial‑mesenchymal transition. In addition, they can regulate the cell cycle via overexpression, inhibit tumor cell proliferation by inactivating PI3K/AKT signaling and have been associated with numerous microRNAs that influence the expression of downstream genes. Therefore, the SPOCK gene family are potential cancer‑regulating genes. The present review summarizes the molecular structure, tissue distribution and biological function of the SPOCK family of proteins, in addition to its association with cancer. Furthermore, the present review documents the progress made in investigations into the role of SPOCK, whilst also discussing prospects for the future of SPOCK‑targeted therapy, to provide novel ideas for clinical application and treatment.

Effects of conceptus proteins on endometrium and blood leukocytes of dairy cattle using transcriptome and meta-analysis

This study investigates the short and long-term effects of IFNT and PAG on the transcriptome of endometrium and blood leukocytes. Holstein heifers received intrauterine infusions of one of the following treatments: 20 mL of a 200 μg/mL bovine serum albumin solution (BSA; vehicle) from day 14 to 16 of the estrous cycle (BSA), vehicle + 10 μg/mL of IFNT from day 14 to 16 (IFNT3), vehicle + 10 μg/mL of IFNT from day 14 to 19 (IFNT6), and vehicle + 10 μg/mL of IFNT from day 14 to 16 followed by vehicle + 10 μg/mL of IFNT + 5 μg/mL of PAG from day 17 to 19 (IFNT+PAG). RNA-seq analysis was performed in endometrial biopsies and blood leukocytes collected after treatments. Acute IFNT signaling in the endometrium (IFNT3 vs BSA), induced differentially expressed genes (DEG) associated with interferon activation, immune response, inflammation, cell death, and inhibited vesicle transport and extracellular matrix remodeling. Prolonged IFNT signaling (IFNT6 vs IFNT3) altered gene expression related to cell invasion, retinoic acid signaling, and embryo implantation. In contrast, PAG induced numerous DEG in blood leukocytes but only 4 DEG in the endometrium. In blood leukocytes, PAG stimulated genes involved in development and TGFB signaling while inhibiting interferon signaling and cell migration. Overall, IFNT is a primary regulator of endometrial gene expression, while PAG predominantly affected the transcriptome of circulating immune cells during early pregnancy. Further research is essential to fully grasp the roles of identified DEG in both the endometrium and blood leukocytes.

CREB1 Is Involved in miR-134-5p-Mediated Endometrial Stromal Cell Proliferation, Apoptosis, and Autophagy

The successful establishment of endometrial receptivity is a key factor in ensuring the fertility of ewes and their economic benefits. Hu sheep have attracted attention due to their high fecundity and year-round estrus. In this study, we found that in the luteal phase, the uterine gland density, uterine coefficient, and number of uterine caruncles of high-fertility Hu sheep were higher than those of low-fertility Hu sheep. Thousands of differentially expressed genes were identified in the endometrium of Hu sheep with different fertility potential using RNA sequencing (RNA-Seq). Several genes involved in endometrial receptivity were screened using bioinformatics analysis. The qRT-PCR analysis further revealed the differential expression of cAMP reactive element binding protein-1 (CREB1) in the Hu sheep endometrium during the estrous cycle. Functionally, our results suggested that CREB1 significantly affected the expression level of endometrial receptivity marker genes, promoted cell proliferation by facilitating the transition from the G1 phase to the S phase, and inhibited cell apoptosis and autophagy. Moreover, we observed a negative linear correlation between miR-134-5p and CREB1 in the endometrium. In addition, CREB1 overexpression prevented the negative effect of miR-134-5p on endometrial stromal cell (ESC) growth. Taken together, these data indicated that CREB1 was regulated by miR-134-5p and may promote the establishment of uterine receptivity by regulating the function of ESCs. Moreover, this study provides new theoretical references for identifying candidate genes associated with fertility.

MSX1 Regulates Goat Endometrial Function by Altering the Plasma Membrane Transformation of Endometrial Epithelium Cells during Early Pregnancy

MSX1 is an important member of the muscle segment homeobox gene (Msh) family and acts as a transcription factor to regulate tissue plasticity, yet its role in goat endometrium remodeling remains elusive. In this study, an immunohistochemical analysis showed that MSX1 was mainly expressed in the luminal and glandular epithelium of goat uterus, and the MSX1 expression was upregulated in pregnancy at days 15 and 18 compared with pregnancy at day 5. In order to explore its function, goat endometrial epithelial cells (gEECs) were treated with 17 β-estrogen (E₂), progesterone (P₄), and/or interferon-tau (IFNτ), which were used to mimic the physiological environment of early pregnancy. The results showed that MSX1 was significantly upregulated with E₂- and P₄-alone treatment, or their combined treatment, and IFNτ further enhanced its expression. The spheroid attachment and PGE2/PGF2α ratio were downregulated by the suppression of MSX1. The combination of E₂, P₄, and IFNτ treatment induced the plasma membrane transformation (PMT) of gEECs, which mainly showed the upregulation of N-cadherin (CDH2) and concomitant downregulation of the polarity-related genes (ZO-1, α-PKC, Par3, Lgl2, and SCRIB). The knockdown of MSX1 partly hindered the PMT induced by E₂, P₄, and IFNτ treatment, while the upregulation of CDH2 and the downregulation of the partly polarity-related genes were significantly enhanced when MSX1 was overexpressed. Moreover, MSX1 regulated the CDH2 expression by activating the endoplasmic reticulum (ER) stress-mediated unfolded protein response (UPR) pathway. Collectively, these results suggest that MSX1 was involved in the PMT of the gEECs through the ER stress-mediated UPR pathway, which affects endometrial adhesion and secretion function.

Screening genes related to embryo implantation in Dazu black goats (Capra Hircus) by morphological and transcriptome analyses

Embryo implantation is a critical step in the establishment of pregnancy. However, the mechanisms of embryo implantation during early pregnancy in goats remain unclear due to the lack of published studies examining the genes involved in embryo implantation. As a popular goat breed in southwest China, Dazu black goats (DBGs) are highly adaptable and exhibit high fertility, making this breed a good model in which to study reproductive performance of goats. Here, morphological analysis showed that compared with the non-pregnant (NP) groups, the endometrial thickness of the goats in the P15 and P19 groups (15 and 19-day pregnant groups, respectively) were increased (P < 0.01). Proliferating Cell Nuclear Antigen (PCNA) staining showed that PCNA was expressed in the NP, P15, and P19 groups. Transcriptome analysis was then conducted to identify gene expression patterns in uterine tissue during DBG embryo implantation. By comparing uterine tissue at different stages of embryonic implantation, 48 in NP_vs._P15, 318 in NP_vs._P19, and 1439 in P15_vs._P19, differentially expressed mRNAs were identified. Gene Ontology (GO) enrichments of the differentially expressed genes were enriched in the extracellular region, extracellular space, transporter activity, extracellular region, immune system process, immune response, and defense response etc. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the biological metabolic pathways with which the differentially expressed genes are associated were explored. Through KEGG analysis, the DBGs were associated with oxidative phosphorylation, complement and coagulation cascades, arginine and proline metabolism, metabolic pathways, arachidonic acid metabolism, and ECM-receptor interaction. These candidate genes (CSF1, C1S, CST6, SLC24A4, HOXA10, HOXA11, MMP9, and ITGA11) and enriched signaling pathways could be valuable references for exploring the molecular mechanisms underlying goat embryo implantation.

RNA sequencing reveals dynamic expression of lncRNAs and mRNAs in caprine endometrial epithelial cells induced by Neospora caninum infection

Background

The effective transmission mode of Neospora caninum, with infection leading to reproductive failure in ruminants, is vertical transmission. The uterus is an important reproductive organ that forms the maternal–fetal interface. Neospora caninum can successfully invade and proliferate in the uterus, but the molecular mechanisms underlying epithelial-pathogen interactions remain unclear. Accumulating evidence suggests that host long noncoding RNAs (lncRNAs) play important roles in cellular molecular regulatory networks, with reports that these RNA molecules are closely related to the pathogenesis of apicomplexan parasites. However, the expression profiles of host lncRNAs during N. caninum infection has not been reported.

Methods

RNA sequencing (RNA-seq) analysis was used to investigate the expression profiles of messenger RNAs (mRNAs) and lncRNAs in caprine endometrial epithelial cells (EECs) infected with N. caninum for 24 h (TZ_24h) and 48 h (TZ_48 h), and the potential functions of differentially expressed (DE) lncRNAs were predicted by using Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of their mRNA targets.

Results

RNA-seq analysis identified 1280.15 M clean reads in 12 RNA samples, including six samples infected with N. caninum for 24 h (TZ1_24h-TZ3_24h) and 48 h (TZ1_48h-TZ3_48h), and six corresponding control samples (C1_24h-C3_24h and C1_48h-C3_48h). Within the categories TZ_24h-vs-C_24h, TZ_48h-vs-C_48h and TZ_48h-vs-TZ_24h, there were 934 (665 upregulated and 269 downregulated), 1238 (785 upregulated and 453 downregulated) and 489 (252 upregulated and 237 downregulated) DEmRNAs, respectively. GO enrichment and KEGG analysis revealed that these DEmRNAs were mainly involved in the regulation of host immune response (e.g. TNF signaling pathway, MAPK signaling pathway, transforming growth factor beta signaling pathway, AMPK signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway), signaling molecules and interaction (e.g. cytokine-cytokine receptor interaction, cell adhesion molecules and ECM-receptor interaction). A total of 88 (59 upregulated and 29 downregulated), 129 (80 upregulated and 49 downregulated) and 32 (20 upregulated and 12 downregulated) DElncRNAs were found within the categories TZ_24h-vs-C_24h, TZ_48h-vs-C_48h and TZ_48h-vs-TZ_24h, respectively. Functional prediction indicated that these DElncRNAs would be involved in signal transduction (e.g. MAPK signaling pathway, PPAR signaling pathway, ErbB signaling pathway, calcium signaling pathway), neural transmission (e.g. GABAergic synapse, serotonergic synapse, cholinergic synapse), metabolism processes (e.g. glycosphingolipid biosynthesis-lacto and neolacto series, glycosaminoglycan biosynthesis-heparan sulfate/heparin) and signaling molecules and interaction (e.g. cytokine-cytokine receptor interaction, cell adhesion molecules and ECM-receptor interaction).

Conclusions

This is the first investigation of global gene expression profiles of lncRNAs during N. caninum infection. The results provide valuable information for further studies of the roles of lncRNAs during N. caninum infection.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05405-5.

Effect of MiR-100-5p on proliferation and apoptosis of goat endometrial stromal cell in vitro and embryo implantation in vivo

The growth of endometrial stromal cells (ESCs) at implantation sites may be a potential factor affecting the success rate of embryo implantation. Incremental proofs demonstrated that ncRNAs (e.g. miRNAs, lncRNAs and circRNAs) were involved in various biological procedures, including proliferation and apoptosis. In this study, the role of miR‐100‐5p on proliferation and apoptosis of goat ESCs in vitro and embryo implantation in vivo was determined. The mRNA expression of miR‐100‐5p was significantly inhibited in the receptive phase (RE) rather than in the pre‐receptive phase (PE). Overexpression of miR‐100‐5p suppressed ESCs proliferation and induced apoptosis. The molecular target of MiR‐100‐5p, HOXA1, was confirmed by 3′‐UTR assays. Meanwhile, the product of HOXA1 mRNA RT‐PCR increased in the RE more than that in the PE. The HOXA1‐siRNA exerted significant negative effects on growth arrest. Instead, incubation of ESCs with miR‐100‐5p inhibitor or overexpressed HOXA1 promoted the cell proliferation. In addition, Circ‐9110 which acted as a sponge for miR‐100‐5p reversed the relevant biological effects of miR‐100‐5p. The intrinsic apoptosis pathway was suppressed in ESCs, revealing a crosstalk between Circ‐9110/miR‐100‐5p/HOXA1 axis, PI3K/AKT/mTOR, and ERK1/2 pathways. To further evaluate the progress in study on embryo implantation regulating mechanism of miR‐100‐5p in vivo, the pinopodes of two phases were observed and analysed, suggesting that, as similar as in situ, miR‐100‐5p was involved in significantly regulating embryo implantation in vivo. Mechanistically, miR‐100‐5p performed its embryo implantation function through regulation of PI3K/AKT/mTOR and ERK1/2 pathways by targeting Circ‐9110/miR‐100‐5p/HOXA1 axis in vivo.

Essential Role of CRIM1 on Endometrial Receptivity in Goat

In domestic ruminants, endometrial receptivity is related to successful pregnancy and economic efficiency. Despite several molecules having been reported in the past regarding endometrial receptivity regulation, much regarding the mechanism of endometrial receptivity regulation remains unknown due to the complex nature of the trait. In this work, we demonstrated that the cysteine-rich transmembrane bone morphogenetic protein (BMP) regulator 1 (CRIM1) served as a novel regulator in the regulation of goat endometrial receptivity in vitro. Our results showed that hormones and IFN-τ increased the expression of CRIM1 in goat endometrial epithelial cells (EECs). Knockdown of CRIM1 via specific shRNA hindered cell proliferation, cell adhesion and prostaglandins (PGs) secretion and thus derailed normal endometrial receptivity. We further confirmed that receptivity defect phenotypes due to CRIM1 interference were restored by ATG7 overexpression in EECs while a loss of ATG7 further impaired receptivity phenotypes. Moreover, our results showed that changing the expression of ATG7 affected the reactive oxygen species (ROS) production. Moreover, mR-143-5p was shown to be a potential upstream factor of CRIM1-regulated endometrial receptivity in EECs. Overall, these results suggest that CRIM1, as the downstream target of miR-143-5p, has effects on ATG7-dependent autophagy, regulating cell proliferation, cell adhesion and PG secretion, and provides a new target for the diagnosis and treatment of early pregnancy failure and for improving the success rates of artificial reproduction.

miR-34a/c induce caprine endometrial epithelial cell apoptosis by regulating circ-8073/CEP55 via the RAS/RAF/MEK/ERK and PI3K/AKT/mTOR pathways

microRNAs (miRNAs) and circular RNAs (circRNAs) are important for endometrial receptivity establishment and embryo implantation in mammals. miR-34a and miR-34c are highly expressed in caprine receptive endometrium (RE). Herein, the functions and mechanisms of miR-34a/c in caprine endometrial epithelial cell (CEEC) apoptosis and RE establishment were investigated. miR-34a/c downregulated the expression level of centrosomal protein 55 (CEP55) and were sponged by circRNA8073 (circ-8073), thereby exhibiting a negative interaction in CEEC. miR-34a/c induced CEEC apoptosis by targeting circ-8073/CEP55 through the regulation of the RAS/RAF/MEK/ERK and phosphoitide 3-kinase (PI3K)/protein kinase B (AKT)/mammalian target of rapamycin (mTOR) pathways. Positive and negative feedback loops and cross-talk were documented between the RAS/RAF/MEK/ERK and PI3K/AKT/mTOR pathways. miR-34a/c regulated the levels of RE marker genes, including forkhead box M1, vascular endothelial growth factor, and osteopontin (OPN). These results suggest that miR-34a/c not only induce CEEC apoptosis by binding to circ-8073 and CEP55 via the RAS/RAF/MEK/ERK and PI3K/AKT/mTOR pathways, but may also regulate RE establishment in dairy goats.

SPOCK1 Involvement in Epithelial-to-Mesenchymal Transition: A New Target in Cancer Therapy?

Background

Cancer metastasis is the main obstacle to increasing the lifespan of cancer patients. Epithelial-to-mesenchymal transition (EMT) plays a significant role in oncogenic processes, including tumor invasion, intravasation, and micrometastasis formation, and is especially critical for cancer invasion and metastasis. The extracellular matrix (ECM) plays a crucial role in the occurrence of EMT corresponding to the change in adhesion between cells and matrices.

Conclusion

SPOCK1 is a critical regulator of the ECM and mediates EMT in cancer cells. This suggests an important role for SPOCK1 in tumorigenesis, migration and invasion. SPOCK1 is a critical regulator of some processes involved in cancer progression, including cancer cell proliferation, apoptosis and migration. Herein, the functions of SPOCK1 in cancer progression are expounded, revealing the association between SPOCK1 and EMT in cancer metastasis. SPOCK1 is a positive downstream regulator of transforming growth factor-β, and SPOCK1-mediated EMT regulates invasion and metastasis through the Wnt/β-catenin pathway and PI3K/Akt signaling pathway. It is of significance that SPOCK1 may be an attractive prognostic biomarker and therapeutic target in cancer treatment.

A comparative analysis of gene expression induced by the embryo in the caprine endometrium

Transcriptomics is an established powerful tool to identify potential mRNAs and ncRNAs (non‐coding RNAs) for endometrial receptivity. In this study, the goat endometrium at estrus day 5 (ED5) and estrus day 15 (ED15) were selected to systematically analyse the differential expressed genes (DEGs) what were induced by the embryo. There were 1,847 genes which were significantly differential expressed in endometrium induced by the embryo at ED5, and 1,346 at ED15 (p‐value < .05). Secreted phosphoprotein 1 (SPP) was the responsive genes for embryo in the goat endometrium during estrus cycle, neurotensis (NTS) and pleiotrophin (PTN) were the responsive genes for embryo in the goat endometrium at ED5, Testin (TES) and Phosphate and Tension Homology Deleted on Chromsome ten (PTEN) at ED15. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) analysis revealed cytoplasm and Endocytosis were indispensable for the endometrium development in dairy goat. In a word, this resulting view of the transcriptome greatly uncovered the global trends in mRNAs expression induced by the embryo in the endometrium of dairy goats.

Transcriptome Changes in the Mink Uterus during Blastocyst Dormancy and Reactivation

Embryo implantation in the mink follows the pattern of many carnivores, in that preimplantation embryo diapause occurs in every gestation. Details of the gene expression and regulatory networks that terminate embryo diapause remain poorly understood. Illumina RNA-Seq was used to analyze global gene expression changes in the mink uterus during embryo diapause and activation leading to implantation. More than 50 million high quality reads were generated, and assembled into 170,984 unigenes. A total of 1684 differential expressed genes (DEGs) in uteri with blastocysts in diapause were compared to the activated embryo group (p < 0.05). Among these transcripts, 1527 were annotated as known genes, including 963 up-regulated and 564 down-regulated genes. The gene ontology terms for the observed DEGs, included cellular communication, phosphatase activity, extracellular matrix and G-protein couple receptor activity. The KEGG pathways, including PI3K-Akt signaling pathway, focal adhesion and extracellular matrix (ECM)-receptor interactions were the most enriched. A protein-protein interaction (PPI) network was constructed, and hub nodes such as VEGFA, EGF, AKT, IGF1, PIK3C and CCND1 with high degrees of connectivity represent gene clusters expected to play an important role in embryo activation. These results provide novel information for understanding the molecular mechanisms of maternal regulation of embryo activation in mink.

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.

The mRNA and lncRNA landscape of the non-pregnant endometrium during the oestrus cycle in dairy goat

Cyclic changes in the endometrium are essential for embryo implantation in mammals; many studies report that such changes constitute a complex process involving numerous molecular mediators. In the present study, goat endometria at oestrus Day 5 and oestrus Day 15 were selected to systematically analyse the transcriptome using strand-specific Ribo-Zero RNA sequencing. Over 120 million high-quality paired-end reads were generated and 440400 transcripts were identified in the endometrial tissue of dairy goats. In total, 489 differentially expressed mRNAs and 854 differentially expressed long non-coding RNAs were identified when comparing the endometrium at goat endometria at oestrus Day 5 and oestrus Day 15. Neurotensin was found to play a potentially important role in the non-pregnant goat endometrium during the oestrus cycle. Furthermore, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses of the cis-target genes of the differentially expressed long non-coding RNAs showed that GO:0005198 (structural molecule activity) and ko04510 (focal adhesion) might be involved in cyclic endometrial changes. Taken together, the resulting transcriptomic profiles elucidate global trends in mRNA and lncRNA expression in non-pregnant endometria during the oestrus cycle in dairy goats.

Endometrial Epithelial Cell Apoptosis Is Inhibited by a ciR8073-miR181a-Neurotensis Pathway during Embryo Implantation

Development of the receptive endometrium (RE) from the pre-receptive endometrium (PE) is essential for embryo implantation, but its molecular mechanisms have not been fully understood. In this study, lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks were constructed to explore the functions of potential competing endogenous RNAs (ceRNA) during the development of RE in dairy goats. We observed that circRNA8073 (ciR8073) decreased the levels of miR-181a by acting as a miRNA sponge. This effect indirectly increased the expression of neurotensin in endometrial epithelial cells (EECs). Neurotensin then inhibited EEC apoptosis by increasing the expression of BCL-2/BAX in favor of BCL-2 via the MAPK pathway and also induced increased expression of leukemia-inhibitory factor, cyclo-oxygenase 2, vascular endothelial growth factor A, and homeobox A10. We have thus identified a ciR8073-miR181a-neurotensin pathway in the endometrium of dairy goats. Through this pathway, ciR8073 functions as a ceRNA that sequesters miR-181a, thereby protecting neurotensin transcripts from miR-181a-mediated suppression in EECs.

Construction of implantation failure related lncRNA-mRNA network and identification of lncRNA biomarkers for predicting endometrial receptivity

Insufficient endometrial receptivity is a major factor leading to implantation failure (IF), and the traditional way of morphological observation of endometrium cannot determine the condition of receptivity sufficiently. Considering that long-noncoding RNAs (lncRNAs) regulate endometrial receptivity and competing endogenous RNA (ceRNA) mechanism works in plenty of biological processes, ceRNA is likely to function in the pathology of IF. In the present study, we aim to construct an implantation failure related lncRNA-mRNA network (IFLMN), and to identify the key lncRNAs as the candidates for predicting endometrial receptivity. The global background network was constructed based on the presumed lncRNA-miRNA and miRNA-mRNA pairs obtained from lncRNASNP and miRTarBase. Differentially expressed genes (DEGs) of IF were calculated using the data of GSE26787, and then re-annotated as differentially expressed mRNAs (DEMs) and lncRNAs (DELs). IFLMN was constructed by hypergeometric test, including 255 lncRNA-mRNA pairs, 10 lncRNAs, and 212 mRNAs. Topological analysis determined the key lncRNAs with the highest centroid. Functional enrichment analyses were performed by unsupervised clustering, GO classification, KEGG pathway, and co-expression module analyses, achieving six key lncRNAs and their ceRNA sub-networks, which were involved in immunological activity, growth factor binding, vascular proliferation, apoptosis, and steroid biosynthesis in uterus and prepared endometrium for embryo implantation. Sixteen endometrial samples were collected during mid-luteal phase, including 8 recurrent implantation failure (RIF) or recurrent miscarriage (RM) women and 8 controls who conceived successfully. Quantitative real-time PCR was performed to compare the expression of the above six lncRNAs, which validated that the expression of all these lncRNAs was significantly elevated in endometrium of RIF/RM patients. Further studies are needed to investigate the underlying mechanism, and the lncRNAs may be developed into predictive biomarkers for endometrial receptivity.

Identification and characterization of long noncoding RNAs and mRNAs expression profiles related to postnatal liver maturation of breeder roosters using Ribo-zero RNA sequencing

BACKGROUND

The liver is mainly hematopoietic in the embryo, and converts into a major metabolic organ in the adult. Therefore, it is intensively remodeled after birth to adapt and perform adult functions. Long non-coding RNAs (lncRNAs) are involved in organ development and cell differentiation, likely they have potential roles in regulating postnatal liver development. Herein, in order to understand the roles of lncRNAs in postnatal liver maturation, we analyzed the lncRNAs and mRNAs expression profiles in immature and mature livers from one-day-old and adult (40 weeks of age) breeder roosters by Ribo-Zero RNA-Sequencing.

RESULTS

Around 21,939 protein-coding genes and 2220 predicted lncRNAs were expressed in livers of breeder roosters. Compared to protein-coding genes, the identified chicken lncRNAs shared fewer exons, shorter transcript length, and significantly lower expression levels. Notably, in comparison between the livers of newborn and adult breeder roosters, a total of 1570 mRNAs and 214 lncRNAs were differentially expressed with the criteria of log2fold change > 1 or < - 1 and P values < 0.05, which were validated by qPCR using randomly selected five mRNAs and five lncRNAs. Further GO and KEGG analyses have revealed that the differentially expressed mRNAs were involved in the hepatic metabolic and immune functional changes, as well as some biological processes and pathways including cell proliferation, apoptotic and cell cycle that are implicated in the development of liver. We also investigated the cis- and trans- regulatory effects of differentially expressed lncRNAs on its target genes. GO and KEGG analyses indicated that these lncRNAs had their neighbor protein coding genes and trans-regulated genes associated with adapting of adult hepatic functions, as well as some pathways involved in liver development, such as cell cycle pathway, Notch signaling pathway, Hedgehog signaling pathway, and Wnt signaling pathway.

CONCLUSIONS

This study provides a catalog of mRNAs and lncRNAs related to postnatal liver maturation of chicken, and will contribute to a fuller understanding of biological processes or signaling pathways involved in significant functional transition during postnatal liver development that differentially expressed genes and lncRNAs could take part in.

Circ-8073 regulates CEP55 by sponging miR-449a to promote caprine endometrial epithelial cells proliferation via the PI3K/AKT/mTOR pathway

Circular RNAs (circRNAs) are a large class of endogenous non-coding RNAs that function as regulators in various cells and tissues. Here, the function and mechanism of circRNA8073 (Circ-8073) on endometrial epithelial cells (EECs) and the development of endometrial receptivity were investigated in dairy goats. Circ-8073 could bind to and inhibit miR-449a activity. Circ-8073 binding to the target site of miR-449a had a negative feedback relationship. Centrosomal protein55 (CEP55) was a direct target gene of miR-449a, and Circ-8073 could increase the expression levels of CEP55 by sponging miR-449a in EECs in vitro. Circ-8073/miR-449a/CEP55 could promote EECs proliferation via the PI3K/AKT/mTOR pathway. In addition, CEP55 could regulate the expression levels of vascular endothelial growth factor (VEGF) and forkhead box M1 (FOXM1) in EECs, which contributed to the development of endometrial receptivity. These findings showed that Circ-8073 regulated CEP55 by sponging miR-449a to promote EEC proliferation via the PI3K/AKT/mTOR pathway, suggesting that it could function as a regulator in the development of endometrial receptivity in dairy goats.

Testin was regulated by circRNA3175-miR182 and inhibited endometrial epithelial cell apoptosis in pre-receptive endometrium of dairy goats

Circular RNAs (circRNAs) in various tissues and cell types from mammalian sources have been studied. However, present knowledge on circRNAs in the development of pre-receptive endometrium (PE) in dairy goats is limited. In the pre-receptive endometrium of dairy goats, higher circRNA3175 (ciR3175) levels, lower miR-182 levels and higher Testin (TES) levels were detected. And ciR3175 could decreased the miR-182 levels by acting as a miRNA sponge, and miR-182 could down-regulated the expression level of TES via the predicted target site in endometrial epithelial cells (EECs) in vitro. Via this way, ciR3175 functioned as a competing endogenous RNAs (ceRNA) that sequestered miR-182, thereby protecting TES transcripts from miR-182-mediated suppression in EECs in vitro. Further, TES inhibited EECs apoptosis by decreasing the expression level of BCL-2/BAX via the MAPK pathway. Thus, a ciR3175-miR182-TES pathway in the endometrium was identified in EECs, and the modulation of which could emerge as a potential target in regulating the pre-receptive endometrium development in dairy goats.

Bovine Milk Comparative Proteome Analysis from Early, Mid, and Late Lactation in the Cattle Breed, Malnad Gidda (Bos indicus)

Bovine milk is important for both veterinary medicine and human nutrition. Understanding the bovine milk proteome at different stages of lactation has therefore broad significance for integrative biology and clinical medicine as well. Indeed, different lactation stages have marked influence on the milk yield, milk constituents, and nourishment of the neonates. We performed a comparative proteome analysis of the bovine milk obtained at different stages of lactation from the Indian indigenous cattle Malnad Gidda (Bos indicus), a widely available breed. The milk differential proteome during the lactation stages in B. indicus has not been investigated to date. Using high-resolution mass spectrometry-based quantitative proteomics of the bovine whey proteins at early, mid, and late lactation stages, we identified a total of 564 proteins, out of which 403 proteins were found to be differentially abundant at different lactation stages. As is expected of any body fluid proteome, 51% of the proteins identified in the milk were found to have signal peptides. Gene ontology analyses were carried out to categorize proteins altered across different lactation stages based on biological process and molecular function, which enabled us to correlate their significance in each lactation stage. We also investigated the potential pathways enriched in different lactation stages using bioinformatics pathway analysis tools. To the best of our knowledge, this study represents the first and largest inventory of milk proteins identified to date for an Indian cattle breed. We believe that the current study broadly informs both veterinary omics research and the emerging field of nutriproteomics during lactation stages.

miR-26a promoted endometrial epithelium cells (EECs) proliferation and induced stromal cells (ESCs) apoptosis via the PTEN-PI3K/AKT pathway in dairy goats

Changes in endometrial cell morphology and function are absolutely necessary for successful embryo implantation. In this study, miR-26a was widely expressed in dairy goats, and was found to be regulated by β-estradiol (E2) and progesterone (P4) in endometrial epithelium cells (EECs) as well as stromal cells (ESCs). Furthermore, miR-26a played a role in the regulation of cells proliferation and apoptosis by directly regulating PTEN and indirectly regulating the PI3K/AKT pathway in EECs but not in ESCs of dairy goats in vitro. In addition, miR-26a regulated the expression of osteopontin (OPN), vascular endothelial growth factor (VEGF), Cyclooxygenase-2 (COX-2), and prolactin (PRL) in endometrial cells. Therefore, we could get a conclusion that miR-26a had very complex and diverse functions in the endometrial cells during the development of endometrial receptivity in dairy goats. This study provided an efficient platform for studying the regulatory effect of miR-26a on endometrial cells during the development of endometrial receptivity in dairy goats.

CircRNA-9119 regulates the expression of prostaglandin-endoperoxide synthase 2 (PTGS2) by sponging miR-26a in the endometrial epithelial cells of dairy goat

Circular RNAs (circRNAs) have been found to play important functional roles in epigenetic regulation under certain physiological and pathological conditions. However, knowledge of circRNAs during the development of receptive endometrium (RE) from pre-RE is limited. In the RE of dairy goats, higher circRNA-9119 levels, with lower miR-26a and higher prostaglandin-endoperoxide synthase 2 (PTGS2) levels, were detected. Further study showed that circRNA-9119 decreased levels of miR-26a by acting as a microRNA sponge, and that miR-26a downregulated the expression of PTGS2 via the predicted target site in endometrial epithelial cells (EECs) of dairy goats in vitro. In this way, circRNA-9119 functioned as a competing endogenous RNAs (ceRNA) that sequestered miR-26a, thereby protecting PTGS2 transcripts from miR-26a-mediated suppression in dairy goat EECs in vitro. Furthermore, PTGS2 participated in the regulation of some protein markers for endometrial receptivity in dairy goat EECs in vitro. Thus, a circRNA-9119–miR-26a–PTGS2 pathway in the endometrium was identified, and modulation of circRNA-9119–miR-26a–PTGS2 expression in EECs may emerge as a potential target to regulate the development of RE.

The muscle development transcriptome landscape of ovariectomized goat

In practical production, almost all rams and about 50% of ewes are used to fatten. Researchers have proved that ewe ovariectomy could improve the productivity significantly, but the specific molecular mechanism is still unknown. In this study, five independent cDNA libraries (three and two from ovariectomized and normal ewe longissimus dorsi samples, respectively) were constructed to thoroughly explore the global transcriptome, further to reveal how the ovariectomized ewes influence muscle development by Illumina2000 sequencing technology. As a result, 205 358 transcripts and 118 264 unigenes were generated. 15 490 simple sequence repeats (SSRs) were revealed and divided into six types, and the short repeat sequence SSR (monomers, dimers, trimers) was the domain type. Single nucleotide polymorphism analysis found that the number of transition was greater than the number of transversion among the five libraries. Furthermore, 1612 differently expressed genes (DEGs) (Log2fold_change > 1 and p < 0.05) were revealed between ovariectomized and normal ewe groups, in which 903 genes were expressed commonly in the two groups, and 288 and 421 genes were uniquely expressed in normal and ovariectomized ewe groups, respectively. Gene Ontology (GO) analysis categorized all unigenes into 555 GO terms and 56 DEGs were significantly categorized into 43 GO terms (p < 0.05). KEGG enrichment analysis annotated 12 976 genes (containing 137 DEGs) to 86 pathways, among them 24 and 11 DEGs involved in development and reproduction associated pathways, respectively. To validate the reliability of the RNA-seq analysis, 22 candidate DEGs were randomly selected to perform quantitative real-time polymerase chain reaction. The result showed that 9 and 1 genes were significantly and approximately significantly expressed in control and treatment group, respectively, and the results of RNA-seq are believable in this study. Overall, these results were helpful for elucidating the molecular mechanism of muscle development of ovariectomized animals and the application of female ovariectomy in fattening.

MiR-449a regulates caprine endometrial stromal cell apoptosis and endometrial receptivity

In this study, an RT-qPCR analysis showed that the expression levels of miR-449a in the pre-receptive endometrium were lower compared to the receptive endometrium, which is consistent with previous sequencing data (previous investigations). To detect the role of miR-449a in endometrial receptivity, we transfected caprine endometrial stromal cells (ESCs) cultured in vitro with miR-449a mimics. The results revealed that miR-449a decreased the mRNA and protein levels of LGR4 by targeting its 3′-untranslated region. The miR-449a mimics significantly reduced the G₁ cell population from 52.56% (mimic NC) to 42.19% with a concordant increase in the G₂ and S cell populations from 47.44% (mimic NC) to 57.81%, suggesting that miR-449a caused ESC cell cycle arrest. In addition, the number of apoptotic cells was significantly higher in ESCs transfected with miR-449a mimics (P < 0.05) than in ESCs transfected with mimic NC. In vivo, rich pinopodes were observed on the endometrial surface in the miR-449a agomir group compared with the miR-449a antagomir group. The results of hematoxylin-eosin staining showed that endometrial thickness was significantly increased in the miR-449a agomir group compared with the miR-449a antagomir group. These results suggest that miR-449a could enhance endometrial receptivity.

Identification of differentially expressed genes through RNA sequencing in goats (Capra hircus) at different postnatal stages

Intramuscular fat (IMF) content and fatty acid composition of longissimus dorsi muscle (LM) change with growth, which partially determines the flavor and nutritional value of goat (Capra hircus) meat. However, unlike cattle, little information is available on the transcriptome-wide changes during different postnatal stages in small ruminants, especially goats. In this study, the sequencing reads of goat LM tissues collected from kid, youth, and adult period were mapped to the goat genome. Results showed that out of total 24 689 Unigenes, 20 435 Unigenes were annotated. Based on expected number of fragments per kilobase of transcript sequence per million base pairs sequenced (FPKM), 111 annotated differentially expressed genes (DEGs) were identified among different postnatal stages, which were subsequently assigned to 16 possible expression patterns by series-cluster analysis. Functional classification by Gene Ontology (GO) analysis was used for selecting the genes showing highest expression related to lipid metabolism. Finally, we identified the node genes for lipid metabolism regulation using co-expression analysis. In conclusion, these data may uncover candidate genes having functional roles in regulation of goat muscle development and lipid metabolism during the various growth stages in goats.

miR-182 aids in receptive endometrium development in dairy goats by down-regulating PTN expression

Increasing evidence has shown that miRNAs play important roles in endometrium development during the menstrual cycle in humans and many other animals. Our previous data indicated that miR-182 levels increase 15.55-fold and pleiotrophin (PTN) levels decrease 20.97-fold in the receptive endometrium (RE, D15) compared with the pre-receptive endometrium (PE, D5) in dairy goats. The present study shows that miR-182 is widely expressed in different tissues of dairy goats and that its expression levels are regulated by E2 and P4 in endometrial epithelium cells (EECs). We confirmed that PTN is a target of miR-182 and that miR-182 regulates the protein levels of AKT, Bcl-2, FAS, MAPK, Caspase-3 and SP1 in EECs. Furthermore, miR-182 up-regulates or maintains the expression levels of osteopontin (OPN), cyclooxygenase-2 (COX-2) and prolactin receptor (PRLR) in EECs, suggesting that miR-182 is an important regulatory factor in the construction of endometrial receptivity in dairy goats. In conclusion, miR-182 participates in the development of endometrial receptivity by down-regulating PTN and affecting the expression of select apoptosis-related genes and increasing or maintaining the expression levels of OPN, COX-2 and PRLR in the EECs of dairy goats.