Endometrial gene expression profile of pregnant sows with extreme phenotypes for reproductive efficiency

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.

Transcriptomic Data Analyses Reveal That Sow Fertility-Related lincRNA NORFA Is Essential for the Normal States and Functions of Granulosa Cells

NORFA, the first lincRNA associated with sow fertility, has been shown to control granulosa cell (GC) functions and follicular atresia. However, the underlying mechanism is not fully understood. In this study, RNA-seq was performed and we noticed that inhibition of NORFA led to dramatic transcriptomic alterations in porcine GCs. A total of 1,272 differentially expressed transcripts were identified, including 1167 DEmRNAs and 105 DEmiRNAs. Furthermore, protein-protein interaction, gene-pathway function, and TF-miRNA-mRNA regulatory networks were established and yielded four regulatory modules with multiple hub genes, such as AR, ATG5, BAK1, CENPE, NR5A1, NFIX, WNT5B, ssc-miR-27b, and ssc-miR-126. Functional assessment showed that these hub DEGs were mainly enriched in TGF-β, PI3K-Akt, FoxO, Wnt, MAPK, and ubiquitin pathways that are essential for GC states (apoptosis and proliferation) and functions (hormone secretion). In vitro, we also found that knockdown of NORFA in porcine GCs significantly induced cell apoptosis, impaired cell viability, and suppressed 17β-estradiol (E2) synthesis. Notably, four candidate genes for sow reproductive traits (INHBA, NCOA1, TGFβ-1, and TGFBR2) were also identified as potential targets of NORFA. These findings present a panoramic view of the transcriptome in NORFA-reduced GCs, highlighting that NORFA, a candidate lincRNA for sow fertility, is crucial for the normal states and functions of GCs.

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

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

Seminal Plasma Modulates miRNA Expression by Sow Genital Tract Lining Explants

The seminal plasma (SP) modulates the female reproductive immune environment after mating, and microRNAs (miRNAs) could participate in the process. Considering that the boar ejaculate is built by fractions differing in SP-composition, this study evaluated whether exposure of mucosal explants of the sow internal genital tract (uterus, utero-tubal junction and isthmus) to different SP-fractions changed the profile of explant-secreted miRNAs. Mucosal explants retrieved from oestrus sows (n = 3) were in vitro exposed to: Medium 199 (M199, Control) or M199 supplemented (1:40 v/v) with SP from the sperm-rich fraction (SRF), the post-SRF or the entire recomposed ejaculate, for 16 h. After, the explants were cultured in M199 for 24 h to finally collect the media for miRNA analyses using GeneChip miRNA 4.0 Array (Affymetrix). Fifteen differentially expressed (False Discovery Rate (FDR) < 0.05 and Fold-change ≥ 2) miRNAs (11 down- versus 4 up-regulated) were identified (the most in the media of uterine explants incubated with SP from post-SRF). Bioinformatics analysis identified that predicted target genes of dysregulated miRNAs, mainly miR-34b, miR-205, miR-4776-3p and miR-574-5p, were involved in functions and pathways related to immune response. In conclusion, SP is able to elicit changes in the miRNAs profile secreted by female genital tract, ultimately depending SP-composition.

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

Simple Summary

Meishan pigs have a lower fetal loss rate during mid-late gestation compared to Duroc pigs. Differentially expressed mRNAs and miRNAs detected in endometrial tissue from Meishan and Duroc sows at mid-late gestation are involved in regulating hormone and oxygen levels, blood vessel development, and developmental processes affecting reproduction. In addition, ssc-miR-503 and ssc-miR-671-5p were shown to target the EGF and ESR1 genes, respectively. These findings provided an important resource for studying embryonic mortality during mid-late gestation in pigs.

Abstract

