Characterization of the endometrial transcriptome in early diestrus influencing pregnancy status in dairy cattle after transfer of in vitro-produced embryos

Endometrial responsiveness to interferon-tau and its association with subsequent reproductive performance in dairy heifers

Our objectives were to evaluate the endometrial responsiveness of dairy heifers to an intrauterine infusion of recombinant bovine interferon-tau (rbIFN-τ) and to associate endometrial responses to rbIFN-τ with subsequent reproductive performance. In Experiments 1 and 2, cyclic heifers were enrolled in a 5-d CIDR Cosynch program for estrous synchronization, and blood sampling and ultrasonography examinations were performed on d 0, 4, 7, 11, and 14 of the estrous cycle. In Experiment 1, heifers were randomly assigned to receive an intrauterine infusion containing 2 µg of rbIFN-τ (rbIFN-τ = 19) or saline (CTRL = 19) into the uterine horn ipsilateral to the corpus luteum (CL) on d 14 of the estrous cycle. Six hours after the infusion, the infused uterine horn was flushed for sampling of the uterine luminal fluid (ULF) for analyses of composition, and the endometrium was biopsied for transcriptomics. In Experiment 2, 100 heifers received an intrauterine infusion of rbIFN-τ, and the same procedures for uterine sample collection described in Experiment 1 were performed. After the intrauterine test, heifers were enrolled in a breeding program and classified as highly fertile (HF; pregnant at first AI) or subfertile (SF; not pregnant at first AI). Statistical analyses were performed using linear regression models, which included the effects of treatment (Experiment 1: CTRL vs. rbIFN-τ) or fertility group (Experiment 2: HF vs. SF) and block of samples. Intrauterine infusion of rbIFN-τ increased the expression of classical interferon-stimulated genes in the endometrium (e.g., ISG15, MX1, OAS2, IRF9, and USP18), and an antiviral response was predicted to be the main downstream effect of the transcriptome changes. In addition, rbIFN-τ increased the abundance of cholesterol, glycerol, and the overall concentration of oxylipins in the ULF. Analysis of endometrial transcriptome between HF and SF heifers revealed important differences in the expression of proteins associated with cell signaling, metabolism, attachment, and migration, with a large representation of genes encoding extracellular matrix proteins. In general, differently expressed genes were expected to be downregulated by IFN-τ but seemed to fail to be downregulated in SF heifers, resulting in higher expression in SF compared with HF heifers. Subfertile heifers had lower concentrations of glycerol and an altered profile of oxylipins in the ULF, with lower abundance of oxylipins derived from arachidonic acid and dihomo-γ-linolenic acid, and greater abundance of oxylipins derived from linoleic acid. Measurements of ovarian function did not differ between groups and, therefore, did not influence the observed results in uterine biology. In conclusion, the endometrial responsiveness to IFN-τ is variable among individuals and associated with subsequent fertility of heifers, indicating that communication between conceptus and endometrium is critical for the uterine receptivity and survival of pregnancy.

Transcriptomics analysis of the bovine endometrium during the perioestrus period

During the oestrous cycle, the bovine endometrium undergoes morphological and functional changes, which are regulated by alterations in the levels of oestrogen and progesterone and consequent changes in gene expression. To clarify these changes before and after oestrus, RNA-seq was used to profile the transcriptome of oestrus-synchronized beef heifers. Endometrial samples were collected from 29 animals, which were slaughtered in six groups beginning 12 h after the withdrawal of intravaginal progesterone releasing devices until seven days post-oestrus onset (luteal phase). The groups represented proestrus, early oestrus, metoestrus and early dioestrus (luteal phase). Changes in gene expression were estimated relative to gene expression at oestrus. Ingenuity Pathway Analysis (IPA) was used to identify canonical pathways and functional processes of biological importance. A total of 5,845 differentially expressed genes (DEGs) were identified. The lowest number of DEGs was observed at the 12 h post-oestrus time point, whereas the greatest number was observed at Day 7 post-oestrus onset (luteal phase). A total of 2,748 DEGs at this time point did not overlap with any other time points. Prior to oestrus, Neurological disease and Organismal injury and abnormalities appeared among the top IPA diseases and functions categories, with upregulation of genes involved in neurogenesis. Lipid metabolism was upregulated before oestrus and downregulated at 48h post-oestrus, at which point an upregulation of immune-related pathways was observed. In contrast, in the luteal phase the Lipid metabolism and Small molecule biochemistry pathways were upregulated.

Mapping Expression Quantitative Trait Loci Targeting Candidate Genes for Pregnancy in Beef Cows

