Integrated Analysis of miRNA-mRNA Network Reveals Different Regulatory Patterns in the Endometrium of Meishan and Duroc Sows during Mid-Late Gestation

Sex and fetal genome influence gene expression in pig endometrium at the end of gestation

BACKGROUND

A fine balance of feto-maternal resource allocation is required to support pregnancy, which depends on interactions between maternal and fetal genetic potential, maternal nutrition and environment, endometrial and placental functions. In particular, some imprinted genes have a role in regulating maternal-fetal nutrient exchange, but few have been documented in the endometrium. The aim of this study is to describe the expression of 42 genes, with parental expression, in the endometrium comparing two extreme breeds: Large White (LW); Meishan (MS) with contrasting neonatal mortality and maturity at two days of gestation (D90-D110). We investigated their potential contribution to fetal maturation exploring genes-fetal phenotypes relationships. Last, we hypothesized that the fetal genome and sex influence their endometrial expression. For this purpose, pure and reciprocally crossbred fetuses were produced using LW and MS breeds. Thus, in the same uterus, endometrial samples were associated with its purebred or crossbred fetuses.

RESULTS

Among the 22 differentially expressed genes (DEGs), 14 DEGs were differentially regulated between the two days of gestation. More gestational changes were described in LW (11 DEGs) than in MS (2 DEGs). Nine DEGs were differentially regulated between the two extreme breeds, highlighting differences in the regulation of endometrial angiogenesis, nutrient transport and energy metabolism. We identified DEGs that showed high correlations with indicators of fetal maturation, such as ponderal index at D90 and fetal blood fructose level and placental weight at D110. We pointed out for the first time the influence of fetal sex and genome on endometrial expression at D90, highlighting AMPD3, CITED1 and H19 genes. We demonstrated that fetal sex affects the expression of five imprinted genes in LW endometrium. Fetal genome influenced the expression of four genes in LW endometrium but not in MS endometrium. Interestingly, both fetal sex and fetal genome interact to influence endometrial gene expression.

CONCLUSIONS

These data provide evidence for some sexual dimorphism in the pregnant endometrium and for the contribution of the fetal genome to feto-maternal interactions at the end of gestation. They suggest that the paternal genome may contribute significantly to piglet survival, especially in crossbreeding production systems.

In-depth proteome characterization of endometrium and extraembryonic membranes during implantation in pig

BACKGROUND

Proteome characterization of the porcine endometrium and extraembryonic membranes is important to understand mother-embryo cross-communication. In this study, the proteome of the endometrium and chorioallantoic membrane was characterized in pregnant sows (PS) during early gestation (d 18 and 24 of gestation) and in the endometrium of non-pregnant sows (NPS) during the same days using LC-MS/MS analysis. The UniProtKB database and ClueGO were used to obtain functional Gene Ontology annotations and biological and functional networks, respectively.

RESULTS

Our analysis yielded 3,254 and 3,457 proteins identified in the endometrium of PS and NPS, respectively; of these, 1,753 being common while 1,501 and 1,704 were exclusive to PS and NPS, respectively. In addition, we identified 3,968 proteins in the extraembryonic membranes of PS. Further analyses of function revealed some proteins had relevance for the immune system process and biological adhesion in endometrium while the embryonic chorion displayed abundance of proteins related to cell adhesion and cytoskeletal organization, suggesting they dominated the moment of endometrial remodeling, implantation and adhesion of the lining epithelia. Data are available via ProteomeXchange with identifier PXD042565.

CONCLUSION

This is the first in-depth proteomic characterization of the endometrium and extraembryonic membranes during weeks 3 to 4 of gestation; data that contribute to the molecular understanding of the dynamic environment during this critical period, associated with the majority of pregnancy losses.

Genome-Wide Selection Signal Analysis to Investigate Wide Genomic Heredity Divergence between Eurasian Wild Boar and Domestic Pig

As important livestock species, pigs provide essential meat resources for humans, so understanding the genetic evolution behind their domestic history could help with the genetic improvement of domestic pigs. This study aimed to investigate the evolution of convergence and divergence under selection in European and Asian domestic pigs by using public genome-wide data. A total of 164 and 108 candidate genes (CDGs) were obtained from the Asian group (wild boar vs. domestic pig) and the European group (wild boar vs. domestic pig), respectively, by taking the top 5% of intersected windows of a pairwise fixation index (FST) and a cross population extended haplotype homozygosity test (XPEHH). GO and KEGG annotated results indicated that most CDGs were related to reproduction and immunity in the Asian group. Conversely, rich CDGs were enriched in muscle development and digestion in the European group. Eight CDGs were subjected to parallel selection of Eurasian domestic pigs from local wild boars during domestication. These CDGs were mainly involved in olfactory transduction, metabolic pathways, and progesterone-mediated oocyte maturation. Moreover, 36 and 18 haplotypes of INPP5B and TRAK2 were identified in this study, respectively. In brief, this study did not only improve the understanding of the genetic evolution of domestication in pigs, but also provides valuable CDGs for future breeding and genetic improvement of pigs.