Embryo loss is a major factor affecting profitability in the pig industry. Embryonic mortality occurs during peri-implantation and mid-late gestation in pigs. Previous investigations have shown that the embryo loss rate in Meishan pigs is significantly lower than in commercial breeds. Most studies have focused on embryonic mortality during early gestation, but little is known about losses during mid-late gestation. In this study, we performed a transcriptome analysis of endometrial tissue in mid-late gestation sows (gestation days 49 and 72) sampled from two breeds (Meishan (MS) and Duroc (DU)) that have different embryo loss rates. We identified 411, 1113, 697, and 327 differentially expressed genes, and 14, 36, 57, and 43 differentially expressed miRNAs in four comparisons (DU49 vs. DU72, DU49 vs. MS49, DU72 vs. MS72, and MS49 vs. MS72), respectively. Subsequently; seven differentially expressed mRNAs and miRNAs were validated using qPCR. Functional analysis suggested the differentially expressed genes and miRNAs target genes mainly involved in regulation of hormone levels, blood vessel development, developmental process involved in reproduction, embryonic placenta development, and the immune system. A network analysis of potential miRNA-gene interactions revealed that differentially expressed miRNAs in Meishan pigs are involved in the response to estradiol and oxygen levels, and affect angiogenesis and blood vessel development. The binding site on ssc-miR-503 for epidermal growth factor (EGF) and the binding site on ssc-miR-671-5p for estrogen receptor α (ESR1) were identified using a dual luciferase assay. The results of this study will enable further exploration of miRNA-mRNA interactions important in pig pregnancy and will help to uncover molecular mechanisms affecting embryonic mortality in pigs during mid-late gestation.

Neonatal lactocrine deficiency affects the adult porcine endometrial transcriptome at pregnancy day 13

Reproductive performance of female pigs that do not receive sufficient colostrum from birth is permanently impaired. Whether lactocrine deficiency, reflected by low serum immunoglobulin immunocrit (iCrit), affects patterns of endometrial gene expression during the periattachment period of early pregnancy is unknown. Here, objectives were to determine effects of low iCrit at birth on the adult endometrial transcriptome on pregnancy day (PxD) 13. On the first day of postnatal life, gilts were assigned to high or low iCrit groups. Adult high (n = 8) and low (n = 7) iCrit gilts were bred (PxD 0), and humanely slaughtered on PxD 13 when tissues and fluids were collected. The endometrial transcriptome was defined for each group using mRNAseq and microRNAseq. Reads were mapped to the Sus scrofa 11.1 genome build. Mature microRNAs were annotated using miRBase 21. Differential expression was defined based on fold change (≥ ±1.5). Lactocrine deficiency did not affect corpora lutea number, uterine horn length, uterine wet weight, conceptus recovery, or uterine luminal fluid estrogen content on PxD 13. However, mRNAseq revealed 1157 differentially expressed endometrial mRNAs in high versus low iCrit gilts. Differentially expressed genes had functions related to solute transport, endometrial receptivity, and immune response. Six differentially expressed endometrial microRNAs included five predicted to target 62 differentially expressed mRNAs, affecting similar biological processes. Thus, lactocrine deficiency on the first day of postnatal life can alter uterine developmental trajectory with lasting effects on endometrial responses to pregnancy as reflected at the level of the transcriptome on PxD 13.

miR-181a promotes porcine granulosa cell apoptosis by targeting TGFBR1 via the activin signaling pathway

Activin/Smad3 signaling plays a pivotal role in follicle development and atresia. However, the precise mechanisms underlying this process are not yet fully understood. Herein, we identified miR-181a as a central component of activin/Smad3-mediated follicle atresia. miR-181a was strikingly upregulated in porcine atretic follicles, which induced the apoptosis of porcine granulosa cells (GCs) in vitro. Furthermore, the transforming growth factor-β type 1 receptor (TGFBR1) was confirmed as a direct target of miR-181a by bioinformatics analysis and luciferase assays. Transfection with an miR-181a agomir repressed the TGFBR1 mRNA and protein levels. In addition, TGFBR1 overexpression repressed GC apoptosis, whereas TGFBR1 inhibition promoted GC apoptosis. miR-181a overexpression downregulated the phosphorylation of Smad3 and blocked the activation of TGF-β signaling. Moreover, activin A downregulated miR-181a expression and upregulated the TGFBR1 and p-Smad3 protein levels. Collectively, these data suggest that miR-181a regulates porcine GC apoptosis by targeting TGFBR1 via the activin signaling pathway.

