Skip to main content
Download PDF

RNA-Sequencing based analysis of bovine endometrium during the maternal recognition of pregnancy

BMC Genomics volume 23, Article number: 494 (2022) Cite this article

Abstract

Background

Maternal recognition is the crucial step for establishing pregnancy in cattle. This study aims to identify endometrial genes and biological pathways involved in the maternal recognition of pregnancy. Caruncular endometrial tissues were collected from Day 15–17 of gestation (pregnant), non-pregnant (absence of conceptus), and cyclic (non-bred) heifers.

Results

Total RNAs were isolated from the caruncular endometrial tissues of pregnant, non-pregnant, and cyclic heifers, and were subjected to high-throughput RNA-sequencing. The genes with at least two-fold change and Benjamini and Hochberg p-value ≤ 0.05 were considered differentially expressed genes and further confirmed with quantitative real-time PCR. A total of 107 genes (pregnant vs cyclic) and 98 genes (pregnant vs non-pregnant) were differentially expressed in the pregnant endometrium. The most highly up-regulated genes in the pregnant endometrium were MRS2, CST6, FOS, VLDLR, ISG15, IFI6, MX2, C15H11ORF34, EIF3M, PRSS22, MS4A8, and TINAGL1. Interferon signaling, immune response, nutrient transporter, synthesis, and secretion of proteins are crucial pathways during the maternal recognition of pregnancy.

Conclusions

The study demonstrated that the presence of conceptus at Day 15–17 of gestation affects the endometrial gene expression related to endometrial remodeling, immune response, nutrients and ion transporters, and relevant signaling pathways in the caruncular region of bovine endometrium during the maternal recognition of pregnancy.

Peer Review reports

Introduction

Beef cattle production is an important source of protein to meet the nutritional needs of a growing population. Improvements to beef cattle reproduction can help increase beef production worldwide to meet the increasing demand [1,2,3]. Early pregnancy failure is one of the critical factors that affect the economic output of the beef industry [4]. In ruminants, the successful establishment of pregnancy requires an intricate dialogue between the uterus and growing conceptus. The majority of pregnancy losses occur in the first month, especially around Day-19 of gestation, mainly due to the inability of the uterus to support conceptus growth and development or poor embryonic development. Understanding uterine changes during early pregnancy provide critical insight into reproductive success.

The endometrium undergoes dynamic changes during the peri-implantation period and provides the biological environment for embryonic growth and development [5]. The bovine endometrium consists of caruncles (aglandular) and intercaruncular tissue (glandular). The caruncle develops the vascular bed and is the site for embryo implantation and metabolic exchange. The endometrial secretions that support conceptus elongation are produced from the luminal and glandular epithelium [6]. In the ruminant, the fertilized oocyte undergoes a series of morphological and biochemical changes as a conceptus in the oviduct and uterus and begins to elongate between Days 12–14 of gestation [7]. By Day 15–17 of gestation, the conceptus develops into a filamentous form and produces interferon tau (IFNτ), which acts as the signal for the maternal recognition of pregnancy [8, 9]. The conceptus derived IFNτ promotes the persistence of the corpus luteum required for adequate progesterone production [10]. It is well-known that progesterone induces endometrial transcriptomes during the peri-implantation periods [10,11,12]. IFNτ also induces endometrial genes and proteins required for immunomodulation, extracellular matrix remodeling, and implantation-specific molecules [13, 14]. Recent studies have suggested that bovine embryos around the peri-implantation period induce endometrial gene expression in the intercaruncular region required to establish gestation [4, 7, 15]. Despite these studies, the transcriptomic changes in the caruncular portion of endometrial tissue during the maternal recognition of pregnancy (Day 15–17 of gestation) are not completely understood. Most previous studies have investigated the conceptus-induced gene expression in the caruncular and intercaruncular endometrial tissues and compared it with cyclic cows around the peri-implantation period [16,17,18], and yet, caruncular endometrial transcriptomes involved in the maternal recognition of pregnancy are not clearly understood. To further improve the conception rate in cattle, the knowledge of specific genes, proteins, and biological pathways during the maternal recognition of pregnancy is required throughout the uterus. As Day-15–17 of gestation is a critical period for the maternal recognition and establishment of pregnancy, we hypothesized that RNA-Sequencing based analysis of bovine caruncular endometrial tissues during the maternal recognition of pregnancy (Day-15–17 of gestation) would reveal important genes and biological pathways required for the maternal recognition of pregnancy. This study analyzed the genes and biological pathways in the caruncular endometrium during the maternal recognition of pregnancy (pregnant vs. cyclic) and (pregnant vs. non-pregnant) using RNA-Sequencing, and Real-time PCR (qPCR).

Methods

Animals and sampling

Animal husbandry, management, and handling procedures were under the Guide for the Care and Use of Agricultural Animals in Research and Teaching (Ag Guide 2020) [19]. Angus heifers (2–3 years old; n = 21) grazing Kikuyu grass (Pennisetum clandestinum) and Pangola (Digitaria eriantha) pastures were used for sampling. The estrous cycles of the heifers (n = 21) were synchronized using 25 mg of prostaglandin F2 alpha (PGF2 α; Lutalyse®, Zoetis, Parsippany, New Jersey, USA) administered intramuscularly on Day-1 (Day-1 designated as the first dose of PGF2 α) and Day-11 (Day-11 designated as the second dose of PGF2 α). Day Fifteen heifers were bred after detecting estrus. Cows were identified as pregnant (presence of conceptus) or non-pregnant (absence of conceptus). After incision of the uterus, the lumen of the uterus was exposed. Caruncles were identified as the small protuberances from the surface of endometrium, and carefully collected the protruded endometrial areas (4–5/heifer) as previously collected [20]. Caruncular endometrial tissues were collected on Day 15–17 of gestation (pregnant; n = 8) or absence of conceptus (non-pregnant; n = 7) or non-bred heifers (cyclic; n = 6) and were stored at -80 °C until further use.

RNA isolation and quality control

Total RNA was isolated from frozen tissues (60–100 mg) using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions as previously described [21]. The total RNA concentration was determined using NanoDrop One (Thermo Fisher Scientific, Madison, WI). RNA quality was determined with the Agilent 2100 Bioanalyzer (Agilent Technologies, Massy, France). The samples with an RNA integrity number (RIN) > 7 were further used for RNA-Sequencing and quantitative real-time PCR. The RNA was stored at -80 °C until further use for RNA-sequencing and quantitative real-time PCR.