A Comprehensive Sequencing Analysis of Testis-Born miRNAs in Immature and Mature Indigenous Wandong Cattle (Bos taurus)

Micro RNAs (miRNAs) have been recognized as important regulators that are indispensable for testicular development and spermatogenesis. miRNAs are endogenous transcriptomic elements and mainly regulate the gene expression at post-transcriptional levels; however, the key role of miRNA in bovine testicular growth is not clearly understood. Thus, supposing to unveil the transcriptomics expression changes in the developmental processes of bovine testes, we selected three immature calves and three sexually mature bulls of the local Wandong breed for testicular-tissue sample collection. The cDNA libraries of experimental animals were established for RNA-sequencing analysis. We detected the miRNA expression in testes by using high-throughput sequencing technology, and bioinformatics analysis followed. The differentially expressed (DE) data showed that 151 miRNAs linked genes were significantly DE between immature and mature bull testes. Further, in detail, 64 were significantly up-regulated and 87 were down-regulated in the immature vs. mature testes (p-value < 0.05). Pathway analyses for miRNA-linked genes were performed and identified JAG2, BCL6, CFAP157, PHC2, TYRO3, SEPTIN6, and BSP3; these genes were involved in biological pathways such as TNF signaling, T cell receptor, PI3KAkt signaling, and functions affecting testes development and spermatogenesis. The DE miRNAs including MIR425, MIR98, MIR34C, MIR184, MIR18A, MIR136, MIR15A, MIR1388 and MIR210 were associated with cattle-bull sexual maturation and sperm production. RT-qPCR validation analysis showed a consistent correlation to the sequencing data findings. The current study provides a good framework for understanding the mechanism of miRNAs in the development of testes and spermatogenesis.

Differential MicroRNA Expression in Porcine Endometrium Related to Spontaneous Embryo Loss during Early Pregnancy

Litter size is an important indicator to measure the production capacity of commercial pigs. Spontaneous embryo loss is an essential factor in determining sow litter size. In early pregnancy, spontaneous embryo loss in porcine is as high as 20–30% during embryo implantation. However, the specific molecular mechanism underlying spontaneous embryo loss at the end of embryo implantation remains unknown. Therefore, we comprehensively used small RNA sequencing technology, bioinformatics analysis, and molecular experiments to determine the microRNA (miRNA) expression profile in the healthy and arresting embryo implantation site of porcine endometrium on day of gestation (DG) 28. A total of 464 miRNAs were identified in arresting endometrium (AE) and healthy endometrium (HE), and 139 differentially expressed miRNAs (DEMs) were screened. We combined the mRNA sequencing dataset from the SRA database to predict the target genes of these miRNAs. A quantitative real-time PCR assay identified the expression levels of miRNAs and mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on differentially expressed target genes of DEMs, mainly enriched in epithelial development and amino acids metabolism-related pathways. We performed fluorescence in situ hybridization (FISH) and the dual-luciferase report gene assay to confirm miRNA and predicted target gene binding. miR-205 may inhibit its expression by combining 3′-untranslated regions (3′ UTR) of tubulointerstitial nephritis antigen-like 1 (TINAGL1). The resulting inhibition of angiogenesis in the maternal endometrium ultimately leads to the formation of arresting embryos during the implantation period. This study provides a reference for the effect of miRNA on the successful implantation of pig embryos in early gestation.

Integrated Analysis of Long Non-Coding RNA and mRNA Expression Profiles in Testes of Calves and Sexually Mature Wandong Bulls (Bos taurus)

The mRNAs and long non-coding RNAs axes are playing a vital role in the regulating of post-transcriptional gene expression. Thereby, elucidating the expression pattern of mRNAs and long non-coding RNAs underlying testis development is crucial. In this study, mRNA and long non-coding RNAs expression profiles were investigated in 3-month-old calves and 3-year-old mature bulls' testes by total RNA sequencing. Additionally, during the gene level analysis, 21,250 mRNAs and 20,533 long non-coding RNAs were identified. As a result, 7908 long non-coding RNAs (p-adjust < 0.05) and 5122 mRNAs (p-adjust < 0.05) were significantly differentially expressed between the distinct age groups. In addition, gene ontology and biological pathway analyses revealed that the predicted target genes are enriched in the lysine degradation, cell cycle, propanoate metabolism, adherens junction and cell adhesion molecules pathways. Correspondingly, the RT-qPCR validation results showed a strong consistency with the sequencing data. The source genes for the mRNAs (CCDC83, DMRTC2, HSPA2, IQCG, PACRG, SPO11, EHHADH, SPP1, NSD2 and ACTN4) and the long non-coding RNAs (COX7A2, COX6B2, TRIM37, PRM2, INHBA, ERBB4, SDHA, ATP6VOA2, FGF9 and TCF21) were found to be actively associated with bull sexual maturity and spermatogenesis. This study provided a comprehensive catalog of long non-coding RNAs in the bovine testes and also offered useful resources for understanding the differences in sexual development caused by the changes in the mRNA and long non-coding RNA interaction expressions between the immature and mature stages.