Despite collective efforts to understand the complex regulation of reproductive traits, no causative genes and/or mutations have been reported yet. By integrating genomics and transcriptomics data, potential regulatory mechanisms may be unveiled, providing opportunities to dissect the genetic factors governing fertility. Herein, we identified regulatory variants from RNA-Seq data associated with gene expression regulation in the uterine luminal epithelial cells of beef cows. We identified 4676 cis and 7682 trans eQTLs (expression quantitative trait loci) affecting the expression of 1120 and 2503 genes, respectively (FDR < 0.05). These variants affected the expression of transcription factor coding genes (71 cis and 193 trans eQTLs) and genes previously reported as differentially expressed between pregnant and nonpregnant cows. Functional over-representation analysis highlighted pathways related to metabolism, immune response, and hormone signaling (estrogen and GnRH) affected by eQTL-regulated genes (p-value ≤ 0.01). Furthermore, eQTLs were enriched in QTL regions for 13 reproduction-related traits from the CattleQTLdb (FDR ≤ 0.05). Our study provides novel insights into the genetic basis of reproductive processes in cattle. The underlying causal mechanisms modulating the expression of uterine genes warrant further investigation.

Harnessing Genomics and Transcriptomics Approaches to Improve Female Fertility in Beef Cattle-A Review

Female fertility is the foundation of the cow-calf industry, impacting both efficiency and profitability. Reproductive failure is the primary reason why beef cows are sold in the U.S. and the cause of an estimated annual gross loss of USD 2.8 billion. In this review, we discuss the status of the genomics, transcriptomics, and systems genomics approaches currently applied to female fertility and the tools available to cow-calf producers to maximize genetic progress. We highlight the opportunities and limitations associated with using genomic and transcriptomic approaches to discover genes and regulatory mechanisms related to beef fertility. Considering the complex nature of fertility, significant advances in precision breeding will rely on holistic, multidisciplinary approaches to further advance our ability to understand, predict, and improve reproductive performance. While these technologies have advanced our knowledge, the next step is to translate research findings from bench to on-farm applications.

Differential expression of pro- and anti-angiogenic factors in the endometrium between repeat breeder and normally fertile cows

In cattle, the establishment of appropriate endometrial vasculature during the estrous cycle is required for preparing a receptive endometrium. This study aimed to investigate 1) mRNA expression of potent pro- and anti-angiogenic factors, 2) protein localization of the anti-angiogenic factor thrombospondin (TSP), and 3) vascularity in the endometrium of repeat breeder (RB) and normally fertile (non-RB) cows. Caruncular and intercaruncular endometrium was collected from RB and non-RB cows during the luteal phase of the estrous cycle. RB cows had greater mRNA expression levels of TSP ligands (TSP1 and TSP2) and receptors (CD36 and CD47) than non-RB cows. Although the mRNA expression levels of most angiogenic factors did not change by repeat breeding, RB cows had greater mRNA expression of fibroblast growth factor receptor 1 (FGFR1), angiopoietin 1 (ANGPT1), and ANGPT2 and a less mRNA expression of vascular endothelial growth factor B (VEGFB) than non-RB cows. By immunohistochemistry, TSP1, TSP2, CD36, and CD47 were detected in the luminal epithelium, glandular epithelium, stromal cells, and blood vessels of the endometrium. Two indexes of vascularity, the number of blood vessels and the percentage of area stained positive for the von Willebrand factor, were lower in the endometrium of RB than in that of non-RB cows. These results demonstrate that RB cows have a greater expression of both ligands and receptors for the anti-angiogenic factor TSP and a reduced vascular distribution in the endometrium compared with non-RB cows, suggesting suppressed endometrial angiogenesis.

mRNA Signatures in Peripheral White Blood Cells Predict Reproductive Potential in Beef Heifers at Weaning

Reproductive failure is a major contributor to inefficiency within the cow-calf industry. Particularly problematic is the inability to diagnose heifer reproductive issues prior to pregnancy diagnosis following their first breeding season. Therefore, we hypothesized that gene expression from the peripheral white blood cells at weaning could predict the future reproductive potential of beef heifers. To investigate this, the gene expression was measured using RNA-Seq in Angus-Simmental crossbred heifers sampled at weaning and retrospectively classified as fertile (FH, n = 8) or subfertile (SFH, n = 7) after pregnancy diagnosis. We identified 92 differentially expressed genes between the groups. Network co-expression analysis identified 14 and 52 hub targets. ENSBTAG00000052659, OLR1, TFF2, and NAIP were exclusive hubs to the FH group, while 42 hubs were exclusive to the SFH group. The differential connectivity between the networks of each group revealed a gain in connectivity due to the rewiring of major regulators in the SFH group. The exclusive hub targets from FH were over-represented for the CXCR chemokine receptor pathway and inflammasome complex, while for the SFH, they were over-represented for immune response and cytokine production pathways. These multiple interactions revealed novel targets and pathways predicting reproductive potential at an early stage of heifer development.