Library preparation and RNA sequence analysis

RNA-Seq libraries were prepared and sequenced at the University of Hawaii Cancer Center Genomics and Bioinformatics Shared Resource (UHCC GBSR) facility. A TruSeq Stranded mRNA kit (Illumina, San Diego, CA) was used to prepare the RNA-Seq libraries from total RNA samples extracted from bovine endometrium, including pregnant (n = 5), non-pregnant (n = 5), and cyclic (n = 5). Libraries were prepared according to the manufacturer’s protocol as previously described [21].

Data analysis of the RNA sequences were done at the University of Hawaii John A. Burns School of Medicine Bioinformatics Core Facility. Single-end reads in the FASTQ format were explored using FastQC (Babraham Institute, Cambridge, UK) and cleaned using Prinseq, a Perl script [22, 23]. The cleaning procedure included trimming low-quality reads from both 3’ and 5’ ends until a base pair of Phred quality score of 30 (99.9% accurate) or greater was found and filtering out reads having a mean quality score less than 30 and length below 30 nucleotides. Cleaned reads were aligned against the bovine reference genome (Bos_taurus.ARS-UCD1.2) using HiSAT2. The resulting SAM files were sorted, and converted to BAM files using SAMtools. Read counts mapped to bovine gene models were generated using htseq-count script from HTSeq package. Finally, bioconductor DESeq2 was used to get the differentially expressed genes among pregnant vs. non-pregnant (P vs. NP), pregnant vs. cyclic (P vs. C), and non-pregnant vs. cyclic (NP vs. C) groups (n = 4/group). In RNA Sequencing, genes having fold change (FC) greater than 2 in the endometrial sample and Benjamini and Hochberg q-value < 0.05 were considered differentially expressed.

Pathways analyses

The ingenuity pathway analysis (IPA) is a human genome-based powerful search tool with several advanced functions that allows insightful data analysis and interpretation. The differentially expressed genes (DEGs) were subjected to the IPA (QIAGEN, Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis) to gain insights into the canonical pathways and network discovery.

Functional annotation and gene ontology enrichment analysis

Functional and pathway analysis was carried out using an open web source named Enrichr (https://maayanlab.cloud/Enrichr/) to gain insight into the various Gene Ontology (GO) terms of the genes in bovine endometrium. The official gene symbol of the up-regulated genes was uploaded to the functional annotation tool in the Enrichr system, and the Bos taurus was selected as the reference genome. The genes that matched up with the genes in Enrichr were annotated into three GO terms: biological process, cellular component, and molecular function. All the GO terms were considered enriched at a modified P-value < 0.05 and a threshold gene count of 2.

Kyoto Encyclopedia of Genes and Genomes (KEGG)

The pathways enrichment for the up-regulated genes in the bovine endometrium using the Kyoto Encyclopedia of Genes and Genomes were analyzed [24]. The official gene symbol of the up-regulated genes was uploaded to the functional annotation tool in the Enrichr system, and the bovine was selected as the reference genome. The enrichment parameters were set to a threshold gene count of 2 and a modified Fisher Exact P-value < 0.05. The over-represented KEGG pathways terms were considered as enriched KEGG pathways.

Quantitative real-time PCR (qPCR)

Among the most highly up-regulated genes, fourteen candidate genes (MRS2, CST6, FOS, VLDLR, ISG15, IF16, MX2, C15H11ORF34, EIF3M, PENK, PRSS22, MS4A8, TINAGL1, and R3HDM1) were selected for validation using qPCR. Primers specific to each gene were designed using the NCBI primer blast tool (Supplementary Table S1). Total RNA (1 μg) was reversed transcribed into cDNA using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, California, USA). The newly synthesized cDNA (20 μL) was diluted (20X), and 3 μl per qPCR reaction was used.

The qPCR assay was performed in a 10 μL reaction mixture containing 3 μL of cDNA and 7 μL of PCR mix using QuantStudio™ 3 System (Applied Biosystems). The PCR mixture was prepared by adding 5 μL of PowerUp SYBR Green Master Mix (Applied Biosystems) and 1 μL each of forward and reverse primers specific to the target gene. The PCR mix and cDNA samples were loaded into a 96-well optical plate and were sealed with clear optical adhesive films (Applied Biosystems) as previously described [21]. The specificity of each primer was validated by running the melting curve, and qPCR products were assessed using gel electrophoresis. To determine the most stable housekeeping gene in the endometrial tissues, the expression of glyceraldehyde 3-phosphate dehydrogenase (GAPDH), beta-actin (β-actin), and TATA-Box Binding Protein (TBP) were analyzed in triplicates across the samples. β-actin was the most uniform housekeeping gene. The target genes were analyzed in triplicates, and the expression level was determined using the cycle threshold (Ct) values after normalization with β-actin. The fold change for each gene was calculated using the comparative CT method (2−ΔΔCt method) [25]. Data for fold change were presented as a mean ± standard error on the bar diagram. Values were subjected to a one-way analysis of variance (ANOVA) followed by the Tukey HSD test for mean separation and comparison to determine differences between the treatments using R Studio. Differences were considered significant at a p-value < 0.05.

Results

RNA sequencing-based differentially expressed genes in the pregnant bovine endometrium during the maternal recognition of pregnancy

Using RNA-Seq, the transcriptomics profile was analyzed in pregnant endometrium (Day 15–17 of gestation) compared to non-pregnant and cyclic cows. Raw sequencing reads in the FASTQ format were obtained from the replicated RNA-Seq libraries and evaluated for their qualities using FastQC. There was an average of 19.5 M, 23.6 M, and 20.6 M original raw reads in pregnant (P), non-pregnant (NP), and cyclic (C) cows, respectively. All the groups (P, NP, and C) had excellent quality sequences (> 96%). Mapping results of the bovine genome database showed that an average of 93.3% of the retained reads from pregnant, 94.3% from non-pregnant, and 94.2% from cyclic were uniquely mapped. A total of 27,270 gene transcripts were annotated using the Ensemble alignment of the bovine genome assembly. Differential expression analysis between pregnant, non-pregnant, and cyclic cows was conducted (DESeq2). A total of 98 genes were differentially expressed between pregnant and non-pregnant cows, 107 genes were differentially expressed (DE) between pregnant and cyclic, and 41 genes were differentially expressed between cyclic and non-pregnant (Fig. 1). Of the 98 genes found in non-pregnant endometrium, 47.9% (47) were uniquely expressed in pregnant cattle, and of the 107 DE genes in pregnant cattle, 46% (50) were uniquely expressed in pregnant cattle (Fig. 1). Only 23% of DE genes were shared between pregnant and non-pregnant cattle (Fig. 2).