Screening of Potential Key Genes Related to Tubal Factor Infertility Based on Competitive Endogenous RNA Network

Background: The molecular biological mechanism of tubal factor infertility (TFI) is still unclear. Long noncoding RNAs (lncRNAs) are considered a major part of the competitive endogenous RNA (ceRNA) network and have attracted growing attention. Our study aimed to explore the regulatory mechanisms of lncRNAs associated with TFI and screen potential key genes related to TFI. Materials and Methods: Differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) were identified by comparing normal and TFI expression patterns of lncRNAs and mRNAs in eutopic endometrial tissues obtained from 3 normal and 3 TFI patients during implantation. These data were used to develop a protein-protein interaction (PPI) network of DEGs using the STRING online software. The identified DELs and DEGs were then used to construct a ceRNA network, and the Network Analyzer Tool Kit in Cytoscape was used to analyze the ceRNA network topology and stability. Finally, the overlapping genes present in both the ceRNA and PPI networks were selected as the potential key genes related to TFI. Results: Ninety-six DEGs (59 up and 37 down) and 45 DELs (28 up and 17 down) were identified. Thirty-four DEGs were mapped in a PPI network. A ceRNA network, including two lncRNAs (LINC00305 and DLX6-AS1), four microRNAs (hsa-miR-20b-5p, hsa-miR-17-5p, hsa-miR-107, and hsa-miR-24-3p), and four mRNAs (MAP3K3, HMGB3, FAM103A1, and TMEM209), was successfully constructed. Importantly, a potential key gene (TMEM209) related to TFI was identified. Conclusion: The construction of a ceRNA network related to TFI may help elucidate the regulatory mechanism by which genes and lncRNAs function as ceRNA networks. Importantly, TMEM209 may be further evaluated as potential therapeutic targets for TFI.

Identification of Differentially Expressed MicroRNAs and Their Potential Target Genes in Adipose Tissue from Pigs with Highly Divergent Backfat Thickness

Simple Summary

The role of microRNA in fat deposition is very important and not clearly understood. We detected 318 pig microRNAs (miRNAs), among high and low backfat tissue samples, by high throughput sequencing. Among them, 18 miRNAs were differentially expressed between the high and low backfat groups. Some of the differentially expressed miRNAs were involved mainly in lipid and carbohydrate metabolism, and glycan biosynthesis and metabolism. In addition, in silico analysis of the mRNA and miRNA transcriptomes, revealed possible regulatory relationships for fat deposition. In particular, three miRNA–mRNA pairs, miR-137–PPARGC1A, miR-141–FASN, and miR-122-5p–PKM, were identified as candidate key regulators of fat deposition. Our findings provide an important insight into miRNA expression patterns in backfat tissue of pig and new insights into the regulatory mechanisms of fat deposition in pig.

Abstract

Fatty traits are very important in pig production. However, the role of microRNAs (miRNAs) in fat deposition is not clearly understood. In this study, we compared adipose miRNAs from three full-sibling pairs of female Landrace pigs, with high and low backfat thickness, to investigate the associated regulatory network. We obtained an average of 17.29 million raw reads from six libraries, 62.27% of which mapped to the pig reference genome. A total of 318 pig miRNAs were detected among the samples. Among them, 18 miRNAs were differentially expressed (p-value < 0.05, |log₂fold change| ≥ 1) between the high and low backfat groups; 6 were up-regulated and 12 were down-regulated. Functional enrichment of the predicted target genes of the differentially expressed miRNAs, indicated that these miRNAs were involved mainly in lipid and carbohydrate metabolism, and glycan biosynthesis and metabolism. Comprehensive analysis of the mRNA and miRNA transcriptomes revealed possible regulatory relationships for fat deposition. Negatively correlated mRNA–miRNA pairs included miR-137–PPARGC1A, miR-141–FASN, and miR-122-5p–PKM, indicating these interactions may be key regulators of fat deposition. Our findings provide important insights into miRNA expression patterns in the backfat tissue of pig and new insights into the regulatory mechanisms of fat deposition in pig.