Machine Learning-Based Co-Expression Network Analysis Unravels Potential Fertility-Related Genes in Beef Cows

Simple Summary

Female reproductive failure is still a challenge for the beef industry. Several biological processes that underlie fertility-related traits, such as the establishment of pregnancy and embryo survival, are still unclear. Increased availability of transcriptomic data has allowed a deep investigation of the potential mechanisms involved in fertility. This study investigated candidate gene biomarkers predictive of pregnancy status and underlying fertility-related networks. To this end, we integrated gene expression profiles through supervised machine learning and gene network modeling. We identified nine biologically relevant endometrial gene biomarkers that could discriminate against pregnancy status in cows. These biomarkers were co-expressed with genes critical for uterine receptivity, including endometrial tissue remodeling, focal adhesion, and embryo development. This study outlined key pathways involved with pregnancy success and provided predictive candidate biomarkers for pregnancy outcome in cows.

Abstract

Reproductive failure is still a challenge for beef producers and a significant cause of economic loss. The increased availability of transcriptomic data has shed light on the mechanisms modulating pregnancy success. Furthermore, new analytical tools, such as machine learning (ML), provide opportunities for data mining and uncovering new biological events that explain or predict reproductive outcomes. Herein, we identified potential biomarkers underlying pregnancy status and fertility-related networks by integrating gene expression profiles through ML and gene network modeling. We used public transcriptomic data from uterine luminal epithelial cells of cows retrospectively classified as pregnant (P, n = 25) and non-pregnant (NP, n = 18). First, we used a feature selection function from BioDiscML and identified SERPINE3, PDCD1, FNDC1, MRTFA, ARHGEF7, MEF2B, NAA16, ENSBTAG00000019474, and ENSBTAG00000054585 as candidate biomarker predictors of pregnancy status. Then, based on co-expression networks, we identified seven genes significantly rewired (gaining or losing connections) between the P and NP networks. These biomarkers were co-expressed with genes critical for uterine receptivity, including endometrial tissue remodeling, focal adhesion, and embryo development. We provided insights into the regulatory networks of fertility-related processes and demonstrated the potential of combining different analytical tools to prioritize candidate genes.

Co-Expression Network and Integrative Analysis of Metabolome and Transcriptome Uncovers Biological Pathways for Fertility in Beef Heifers

Reproductive failure remains a significant challenge to the beef industry. The omics technologies have provided opportunities to improve reproductive efficiency. We used a multistaged analysis from blood profiles to integrate metabolome (plasma) and transcriptome (peripheral white blood cells) in beef heifers. We used untargeted metabolomics and RNA-Seq paired data from six AI-pregnant (AI-P) and six nonpregnant (NP) Angus-Simmental crossbred heifers at artificial insemination (AI). Based on network co-expression analysis, we identified 17 and 37 hub genes in the AI-P and NP groups, respectively. Further, we identified TGM2, TMEM51, TAC3, NDRG4, and PDGFB as more connected in the NP heifers’ network. The NP gene network showed a connectivity gain due to the rewiring of major regulators. The metabolomic analysis identified 18 and 15 hub metabolites in the AI-P and NP networks. Tryptophan and allantoic acid exhibited a connectivity gain in the NP and AI-P networks, respectively. The gene–metabolite integration identified tocopherol-a as positively correlated with ENSBTAG00000009943 in the AI-P group. Conversely, tocopherol-a was negatively correlated in the NP group with EXOSC2, TRNAUIAP, and SNX12. In the NP group, α-ketoglutarate-SMG8 and putrescine-HSD17B13 were positively correlated, whereas a-ketoglutarate-ALAS2 and tryptophan-MTMR1 were negatively correlated. These multiple interactions identified novel targets and pathways underlying fertility in bovines.

Progesterone-dependent and progesterone-independent modulation of luminal epithelial transcription to support pregnancy in cattle