Fig. 1
figure 1

Venn diagram comparing the number of identified genes in the endometrium of pregnant (P) vs. cyclic (C), pregnant (P) vs. non-pregnant (NP), and non-pregnant (NP) vs. cyclic (C) cows along with the list of the name of the genes. Number of differentially expressed genes (DEGs) in P vs. C (n = 107), P vs. NP (n = 98), and NP vs. C (n = 41)

Full size image
Fig. 2
figure 2

Heat-map of the 35 up-regulated genes in the bovine endometrium compared with pregnant versus non-pregnant (P vs. NP) and pregnant versus cyclic (P vs. C). Extremes of green color indicate the gene's higher expression, while blue color indicates a lesser expression of genes

Full size image

To characterize the conceptus-induced transcriptomic changes in the pregnant endometrium, we shorted out the most highly up-regulated or down-regulated genes in the pregnant endometrium compared to cyclic and non-pregnant. The highly up-regulated genes in the pregnant endometrium (vs. cyclic) are MRS2, CST6, FOS, VLDLR, ISG15, IF16, MX1, MX2, C15H11ORF34, EIF3M, and TINAGL1 (Table 1). Among these genes, MRS2 was most highly expressed followed by CST6, FOS, VLDLR, and ISG15. The highly down-regulated genes in the pregnant endometrium (vs. cyclic) are LRFN4, PKMYT1, ZNHIT2, VARS2, ARL2BP, LINGO2, SLC16A11, SLC16A4, TMEM151B, and NME7 (Table 2). Among these genes, LRFN4 was highly down-regulated followed by PKMYT1, ZNHIT2, VARS2, and ARL2BP. The highly up-regulated genes in the pregnant endometrium (vs. non-pregnant) are ISG15, IFI6, PENK, PRSS22, MS4A8, CLDN4, C15H11ORF34, MRS2, TINAGL1, and R3HDM1 (Table 3). Among these genes, ISG15 was highly expressed followed by genes IFI6, PENK, PRSS22, and MS4A8. The highly down-regulated genes in the pregnant endometrium (vs. non-pregnant) are SNX20, PACSIN1, NTN3, COL1A1, SLC16A4, TERF1, NOC2L, COL3A1, COL6A3, and SPARCL1 (Table 4). Among these genes, SNX20 was highly down-regulated followed by PACSIN1, NTN3, COL1A1, and SLC16A4. However, some of the genes such as ISG15, IFI6, MRS2, MX2, C15H11ORF34, and TINAGL1 were common in the pregnant endometrium as compared to both cyclic and pregnant endometrium indicating the crucial involvement of these genes in the maternal recognition of pregnancy.

Table 1 Up-regulated genes in the bovine pregnant endometrium (P vs. C)
Full size table
Table 2 Down-regulated genes in the bovine pregnant endometrium (P vs. C)
Full size table
Table 3 Up-regulated genes in the bovine pregnant endometrium (P vs. NP)
Full size table
Table 4 Down-regulated genes in the bovine pregnant endometrium (P vs. NP)
Full size table

Functional annotation and pathways enrichment analysis of DEGs (2)

The gene ontology analysis demonstrates Type-1 interferon signaling, immune response, and extracellular matrix organization were important functional pathways observed in the biological process. Similarly, ion transporters such as SLC34A2, SLC2A1, and SLC16A11 were important in the molecular functions. The cellular component functions on the endoplasmic reticulum lumen were governed by genes such as WNT5B, IL23A1, and PENK (Table 5).

Table 5 Gene Ontology analysis in the bovine pregnant endometrium versus cyclic and nonpregnant
Full size table

The pathways enrichment for the up-regulated genes in the bovine endometrium using the Kyoto Encyclopedia of Genes and Genomes (KEGG) were analyzed [24]. The official gene symbol of the up-regulated genes was uploaded to the functional annotation tool in the Enrichr system, and the bovine was selected as the reference genome. The enrichment parameters were set to a threshold gene count of 2 and a modified Fisher Exact p-value < 0.05. The over-represented KEGG pathways terms were considered as enriched KEGG pathways.

In the KEGG pathway, both groups (P vs. C) and (P vs. NP) had some common pathways, i.e., the mineral absorption pathway. On the other hand, Th17 cell differentiation, endocrine, and other factor-regulated calcium reabsorption and progesterone-mediated oocyte maturation pathways were highly enriched in P vs. C (Table 6), whereas extracellular matrix (ECM) receptor interactions, C-type lectin receptor signaling pathway, and IL-17 signaling pathway were enriched in P vs. NP females (Table 7). Collagen genes were found abundantly in bovine pregnant endometrium.

Table 6 KEGG Pathway (P vs. C)
Full size table
Table 7 KEGG Pathway (P vs. NP)
Full size table

In the ingenuity canonical pathways, two pathways were common in both groups: interferon signaling and IL-12 signaling and production in macrophages. Other uncommon pathways to the group were the Th17 activation pathway, MIF regulation of innate immunity, and IL-23 signaling pathway enriched in P vs. C (Table 8). In contrast, oxidative phosphorylation, GP6 signaling pathway, and inhibition of matrix metalloproteases were pathways that were more highly expressed in the P vs. NP females. (Table 9). In the IPA network, lipid metabolism molecules such as Alp, CDH1, COLQ, EIF3M, EIF4A1, EIF4A3, ELOA, EPAS1, FARP1, Fgf, HELZ, HISTONE, Histoneh3, HNF1A, Insulin, KMT2E, mediator; MMP2, NFIA, NOC2L, PTEFb, POLR2B, Proinsulin, RBBP4, RNA polymerase II, Rnr, RPH3AL, RPSA, SKIDA1, SLC16A4, TCF/LEF, TMEM132A, TMPRSS2, ZFC3H1, and ZNHIT2 were significantly enriched (Fig. 3).