Morphological changes in the porcine cervix: A comparison between nulliparous and multiparous sows with regard to post-cervical artificial insemination

In recent decades, new artificial insemination (AI) methods, such as post-cervical AI (PCAI), have been developed in pig. PCAI involves crossing the cervix to deposit the sperm in the uterine body. Although PCAI application in sows is frequent, its application in nulliparous (gilts) females it is still limited due to the difficulty of passing through the cranial part of the cervical lumen. We hypothesized that ageing and parity would modify the cervical canal, facilitating the introduction of AI devices through the cervix. The aim was to compare the morphology of the uterus at different levels between multiparous and nulliparous females. Morphological analysis of the uterus pointed to a longer cervix (25.9 ± 4.6 vs. 21.6 ± 3.3 cm, p < 0.001) and greater length of the part of the reproductive tract involved in PCAI (from rima vulvae to the last cervical cushion) (56.2 ± 6.0 vs. 50.3 ± 5.2 cm, p < 0.001) in multiparous sows compared with nulliparous animals. As regards the structure of the vaginal and uterine parts of the cervix (the part in contact with the vagina and uterine body, respectively), the cross-sectional area, perimeter and total thickness were greater in the uterine part of multiparous than of nulliparous animals (area: 4.07 ± 1.46 vs. 2.46 ± 0.56 cm², p < 0.01; perimeter: 8.50 ± 1.44 cm vs. 6.28 ± 0.92 cm, p < 0.001; thickness: 10.79 ± 0.96 vs. 8.35 ± 0.62 mm, p < 0.05), but not in the vaginal part. The tissue content analysed in histological cross-sections also showed differences between female groups, a greater content of connective tissue (58.86 ± 10.78 vs. 67.60 ± 13.38%, p < 0.001) and a lower amount of muscle fibres (39.79 ± 10.24 vs. 30.66 ± 13.69%, p < 0.001) being observed in multiparous sows. Finally, silicone casts of the cervical lumen revealed differences between the two groups in the size and shape of the ridges in the lumen trajectory. Parity, which is also influenced by ageing, determines important changes in the size, structure and tissue content of the cervix wall, as well as in the morphology of the cervical canal, which may be responsible for the different levels of performance of PCAI in the female populations. Therefore, the future design of AI strategies and catheters should take into consideration the morphological variations of the cervix lumen, which will depend on age and parity of the females.

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein-coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.

Defining age- and lactocrine-sensitive elements of the neonatal porcine uterine microRNA-mRNA interactome

Factors delivered to offspring in colostrum within 2 days of birth support neonatal porcine uterine development. The uterine mRNA transcriptome is affected by age and nursing during this period. Whether uterine microRNA (miRNA) expression is affected similarly is unknown. Objectives were to (1) determine effects of age and nursing on porcine uterine miRNA expression between birth and postnatal day (PND) 2 using miRNA sequencing (miRNAseq) and; (2) define affected miRNA–mRNA interactions and associated biological processes using integrated target prediction analysis. At birth (PND 0), gilts were euthanized, nursed ad libitum, or gavage-fed milk replacer for 48 h. Uteri were collected at birth or 50 h postnatal. MicroRNAseq data were validated using quantitative real-time PCR. Targets were predicted using an established mRNA database generated from the same tissues. For PND 2 versus PND 0 comparisons, 31 differentially expressed (DE) miRNAs were identified for nursed, and 42 DE miRNAs were identified for replacer-fed gilts. Six DE miRNAs were identified for nursed versus replacer-fed gilts on PND 2. Target prediction for inversely correlated DE miRNA–mRNA pairings indicated 20 miRNAs targeting 251 mRNAs in nursed, versus 29 miRNAs targeting 585 mRNAs in replacer-fed gilts for PND 2 versus PND 0 comparisons, and 5 miRNAs targeting 81 mRNAs for nursed versus replacer-fed gilts on PND 2. Biological processes predicted to be affected by age and nursing included cell-to-cell signaling, cell morphology, and tissue morphology. Results indicate novel age- and lactocrine-sensitive miRNA–mRNA relationships associated with porcine neonatal uterine development between birth and PND 2.