In cattle, starting 4-5 days after estrus, preimplantation embryonic development occurs in the confinement of the uterine lumen. Cells in the endometrial epithelial layer control the molecular traffic to and from the lumen and, thereby determine luminal composition. Starting early postestrus, endometrial function is regulated by sex steroids, but the effects of progesterone on luminal cells transcription have not been measured in vivo. The first objective was to determine the extent to which progesterone controls transcription in luminal epithelial cells 4 days (D4) after estrus. The second objective was to discover luminal transcripts that predict pregnancy outcomes when the effect of progesterone is controlled. Endometrial luminal epithelial cells were collected from embryo transfer recipients on D4 using a cytological brush and their transcriptome was determined by RNASeq. Pregnancy by embryo transfer was measured on D30 (25 pregnant and 18 nonpregnant). Progesterone concentration on D4 was associated positively (n = 182) and negatively (n = 58) with gene expression. Progesterone-modulated transcription indicated an increase in oxidative phosphorylation, biosynthetic activity, and proliferation of epithelial cells. When these effects of progesterone were controlled, different genes affected positively (n = 22) and negatively (n = 292) odds of pregnancy. These set of genes indicated that a receptive uterine environment was characterized by the inhibition of phosphoinositide signaling and innate immune system responses. A panel of 25 genes predicted the pregnancy outcome with sensitivity and specificity ranging from 64%-96% and 44%-83%, respectively. In conclusion, in the early diestrus, both progesterone-dependent and progesterone-independent mechanisms regulate luminal epithelial transcription associated with pregnancy outcomes in cattle.

Metabolites Secreted by Bovine Embryos In Vitro Predict Pregnancies That the Recipient Plasma Metabolome Cannot, and Vice Versa

This work describes the use of mass spectrometry-based metabolomics as a non-invasive approach to accurately predict birth prior to embryo transfer (ET) starting from embryo culture media and plasma recipient. Metabolomics was used here as a predictive platform. Day-6 in vitro produced embryos developed singly in modified synthetic oviduct fluid culture medium (CM) drops for 24 h were vitrified as Day-7 blastocysts and transferred to recipients. Day-0 and Day-7 recipient plasma (N = 36 × 2) and CM (N = 36) were analyzed by gas chromatography coupled to the quadrupole time of flight mass spectrometry (GC-qTOF). Metabolites quantified in CM and plasma were analyzed as a function to predict pregnancy at Day-40, Day-62, and birth (univariate and multivariate statistics). Subsequently, a Boolean matrix (F1 score) was constructed with metabolite pairs (one from the embryo, and one from the recipient) to combine the predictive power of embryos and recipients. Validation was performed in independent cohorts of ETs analyzed. Embryos that did not reach birth released more stearic acid, capric acid, palmitic acid, and glyceryl monostearate in CM (i.e., (p < 0.05, FDR < 0.05, Receiver Operator Characteristic—area under curve (ROC-AUC) > 0.669)). Within Holstein recipients, hydrocinnamic acid, alanine, and lysine predicted birth (ROC-AUC > 0.778). Asturiana de los Valles recipients that reached birth showed lower concentrations of 6-methyl-5-hepten-2-one, stearic acid, palmitic acid, and hippuric acid (ROC-AUC > 0.832). Embryonal capric acid and glyceryl-monostearate formed F1 scores generally >0.900, with metabolites found both to differ (e.g., hippuric acid, hydrocinnamic acid) or not (e.g., heptadecanoic acid, citric acid) with pregnancy in plasmas, as hypothesized. Efficient lipid metabolism in the embryo and the recipient can allow pregnancy to proceed. Changes in phenolics from plasma suggest that microbiota and liver metabolism influence the pregnancy establishment in cattle.

The incompletely fulfilled promise of embryo transfer in cattle-why aren't pregnancy rates greater and what can we do about it?

Typically, bovine embryos are transferred into recipient females about day 7 after estrus or anticipated ovulation, when the embryo has reached the blastocyst stage of development. All the biological and technical causes for failure of a female to produce a blastocyst 7 d after natural or artificial insemination (AI) are avoided when a blastocyst-stage embryo is transferred into the female. It is reasonable to expect, therefore, that pregnancy success would be higher for embryo transfer (ET) recipients than for inseminated females. This expectation is not usually met unless the recipient is exposed to heat stress or is classified as a repeat-breeder female. Rather, pregnancy success is generally similar for ET and AI. The implication is that either one or more of the technical aspects of ET have not yet been optimized or that underlying female fertility that causes an embryo to die before day 7 also causes it to die later in pregnancy. Improvements in pregnancy success after ET will depend upon making a better embryo, improving uterine receptivity, and forging new tools for production and transfer of embryos. Key to accelerating progress in improving pregnancy rates will be the identification of phenotypes or phenomes that allow the prediction of embryo competence for survival and maternal capacity to support embryonic development.