Table 8 Ingenuity Canonical Pathway (P vs. C)
Full size table
Table 9 Ingenuity Canonical Pathway (P vs. NP)
Full size table
Fig. 3
figure 3

Lipid Metabolism in the bovine endometrium between A pregnant vs. cyclic, and B pregnant vs. non-pregnant. The network is displayed graphically as nodes (genes). The node color intensity indicates the expression of genes, with red representing up-regulation and green with down-regulation in the pregnant endometrium

Full size image

The RNA-Seq data identified the differentially expressed genes in the pregnant bovine endometrium. Among the most highly up-regulated genes, fourteen candidate genes (MRS2, CST6, FOS, VLDLR, ISG15, IF16, MX2, C15H11ORF34, EIF3M, PENK, PRSS22, MS4A8, TINAGL1, and R3HDM1) were selected for validation using qPCR. The results of relative fold change for candidate genes obtained from qPCR are shown in Fig. 4. MRS2, MS4A8, PRSS22, CST6, VLDLR, IFI6, C15H11ORF34, ISG15, TINAGL1, and MX2 were significantly higher (p < 0.05) in the pregnant endometrium compared to NP and C, whereas R3DHM1, EIF3M, FOS, and PENK remained unchanged.

Fig. 4
figure 4

Validation of the gene expression in the endometrium. The fold changes were normalized with the β-actin. Data are represented as the mean ± standard error. The x-axis represents the different physiological statuses of the cows (C; cyclic, NP; non-pregnant, and P; pregnant cattle). Y-axis represents relative fold change for gene expression. All the fold change values were subjected to a one-way analysis of variance followed by the Tukey HSD test for mean separation and comparison to determine differences among the groups. Differences were considered significant at a *p-value < 0.05”

Full size image

Discussion

In ruminants, the successful establishment of pregnancy requires the intricate dialogue between the uterus and growing conceptus. During the maternal recognition of pregnancy, a conceptus-derived signal (IFNτ) leads to the persistence of the corpus luteum and induces the endometrial transcript to establish gestation [7]. Previous studies have identified several genes induced by IFNτ in the caruncular and intercaruncular endometrium during the peri-implantation period [16,17,18]. In the present study, both RNA-Seq and qPCR analysis confirmed the differential expression of several pre-discovered and novel genes and their biological pathways in the pregnant endometrium compared to cyclic and non-pregnant endometrium. Interferon signaling, immune response, nutrient transporter, synthesis, and secretion of proteins are crucial pathways during the maternal recognition of pregnancy. This study found some important molecules such as Type-1 interferon signaling (MX1, MX2, IF16, IRF1, and ISG15), ion transporters (SLC34A2, SLC2A1, SLC16A11, and SLC16A4), ECM organization (COL1A1, COL1A2, and COL3A1), and novel genes (MRS2, C15H11ORF34).

Interferon signaling pathway

The interferon signaling pathway is important in the pregnant endometrium around the peri-implantation period [7, 21]. In this study, the interferon signaling pathway was highly enriched in the pregnant endometrium. Under interferon signaling, we identified several genes such as MX2, IFI6, IFIT1, IRF1, IRF9, ISG15, MX1, STAT1, and TAP1. ISG15 is among the highly up-regulated genes in the pregnant endometrium. Gene Ontology and Ingenuity pathway analysis from our study revealed that ISG15 functions in the interferon signaling activation pathway in the pregnant endometrium, which is very important for the maternal recognition of the pregnancy. Previous studies have detected the ISG15 expression in the stromal cells and glandular epithelial cells [6, 12]. ISG15 expression is linked with several activities such as gene transcription, DNA repair, signal transduction, apoptosis, and cell cycle by conjugating itself to those target proteins [26]. ISG15 can regulate the innate immunity of embryonic cells. IFI6 was highly up-regulated in the pregnant endometrium and is part of the interferon signaling pathway. IFI6 is not mapped yet in cattle, but the human IFI6 gene is located on HSAP1p35, homologous to a chromosomal region between cattle and humans [27]. The timing of the up-regulation of ISGs, such as IFI6 in pregnant heifers was observed in previous studies [12]. Our study and previously reported suggest that IFI6 signaling is required to maintain the endometrial immune status required for embryonic survival and maternal recognition of pregnancy. In the present study, MX2 was among the highly up-regulated genes. MX2, which is recognized as an intracellular antiviral protein, belongs to a large GTPase family. It takes part in the interferon signaling pathway and innate immune system. In our study, GO annotations related to this gene include GTP binding and GTPase activity. MX2 is up-regulated in response to IFNτ from the elongating conceptus [7]. MX2 expression in response to elongating conceptus is consistent across mammalian species, including cattle, sheep, and humans. The MX2 expression was increased in the pregnant endometrium of ewes [28] and cows [29] in response to IFNτ. It is reported that MX2 mRNA is detectable in the peripheral blood lymphocytes with higher intensity at Day 16 of gestation than in non-pregnant cows [29]. Further, gene ontology analysis suggested that MX2 is an important antiviral gene in the uterus during the preimplantation period.

Extracellular matrix signaling

Extensive extracellular matrix and cellular remodeling occur in the endometrium during the estrous cycle, peri-implantation period, and different gestation stages [20, 30,31,32,33,34]. Collagen, the most abundant extracellular protein in mammals, is the main structural protein in the extracellular matrix in the various connective tissues [33]. In this study, GO analysis detected several collagen genes (COL1A1, COL1A2, and COL3A1) involved in the extracellular matrix organization in the pregnant bovine endometrium. We also found that endometrial collagen genes are predicted in platelet-derived growth factor(s). The KEGG pathway revealed many collagen genes in the pregnant endometrium having different functions. The genes involved in protein digestion and absorption (COL1A1, COL3A1, COL1A2, and COL6A3), ECM-receptor interaction (COL1A1, COL1A2, and COL6A3), and platelet activation (COL1A1, COL3A1, and COL1A2). Collagen genes have an essential role in cell adhesion. Some of the important collagen genes that help in adhesion are COL1A1, COL1A2, and COL6A3. According to our findings from the Ingenuity Canonical Pathways analysis, COL12A1, COL1A1, COL1A2, COL3A1, COL4A2, COL5A1, COL6A1, and COL6A3 are involved in the Glycoprotein 6 (GP 6) signaling pathway.

