- Research article
- Open Access
Global gene expression in endometrium of high and low fertility heifers during the mid-luteal phase of the estrous cycle
BMC Genomics volume 15, Article number: 234 (2014)
In both beef and dairy cattle, the majority of early embryo loss occurs within the first 14 days following insemination. During this time-period, embryos are completely dependent on their maternal uterine environment for development, growth and ultimately survival, therefore an optimum uterine environment is critical to their survival. The objective of this study was to investigate whether differences in endometrial gene expression during the mid-luteal phase of the estrous cycle exist between crossbred beef heifers ranked as either high (HF) or low fertility (LF) (following four rounds of artificial insemination (AI)) using the Affymetrix® 23 K Bovine Gene Chip.
Conception rates for each of the four rounds of AI were within a normal range: 70–73.3%. Microarray analysis of endometrial tissue collected on day 7 of the estrous cycle detected 419 differentially expressed genes (DEG) between HF (n = 6) and LF (n = 6) animals. The main gene pathways affected were, cellular growth and proliferation, angiogenesis, lipid metabolism, cellular and tissue morphology and development, inflammation and metabolic exchange. DEG included, FST, SLC45A2, MMP19, FADS1 and GALNT6.
This study highlights, some of the molecular mechanisms potentially controlling uterine endometrial function during the mid-luteal phase of the estrous cycle, which may contribute to uterine endometrial mediated impaired fertility in cattle. Differentially expressed genes are potential candidate genes for the identification of genetic variation influencing cow fertility, which may be incorporated into future breeding programmes.
The failure of breeding females to become pregnant, in both dairy and beef cattle production systems, directly impacts the economic viability of these enterprises, and ultimately hinders genetic progress. Significant decreases in dairy cow fertility, ranging from 0.45% to 1% per annum, have been reported in cattle populations across the globe [1–3]. Following insemination the greatest increment of cow reproductive wastage occurs in the form of early embryo mortality with approximately 80% of this occurring within 14–16 days [4–6]. More specifically, previous studies have highlighted that the majority of early embryo loss typically commences around the mid-luteal phase of an estrous cycle i.e. day 7 of pregnancy [7, 8] concurrent with the critical blastulation stage of embryo development .
There is evidence of repeatable differences between cows in their ability to become pregnant. McMillan  reported a 65% difference in pregnancy rate at 60 days of gestation, following 6 consecutive in vitro embryo transfer events, between two groups of cows. Differences in follicle wave dynamics, duration of estrus, site of ovulation, or subsequent progesterone profiles were not found to contribute to the observed difference in pregnancy rate. Indeed, the authors suggested that “uterine” rather than “ovarian” factors may be responsible for the variation observed. This uterine effect was also hypothesized in similar studies examining phenotypic differences between high and low fertility animals [11–13]. Furthermore, data from our laboratory suggest a repeatability estimate of 0.18 for embryo survival in beef heifers  and heritability estimates for conception rate have been reported to exceed 0.20 [14, 15].
The prerequisites to the establishment and maintenance of a successful pregnancy include a viable embryo, an appropriate steroidal environment and an optimally functioning and receptive endometrium [16–18]. The endometrium plays a pivotal role in orchestrating the events that lead to fertilization, implantation and pregnancy. Throughout the estrous cycle and pregnancy, the endometrium is subjected to a host of functional and morphological changes, regulated by the hormones progesterone, estradiol and oxytocin . The endometrium also functions to secrete a multitude of growth factors, proteins and cytokines, all of which constitute the histotroph, an important source of energy and nutrition to a growing embryo in vivo[20–22].
Using conventional candidate approaches many studies have examined bovine endometrial gene expression under various conditions; during early pregnancy in animals that produced viable and non-viable embryos , in pregnant and cycling animals with artificially induced high, and normal systemic progesterone concentrations [24–26] and during the various phases of the estrous cycle . Furthermore global endometrial gene expression analyses have been conducted and include comparisons between cycling and pregnant animals [28, 29], fertile and sub-fertile animal strains [30, 31], progesterone supplementation treatments , and specific estrous cycle phases [33, 34]. Despite these efforts, endometrial gene expression of animals characterized as either high or low fertility has not been investigated. Given the critical importance of day 7 [7, 8], we hypothesise that uterine endometrial gene expression patterns will be different between high and low fertility heifers on day 7 of the estrous cycle. Thus, the objective of this study was to characterize differential gene expression profiles in endometrial tissue harvested on day 7 of the estrous cycle from heifers ranked as either HF or LF fertility based on four successive inseminations and pregnancy diagnoses. Intercaruncular endometrial tissue was examined due to the fact that caruncular endometrium lacks uterine glands which are essential to the exchange, transport and secretion of pertinent metabolites which constitute the uterine histoptroph and are required to support pregnancy [35, 36].
All experimental procedures involving heifers were licensed by the Department of Health and Children, Ireland (licence number B100/846). Protocols were in accordance with the Cruelty to Animals Act (Ireland 1876, as amended by European Communities regulations 2002 and 2005) and the European Community Directive 86/609/EC and were sanctioned by the Institutional Animal Research Ethics Committee.
Estrous cycles of reproductively normal nulliparous crossbred beef heifers (Bos taurus n = 120) were synchronized using two intramuscular administrations of 500 μg of the prostaglandin F2α analogue (PG), cloprostenol (Estrumate®, Schering-Plough Ltd., Shire Park, Welwyn Garden City, Hertfordshire, UK). Animals were visually observed for signs of estrous activity 3- to 5- times daily as described by Lynch et al. . Only heifers observed to be in standing estrus were inseminated 6–18 hrs after onset of heat . Inseminations were carried out artificially by one trained technician. Heifers were given a single insemination of frozen-thawed semen, collected from a single ejaculate of one high fertility bull. Sire breed was Limousin and named Bolide (FL17). At the time of the 1st insemination, heifers were on average 20 months of age and weighed 440 ± 9.0 kg (Mean ± SEM).
Using an Aloka SSD-500 V ultrasound scanner, fitted with a 7.5 MHZ transducer (Aloka Co. Ltd., Tokyo, Japan), pregnancy was diagnosed 28 days after insemination using the criteria set out by Kastelic et al. . Following diagnosis, all pregnant heifers received PG on day 28 to induce embryo loss. Six weeks after induced embryo loss all heifers were subjected to estrous reprogramming using a two-injection PG-regimen (11 days apart), inseminated and pregnancy scanned as described above.
For the purpose of establishing an accurate high versus low heifer fertility model, this schedule was followed for a further two occasions. Thus, following four inseminations, animals that established a pregnancy on all four occasions were categorized as “HF” heifers while those achieving pregnancy on only one occasion were categorized as “LF” heifers. To eliminate the possibility of a physical or anatomical abnormality that may have impeded heifers from becoming pregnant, animals with zero recorded pregnancies were omitted from the study.
After the fourth insemination, and subsequent pregnancy diagnosis, pregnant heifers were returned to estrous. Approximately three months later, estrous cycles of animals were synchronized again in preparation for endometrial harvesting on D7. Figure 1 illustrates the timeline of events during the experimental period.
Throughout the experimental period, animals had ad libitum access to grass silage supplemented with 2 kg of concentrates per heifer per day. Heifers were housed on concrete slats in groups of 15, at 2.5 m2 per heifer, for the duration of the study (15 months). Slaughter liveweight averaged 625 kg, with BCS of 4.0. Heifers were gaining weight during the course of inseminations in the region of 0.60 kg/day.
Animals from HF (n = 6) and LF groups (n = 6) were slaughtered on D7 in a licensed abattoir (KEPAK, Athleague, Co. Roscommon, Ireland). Following slaughter the reproductive tract and ovaries were checked for gross abnormalities but none were recorded. Uterine tissues were opened longitudinally along the mesenteric border. Intercaruncular endometrial cross-sections approximately 4 sq cm, and weighing 2.5 g, were harvested from the animals by peeling from the underlying uterine myometrium from the middle-third of the uterine horn ipsilateral to the corpus luteum (CL) within 20 min of slaughter.
Samples were washed in sterile PBS, and stored in RNAlater® at 4°C for 24 h before being transferred for long-term storage at −20°C. All surgical instruments used for tissue collection were sterilized and treated with RNA Zap (Ambion, Applera Ireland, Dublin, Ireland). In addition, on the day of slaughter CL diameter for each heifer was determined using vernier calipers.
Heifers were blood sampled via jugular venipuncture for subsequent measurement of progesterone at 0900 and 2100 h commencing 24 h after PG for a cycle length. All blood samples were collected into 10 ml ethylenediamine tetraacetic acid (EDTA) heparinized Vacutainers (Becton Dickson Vacutainer Systems, Plymouth, UK). Samples were held in iced water until centrifuged at 1500 × g at 4°C for 15 mins after which plasma was extracted and stored in sterile 7 ml vials at −20°C until assayed.
Progesterone profiles for each of the six heifers within HF and LF groups were established. Concentration of progesterone was measured in plasma as the mean of the two samples taken on each cycle day of the previous cycle and on 7 days prior to slaughter using the Coat-a-Count assay procedure (Coat-a-Count Diagnostic Products Corporation, Los Angeles, CA, USA) with each sample tested in duplicate. The inter-assay and intra-assay coefficients of variation for low, medium and high control samples were 17.4% and 4.4%, 5.6% and 28.4%, and 4.2% and 4.9% with mean concentrations of 0.24, 2.54 and 7.21 ng/mL, respectively. The minimum detectable limit for this assay was 0.06 ng/mL.
RNA extraction and quality analysis
Total RNA was prepared from 100–200 mg of endometrial tissue using the TRIzol reagent (Sigma-Aldrich Ireland Ltd., Dublin, Ireland). Tissue samples were homogenized in 3 ml of TRIzol reagent and chloroform, and subsequently precipitated using isopropanol (Sigma-Aldrich Ireland Ltd., Dublin, Ireland). RNA samples were stored at −80°C. Samples of RNA, (20 μg), were purified and treated for contaminating genomic DNA using RNeasy clean-up kits in accordance with manufacturer’s guidelines supplied (QIAGEN, Crawley, West Sussex, UK). This protocol included an on-column DNase treatment step. RNA quality and quantity were assessed using automated capillary gel electrophoresis on a Bioanalyzer 2100 with RNA 6000 Nano Lab-chips according to manufacturer’s instructions (Agilent Technologies Ireland, Dublin, Ireland). Absorbance ratios (28S/18S) and RNA integrity values recorded for all RNA samples extracted post clean-up ranged between 1.8 and 2.0, and 7.5 and 9.8, respectively.
Gene expression was determined using a 24,027 probe set bovine oligonucleotide array (Affymetrix®), representing ~23,000 bovine transcripts based on the original mapping using Unigene build 57 (March 24, 2004). RNA from each heifer was hybridized to a separate array. All 12 RNA samples were hybridized and scanned by the German Resource Centre for Genomic Research (RZPD), Germany, according to the manufacturer’s instructions.
All microarray analyses including preprocessing, normalization and statistical analysis were carried out using R (R, 2007) version 2.6 and Bioconductor  version 2.1 as previously described by . Data were quality assessed before and after normalization using a number of in-built quality control methods implemented in the Bioconductor affycoretools and associated packages to identify problems if they existed with array hybridization, RNA degradation and data normalization. Microarray data were preprocessed using the mmgMOS normalization method [41, 42] using the default settings and differential expression (DE) was calculated using the pumaDE method both implemented in the Bioconductor package “puma” [42–45]. The puma method uses a Bayesian hierarchical model to calculate the probability of positive likelihood ratio (PPLR). The PPLR associates probability values of genes being differentially expressed, which is a measure of false positive detection of DE, to each ratio and generates lists of genes ranked by the probability of DE. This PPLR statistic was converted into “P-like values” using the recommended formula in the puma method prior to subsequent analysis.
As many of the original annotations for the Affymetrix bovine chip are erroneous [6, 46], remapped annotations were determined using the “bovinedaiplusv6cdf” chip definition file (CDF). This annotation is based on the CDF-Merger procedure as described by De Leeuw et al. , which generates a hybrid CDF based on the standard Affymetrix CDF (version 26) and the custom Brainarray (version 11.0.1) CDF. This re-mapped annotation includes mapping to all RefSeq (mature RNA protein coding transcripts and validated complete coding sequences in GenBank). Annotations were also supplemented by interrogating the Ensembl Bos taurus database version 46 using the BioMart package in Bioconductor and manual annotation where possible with recent entries in Entrez Gene.
To examine the molecular functions and genetic networks, the microarray data were further analyzed using Ingenuity Pathway Analysis (v. 8.8, Ingenuity Systems, Mountain View, CA; http://www.ingenuity.com), a web-based software application that enables identification of over-represented biological mechanisms, pathways, and functions most relevant to experimental datasets or genes of interest [40, 48–50].
A dataset containing gene identifiers and corresponding expression and P-like values was uploaded into IPA. Briefly, each identifier was mapped to its corresponding gene object in the Ingenuity knowledge base. A P-like value of P < 0.05 from the puma analysis was set to identify genes whose expression was statistically significantly up- or down-regulated. These genes, called “focus” genes, were overlaid onto a global molecular network developed from information contained within the Ingenuity knowledge base. Networks of these focus genes were then algorithmically generated based on their connectivity. Network analysis returns a score that ranks networks according to their degree of relevance to the network eligible molecules in the dataset. The score takes into account the number of network eligible molecules in the network and its size, as well as the total number of network eligible molecules analyzed and the total number of molecules in the knowledge base that could potentially be included in networks.
The microarray results were validated by carrying out RT-qPCR on 18 genes. Candidate genes were chosen based on the following criteria; those that were top ranking in our microarray DEG list, genes with known functional importance in uterine mediated sub-fertility which were either up- or down-regulated and genes which were not differentially expressed between the two treatment groups.
Using the same RNA samples that were analyzed in the microarray studies, first strand cDNA was synthesized using the High Capacity cDNA Reverse Transcription kit according to manufacturer’s instructions (Applied Biosciences, Ireland). Purified total RNA (1 μg) was reverse transcribed using random hexamers. The converted cDNA was quantified by absorbance at 260 nm, diluted to 50 ng/μl working stocks and stored at −20°C, for subsequent analyses.
Analysis of putative reference genes for RT-qPCR studies was carried out using GeNorm version 3.5 Microsoft Excel Add in (Microsoft, Redmond, WA) . The stability of the expression of several cited reference genes including, ribosomal protein L15 , 18 s ribosomal RNA , ubiquitin , glyceraldehyde phosphate dehydrogenase and β-actin [55, 56], was investigated across all samples in this study. Similar to Coyne et al. , ubiquitin (at an optimal concentration of 2.5 μM) exhibited the greatest stability during qPCR analysis of endometrial mRNA samples analyzed, with an M value of 0.022. Based on a recommended cut-off V value of 0.15; ubiquitin was selected as a single standard reference gene for these experiments as the use of additional reference genes did not contribute to a more accurate normalization factor.
Primers were designed, to span exon boundaries where possible, using the Primer3 software programme  and oligos were aligned by Basic Local Alignment Search Tool (BLASTN) on the National Centre for Biotechnology Information (NCBI) web page, to verify their identity and homology to the bovine genome (http://www.ncbi.nlm.nih.gov/BLAST/). All oligonucleotides were commercially synthesized as highly purified salt-free products by Sigma Aldrich Ireland Ltd. Primers were first tested using end point PCR to optimize amplification conditions. All amplified PCR products generated in this study were purified using the PCR purification kit (Roche, Basel, Switzerland) and sequenced (Macrogen; Nucleics Pty Ltd, Bendigo, Australia) to verify their identity. Primer sequences used in this study are listed in Table 1.
Primer concentrations were optimized for each gene by titrating 5, 10, and 20 μM per primer. The most suitable primer concentration was chosen based on four criteria in order of decreasing importance: i) a clear distinct melt curve absent of any additional peak(s) caused by non-specific binding, ii) a curve within the temperature range 75–85°C, iii) the primer concentration producing the lowest threshold cycle number (Ct) and lastly, iv) replication amongst Ct values and melting temperatures (Tm). Subsequently, efficiencies of chosen primer concentrations were determined over a 5-fold dilution series, whereby cDNA was diluted into working solutions: stock, 1:2, 1:4, 1:8, 1:16, and RT-qPCR assays carried out. This was repeated for every gene. The r2 and amplification efficiency (E) values for RT-qPCR were calculated from linear regression analysis of log (input cDNA) versus Ct plot. The slope for each set of standards was used to determine E = 10(−1/slope) – 1. Slopes, amplification efficiencies and R2 estimates for individual genes are reported in Table 2. Only primers with PCR efficiencies between 90% and 110% were used.
Each RT-qPCR reaction was carried out in a 96-well plate format with a total volume of 20 μl, containing 1 μl cDNA, (10 ng/μl), 10 μl Fast SYBR® Green Master Mix (Applied Biosystems, Ireland), 1 μl forward and reverse primers and 8 μl nuclease-free H2O. Performance of RT-qPCR was carried out using the Applied Biosystems Fast 7500 v2.0.1 with the following cycling parameters: 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 60°C for 60 s, followed by amplicon dissociation (95°C for 15 s, 60°C for 60 s, 95°C for 15 s and 60°C for 15 s). Dissociation curves were examined for the presence of a single PCR product. The software package GenEx 220.127.116.11 (MultiD Analyses AB, Gothenburg, Sweden) was used for efficiency correction of the raw cycle threshold (Ct) values, interplate calibration based on a calibrator sample included on all plates, averaging of replicates, normalization to the reference gene and the calculation of quantities relative to the greatest Ct. Expression of each target gene was normalised to the reference gene and relative differences in gene expression were calculated using the 2-ΔΔCT method .
All data were analyzed using the Statistical Analysis Systems software package (SAS Inst. Inc., Cary, NC) version 9.1. Data from RT-qPCR studies were tested for adherence to normality using PROC UNIVARIATE (SAS, 2003). Non-normal data were subsequently transformed using the best fit function as described by PROC TRANSREG (SAS, 2003). Differences in mean values between the two groups (HF and LF) were tested using ANOVA (PROC MIXED). Animal within treatment was used as the error term. The Tukey critical difference test was used to determine statistical difference between LF and HF mean values. The CORR procedure of SAS (PROC CORR, SAS 2003) was used to determine correlations between microarray and RT-qPCR data. Pearson correlation coefficients were estimated for each individual gene across all animals (n = 12). A P value of P < 0.05 was considered to be statistically significant. Data collected from CL diameter measurements were tested for adherence to normality using PROC UNIVARIATE (SAS, 2003). CL differences in mean values between the two groups (HF and LF) were tested using ANOVA (PROC MIXED). Animal within treatment was used as the error term. For the analysis of progesterone profiles individual profiles were normalized relative to day of estrus (Day 0). The effect of fertility status “HF” versus “LF” was established using a repeated measured analysis (PROC MIXED; SAS).
Embryo survival rates were 73.3%, 71.7%, 73.3% and 70.0% for A.I. rounds 1, 2, 3 and 4 respectively. A total of 31 heifers qualified as HF or LF; 15 HF and 16 LF, of which three of these were eliminated from the study due to the presence of ovarian abnormalities detected at ultrasound scanning. Pregnancy rate for LF heifers was consistent across all four replicates. Six HF and 6 LF heifers were randomly chosen within their respective fertility groups for slaughter on D7. The mean inter-estrous intervals in a previous recorded estrous cycle were 20.17 ± 0.96 and 20.83 ± 0.96 days (P > 0.10) for the HF and LF heifers, respectively. At day of slaughter, mean CL diameters were 22.58 ± 3.48 (SD) mm and 23.55 ± 4.4 (SD) mm for HF and LF heifers, respectively, i.e., there was no significant difference in CL diameter between fertility groups (P > 0.10).
There was no effect of fertility status, or interaction effect of fertility status and day of cycle (P > 0.10), on the concentration of progesterone. On the day of slaughter plasma concentrations did not differ between the high and low fertility groups (HF 5.96 ng ml−1; LF 5.65 ng ml−1, P = 0.589).
Microarray differential gene expression
A total of 419 genes were found to be differentially expressed between LF and HF (n = 6 vs. 6). Of these, 171 were up-regulated and 248 down-regulated in the LF compared with HF heifers, respectively. Transcript abundance differences between LF and HF groups resulted in fold changes ranging from 6.6-fold down to 8-fold up-regulated in LF animals. The microarray data have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE29853. Hierarchical clustering of differentailly expressed genes is presented as a heatmap and dendogram in Additional file 1: Figure S1.
Of the 419 DEG, a total of 227 genes were successfully mapped to a molecular/biological pathway and/or category in the IPA database, while 202 of these were network eligible using IPA. Among the mapped DEG, 73 were up-regulated (Additional file 2: Table S1) and 154 down regulated (Additional file 2: Table S2).
Biological categories with the largest number of up regulated genes included DNA replication, recombination and repair, nucleic acid metabolism and carbohydrate metabolism. Categories with the largest number of down-regulated genes were organ morphology, and connective tissue development and function. Of the top 20 most statistically significantly over-represented biological categories, DNA replication, recombination and repair had the greatest ratio of up- to down-regulated genes (Figure 2). Pathways with the greatest number of DEG, including their respective number of DEG, were cellular growth and proliferation (n = 57), inflammatory disease (n = 55), cell death (n = 49), cellular development (n = 43), small molecule biochemistry (n = 37), cellular morphology (n = 36) and tissue development (n = 36) as shown in Table 3.
Canonical signaling pathway analysis uncovered genes with functions in ILK-signaling, TR/RXR activation, regulation of actin based motility by Rho and Integrin signaling (Table 4). Genes associated with canonical signaling pathways were down-regulated in LF animals for all statistically significant pathways mapped with the exception of TR/RXR activation where the ratio of up- to down-regulated genes was uniform. Canonical metabolic pathways over-represented within the microarray data included fatty acid biosynthesis, o-glycan biosynthesis and purine metabolism. There were more genes up-regulated in canonical metabolic than canonical signaling pathways with the greatest ratio of up- to down-regulated genes expressed in the metabolic pathway: o-glycan biosynthesis (Table 4).
Using IPA a total of 19 gene networks were identified, 12 of which had 13 to 25 focus genes among DEG (Additional file 2: Tables S1 and S2). The 12 top networks are listed in Table 5. Lipid metabolism featured in three of the top 12 networks. In addition, organ/tissue/cell morphology and development appeared a central biological theme over-represented among DEG. Illustrations of gene interactions among DEG contained within the top two scoring networks can be seen in Figures 3 and 4. Biological pathways; lipid metabolism, cell growth and proliferation, and tissue development and function, were repeatedly featured pathways that constituted these top scoring networks.
Eighteen genes were validated by real-time RT-qPCR (Table 1). There was moderate to good consistency between methodologies for direction and magnitude of differential gene expression among genes analyzed. Correlation coefficients exceeded 0.60 in fourteen of the eighteen genes validated (Figure 5, Additional file 2: Table S3).
The animal model generated in this study, is the first of its kind. Two groups of heifers consistently divergent in conception rate; HF and LF were successfully generated and endometrial gene expression examined. We identified key genes and pathways potentially contributing to endometrial related conception rate variance, the most extreme of which had no previously known involvement in endometrial function, including cellular growth and proliferation NPPC and GJA1; angiogenesis MMP19 and HMGB1; lipid metabolism FASN and PPARA; cellular and tissue morphology and development FST and TGFB1I1; inflammation IL-33; and metabolic exchange SLC1A3 and SLC25A24.
Several studies have highlighted the vital role progesterone plays in early embryo development to the extent that decreased conception rates were observed in heifers with a delayed postovulatory progesterone peak . Furthermore, it has been well documented that progesterone influences endometrial and oviductal function [61, 62]. In the present study, progesterone concentrations were within the normal range for both HF and LF heifers, and did not vary between groups. In addition, CL diameter measurements were not different between HF and LF animals and were consistent with observations from other studies examining CL diameters during this period of the estrous cycle . The high conception rates achieved across successive breedings was indicative of reproductively healthy animals, with good heat detection and insemination technique providing confidence in retrospective fertility status. However, it is important to note, other factors potentially contributing to the conception rate differences observed between HF and LF heifers, including oocyte quality and oviductal environment, were not analysed in this study.
Endometrial function plays a critical role in pre-implantation embryo survival. Consequently, much work has focused on the biochemical and molecular phenomena surrounding the progression of an estrous cycle . The present study is novel as it provides information on gene expression during an important period of the estrous cycle: the mid-luteal phase, otherwise recognized as a critical period of embryo loss during pregnancy [5, 7] between animals of high and low reproductive capacity. Reiterating the importance of examining transcription during this phase, Salilew-Wondim et al.  recently found more extensive differential gene expression in endometrium harvested from heifers on D7 (an estrous cycle prior to embryo transfer) between heifers that conceived and those that returned to estrus before day 21, when compared with D14.
GALNT6, encoding enzyme UDP-N-acetyl-alpha-D-galactosamine: polypeptide N-acetylgalactosaminyltransferase 6, was the most abundantly expressed gene in LF heifers. It was 6.7 fold up-regulated in LF compared with HF heifers. This is the first report of expression of this gene in Bos taurus. The GALNT6 gene is located on chromosome 5 in the bovine genome and shares a coding region with SLC4A8, a sodium bicarbonate co-transporter. Expression of this gene in humans is implicated in the synthesis of oncofetal fibronectin (onfFN) , a protein found in plasma and cervicovaginal secretions; increased concentrations of which has been associated with abnormal pregnancy [65, 66]. However, Feinberg et al.  reported increased protein levels of onfFN at the trophoblast–endometrial ECM interface in human pregnancy tissues from gestational day 20 to full term in healthy pregnancies. These observations suggest that differential expression of the enzyme GALNT6 may have consequences for embryo survival and that this may be time specific however its role is currently unclear.
Pathway analysis is widely used to analyze gene expression data and serves as an effective tool for delineating the underlying biological processes involved in mRNA aberrations [68–71]. Biological pathways altered in the current study included: cellular growth and proliferation, lipid metabolism, tissue remodeling, ECM mineralization, inflammation, angiogenesis, and metabolic exchange.
Cellular growth and proliferation
Owing to its regenerative nature, the endometrium undergoes highly complex but tightly regulated cellular proliferation and differentiation throughout the estrous cycle . There is little published information on the molecular mechanism of bovine endometrial proliferation throughout the estrous cycle however, studies examining uterine tissue of non-pregnant ewes during cycle days 0 to 15 showed an increased rate of cellular proliferation between days 0 and 4, decreasing by day 15, suggesting a proliferative disposition is normal earlier in the estrous cycle . Results from our study indicate that LF animals could be experiencing an abnormal decline in cellular growth/proliferation i.e. 21 genes implicated in cellular proliferation inhibition, including FST, NPPC, GJA1, SOX6, were up-regulated in the LF animals. Of these genes FST, NPPC and GJA1 were previously found to be expressed in bovine endometrial tissue [78, 79]. Substantial inhibition of endometrial cellular proliferation would retard the development of a secretory endometrium and suppress endometrial maturation , thus making successful implantation unlikely.
A critical element of tissue growth and development is the growth of new blood vessels, also known as angiogenesis . Generally inactive in healthy individuals and animals, angiogenesis plays an active role in endometrial function, as well as growth of ovarian follicles and CL during the reproductive cycle [82, 83]. In a highly proliferating tissue such as endometrium, and particularly during the hypothesized window of proliferation day 0 to 14/15, angiogenesis is necessary for the provision of nutrients. Factors controlling angiogenesis include growth factors, nitric oxide and matrix metalloproteinases (MMPs), of which MMP19 was down-regulated in the LF animals . Also down-regulated, high-mobility group box 1 (HMGB1) which codes for a protein which has previously been identified in uterine fluid of dairy heifers on day 7 post estrus . A role for members of the HMBG family in angiogenesis is supported by their expression during mouse embryogenesis  with higher expression levels found in proliferating cells  and lower expression in fibroblasts from old-age humans . Down-regulation of these and other angiogenic genes, which was the case in LF animals, could prevent the necessary angiogenic cascades synergistic with cellular proliferation that dominate the mid-luteal phase .
Lipid metabolism appears in three of the top 5 networks, suggesting its importance as a metabolic process in uterine physiology. Genes involved included ACAT1, CCAT, LGALS1, PCCB, SRD5A1, FASN and PPARA. In particular, increased PPARA transcript abundance, as observed in LF heifers, coincides with increased fatty acid catabolism . Fatty acids are essential precursors to steroids and eicosanoids, metabolites necessary for normal ovarian and uterine function . Furthermore, studies have shown fatty acid supplementation positively influences reproductive performance [9, 54].
Fatty acid synthase (FASN) exhibits its anabolic capacity by aiding in the conversion of dietary carbohydrate to fat, which is subsequently organized into hepatic adipocytes and lactating mammary tissue as triglyceride and milk lipids, respectively [89, 90]. It has also been found that expression of FASN peaks during the proliferative phases of the menstrual cycle . Metabolic demands are particularly high during this phase as a result of the extensive endometrial remodeling and reconstruction, a central theme to both the estrous and menstrual cycles. Increased FASN would be favorable in such a demanding situation to deliver the required fatty acid for the assembly of new cell membranes, modification of DNA transcriptional machinery and hormone construction. Interestingly, expression of FASN was down-regulated in the LF heifers suggesting the aforementioned processes were compromised in these animals, which potentially affecting their ability to conceive.
Steroid 5α-reductase type 1 enzyme is involved in the metabolism of progesterone that is found in uterine and cervical cavities. Murine gene knock-out studies have shown that parturition is adversely affected by aberrant expression of this gene, impeding cervical ripening and fetal delivery as a result of elevated progesterone levels in the cervix . Expression of the gene coding for this enzyme was up-regulated in LF heifers, thus progesterone catabolism is likely to be active in these animals. As high progesterone levels are positively associated with embryo survival [60, 93], it is therefore possible that the LF animals are experiencing low local progesterone concentrations and ultimately, this could be contributing to their low conception rates.
Cellular and tissue morphology and development
The ability of cells to generate alternate cell types whose phenotype is different from that of the source tissue is known as plasticity. Endometrial epithelial and stromal cell proliferation, as discussed previously, is a complex multi-component process involving cues from extra-cellular growth factors and ovarian hormones [72, 94]. However, in their absence, isolated bovine endometrial stromal cells exhibit the ability to develop into bone . Results from our microarray study showed a large representation from this biological category, with 36 DEG enriched. Genes implicated in cell and tissue morphology and development which were down-regulated in low fertility heifers included, PPARA, IL6ST, GJA1, SFRP1 and IL-33.
One particular biochemical pathway which facilitates cellular transformation includes extracellular matrix mineralization (ECM) . A well known regulator of ECM mineralization is the activin a-FST system. Activin A inhibits ECM mineralization whereas FST, an activin antagonist which prevents activin-receptor interaction , increases mineralization in cell cultures . Transgenic female mice with gain-of-function FST, in which mouse follistatin was over-expressed, developed thin uteri and small ovaries, resulting in infertility . FST was differentially expressed between HF and LF heifers, indicating a role for this gene pathway in mid-luteal endometrial homeostasis and early embryo survival.
ECM remodeling, occurring during both pregnancy and the estrous cycle, facilitated by the matrix-metalloproteinases, ensures the provision of a suitable structural microenvironment where the embryo can grow [100, 101]. Matrix-metalloproteinase-19 (MMP-19), an important molecule in this pathway and which was down-regulated in LF heifers, plays a significant role in ECM remodelling . Interestingly, Wathes et al.  reported that differential expression of genes MMP - 1, 2, 3, and 13 two week post partum in the bovine endometrium, was highly correlated with differential expression of IGF binding protein 4, a known antagonist of IGF1 expression . The IGF system, in particular IGF1, is associated with several reproductive processes in cattle including preimplantation embryo development [105–107].
The transforming growth factor βs (TGF-β) are multifunctional cytokines that also regulate tissue remodelling and repair [108, 109]. High expression of TGF-β has been observed during pro-estrus and diestrus [110, 111] thereby highlighting the role for TGF-βs in endometrial remodelling, an important process impeding estrous cycle transition . Transforming growth factor beta 1 induced transcript (TGFB1I1) was down-regulated in the LF animals, suggesting altered or irregular endometrial remodelling in these animals which may be contributing to the conception rate differences observed between the two divergent fertility groups.
Inflammation is an innate cyclical physiological process facilitating progression of reproductive cycles in the endometrium. The animal model in this study isparticularly useful for the identification of inflammatory pathways associated with uterine low fertility for numerous reasons. Firstly, there has been no mitogenic challenge. This study strictly examines gene expression between high and low conception rate animals without influence from any exogenous metabolites, either dietary or pharmaceutical. Secondly, tissue sampling occurred during an estrous cycle where no embryo was present. Lastly, the study employed nulliparous heifers where the likelihood of uterine infection is low, as was demonstrated by the lack of clinical evidence of metritis, endometritis, pyometra or metaplasia across all heifers examined.
In total 55 DEG featured in inflammatory linked pathways. It is clear from the high proportion of DEG that inflammation is a central theme in estrous cycle and uterine sub-fertility physiology. IL-33, a cytokine which influences the production of other pro-inflammatory cytokines IL-5, IL-13 and chemokine GM-CSF was more highly expressed in LF animals. In addition, IL-33 regulates transcription of endothelial cells in inflamed rheumatic tissues . As mentioned previously, cell plasticity is altered in a state of chronic inflammation or trauma. Inflammation due to up-regulated IL-33 could be altering the constitution of the endometrium in the LF animals, and thus impeding embryo implantation. Hence, low conception rates could be directly linked to inflammation induced, altered cellular plasticity, in uterine endometrial tissues.
Similar to Forde et al. , Bauersachs et al.  and Salilew-Wondim et al. , genes coding metabolite transporters, specifically the solute carrier (SLC) family members were found to be differentially expressed between HF and LF animals. The five SLC genes identified were; SLC1A3, SLC17A5, SLC25A12, SLC25A24, SLC45A2. The most abundantly expressed gene of the entire DEG list, SLC45A2, was 8-fold more highly expressed in the uterus of LF relative to HF heifers. As the name suggests SLC genes are involved in the transfer of solutes across the cell membrane, particularly amino acids [115–117]. Amino acids are fundamental for the normal growth and development of the early embryo, acting as precursors of nucleic acids and proteins, osmolytes and signaling molecules. Concentrations of amino acids in oviductal and uterine fluid during the estrous cycle have been reported to modulate with stage of cycle, systemic progesterone environment and differ compared with plasma, demonstrating their active transport in these tissues [118–121]. The endometrium functions as a secretory layer, suggesting the importance of metabolite exchange in this specific tissue. Animals with less efficient metabolic exchange in the uterus may be unable to sustain embryo development during early pregnancy, and thus be experiencing recurring early embryo loss.
Microarray analysis was carried out on endometrial tissue, an amalgam of varying cell types. Examining tissue mRNA gene expression provides an insight into the genetic regulation of multiple cell types from the host. It was essential to use RNA from all endometrial cell types as it is not apparent, as of yet, whether or which individual endometrial cell types are contributing to low conception rates in cattle. Investigations into the types and locations of contributing cell types via in situ hybridisation or immunofluorescence would assist in the development of proposed hypotheses.
Global endometrial gene expression profiles during the mid-luteal phase of the estrous cycle, in HF and LF heifers was investigated, and the most significant biological pathways likely to be involved in uterine function and embryo survival identified. The new knowledge generated offers substantial insight into some of the molecular mechanisms underlying uterine endometrial function and uterine mediated low-fertility, during the early to mid luteal phase of the estrous cycle in cattle. Furthermore, expression analysis provides invaluable data on key differentially expressed genes which may be selected for future SNP discovery analysis which following validation may be used as genetic markers for fertility and incorporated into breeding programmes.
Availability of supporting data
The data sets supporting the results of this article are available in the NCBI’s Gene Expression Omnibus repository, GSE29853 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE29853.
Butler WR, Smith RD: Interrelationships between energy balance and postpartum reproductive function in dairy cattle. J Dairy Sci. 1989, 72 (3): 767-783. 10.3168/jds.S0022-0302(89)79169-4.
Royal M, Mann GE, Flint AP: Strategies for reversing the trend towards subfertility in dairy cattle. Vet J. 2000, 160 (1): 53-60. 10.1053/tvjl.1999.0450.
Royal MD, Darwash AO, Flint APF, Webb R, Woolliams JA, Lamming GE: Declining fertility in dairy cattle: changes in traditional and endocrine parameters of fertility. Anim Sci. 2000, 70: 487-501.
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.
Berg DK, van Leeuwen J, Beaumont S, Berg M, Pfeffer PL: Embryo loss in cattle between days 7 and 16 of pregnancy. Theriogenology. 2010, 73 (2): 250-260. 10.1016/j.theriogenology.2009.09.005.
Diskin MG, Sreenan JM: Fertilization and embryonic mortality rates in beef heifers after artificial insemination. J Reprod Fertil. 1980, 59 (2): 463-468. 10.1530/jrf.0.0590463.
Diskin MG, Morris DG: Embryonic and early foetal losses in cattle and other ruminants. Reprod Domest Anim. 2008, 43 (Suppl 2): 260-267.
Sartori R, Bastos MR, Wiltbank MC: Factors affecting fertilisation and early embryo quality in single and super-ovulated dairy cattle. Reprod Fer Dev. 2010, 22 (1): 151-158. 10.1071/RD09221.
Bazer FW, Ott TL, Spencer TE: Pregnancy recognition in ruminants, pigs and horses: Signals from the trophoblast. Theriogenology. 1994, 41 (1): 79-94. 10.1016/S0093-691X(05)80052-4.
McMillan WH, Donnison MJ: Understanding maternal contributions to fertility in recipient cattle: development of herds with contrasting pregnancy rates. Anim Reprod Sci. 1999, 57 (3): 127-140.
Cox SF, McMillan WH, Donnison MJ: Establishment of a herd of cattle with divergent potential pregnancy rate. Theriogenology. 1998, 49 (1): 242-10.1016/S0093-691X(98)90595-7.
Peterson AJ, Donnison MJ, Pearson S, McMillan WH: Contrasting early embryo development in a herd of recipient cattle with previously high or low pregnancy rates. Theriogenology. 1999, 51 (1): 229-229. 10.1016/S0093-691X(99)91788-0.
Diskin MG, Sreenan JM: Repeatability estimate for embryo survival following insemination at PG-induced heats in beef heifers. J Anim Sci. 2005, 83 (Supplement 1): 38-
Doyle SP, Golden BL, Green RD, Brinks JS: Additive genetic parameter estimates for heifer pregnancy and subsequent reproduction in Angus females. J Anim Sci. 2000, 78 (8): 2091-2098.
Thallman RM, Cundiff LV, Gregory KE, Koch RM: Germplasm evaluation in beef cattle–Cycle IV: postweaning growth and puberty of heifers. J Anim Sci. 1999, 77 (10): 2651-2659.
Bersinger NA, Wunder DM, Birkhauser MH, Mueller MD: Gene expression in cultured endometrium from women with different outcomes following IVF. Mol Hum Reprod. 2008, 14 (8): 475-484. 10.1093/molehr/gan036.
Morris D, Diskin M: Effect of progesterone on embryo survival. Animal. 2008, 2 (8): 1112-1119.
Diskin MG, Kenny DA, Dunne L, Sreenan JM: Systemic Progesterone pre and Post AI and Early Embryo Survival in Cattle. 2002, Tullamore, Ireland: Proceedings of the Agricultural Research Forum, 27-
Spencer TE, Johnson GA, Burghardt RC, Bazer FW: Progesterone and placental hormone actions on the uterus: insights from domestic animals. Biol Reprod. 2004, 71 (1): 2-10. 10.1095/biolreprod.103.024133.
Igwebuike UM: A review of uterine structural modifications that influence conceptus implantation and development in sheep and goats. Anim Reprod Sci. 2009, 112 (1–2): 1-7.
Faulkner S, Elia G, Mullen MP, O’Boyle P, Dunn MJ, Morris D: A comparison of the bovine uterine and plasma proteome using iTRAQ proteomics. Proteomics. 2012, 12 (12): 2014-2023. 10.1002/pmic.201100609.
Mullen MP, Elia G, Hilliard M, Parr MH, Diskin MG, Evans AC, Crowe MA: Proteomic characterization of histotroph during the preimplantation phase of the estrous cycle in cattle. J Proteome Res. 2012, 11 (5): 3004-3018. 10.1021/pr300144q.
Beltman ME, Forde N, Furney P, Carter F, Roche JF, Lonergan P, Crowe MA: Characterisation of endometrial gene expression and metabolic parameters in beef heifers yielding viable or non-viable embryos on Day 7 after insemination. Reprod Fertil Dev. 2010, 22 (6): 987-999. 10.1071/RD09302.
McCarthy SD, Roche JF, Forde N: Temporal changes in endometrial gene expression and protein localization of members of the IGF family in cattle: Effects of progesterone and pregnancy. Physiol Genomics. 2012, 44 (2): 130-140. 10.1152/physiolgenomics.00106.2011.
Forde N, Spencer TE, Bazer FW, Song G, Roche JF, Lonergan P: Effect of pregnancy and progesterone concentration on expression of genes encoding for transporters or secreted proteins in the bovine endometrium. Physiol Genomics. 2010, 41 (1): 53-62. 10.1152/physiolgenomics.00162.2009.
Mullen MP, Forde N, Parr MH, Diskin MG, Morris DG, Nally JE, Evans AC, Crowe MA: Alterations in systemic concentrations of progesterone during the early luteal phase affect RBP4 expression in the bovine uterus. Reprod Fertil Dev. 2012, 24 (5): 715-722. 10.1071/RD11246.
Wolf E, Bauersachs S: Functional genome research in reproductive biology and biotechnology - a minireview. Anim Sci Paper Rep. 2010, 28 (2): 123-130.
Forde N, Carter F, Spencer TE, Bazer FW, Sandra O, Mansouri-Attia N, Okumu LA, McGettigan PA, Mehta JP, McBride R, O'Gaora P, Roche JF, Lonergan P: Conceptus-induced changes in the endometrial transcriptome: how soon does the cow know she is pregnant?. Biol Reprod. 2011, 85 (1): 144-156. 10.1095/biolreprod.110.090019.
Bauersachs S, Ulbrich SE, Reichenbach H-D, Reichenbach M, Büttner M, Meyer HHD, Spencer TE, Minten M, Sax G, Winter G, Wolf E: Comparison of the effects of early pregnancy with human interferon, alpha 2 (IFNA2), on gene expression in bovine endometrium. Biol Reprod. 2012, 86 (2): 46-10.1095/biolreprod.111.094771.
Walker CG, Littlejohn MD, Mitchell MD, Roche JR, Meier S: Endometrial gene expression during early pregnancy differs between fertile and sub-fertile dairy cow strains. Physiol Genomics. 2011, 44 (1): 47-58.
Salilew-Wondim D, Holker M, Rings F, Ghanem N, Ulas-Cinar M, Peippo J, Tholen E, Looft C, Schellander K, Tesfaye D: Bovine pretransfer endometrium and embryo transcriptome fingerprints as predictors of pregnancy success after embryo transfer. Physiol Genomics. 2010, 42 (2): 201-218. 10.1152/physiolgenomics.00047.2010.
Forde N, Carter F, Fair T, Crowe MA, Evans ACO, Spencer TE, Bazer FW, McBride R, Boland MP, O’Gaora P, Lonergan P, Roche JF: Progesterone-regulated changes in endometrial gene expression contribute to advanced conceptus development in cattle. Biol Reprod. 2009, 81 (4): 784-794. 10.1095/biolreprod.108.074336.
Forde N, Beltman ME, Duffy GB, Duffy P, Mehta JP, O’Gaora P, Roche JF, Lonergan P, Crowe MA: Changes in the endometrial transcriptome during the bovine estrous cycle: effect of low circulating progesterone and consequences for conceptus elongation. Biol Reprod. 2010, 84 (2): 266-278.
Mitko K, Ulbrich SE, Wenigerkind H, Sinowatz F, Blum H, Wolf E, Bauersachs S: Dynamic changes in messenger RNA profiles of bovine endometrium during the oestrous cycle. Reproduction. 2008, 135 (2): 225-240. 10.1530/REP-07-0415.
Cooke PS, Spencer TE, Bartol FF, Hayashi K: Uterine glands: development, function and experimental model systems. Mol Hum Reprod. 2013, 19 (9): 547-558. 10.1093/molehr/gat031.
Cooke PS, Ekman GC, Kaur J, Davila J, Bagchi IC, Clark SG, Dziuk PJ, Hayashi K, Bartol FF: Brief exposure to progesterone during a critical neonatal window prevents uterine gland formation in mice. Biol Reprod. 2012, 86 (3): 63-10.1095/biolreprod.111.097188.
Morris DG, Waters SM, McCarthy SD, Patton J, Earley B, Fitzpatrick R, Murphy JJ, Diskin MG, Kenny DA, Brass A, Wathes DC: Pleiotropic effects of negative energy balance in the postpartum dairy cow on splenic gene expression: repercussions for innate and adaptive immunity. Physiol Genomics. 2009, 39 (1): 28-37. 10.1152/physiolgenomics.90394.2008.
Kadota K, Nakai Y, Shimizu K: Ranking differentially expressed genes from Affymetrix gene expression data: methods with reproducibility, sensitivity, and specificity. Algorithms Mol Biol. 2009, 4: 7-10.1186/1748-7188-4-7.
Pearson RD, Liu X, Sanguinetti G, Milo M, Lawrence ND, Rattray M: Puma: a Bioconductor package for propagating uncertainty in microarray analysis. BMC Bioinforma. 2009, 10: 211-10.1186/1471-2105-10-211.
Liu X, Milo M, Lawrence ND, Rattray M: Probe-level measurement error improves accuracy in detecting differential gene expression. Bioinformatics. 2006, 22 (17): 2107-2113. 10.1093/bioinformatics/btl361.
Pearson RD: A comprehensive re-analysis of the Golden Spike data: towards a benchmark for differential expression methods. BMC Bioinforma. 2008, 9: 164-10.1186/1471-2105-9-164.
Sanchez-Calderon H, Rodriguez-De La Rosa L, Milo M, Pichel JG, Holley M, Varela-Nieto I: RNA microarray analysis in prenatal mouse cochlea reveals novel IGF-I target genes: implication of MEF2 and FOXM1 transcription factors. PLoS One. 2010, 5 (1): e8699-10.1371/journal.pone.0008699.
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucl Acids Res. 2005, 33 (20): e175-10.1093/nar/gni179.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20 (3): 307-315. 10.1093/bioinformatics/btg405.
de Leeuw W, Rauwerda H, Jonker M, Breit T: Salvaging Affymetrix probes after probe-level re-annotation. BMC Res Notes. 2008, 1 (1): 66-10.1186/1756-0500-1-66.
McCarthy SD, Waters SM, Kenny DA, Diskin MG, Fitzpatrick R, Patton J, Wathes DC, Morris DG: Negative energy balance and hepatic gene expression patterns in high-yielding dairy cows during the early postpartum period: a global approach. Physiol Genomics. 2010, 42A (3): 188-199. 10.1152/physiolgenomics.00118.2010.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): RESEARCH0034-
Loor JJ, Everts RE, Bionaz M, Dann HM, Morin DE, Oliveira R, Rodriguez-Zas SL, Drackley JK, Lewin HA: Nutrition-induced ketosis alters metabolic and signaling gene networks in liver of periparturient dairy cows. Physiol Genomics. 2007, 32 (1): 105-116. 10.1152/physiolgenomics.00188.2007.
Barilli A, Rotoli BM, Visigalli R, Bussolati O, Gazzola GC, Kadija Z, Rodi G, Mariani F, Ruzza ML, Luisetti M, Dall'Asta V: Lysinuric Protein Intolerance system y + L activity is defective in monocytes and in GM-CSF-differentiated macrophages. Orphanet J Rare Dis. 2010, 5: 32-10.1186/1750-1172-5-32.
Selvey S, Thompson EW, Matthaei K, Lea RA, Irving MG, Griffiths LR: [beta]-Actin–an unsuitable internal control for RT-PCR. Mol Cell Probes. 2001, 15 (5): 307-311. 10.1006/mcpr.2001.0376.
Coyne GS, Kenny DA, Childs S, Sreenan JM, Waters SM: Dietary n-3 polyunsaturated fatty acids alter the expression of genes involved in prostaglandin biosynthesis in the bovine uterus. Theriogenology. 2008, 70 (5): 772-782. 10.1016/j.theriogenology.2008.05.048.
Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004, 313 (4): 856-862. 10.1016/j.bbrc.2003.11.177.
Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29 (2): 332-337.
Rozen S, Skaletsky H: Primer3 on the WWW for General Users and for Biologist Programmers. Methods Mol Biol. 2000, 132: 365-386.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.
Diskin MG KD, Dunne LD, Sreenan JM: Systemic Progesterone pre- and Post- AI and Embryo Survival in Heifers. 2002, Tullamore, Ireland: ‘Irish Agricultural Research Forum’, 27-
Childs S, Lynch CO, Hennessy AA, Stanton C, Wathes DC, Sreenan JM, Diskin MG, Kenny DA: Effect of dietary enrichment with either n-3 or n-6 fatty acids on systemic metabolite and hormone concentration and ovarian function in heifers. Animal. 2008, 2 (06): 883-893.
Chen X, Wang L, Smith JD, Zhang B: Supervised principal component analysis for gene set enrichment of microarray data with continuous or survival outcomes. Bioinformatics. 2008, 24 (21): 2474-2481. 10.1093/bioinformatics/btn458.
Goeman JJ, Buhlmann P: Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007, 23 (8): 980-987. 10.1093/bioinformatics/btm051.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Wang K, Zhang H, Kugathasan S, Annese V, Bradfield JP, Russell RK, Sleiman PM, Imielinski M, Glessner J, Hou C, Wilson DC, Walters T, Kim C, Frackelton EC, Lionetti P, Barabino A, Van Limbergen J, Guthery S, Denson L, Piccoli D, Li M, Dubinsky M, Silverberg M, Griffiths A, Grant SF, Satsangi J, Baldassano R, Hakonarson H: Diverse genome-wide association studies associate the IL12/IL23 pathway with Crohn Disease. Am J Hum Genet. 2009, 84 (3): 399-405. 10.1016/j.ajhg.2009.01.026.
Waagmeester A, Kelder T, Evelo C: The role of bioinformatics in pathway curation. Genes Nutr. 2008, 3 (3): 139-142.
Kurita T, Young P, Brody JR, Lydon JP, O’Malley BW, Cunha GR: Stromal progesterone receptors mediate the inhibitory effects of progesterone on estrogen-induced uterine epithelial cell deoxyribonucleic acid synthesis. Endocrinology. 1998, 139 (11): 4708-4713.
Johnson ML, Redmer DA, Reynolds LP: Uterine growth, cell proliferation, and c-fos proto-oncogene expression throughout the estrous cycle in ewes. Biol Reprod. 1997, 56 (2): 393-401. 10.1095/biolreprod56.2.393.
Di Simone N, Crowley WF, Wang QF, Sluss PM, Schneyer AL: Characterization of inhibin/activin subunit, follistatin, and activin type II receptors in human ovarian cancer cell lines: a potential role in autocrine growth regulation. Endocrinology. 1996, 137 (2): 486-494.
Tao J, Mallat A, Gallois C, Belmadani S, Mery PF, Nhieu JT, Pavoine C, Lotersztajn S: Biological effects of C-type natriuretic peptide in human myofibroblastic hepatic stellate cells. J Biol Chem. 1999, 274 (34): 23761-23769. 10.1074/jbc.274.34.23761.
Zhang YW, Nakayama K, Morita I: A novel route for connexin 43 to inhibit cell proliferation: negative regulation of S-phase kinase-associated protein (Skp 2). Cancer Res. 2003, 63 (7): 1623-1630.
Iguchi H, Urashima Y, Inagaki Y, Ikeda Y, Okamura M, Tanaka T, Uchida A, Yamamoto TT, Kodama T, Sakai J: SOX6 suppresses cyclin D1 promoter activity by interacting with beta-catenin and histone deacetylase 1, and its down-regulation induces pancreatic beta-cell proliferation. J Biol Chem. 2007, 282 (26): 19052-19061. 10.1074/jbc.M700460200.
Shimizu T, Krebs S, Bauersachs S, Blum H, Wolf E, Miyamoto A: Actions and interactions of progesterone and estrogen on transcriptome profiles of the bovine endometrium. Physiol Genomics. 2010, 42A (4): 290-300. 10.1152/physiolgenomics.00107.2010.
Mansouri-Attia N, Aubert J, Reinaud P, Giraud-Delville C, Taghouti G, Galio L, Everts RE, Degrelle S, Richard C, Hue I, Yang X, Tian XC, Lewin HA, Renard JP, Sandra O: Gene expression profiles of bovine caruncular and intercaruncular endometrium at implantation. Physiol Genomics. 2009, 39 (1): 14-27. 10.1152/physiolgenomics.90404.2008.
Brenner RM, Slayden OD, Critchley HO: Anti-proliferative effects of progesterone antagonists in the primate endometrium: a potential role for the androgen receptor. Reproduction. 2002, 124 (2): 167-172. 10.1530/rep.0.1240167.
Kooy J, Taylor NH, Healy DL, Rogers PA: Endothelial cell proliferation in the endometrium of women with menorrhagia and in women following endometrial ablation. Hum Reprod. 1996, 11 (5): 1067-1072. 10.1093/oxfordjournals.humrep.a019298.
Klagsbrun M, D’Amore PA: Regulators of angiogenesis. Annu Rev Physiol. 1991, 53: 217-239. 10.1146/annurev.ph.53.030191.001245.
Peeters LL, Vigne JL, Tee MK, Zhao D, Waite LL, Taylor RN: PPAR gamma represses VEGF expression in human endometrial cells: implications for uterine angiogenesis. Angiogenesis. 2005, 8 (4): 373-379.
Fraser HM, Duncan WC: SRB reproduction, fertility and development award lecture 2008. Regulation and manipulation of angiogenesis in the ovary and endometrium. Reprod Fertil Dev. 2009, 21 (3): 377-392. 10.1071/RD08272.
Egeblad M, Werb Z: New functions for the matrix metalloproteinases in cancer progression. Nat Rev Cancer. 2002, 2 (3): 161-174. 10.1038/nrc745.
Ronfani L, Ferraguti M, Croci L, Ovitt CE, Scholer HR, Consalez GG, Bianchi ME: Reduced fertility and spermatogenesis defects in mice lacking chromosomal protein Hmgb2. Development. 2001, 128 (8): 1265-1273.
Seyedin SM, Kistler WS: Levels of chromosomal protein high mobility group 2 parallel the proliferative activity of testis, skeletal muscle, and other organs. J Biol Chem. 1979, 254 (22): 11264-11271.
Ly DH, Lockhart DJ, Lerner RA, Schultz PG: Mitotic misregulation and human aging. Science. 2000, 287 (5462): 2486-2492. 10.1126/science.287.5462.2486.
Klune JR, Dhupar R, Cardinal J, Billiar TR, Tsung A: HMGB1: endogenous danger signaling. Mol Med. 2008, 14 (7–8): 476-484.
Lee H, Shin N, Song M, Kang UB, Yeom J, Lee C, Ahn YH, Yoo JS, Paik YK, Kim H: Analysis of nuclear high mobility group box 1 (HMGB1)-binding proteins in colon cancer cells: clustering with proteins involved in secretion and extranuclear function. J Proteome Res. 2010, 9 (9): 4661-4670. 10.1021/pr100386r.
Zicari A, Centonze C, Realacci M, Buchetti B, Pietropolli A, Ticconi C: Estradiol 17-beta and progesterone modulate inducible nitric oxide synthase and high mobility group box 1 expression in human endometrium. Reprod Sci. 2008, 15 (6): 559-566. 10.1177/1933719107312560.
Francis GA, Fayard E, Picard F, Auwerx J: Nuclear receptors and the control of metabolism. Ann Rev Physiol. 2003, 65 (1): 261-311. 10.1146/annurev.physiol.65.092101.142528.
Mattos R, Staples CR, Thatcher WW: Effects of dietary fatty acids on reproduction in ruminants. Rev Reprod. 2000, 5 (1): 38-45. 10.1530/ror.0.0050038.
Thompson BJ, Smith S: Biosynthesis of fatty acids by lactating human breast epithelial cells: an evaluation of the contribution to the overall composition of human milk fat. Pediatr Res. 1985, 19 (1): 139-143. 10.1203/00006450-198501000-00036.
Wakil SJ: Fatty acid synthase, a proficient multifunctional enzyme. Biochemistry. 1989, 28 (11): 4523-4530. 10.1021/bi00437a001.
Pizer ES, Kurman RJ, Pasternack GR, Kuhajda FP: Expression of fatty acid synthase is closely linked to proliferation and stromal decidualization in cycling endometrium. Int J Gynecol Pathol. 1997, 16 (1): 45-51. 10.1097/00004347-199701000-00008.
Mahendroo MS, Porter A, Russell DW, Word RA: The parturition defect in steroid 5alpha-reductase type 1 knockout mice is due to impaired cervical ripening. Mol Endocrinol. 1999, 13 (6): 981-992.
Gargett CE, Chan RW, Schwab KE: Hormone and growth factor signaling in endometrial renewal: role of stem/progenitor cells. Mol Cell Endocrinol. 2008, 288 (1–2): 22-29.
Donofrio G, Franceschi V, Capocefalo A, Cavirani S, Sheldon IM: Bovine endometrial stromal cells display osteogenic properties. Reprod Biol Endocrinol. 2008, 6: 65-10.1186/1477-7827-6-65.
Walter I, Helmreich M, Handler J, Aurich C: Mineralised deposits in the uterine glands of mares with chronic endometrial degeneration. Vet Rec. 2003, 153 (23): 708-710.
Lainas T, Zorzovilis I, Petsas G, Alexopoulou E, Lainas G, Ioakimidis T: Osseous metaplasia: case report and review. Fertil Steril. 2004, 82 (5): 1433-1435. 10.1016/j.fertnstert.2004.04.055.
Roth E, Taylor HB: Heterotopic cartilage in the uterus. Obstet Gynecol. 1966, 27 (6): 838-844.
Murshed M, Harmey D, Millán JL, McKee MD, Karsenty G: Unique coexpression in osteoblasts of broadly expressed genes accounts for the spatial restriction of ECM mineralization to bone. Genes Dev. 2005, 19 (9): 1093-1104. 10.1101/gad.1276205.
Abedin M, Tintut Y, Demer LL: Vascular calcification: mechanisms and clinical ramifications. Arterioscler Thromb Vasc Biol. 2004, 24 (7): 1161-1170. 10.1161/01.ATV.0000133194.94939.42.
Nakamura T, Takio K, Eto Y, Shibai H, Titani K, Sugino H: Activin-binding protein from rat ovary is follistatin. Science. 1990, 247 (4944): 836-838. 10.1126/science.2106159.
Eijken M, Swagemakers S, Koedam M, Steenbergen C, Derkx P, Uitterlinden AG, van der Spek PJ, Visser JA, de Jong FH, Pols HAP, van Leeuwen JP: The activin A-follistatin system: potent regulator of human extracellular matrix mineralization. FASEB J. 2007, 21 (11): 2949-2960. 10.1096/fj.07-8080com.
Guo Q, Kumar TR, Woodruff T, Hadsell LA, DeMayo FJ, Matzuk MM: Overexpression of mouse follistatin causes reproductive defects in transgenic mice. Mol Endocrinol. 1998, 12 (1): 96-106. 10.1210/mend.12.1.0053.
Hashizume K, Takahashi T, Shimizu M, Todoroki J, Shimada A, Hirata M, Sato T, Ito A: Matrix-metalloproteinases-2 and −9 production in bovine endometrial cell culture. J Reprod Dev. 2003, 49 (1): 45-53. 10.1262/jrd.49.45.
Salamonsen LA: Role of proteases in implantation. Rev Reprod. 1999, 4 (1): 11-22. 10.1530/ror.0.0040011.
Vu TH, Werb Z: Matrix metalloproteinases: effectors of development and normal physiology. Genes Dev. 2000, 14 (17): 2123-2133. 10.1101/gad.815400.
Wathes DC, Cheng Z, Fenwick MA, Fitzpatrick R, Patton J: Influence of energy balance on the somatotrophic axis and matrix metalloproteinase expression in the endometrium of the postpartum dairy cow. Reproduction. 2011, 141 (2): 269-281. 10.1530/REP-10-0177.
Contois LW, Nugent DP, Caron JM, Cretu A, Tweedie E, Akalu A, Liebes L, Friesel R, Rosen C, Vary C, Brooks PC: Insulin-like growth factor binding protein-4 (IGFBP-4) differentially inhibits growth factor-induced angiogenesis. J Biol Chem. 2011, 287 (3): 1779-1789.
Coyne GS, Kenny DA, Waters SM: Effect of dietary n-3 polyunsaturated fatty acid supplementation on bovine uterine endometrial and hepatic gene expression of the insulin-like growth factor system. Theriogenology. 2011, 75 (3): 500-512. 10.1016/j.theriogenology.2010.09.018.
Mullen MP, Berry DP, Howard DJ, Diskin MG, Lynch CO, Giblin L, Kenny DA, Magee DA, Meade KG, Waters SM: Single nucleotide polymorphisms in the insulin-like growth factor 1 (IGF-1) gene are associated with performance in Holstein-Friesian dairy cattle. Front Gen. 2011, 2: 3-
Velazquez MA, Newman M, Christie MF, Cripps PJ, Crowe MA, Smith RF, Dobson H: The usefulness of a single measurement of insulin-like growth factor-1 as a predictor of embryo yield and pregnancy rates in a bovine MOET. program. Theriogenology. 2005, 64 (9): 1977-1994. 10.1016/j.theriogenology.2005.05.001.
Massague J: The transforming growth factor-beta family. Annu Rev Cell Biol. 1990, 6 (1): 597-641. 10.1146/annurev.cb.06.110190.003121.
Rizzino A: Transforming growth factor-[beta]: Multiple effects on cell differentiation and extracellular matrices. Dev Biol. 1988, 130 (2): 411-422. 10.1016/0012-1606(88)90337-5.
Doré JJ, Wilkinson JE, Godkin JD: Ovine endometrial expression of transforming growth factor beta isoforms during the peri-implantation period. Biol Reprod. 1996, 54 (5): 1080-1087. 10.1095/biolreprod54.5.1080.
Godkin JD, Dore JJ: Transforming growth factor beta and the endometrium. Rev Reprod. 1998, 3 (1): 1-6. 10.1530/ror.0.0030001.
Chegini N, Gold LI, Williams RS: Localization of transforming growth factor beta isoforms TGF-beta 1, TGF-beta 2, and TGF-beta 3 in surgically induced endometriosis in the rat. Obstet Gynecol. 1994, 83 (3): 455-461.
Allakhverdi Z, Smith DE, Comeau MR, Delespesse G: Cutting edge: the ST2 ligand IL-33 potently activates and drives maturation of human mast cells. J Immunol. 2007, 179 (4): 2051-2054.
Barksby HE, Lea SR, Preshaw PM, Taylor JJ: The expanding family of interleukin-1 cytokines and their role in destructive inflammatory disorders. Clin Exp Immunol. 2007, 149 (2): 217-225. 10.1111/j.1365-2249.2007.03441.x.
Hediger MA, Romero MF, Peng JB, Rolfs A, Takanaga H, Bruford EA: The ABCs of solute carriers: physiological, pathological and therapeutic implications of human membrane transport proteinsIntroduction. Pflugers Arch. 2004, 447 (5): 465-468. 10.1007/s00424-003-1192-y.
Kanai Y, Hediger M: The glutamate/neutral amino acid transporter family SLC1: molecular, physiological and pharmacological aspects. Pflügers Archiv Eur J Physiol. 2004, 447 (5): 469-479. 10.1007/s00424-003-1146-4.
Palmieri F: The mitochondrial transporter family (SLC25): physiological and pathological implications. Pflügers Archiv Eur J Physiol. 2004, 447 (5): 689-709. 10.1007/s00424-003-1099-7.
Hugentobler SA, Diskin MG, Leese HJ, Humpherson PG, Watson T, Sreenan JM, Morris DG: Amino acids in oviduct and uterine fluid and blood plasma during the estrous cycle in the bovine. Mol Reprod Dev. 2007, 74 (4): 445-454. 10.1002/mrd.20607.
Leese HJ, Hugentobler SA, Gray SM, Morris DG, Sturmey RG, Whitear SL, Sreenan JM: Female reproductive tract fluids: composition, mechanism of formation and potential role in the developmental origins of health and disease. Reprod Fertil Dev. 2008, 20 (1): 1-8. 10.1071/RD07153.
Mullen MP, Bazer FW, Wu G, Parr MH, Evans ACO, Crowe MA, Diskin MG: Effects of systemic progesterone during the early luteal phase on the availabilities of amino acids and glucose in the bovine uterine lumen. Reprod Fertil Dev. 2013, 26 (2): 282-292.
Hugentobler SA, Sreenan JM, Humpherson PG, Leese HJ, Diskin MG, Morris DG: Effects of changes in the concentration of systemic progesterone on ions, amino acids and energy substrates in cattle oviduct and uterine fluid and blood. Reprod Fertil Dev. 2010, 22 (4): 684-694. 10.1071/RD09129.
Bennett EP, Hassan H, Mandel U, Hollingsworth MA, Akisawa N, Ikematsu Y, Merkx G, van Kessel AG, Olofsson S, Clausen H: Cloning and characterization of a close homologue of human UDP-N-acetyl-alpha-D-galactosamine:Polypeptide N-acetylgalactosaminyltransferase-T3, designated GalNAc-T6. Evidence for genetic but not functional redundancy. J Biol Chem. 1999, 274 (36): 25362-25370. 10.1074/jbc.274.36.25362.
The authors would like to thank staff at Teagasc Athenry for the management and care of the cows used in this study. A. Killeen was in receipt of a Teagasc Walsh Fellowship. The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific work.
The authors declare that they have no competing interests.
AK and SW carried out molecular gene expression studies, participated in the microarray data analysis, carried out the pathway analysis, statistical analysis and drafted the manuscript. SW, MD and DK conceived of the study, and participated in its design and coordination and helped to draft the manuscript. DM performed the microarray data analysis, participated in the pathway analysis and assisted in the manuscript editing. MD and MM participated in the metabolite measurement assays. MD and DK conducted statistical analyses. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S1: Up-regulated DEG (P < 0.05): Entrez ID, Symbol, Entrez Gene Name, Fold Change. Table S2. Down-regulated DEG (P < 0.05): Entrez ID, Symbol, Entrez Gene Name, Fold Change. Table S3. Genes validated between RT-qPCR and microarray methodologies, including Fold changes, P-values and correlation coefficients. (DOC 290 KB)
About this article
Cite this article
Killeen, A.P., Morris, D.G., Kenny, D.A. et al. Global gene expression in endometrium of high and low fertility heifers during the mid-luteal phase of the estrous cycle. BMC Genomics 15, 234 (2014). https://doi.org/10.1186/1471-2164-15-234
- Estrous Cycle
- Differentially Express Gene
- Endometrial Tissue
- Conception Rate
- Embryo Survival