Machine learning approach to integrated endometrial transcriptomic datasets reveals biomarkers predicting uterine receptivity in cattle at seven days after estrous

The main goal was to apply machine learning (ML) methods on integrated multi-transcriptomic data, to identify endometrial genes capable of predicting uterine receptivity according to their expression patterns in the cow. Public data from five studies were re-analyzed. In all of them, endometrial samples were obtained at day 6–7 of the estrous cycle, from cows or heifers of four different European breeds, classified as pregnant (n = 26) or not (n = 26). First, gene selection was performed through supervised and unsupervised ML algorithms. Then, the predictive ability of potential key genes was evaluated through support vector machine as classifier, using the expression levels of the samples from all the breeds but one, to train the model, and the samples from that one breed, to test it. Finally, the biological meaning of the key genes was explored. Fifty genes were identified, and they could predict uterine receptivity with an overall 96.1% accuracy, despite the animal’s breed and category. Genes with higher expression in the pregnant cows were related to circadian rhythm, Wnt receptor signaling pathway, and embryonic development. This novel and robust combination of computational tools allowed the identification of a group of biologically relevant endometrial genes that could support pregnancy in the cattle.

Bovine in-vitro produced embryos: Development of embryo proper and associated membranes from day 26 to 47 of gestation

Based on in-vitro produced (IVP) bovine embryos, embryo proper and embryonic/fetal membranes were studied in 12 pregnancies from day 26 to 47. The embryos/fetuses displayed external as well as internal development of organs and structures according to the expectations from comparable in-vivo studies. However, the embryonic/fetal membranes were shorter than those reported for in-vivo-derived embryos/fetuses on days 26-35 of calculated age, whereas on days 41-47 they were of comparable lengths.

Inflammatory diseases in dairy cows: Risk factors and associations with pregnancy after embryo transfer

The objectives of the present prospective cohort study were to identify risk factors for inflammatory diseases in Holstein-Gyr crossbred dairy cows and characterize the associations of those diseases with pregnancy per embryo transfer (ET). Diseases were diagnosed in the first 60 d postpartum in 252 primiparous and 481 multiparous cows. Uterine diseases (UTD) included retained placenta, metritis, clinical endometritis, and subclinical endometritis. Nonuterine diseases (NUTD) included mastitis, lameness, pneumonia, and displaced abomasum. Blood was sampled on d 0, 1, and 2 postpartum and analyzed for concentrations of haptoglobin, fatty acids, total Ca (tCa), P, and Mg, and again on d 8 postpartum and analyzed for concentration of β-hydroxybutyrate. The association between concentrations of metabolites in serum and inflammatory diseases was determined. Cows received a timed ET program starting 28 ± 3 d postpartum with first ET at 46 ± 3 d postpartum using fresh in vitro-produced embryos. Pregnancy was diagnosed on d 31 and 59 of presumptive gestation. Overall, 63.3% of the cows were diagnosed with UTD and 20.6% with NUTD. The risk factors for UTD included season of calving, parity group, calving problems, days with subclinical hypocalcemia, and serum concentrations of haptoglobin and Mg, whereas the risk factors for NUTD were parity group and serum Mg concentration. Cows that developed UTD had increased concentrations of haptoglobin on d 2 and fatty acids on d 1 and 2, and reduced concentrations of tCa on d 1 and 2 and of P and Mg on d 2 postpartum compared with cows without UTD. Cows that developed NUTD had increased concentrations of fatty acids on d 0 to 2 postpartum, and decreased concentrations of tCa and P on d 0 and 1, and of Mg on d 1 and 2 postpartum compared with cows without NUTD. Cows that developed NUTD had a 340-kg reduction in milk yield in the first 60 d postpartum. Inflammatory diseases were associated with lesser body condition score and increased loss of body condition in the first 70 d postpartum. Maintenance of pregnancy after ET was reduced in UTD cows following the first (41.7 vs. 25.4%) or all ET (46.4 vs. 36.2%), whereas maintenance of pregnancy was reduced in NUTD cows only at the second ET (39.0 vs 25.9%). The reduced pregnancy maintenance in UTD cows combined with a reduced 21-d service rate (61.9 vs. 54.8%) decreased the 21-d cycle pregnancy rate (28.6 vs. 19.9%) and the hazard of pregnancy to 300 d postpartum by 35%, resulting in an extra 32 d open. In conclusion, inflammatory diseases depressed fertility in dairy cows receiving ET, with the greatest impact observed in UTD cows. This suggests that local inflammation of the uterus impairs maintenance of pregnancy in dairy cows following ET.