Matrix metalloproteinases (MMP) are known to degrade the ECM for cellular proliferation, differentiation, migration, and apoptosis [30,31,32,33,34]. In our study, ingenuity canonical pathway identified the genes (MMP-2, MMP-14, SDC1, TIMP2). It is well-known that MMP-14 regulates the MMP-2 through binding TIMP-2. This MMP cascade regulated the endometrial cell functions required for embryo implantation [30,31,32,33,34]. Syndecan-1 (SDC1/CD138)) is an integral membrane protein and takes part in cell proliferation, cell migration, and cell–matrix interactions [35]. Although cyclic and non-pregnant cows lacked a conceptus, interestingly, in our study, we found some molecules such as COL1A1, TRADD, and COL3A1 were up-regulated in both non-pregnant and pregnant endometrium. One reason could be that those collagen genes (COL1A1 and COL3A1) are abundantly present in platelet-derived growth factor binding [36]. In contrast, TRADD is essential for death-inducing signaling, which is likely a pathway important for pregnancy and nonpregnancy in cattle [37].

Cystatin 6 or Cystatin E/M (CST6) is among the highly up-regulated genes in the pregnant endometrium. GO annotations related to this gene include cysteine-type endopeptidase inhibitor activity. The functions of the CST6 include uterine endometrial and placental tissue remodeling and facilitating transplacental transport of nutrients. CST6 expression was detected in the endometrium during the estrous cycle and pregnancy. The expression of CST6 in chorionic epithelia of the placental membrane was increasing during late pregnancy, suggesting that cell type-specific expression and function of CST6 are critical for appropriate maternal–fetal interaction [38].

Ions transporter signaling pathways

Solute carrier (SLC) genes, consisting of 52 families, are mostly located in the cell membrane and code for membrane transport proteins [39]. The function of the SLC gene includes the transport of glucose, electrolytes, and amino acids. Since numerous nutrients and electrolytes must be transported from the blood to the uterine environment, SLC genes play a critical role in the bovine endometrium. SLC2A1 was among the top 20 most up-regulated SLC genes in the pregnant endometrium. Conversely, SLC16A11 (monocarboxylate transporter), SLC16A4 (monocarboxylate transporter), and SLC7A4 (cationic amino acid transporter/glycoprotein- associated amino acid) were found to be down-regulated in the pregnant endometrium. SLC46A1 was up-regulated in the non-pregnant endometrium. GO analysis showed SLC34A2 (Type-II Na + /HPO42- co-transporter), SLC2A1(glucose transporter), SLC16A11 (monocarboxylate transporter), and SLC16A4 (monocarboxylate transporter) involved in ion transportation [40]. SLC34A2 was included in parathyroid hormone synthesis, secretion, and action. This gene was also found in the mineral absorption pathway. The SLC46A1 (folic acid transporter) was also in the pregnant bovine endometrium. These results suggest that transporter molecules transport nutrients from blood circulation to the endometrial cells for their growth and development and then transported to the uterine lumen to nourish the embryo. Magnesium Transporter MRS2 (MRS2) is one of the highly up-regulated genes in the pregnant endometrium. MRS2 is a protein-coding gene located in mitochondria [41]. The pathway analysis from our result showed that MRS2 plays a significant role in the cell cycle, cellular assembly, and organization while having a significant role in DNA replication, recombination, and repair. According to GO analysis, MRS2 is associated with magnesium ion transportation. For the first time, this study reported MRS2 in the bovine endometrium during the maternal recognition of pregnancy. However, the spatiotemporal expression of MRS2 is completely unknown in the bovine endometrium.

Immunity and inflammation signaling pathways

Immunity and inflammation play a vital role in the endometrium during implantation, placentation, and fetal development [42]. The developing embryo modulates the maternal immune system and promotes maternal tolerance to the embryo during the peri-implantation period [43]. T helper 17 (Th17) cells produce IL-17 which recruits neutrophils via granulocyte colony-stimulating factor and IL-8 [44]. Th17 cells are known to regulate the inflammatory process in the endometrium required for the establishment and maintenance of pregnancy [45]. Non-pregnant women with unexplained recurrent pregnancy loss had a higher proportion of Th17 cell levels in the circulating blood than parous controls [45]. FOS is among the highly up-regulated gene in pregnant bovine endometrium. Ingenuity canonical pathway analysis from this study revealed that FOS impacts IL-12 signaling, and macrophage production has an important immune function and regulates innate immunity. It acts on the cAMP signaling pathway. Furthermore, the KEGG pathway shows that FOS takes part in endocrine activities such as parathyroid hormone synthesis, secretion, and action, and it has a significant role in Th17 cell differentiation. GO annotations related to this gene include DNA-binding transcription factor activity. The altered expression and distribution of FOS protein prompted endometriosis in the baboon [46].

Early-pregnancy associated novel genes in the bovine endometrium

Very Low-Density Lipoprotein Receptor (VLDLR) is highly up-regulated in the pregnant endometrium. VLDLR belongs to the low-density-lipoprotein (LDL) transmembrane receptor family that localizes to the plasma membrane and is located at chromosome 9 [47]. VLDLR consists of cell surface proteins involved in receptor-mediated endocytosis of specific ligands. This gene encodes a lipoprotein receptor, a member of the LDLR family, and plays a vital role in VLDL-triglyceride metabolism. VLDLR is considered as a potential mediator of P4-dependent signaling through membrane progesterone receptors (mPR) to drive oocyte maturation and meiosis progression. It was found that the knocking down of VLDLR inhibited the oocyte maturation and meiosis progression. In contrast, overexpression of the VLDLR showed the exact opposite action of what it did when it was knocked down, confirming its importance on P4-dependent oocyte maturation [48]. VLDLR is known to permit cholesterol to reach tissues from the bloodstream, and it may be used as an energy source. Endometrial cells secrete large amounts of cytokines, growth factors, and other molecules for embryonic growth and development during pregnancy. Therefore, VLDLR might play an important role in endometrial cell function for the establishment of pregnancy.