EPS-LASSO: test for high-dimensional regression under extreme phenotype sampling of continuous traits

Motivation

Extreme phenotype sampling (EPS) is a broadly-used design to identify candidate genetic factors contributing to the variation of quantitative traits. By enriching the signals in extreme phenotypic samples, EPS can boost the association power compared to random sampling. Most existing statistical methods for EPS examine the genetic factors individually, despite many quantitative traits have multiple genetic factors underlying their variation. It is desirable to model the joint effects of genetic factors, which may increase the power and identify novel quantitative trait loci under EPS. The joint analysis of genetic data in high-dimensional situations requires specialized techniques, e.g. the least absolute shrinkage and selection operator (LASSO). Although there are extensive research and application related to LASSO, the statistical inference and testing for the sparse model under EPS remain unknown.

Results

We propose a novel sparse model (EPS-LASSO) with hypothesis test for high-dimensional regression under EPS based on a decorrelated score function. The comprehensive simulation shows EPS-LASSO outperforms existing methods with stable type I error and FDR control. EPS-LASSO can provide a consistent power for both low- and high-dimensional situations compared with the other methods dealing with high-dimensional situations. The power of EPS-LASSO is close to other low-dimensional methods when the causal effect sizes are small and is superior when the effects are large. Applying EPS-LASSO to a transcriptome-wide gene expression study for obesity reveals 10 significant body mass index associated genes. Our results indicate that EPS-LASSO is an effective method for EPS data analysis, which can account for correlated predictors.

Availability and implementation

The source code is available at https://github.com/xu1912/EPSLASSO.

Supplementary information

Supplementary data are available at Bioinformatics online.

Insertion/Deletion Within the KDM6A Gene Is Significantly Associated With Litter Size in Goat

A previous whole-genome association analysis identified lysine demethylase 6A (KDM6A), which encodes a type of histone demethylase, as a candidate gene associated to goat fecundity. KDM6A gene knockout mouse disrupts gametophyte development, suggesting that it has a critical role in reproduction. In this study, goat KDM6A mRNA expression profiles were determined, insertion/deletion (indel) variants in the gene identified, indel variants effect on KDM6A gene expression assessed, and their association with first-born litter size analyzed in 2326 healthy female Shaanbei white cashmere goats. KDM6A mRNA was expressed in all tissues tested (heart, liver, spleen, lung, kidney, muscle, brain, skin and testis); the expression levels in testes at different developmental stages [1-week-old (wk), 2, 3 wk, 1-month-old (mo), 1.5 and 2 mo] indicated a potential association with the mitosis-to-meiosis transition, implying that KDM6A may have an essential role in goat fertility. Meanwhile, two novel intronic indels of 16 bp and 5 bp were identified. Statistical analysis revealed that only the 16 bp indel was associated with first-born litter size (P < 0.01), and the average first-born litter size of individuals with an insertion/insertion genotype higher than that of those with the deletion/deletion genotype (P < 0.05). There was also a significant difference in genotype distributions of the 16 bp indel between mothers of single-lamb and multi-lamb litters in the studied goat population (P = 0.001). Consistently, the 16 bp indel also had a significant effect on KDM6A gene expression. Additionally, there was no significant linkage disequilibrium (LD) between these two indel loci, consistent with the association analysis results. Together, these findings suggest that the 16 bp indel in KDM6A may be useful for marker-assisted selection (MAS) of goats.

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.