C15H11ORF34, also known as Placenta Expressed Transcript 1 (PLET1), is highly up-regulated in the pregnant endometrium. RNA-Seq data identified unannotated genes such as C15H11ORF34 up-regulated in embryos derived from T-cells [49]. This gene was among the top 10 up-regulated genes in our study, and the qPCR validation result showed that it was the most highly expressed gene in the pregnant bovine endometrium. The spatiotemporal expression and function of this gene are entirely unknown. Eukaryotic Translation Initiation Factor 3 Subunit M (EIF3M) is among the up-regulated genes in the pregnant endometrium. EIF3M is found in the cytosol [23, 50]. EIF3M functions in lipid metabolism, molecular transport, and protein synthesis. Furthermore, it is also involved in cell cycling, cell morphology, and apoptosis. IPA analysis from the study shows it is associated with the regulation of eIF4 (Eukaryotic Initiation Factor-4) signaling.

TINAGL1 is among the highly up-regulated gene in the pregnant endometrium. TINAGL1 is important in antimicrobial response, cell signaling, and inflammatory response. A previous study has shown an increased expression of TINAGL1 on Day-13 in pregnant heifers [26]

Proenkephalin (PENK) is a highly up-regulated gene in the pregnant endometrium. PENK is a member of the opioid polypeptide hormone found in various mammals, rodents, and avian species [51, 52]. It is known to play a role in many physiologic functions, including pain perception and stress responses. The previous study showed an increased PENK expression on Day 13 in heifers’ endometrium [26]. Besides its dominance in the CNS, it is also expressed in the oviducts in chicken and is associated with eggshell calcification [51, 52]. However, its mechanistic role in the bovine uterus has not been established. PENK was detected in the myometrial region of the pregnant mouse's uterus until Day 18 of pregnancy and helped in maternal adaptation to pregnancy and supporting embryo growth. PENK detected in the uterus was suggested to have multiple material adaptation roles to pregnancy and support embryo growth [53].

PRSS22 is a highly up-regulated gene in the pregnant endometrium. PRSS22 is located on chromosome 16 and has a predicted function related to serine-type endopeptidase activity. PRSS22 has been reported in the endometrial region of mice and humans [54] but has not previously been reported in cattle. MS4A8 is highly up-regulated in the pregnant bovine endometrium. It has membranous and cytoplasmic gene expression, especially in Fallopian tubes and respiratory epithelium in humans [55]. R3HDM1 is among the up-regulated gene in the pregnant endometrium. According to Entrez Gene, R3HDM1 maps on chromosome 2 at 2q21.3. GO annotations related to this gene include nucleic acid-binding.

In conclusion, our study showed a significant difference in gene expression between pregnant versus cyclic or non-pregnant cows. The presence of an embryo induces the endometrial transcripts related to endometrial remodeling, immune response, nutrient, ion transporters, and relevant signaling pathways in the caruncular region of bovine endometrial tissue. Further, in the absence of the embryo, these transcripts are downregulated in the endometrium. In this study, using RNA-sequencing, we found some novel genes (MRS2, C15H11ORF34, VLDLR, and PRSS22). This study provides a comprehensive dataset of transcript changes associated with maternal recognition of pregnancy, which can further be linked with the specific functions of the identified genes for future experimentation.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. RNA-Sequencing data were submitted to the NCBI GEO repository (GSE196789) https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE196789.

References

  1. Steinfeld H. The livestock revolution–a global veterinary mission. Vet Parasitol. 2004;125(1–2):19–41.

    Article  PubMed  Google Scholar 

  2. Dangour AD, Green R, Hasler B, Rushton J, Shankar B, Waage J. Linking agriculture and health in low- and middle-income countries: an interdisciplinary research agenda. Proc Nutr Soc. 2012;71(2):222–8.

    Article  PubMed  Google Scholar 

  3. Turk J. Meeting projected food demands by 2050: understanding and enhancing the role of grazing ruminants. J Anim Sci. 2016;94(suppl_6):53–62.

    Article  CAS  Google Scholar 

  4. Davoodi S, Cooke RF, Fernandes AC, Cappellozza BI, Vasconcelos JL, Cerri RL. Expression of estrus modifies the gene expression profile in reproductive tissues on Day 19 of gestation in beef cows. Theriogenology. 2016;85(4):645–55.

    Article  CAS  PubMed  Google Scholar 

  5. Bazer FW, Burghardt RC, Johnson GA, Spencer TE, Wu G. Mechanisms for the establishment and maintenance of pregnancy: synergies from scientific collaborations. Biol Reprod. 2018;99(1):225–41.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Bazer FW, Song G, Thatcher WW. Roles of conceptus secretory proteins in establishment and maintenance of pregnancy in ruminants. Asian-Australas J Anim Sci. 2012;25(1):1–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Mansouri-Attia N, Aubert J, Reinaud P, Giraud-Delville C, Taghouti G, Galio L, Everts RE, Degrelle S, Richard C, Hue I, et al. Gene expression profiles of bovine caruncular and intercaruncular endometrium at implantation. Physiol Genomics. 2009;39(1):14–27.

    Article  CAS  PubMed  Google Scholar 

  8. Imakawa K, Anthony RV, Kazemi M, Marotti KR, Polites HG, Roberts RM. Interferon-like sequence of ovine trophoblast protein secreted by embryonic trophectoderm. Nature. 1987;330(6146):377–9.

    Article  CAS  PubMed  Google Scholar 

  9. Moraes JGN, Behura SK, Geary TW, Hansen PJ, Neibergs HL, Spencer TE. Uterine influences on conceptus development in fertility-classified animals. Proc Natl Acad Sci. 2018;115(8):E1749–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Klein C, Bauersachs S, Ulbrich SE, Einspanier R, Meyer HH, Schmidt SE, Reichenbach HD, Vermehren M, Sinowatz F, Blum H, et al. Monozygotic twin model reveals novel embryo-induced transcriptome changes of bovine endometrium in the preattachment period. Biol Reprod. 2006;74(2):253–64.

    Article  CAS  PubMed  Google Scholar 

  11. Walker CG, Littlejohn MD, Meier S, Roche JR, Mitchell MD. DNA methylation is correlated with gene expression during early pregnancy in Bos taurus. Physiol Genomics. 2013;45(7):276–86.

    Article  CAS  PubMed  Google Scholar 

  12. Forde N, Beltman ME, Lonergan P, Diskin M, Roche JF, Crowe MA. Oestrous cycles in Bos taurus cattle. Anim Reprod Sci. 2011;124(3–4):163–9.

    Article  CAS  PubMed  Google Scholar 

  13. Bazer FW, Spencer TE, Ott TL. Interferon tau: a novel pregnancy recognition signal. Am J Reprod Immunol. 1997;37(6):412–20.

    Article  CAS  PubMed  Google Scholar 

  14. Zhao G, Jiang K, Zhang T, Wu H, Qiu C, Deng G. Specific interferon tau gene-regulation networks in bovine endometrial luminal epithelial cells. Theriogenology. 2018;105:51–60.

    Article  CAS  PubMed  Google Scholar 

  15. Dunne LD, Diskin MG, Sreenan JM. Embryo and foetal loss in beef heifers between day 14 of gestation and full term. Anim Reprod Sci. 2000;58(1–2):39–44.

    Article  CAS  PubMed  Google Scholar 

  16. Sandra O, Mansouri-Attia N, Lea RG. Novel aspects of endometrial function: a biological sensor of embryo quality and driver of pregnancy success. Reprod Fertil Dev. 2011;24(1):68–79.

    Article  PubMed  Google Scholar 

  17. Biase FH, Hue I, Dickinson SE, Jaffrezic F, Laloe D, Lewin HA, Sandra O. Fine-tuned adaptation of embryo–endometrium pairs at implantation revealed by transcriptome analyses in Bos taurus. Plos Biol. 2019;17(4): e3000046.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Roberts RM, Chen Y, Ezashi T, Walker AM. Interferons and the maternal-conceptus dialog in mammals. Semin Cell Dev Biol. 2008;19(2):170–7.

    Article  CAS  PubMed  Google Scholar 

  19. Ag Guide. Guide for the Care and Use of Agricultural Animals in Research and Teaching, 4th ed.; American Dairy Science Association, American Society of Animal Science and Poultry Science Association: Champaign, IL, USA, 2020; 227, Available online: https://www.asas.org/services/ag-guide.

  20. Mishra B, Kizaki K, Koshi K, Ushizawa K, Takahashi T, Hosoe M, Sato T, Ito A, Hashizume K. Expression of extracellular matrix metalloproteinase inducer (EMMPRIN) and its related extracellular matrix degrading enzymes in the endometrium during estrous cycle and early gestation in cattle. Reprod Biol Endocrinol. 2010;8(1):60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Sah N, Kuehu DL, Khadka VS, Deng Y, Peplowska K, Jha R, Mishra B. RNA sequencing-based analysis of the laying hen uterus revealed the novel genes and biological pathways involved in the eggshell biomineralization. Sci Rep. 2018;8(1):16853.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein TI, Nudel R, Lieder I, Mazor Y, et al. The genecards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 2016;54(1):1 30 31-31 30 33.

    Article  Google Scholar 

  24. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

  26. Forde N, Mehta JP, McGettigan PA, Mamo S, Bazer FW, Spencer TE, Lonergan P. Alterations in expression of endometrial genes coding for proteins secreted into the uterine lumen during conceptus elongation in cattle. BMC Genomics. 2013;14(1):321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Kerscher O, Felberbaum R, Hochstrasser M. Modification of proteins by ubiquitin and ubiquitin-like proteins. Annu Rev Cell Dev Biol. 2006;22(1):159–80.

    Article  CAS  PubMed  Google Scholar 

  28. Ott TL, Yin J, Wiley AA, Kim HT, Gerami-Naini B, Spencer TE, Bartol FF, Burghardt RC, Bazer FW. Effects of the estrous cycle and early pregnancy on uterine expression of Mx protein in sheep (Ovis aries). Biol Reprod. 1998;59(4):784–94.

    Article  CAS  PubMed  Google Scholar 

  29. Shirozu T, Sasaki K, Kawahara M, Yanagawa Y, Nagano M, Yamauchi N, Takahashi M. Expression dynamics of bovine MX genes in the endometrium and placenta during early to mid pregnancy. J Reprod Dev. 2016;62(1):29–35.

    Article  CAS  PubMed  Google Scholar 

  30. Mishra B, Kizaki K, Sato T, Ito A, Hashizume K. The role of extracellular matrix metalloproteinase inducer (EMMPRIN) in the regulation of bovine endometrial cell functions. Biol Reprod. 2012;87(6):149.

    Article  PubMed  CAS  Google Scholar 

  31. Mishra B, Kizaki K, Koshi K, Ushizawa K, Takahashi T, Hosoe M, Sato T, Ito A, Hashizume K. Expression of extracellular matrix metalloproteinase inducer (EMMPRIN) and its expected roles in the bovine endometrium during gestation. Domest Anim Endocrinol. 2012;42(2):63–73.

    Article  CAS  PubMed  Google Scholar 

  32. Mishra B, Koshi K, Kizaki K, Ushizawa K, Takahashi T, Hosoe M, Sato T, Ito A, Hashizume K. Expression of ADAMTS1 mRNA in bovine endometrium and placenta during gestation. Domest Anim Endocrinol. 2013;45(1):43–8.

    Article  CAS  PubMed  Google Scholar 

  33. Di Lullo GA, Sweeney SM, Körkkö J, Ala-Kokko L, San Antonio JD. Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen. J Biol Chem. 2002;277(6):4223–31.

    Article  PubMed  CAS  Google Scholar 

  34. Salamonsen LA. Role of proteases in implantation. Rev Reprod. 1999;4(1):11–22.

    Article  CAS  PubMed  Google Scholar 

  35. Altemeier WA, Schlesinger SY, Buell CA, Parks WC, Chen P. Syndecan-1 controls cell migration by activating Rap1 to regulate focal adhesion disassembly. J Cell Sci. 2012;125(Pt 21):5188–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Kramer F, Dernedde J, Mezheyeuski A, Tauber R, Micke P, Kappert K. Platelet-derived growth factor receptor beta activation and regulation in murine myelofibrosis. Haematologica. 2020;105(8):2083–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Pobezinskaya YL, Kim YS, Choksi S, Morgan MJ, Li T, Liu C, Liu Z. The function of TRADD in signaling through tumor necrosis factor receptor 1 and TRIF-dependent Toll-like receptors. Nat Immunol. 2008;9(9):1047–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yockey LJ, Iwasaki A. Interferons and proinflammatory cytokines in pregnancy and fetal development. Immunity. 2018;49(3):397–412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Cedernaes J, Olszewski PK, Almen MS, Stephansson O, Levine AS, Fredriksson R, Nylander O, Schioth HB. Comprehensive analysis of localization of 78 solute carrier genes throughout the subsections of the rat gastrointestinal tract. Biochem Biophys Res Commun. 2011;411(4):702–7.

    Article  CAS  PubMed  Google Scholar 

  40. He L, Vasiliou K, Nebert DW. Analysis and update of the human solute carrier (SLC) gene superfamily. Hum Genomics. 2009;3(2):195–206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Li Y, Sun J, Ling Y, Ming H, Chen Z, Fang F, Liu Y, Cao H, Ding J, Cao Z, et al. Transcription profiles of oocytes during maturation and embryos during preimplantation development in vivo in the goat. Reprod Fertil Dev. 2020;32(7):714–25.

    Article  CAS  PubMed  Google Scholar 

  42. Mor G, Cardenas I, Abrahams V, Guller S. Inflammation and pregnancy: the role of the immune system at the implantation site. Ann N Y Acad Sci. 2011;1221(1):80–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Maeda Y, Ohtsuka H, Tomioka M, Oikawa M. Effect of progesterone on Th1/Th2/Th17 and regulatory T cell-related genes in peripheral blood mononuclear cells during pregnancy in cows. Vet Res Commun. 2013;37(1):43–9.

    Article  PubMed  Google Scholar 

  44. McKenzie BS, Kastelein RA, Cua DJ. Understanding the IL-23-IL-17 immune pathway. Trends Immunol. 2006;27(1):17–23.

    Article  CAS  PubMed  Google Scholar 

  45. Lee SK, Kim JY, Lee M, Gilman-Sachs A, Kwak-Kim J. Th17 and Regulatory T cells in Women with Recurrent Pregnancy Loss. Am J Reprod Immunol. 2012;67(4):311–8.

    Article  CAS  PubMed  Google Scholar 

  46. Hastings JM, Jackson KS, Mavrogianis PA, Fazleabas AT. The estrogen early response gene fos is altered in a baboon model of endometriosis1. Biol Reprod. 2006;75(2):176–82.

    Article  CAS  PubMed  Google Scholar 

  47. Go GW, Mani A. Low-density lipoprotein receptor (LDLR) family orchestrates cholesterol homeostasis. Yale J Biol Med. 2012;85(1):19–28.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Nader N, Dib M, Courjaret R, Hodeify R, Machaca R, Graumann J, Machaca K. The VLDL receptor regulates membrane progesterone receptor trafficking and non-genomic signaling. J Cell Sci. 2018;131(10):jcs212522.

  49. Srirattana K, St John JC. Manipulating the Mitochondrial Genome To Enhance Cattle Embryo Development. G3 (Bethesda). 2017;7(7):2065–80.

    Article  CAS  Google Scholar 

  50. Proshkin SA, Shematorova EK, Souslova EA, Proshkina GM, Shpakovski GV. A minor isoform of the human RNA polymerase II subunit hRPB11 (POLR2J) interacts with several components of the translation initiation factor eIF3. Biochem Mosc. 2011;76(8):976.

    Article  CAS  Google Scholar 

  51. Jeong W, Kim J, Ahn SE, Lee SI, Bazer FW, Han JY, Song G. AHCYL1 is mediated by estrogen-induced ERK1/2 MAPK cell signaling and microRNA regulation to effect functional aspects of the avian oviduct. Plos One. 2012;7(11): e49204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Brionne A, Nys Y, Hennequet-Antier C, Gautron J. Hen uterine gene expression profiling during eggshell formation reveals putative proteins involved in the supply of minerals or in the shell mineralization process. BMC Genomics. 2014;15(1):220.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Zhu Y, Pintar JE. Expression of opioid receptors and ligands in pregnant mouse uterus and placenta. Biol Reprod. 1998;59(4):925–32.

    Article  CAS  PubMed  Google Scholar 

  54. Pelch KE, Schroder AL, Kimball PA, Sharpe-Timms KL, Davis JW, Nagel SC. Aberrant gene expression profile in a mouse model of endometriosis mirrors that observed in women. Fertil Steril. 2010;93(5):1615-1627 e1618.

    Article  CAS  PubMed  Google Scholar 

  55. Huang T, Huang X, Shi B, Liang X, Luo J, Yao M. Relationship among MS4A8 expression, its variants, and the immune response in a porcine model of Salmonella. Can J Anim Sci. 2018;98(4):778–86.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank Marla Fergerstrom of CTAHR Mealani Research Station for taking care of heifers, and Hawaii Beef Producers LLC, Big Island, Hawaii, for their help in sampling.

Funding

This work was supported by the Hawaii Department of Agriculture (HDOA #38993) to B.M.

Author information

Authors and Affiliations

  1. Department of Human Nutrition, Food and Animal Sciences, University of Hawaii at Manoa, Honolulu, HI, 96822, USA

    Bindu Adhikari, Chin N. Lee, Glen Fukumoto, Mark Thorne, Kyle Caires, Jenee Odani & Birendra Mishra

  2. Department of Quantitative Health Sciences, University of Hawaii at Manoa, Honolulu, HI, 96822, USA

    Vedbar S. Khadka & Youping Deng

Authors
  1. Bindu Adhikari
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chin N. Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Vedbar S. Khadka
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Youping Deng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Glen Fukumoto
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Mark Thorne
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Kyle Caires
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jenee Odani
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Birendra Mishra
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

B.A. collected samples, conducted the experiment, analyzed the data, and drafted the manuscript. C.N.L., G.F., M.T., K.C., and J.O. helped in animal preparation, sampling, and manuscript preparation. V.S.K. and Y.D. helped in RNA-seq data analysis. B.M. designed the study, assisted in sampling, data analysis, and drafted the manuscript. All authors read and approved the manuscript for publication.

Corresponding author

Correspondence to Birendra Mishra.

Ethics declarations

Ethics approval and consent to participate

Ethics approval was obtained from the Institutional Animal Care and Use Committee (IACUC) of the University of Hawaii at Manoa (Approval no. 18–3008). The study was conducted following IACUC. All methods involving animals were performed in accordance with the relevant guidelines and regulations of the University of Hawaii at Manoa. The study was in accordance with ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table S1.

 Primers used to quantify the expression of target genes by qPCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Adhikari, B., Lee, C.N., Khadka, V.S. et al. RNA-Sequencing based analysis of bovine endometrium during the maternal recognition of pregnancy. BMC Genomics 23, 494 (2022). https://doi.org/10.1186/s12864-022-08720-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08720-4

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us