Skip to main content

Evaluation of real-time PCR endogenous control genes for analysis of gene expression in bovine endometrium



Quantitative real-time PCR gene expression results are generally normalised using endogenous control genes. These reference genes should be expressed at a constant level across all sample groups in a study, and should not be influenced by study treatments or conditions. There has been no systematic investigation of endogenous control genes for bovine endometrium to date. The suitability of both commonly used and novel endogenous control genes was evaluated in this study, with the latter being selected from stably expressed transcripts identified through microarray analysis of bovine endometrium. Fifteen candidate endogenous control genes were assessed across different tissue subtypes in pregnant and cycling Holstein-Friesian dairy cows from two divergent genetic backgrounds.


The expression profiles of five commonly used endogenous control genes (GAPDH, PPIA, RPS9, RPS15A, and UXT) and 10 experimentally derived candidate endogenous control genes (SUZ12, C2ORF29, ZNF131, ACTR1A, HDAC1, SLC30A6, CNOT7, DNAJC17, BBS2, and RANBP10) were analysed across 44 samples to determine the most stably expressed gene. Gene stability was assessed using the statistical algorithms GeNorm and Normfinder. All genes presented with low overall variability (0.87 to 1.48% CV of Cq). However, when used to normalise a differentially expressed gene (oxytocin receptor - OXTR) in the samples, the reported relative gene expression levels were significantly affected by the control gene chosen. Based on the results of this analysis, SUZ12 is proposed as the most appropriate control gene for use in bovine endometrium during early pregnancy or the oestrus cycle.


This study establishes the suitability of novel endogenous control genes for comparing expression levels in endometrial tissues of pregnant and cycling bovines, and demonstrates the utility of microarray analysis as a method for identifying endogenous control gene candidates.


Quantitative real-time reverse transcription PCR (RT-PCR) is an extremely sensitive technique that allows the precise measurement of gene expression across more than seven orders of magnitude[1, 2]. RT-PCR is often considered the gold standard for quantifying gene expression, and is commonly used to validate techniques with greater throughput but less overall sensitivity, such as microarray analysis [35]. RT-PCR relies on the use of fluorescent dyes to quantify transcript amplification, with the amplification cycle number at which these dyes/transcripts are detected (above background) giving an indication as to the relative abundance of the target molecules. The sensitivity of RT-PCR makes it a powerful tool for gene expression measurement, especially when sample quantities are limited or a transcript is expressed at a low level. However, this sensitivity also means that a great deal of care must be taken with regards to experimental design and implementation of the procedure.

When designing an experiment to evaluate gene expression in a group of samples, a number of critical factors must be kept constant. These include RNA extraction, DNase treatments, and cDNA synthesis. Normalisation of RT-PCR results is required to control inter-sample differences that may arise as a result of these sample processing steps, and ensure the gene expression of target transcripts are robustly quantifiable [6, 7].

The most common method for normalising RT-PCR data involves the use of one or more endogenous control genes. An ideal endogenous control gene is one that is stably expressed within the samples to be compared, regardless of tissue differences, experimental conditions, or treatments. Choosing an endogenous control gene to normalise gene expression data is one of the most crucial steps in the experimental design. Genes used as endogenous controls in RT-PCR experiments are often chosen with little prior knowledge of their expression over the experimental conditions examined, and are often selected arbitrarily from a pool of commonly used endogenous control genes such as GAPDH, and β-actin [8]. The most widely used endogenous control gene in studies of endometrial gene expression is GAPDH [912]. However, the suitability of GAPDH as an endogenous control gene has recently come into question, especially due to its potential regulation in a wide variety of physiological states [13], making it a questionable choice for RT-PCR normalisation [14].

Over the past three decades, genetic selection for milk production has resulted in a significant decline in dairy cattle fertility [15]. The fertilisation rate in dairy cattle is around 90% and does not differ between low-moderate and high-producing animals. However, the calving rate in lower producing animals is approximately 55%, whereas in high-producing animals this rate is approximately 35%. Pregnancy losses are thought to occur primarily during the pregnancy recognition/pre-implantation period [16], making studies of endometrial gene expression critical to further understanding of pregnancy establishment, recognition and maintenance within the bovine reproductive cycle.

The primary aim of this study was to identify suitable endogenous control genes for analysis of endometrial tissues from pregnant and cycling bovines. This study also aimed to investigate the potential use of microarray data analysis for identification of novel endogenous control genes, and the effect of endogenous control gene selection on the calculated expression of a target gene.

A total of 15 candidate endogenous control genes were analysed in 44 samples representing two different tissues (intercaruncular and caruncular) from 22 animals. These animals were either pregnant or cycling at day 17 of the reproductive cycle, and represented Holstein-Friesian cows from two divergent genetic backgrounds (North American (NA), and New Zealand (NZ)).

Two strategies were employed to identify the candidates. Five genes were selected on the basis that they had been previously used as housekeeping genes [1719], and an additional 10 novel genes were derived from a microarray experiment based around the same 44 samples used in the current analysis. Genespring GX software was used to generate a list through filtering on expression stability across the 44 samples. This list was subjected to GeNorm [20] and Normfinder [21] analysis to identify the 10 most suitable genes. The suitability of all 15 genes was then tested through statistical analyses, including a comparison of expression stability as determined by GeNorm and Normfinder algorithms. The effect of using these endogenous control genes was then evaluated using relative quantification of a gene known to be differentially expressed in the study.

Results and Discussion

Microarray analysis of the 44 endometrial samples revealed 27 transcripts with a high degree of expression stability (filtering on expression level - upper limit 1.2, lower limit 0.833). GeNorm and Normfinder were utilised to identify the 10 most stably expressed transcripts for further analysis. For RT-PCR design, full length transcripts were identified by querying microarray probe sequences against the bovine genome (Btau3,1) using NCBI BLAST

Gene expression levels of the candidate endogenous control genes (expressed in Cq values) are displayed in Table 1 and Figure 1. Cq values for sample replicates had very low variability with a mean intra-assay coefficient of variation (CV) of 0.41%. All genes had low overall variability, with the Cq range between 1.06 and 2.04 cycles, standard deviations ranging from 0.25 to 0.53 cycles, and CV values ranging from 0.87 to 1.48% Table 1.

Table 1 Quantification cycle (Cq) values and statistics for 15 candidate endogenous control genes assayed across 44 bovine endometrial samples.
Figure 1

Expression levels of candidate endogenous control genes in pregnant (P) and cycling (Cy) endometrial tissue samples. Values are given as quantification cycle (Cq) numbers. Boxes represent the lower and upper quartiles with medians; whiskers illustrate the maximum and minimums of the samples. There were no significant differences (P > 0.05) between the Cq means (T-test) or variances (F-test) of pregnant and cycling animals for any genes tested (except the variance of BBS2, P = 0.039).

Significance calculations between gene expression data for pregnant and cycling animals were performed on Cq and relative concentration values, as estimated through absolute quantification using the Roche LC480 software. No significant differences (P > 0.05) in means or variances (Cq or concentration) between pregnant and cycling endometrial tissues were apparent for any candidate genes (except BBS2, which had a p-value of 0.04 (F-test) for variance between pregnant and cycling animals). The means and variances between the two different strains of Holstein-Friesians also lacked significance (data not shown).

Cqs and relative concentration variances between the two different tissue types showed no significant differences, however means were significantly different for all genes, except RPS15A (P = 0.39 (T-test), P = 0.032 (ANOVA)).

The differences in means calculated between tissue subtypes likely reflect the distinct morphological and functional differences between caruncular and intercaruncular endometrium, which relate to their respective roles in reproduction. The 'caruncules' of the endometrium are specialised projections that are the site of embryo attachment. Caruncules become highly vascularised, and are the major site for small molecule and gaseous exchange. In comparison, intercaruncular tissue is highly glandular and responsible for early nourishment of the embryo through secretions of large molecules into the uterus [2224]. The intercaruncular tissue is often thought to be more important in early pregnancy and, therefore, the majority of expression studies in pre-implantation bovine endometrium focus solely on intercaruncular tissue gene expression[2527]. There is very little reported expression analysis of caruncular endometrium and no known studies comparing expression profiles of the two tissues.

Expression stability testing of candidate genes

To further analyse the suitability of the candidate genes for use as endogenous controls for bovine endometrial tissues, expression stability was assessed using the GeNorm [20] and Normfinder [21] algorithms.

GeNorm rankings for the 15 genes tested are presented in Table 2 and Figure 2A. GeNorm identified SUZ12 and ZNF131 as the most stably expressed genes of the 15 candidates. Four of the five most stable genes were those derived from microarray data; the other was RPS15A, which was selected from the literature. By contrast, four of the five least stable genes were chosen from the literature, including GAPDH. The most stably expressed gene identified by Normfinder was SUZ12 (Table 2), which was also one of the two most stable genes identified in the GeNorm analysis. The best combination of genes identified by Normfinder was SUZ12 and C2ORF29. The five most stably expressed genes according to Normfinder consisted entirely of microarray-derived genes, while the least stable five contained three out of five genes selected from the literature, again including GAPDH. The comparative ranking of all genes for both GeNorm and Normfinder algorithm analyses is displayed in Table 2.

Table 2 Stability ranking of candidate endogenous control genes from Normfinder and GeNorm analyses.
Figure 2

GeNorm output. Average expression stability value (M) for candidate endogenous control genes in bovine endometrial tissue samples. The M-value threshold for stability of a gene according to GeNorm is 1.5. The most stable genes as determined through analysis of the pairwise variation of each gene with all other genes were SUZ12 and ZNF131.

GeNorm calculates an expression stability value (M) for each candidate gene based on pairwise comparisons of variability. Each gene is compared to every other gene to determine the two genes that contain the least variation. The stability value calculated for each gene is used to rank genes from least to most stable. The authors of the method give an M-value of 1.5 as a cut-off for suitability as an endogenous control gene. The principal behind the pairwise stability measure ranking is that two ideal candidate normalisation genes should have an equal expression ratio in all samples [28]. In the present study, all genes are well below the stipulated 1.5 M-value. The program then calculates a normalisation factor in each sample for the most stable genes (data not presented). GeNorm also calculates the optimal number of endogenous control genes to be used in the analysis of gene expression (Table 3). This value is determined by locating the point where addition of the next most stable gene does not significantly affect the normalisation factor, using a cut-off value of 0.15. In this study, the value for 2 genes was 0.042, suggesting that 2 genes should be sufficient to normalise the experimental data.

Table 3 GeNorm output used to determine the optimal number of endogenous control genes for normalisation in bovine endometrial tissue samples.

Normfinder is another freely available tool for the identification of stable endogenous control genes. The main point of difference between the two methods is that Normfinder takes into account both inter- and intra-group variability. The program not only identifies the most stable pair of endogenous control genes but also identifies the best overall endogenous control gene. The calculation of variability between groups is especially important in the present study considering the significant expression differences between the two tissue subtypes. The use of the two most stably expressed genes, in this case SUZ12 and C2ORF29 (Normfinder), should provide sufficient normalisation for tissue comparisons as Normfinder selects the best combination of genes whilst taking into account any grouping effects such as tissue type.

The differences in rankings of gene stability using the two algorithms could be due to the fact that they use very different methods to assess gene stability. GeNorm selects genes based on the pairwise variation between genes. The two most stably expressed genes are therefore those genes that share an expression profile. In contrast, Normfinder uses a model based algorithm that takes into account overall stability as well as any groups that may be present in the sample set. For example, if there are any grouping effects on gene expression a gene would be ranked lower than one that demonstrated variability not associated with any particular group.

Effect of Endogenous control gene on target gene relative quantification

Temporal down-regulation of endometrial oxytocin receptor (OXTR) expression is a hallmark of early pregnancy, with embryonic interferon tau (IFNτ) thought to elicit this response [29, 30]. Given the expectation of differential OXTR expression between pregnant and cycling animals, the effect of control gene stability on gene expression values of OXTR in the 44 endometrial samples was tested.

Figure 3 presents relative OXTR expression levels when normalised with endogenous control genes of varying stability - the most stable gene identified by Normfinder (SUZ12), the two most stable genes identified by GeNorm (SUZ12 and ZNF131), and two of the least stably expressed genes identified by both Normfinder and GeNorm (GAPDH and UXT). Normfinder identified SUZ12 and C2ORF29 as the best combination of genes, but the relative expression was not significantly different from that calculated using only SUZ12 or a combination of SUZ12 and ZNF131 (data not presented). Correlation of normalised RT-PCR data to OXTR microarray expression data was not affected by choice of normalisation strategy. When compared to microarray reported expression, all calculated expression values had correlation coefficients of 0.79.

Figure 3

Relative gene expression results for OXTR when normalised with endogenous control genes of different stability in pregnant (A) and cycling (B) animals. These graphs show calculated relative expression values for OXTR when normalised to the most stable gene (SUZ12), two of the least stable genes (UXT and GAPDH), and to the two best genes according to GeNorm analysis (SUZ12 and ZNF131), with the latter using a normalisation factor calculated from the GeNorm analysis. (Pregnant = P, cycling = Cy, caruncular = a, intercaruncular = b, North American = NA, New Zealand = NZ)

OXTR expression was significantly greater in cycling than in pregnant cows regardless of the endogenous control gene used (Figure 3, ANOVA, P < 0.01). Notably, OXTR expression in pregnant animals was greater on average and more variable in NA animals (which have lower fertility in general [31]) than for NZ animals (Figure 3A and 3B), and was also related to embryo size (data not presented). The use of different endogenous control genes had no effect on OXTR expression differences in these group comparisons.

However, there were differences in the normalised OXTR expression values that were inconsistent across the 44 samples, depending on which endogenous control genes were employed. For example, in sample Pregnant NZ-5a the calculated relative OXTR expression value was 12.46 when normalised to SUZ12, 6.32 when normalised to GAPDH, 5.10 when normalised to UXT, and 12.04 when normalised to both SUZ12 and ZNF131. Another sample (Pregnant NZ-2a), the values were 10.15, 9.60, 5.35, and 11.26 when normalised to SUZ12, GAPDH, UXT and both SUZ12 and ZNF131, respectively.

The average calculated fold change difference between pregnant and cycling animals was not affected by the choice of normalisation gene, possibly due to the large difference in expression level for this comparison (10-fold average difference between pregnant and cycling animals). However, the choice of reference gene could be important when normalising genes that exhibit more subtle variation between experimental groups, given the considerable variation in expression shown between individual samples.


This study provides the first reported assessment of endogenous control genes for use in expression studies in bovine endometrium. Normalisation is a critical factor in reporting RT-PCR expression data, providing a necessary control for error associated with sample preparation. Normalisation using endogenous control genes provides a means of controlling this error, provided the gene(s) used are stably expressed across all samples under investigation. The study described here tested 15 candidate reference genes across 44 bovine endometrial samples representing a range of physiological states and tissue subtypes. This study evaluated the suitability of both commonly used and novel experimentally derived reference genes for use in normalisation of RT-PCR data. Candidates derived via microarray analysis were superior to existing, commonly used endogenous control genes, demonstrating the suitability of using microarray data for deriving novel endogenous control genes. This study also highlighted the importance of accurate normalisation with a stable endogenous control gene, by demonstrating relative expression of a differentially expressed gene when normalised using control genes of varying stability. SUZ12 was ranked first for stability across samples as determined by the statistical algorithms used in GeNorm and Normfinder, and is therefore proposed as the best gene for normalisation of RT-PCR data in the current study.


All animal manipulations were carried out with the approval of the Ruakura Animal Ethics Committee (Hamilton, New Zealand). This work was conducted at No 5 Dairy, DairyNZ Ltd (Hamilton, New Zealand).

Sample information

Endometrial tissue samples (intercaruncular and caruncular) were obtained immediately post-mortem from 22 Holstein-Friesian dairy cows. The exact details for sample information and collection are described in Meier et al 2009 [32]. Briefly, animals had estrous cycles synchronised, with half of the animals receiving a blastocyst stage embryo at day 7 of the estrous cycle. Endometrial samples were obtained post-mortem from each animal at day 17 of the reproductive cycle. The animals consisted of 12 pregnant and 10 cycling, further divided into North American (NA) and New Zealand (NZ) genetic ancestry (Table 4). This genetic strain model was chosen because NA cows have been reported to have poorer reproductive performance than NZ cows [31].

Table 4 Tissue categories: Two different tissue types each represented by four different groups of Holstein Friesians.

RNA Extraction

Endometrial tissue was homogenised in Qiagen buffer RLT (QIAGEN GmbH, QIAGEN, Hilden, Germany) using FastPrep Lysing Matrix D tubes in a FastPrep instrument (MP Biomedicals, Solon, OH).

Total RNA was extracted using a Qiagen RNeasy kit (QIAGEN). All samples were DNase treated using the Ambion DNA-free kit (Ambion, Austin, TX) according to the manufacturer's instructions. RNA quantity was determined by spectrophotometry in a Nanodrop ND-1000 (Nanodrop Technologies, Wilmington, DE). RNA integrity was checked using the Agilent 2100 Bioanalyzer with a RNA 6000 Nano LabChip kit (Agilent Technologies, Palo Alto, CA).

cDNA Synthesis

One μg of an endometrial RNA sample was used for cDNA synthesis using the Invitrogen Superscript III Supermix kit (Invitrogen Corporation, Carlsbad, California). Total RNA was reverse transcribed according to the manufacturer's instructions using a final concentration of 27 μM of random pentadecamers primers. Briefly, RNA and random primers were mixed and denatured at 65 C for 5 minutes, followed by 1 minute on ice. Annealing buffer and Superscript/RNase was added to samples; these were then incubated for 10 minutes at 25 C (primer annealing), followed by 50 minutes at 50 C, and finally 5 minutes at 85 C to inactivate the enzyme. Reverse transcription (RT) negative controls were performed to test for the presence of genomic DNA contamination in RNA samples. Duplicate experimental samples were processed for cDNA synthesis as described above, but without the inclusion of the reverse transcriptase enzyme. Amplification was then tested for all genes using RT-PCR followed by assessment on a 3% agarose gel. No amplification was found in any of these samples.

Candidate endogenous control genes

Fifteen potential endogenous control genes were selected either from a literature search or were identified from a microarray study through statistical analysis using Genespring GX (Agilent Technologies) software in combination with the GeNorm [20]and Normfinder [21] algorithms.

Briefly, Genespring GX software was used to analyse an array data set (Agilent 44k 60-mer oligonucleotide bovine array) representing 44 bovine endometrial samples collected as described above. All array data had undergone standard quality control and statistical analysis, including filtering on flags (present or marginal in all samples) and filtering on raw expression level of 200 to obtain reliable data. A list of 27 genes were derived from this dataset by filtering on normalised expression level (upper limit 1.2 and lower limit 0.833) to obtain genes that had stable expression across all 44 samples. This dataset was further analysed using the Microsoft Excel applets GeNorm and Normfinder. These programs were used to determine the 10 most stably expressed genes in the array data set.

Quantitative Real-Time PCR

Real-time PCR using the Roche Lightcycler 480 (Roche) was performed on 15 candidate endogenous control genes for each of the 44 bovine endometrial samples using the Roche real-time PCR master mix (Lightcycler 480 Probes Master) in combination with Roche Universal Probe Library (UPL) assays. Assays were designed to publicly available bovine gene sequences (NCBI) using Roche UPL design software (ProbeFinder, v.2.45). All assays were designed to span an intron-exon boundary to prevent amplification of DNA. The primer and probe sequences are presented in Table 5.

Table 5 Characteristics of gene specific real-time PCR assays.

The PCR reaction volume was 10 μL consisting of 0.5 μM of each primer and 0.1 μM of probe. Standard cycling conditions were used [95 C for 10 minutes, (95 C for 10 seconds, 60 C for 30 seconds) × 50 cycles, 40 C for 40 seconds].

To quantify gene expression, cDNA was diluted 100-fold for all genes except ACTR1A, DNAJC17, HDAC1, RANBP10, and the target gene OXTR where cDNA was diluted 10-fold.

Each PCR run included a no-template control with water added instead of cDNA, as well as a RT negative control for each gene. Triplicate measurements were performed for all samples and standard curves. All samples for each gene were run on the same plate.

The Roche Lightcycler 480 software was used to perform quantification analysis of gene expression using the relative standard curve second derivative maximum analysis method, a non-linear regression line method. A six point relative standard curve of serial dilutions of cDNA was used with an estimated starting concentration of 1.0 and final concentration of 1.6E-03.

Relative quantification of OXTR

The Roche Lightcycler 480 Software was used to perform advanced relative quantification analysis of OXTR gene expression. SUZ12, UXT and GAPDH were each used as a reference for the quantification of OXTR expression. Relative quantification was also performed using the normalisation factor of the two most stably expressed genes identified through GeNorm analysis (SUZ12 and ZNF131). Normalisation factors were calculated by taking the geometric mean of the two genes for each sample.

Data analysis

Statistical analyses were performed using the Excel applets GeNorm (version 3.5)[20]and Normfinder [21] to estimate expression stability of the candidate endogenous control genes. GeNorm requires expression data to be input as concentrations determined via quantification, taking into consideration PCR efficiencies (Table 5). The program then estimates the most stable genes based upon pairwise comparisons of sample variability. The two most stable genes are identified and a normalisation factor calculated. Normfinder analyses the stability of the candidate genes taking into consideration inter-group variability. The program then ranks genes based on a stability value, with the lowest value indicating the most stably expressed gene. Significance tests were also performed on the data. Microsoft Excel was used to perform T-tests (with Bonferroni multiple testing correction applied) to test for significant differences between the means, and F-tests were used to assess differences in variance between experimental groups for each gene. ANOVA was also used to test for any significant difference in the means (GenSTAT). Pearson correlation calculation was used to assess the correlation of the microarray reported expression for the target gene (OXTR) and the RT-PCR reported expression.


  1. 1.

    Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochemical and Biophysical Research Communications. 2006, 345: 646-651. 10.1016/j.bbrc.2006.04.140..

    Article  CAS  PubMed  Google Scholar 

  2. 2.

    Valasek MA, Repa JJ: The power of real-time PCR. Advan Physiol Edu. 2005, 29: 151-159. 10.1152/advan.00019.2005..

    Article  Google Scholar 

  3. 3.

    Wang Y, Barbacioru C, Hyland F, Xiao W, Hunkapiller K, Blake J, Chan F, Gonzalez C, Zhang L, Samaha R: Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC Genomics. 2006, 7: 59- 10.1186/1471-2164-7-59

    PubMed Central  Article  PubMed  Google Scholar 

  4. 4.

    Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nature Protocols. 2006, 1: 1559-1582. 10.1038/nprot.2006.236

    Article  CAS  PubMed  Google Scholar 

  5. 5.

    Morey JS, Ryan JC, Van Dolah FM: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biological Procedures Online. 2006, 8: 175-193. 10.1251/bpo126

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  6. 6.

    Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR - a perspective. J Mol Endocrinol. 2005, 34: 597-601. 10.1677/jme.1.01755

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. Journal of Molecular Endocrinology. 2002, 29: 23-39. 10.1677/jme.0.0290023

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. BioTechniques. 2004, 37: 112-114. 116, 118-119

    CAS  PubMed  Google Scholar 

  9. 9.

    Kizaki K, Ushizawa K, Takahashi T, Yamada O, Todoroki J, Sato T, Ito A, Hashizume K: Gelatinase (MMP-2 and -9) expression profiles during gestation in the bovine endometrium. Reproductive Biology and Endocrinology: RB&E. 2008, 6: 66- 10.1186/1477-7827-6-66

    Article  Google Scholar 

  10. 10.

    Davies D, Meade KG, Herath S, Eckersall PD, Gonzalez D, White JO, Conlan RS, O'Farrelly C, Sheldon IM: Toll-like receptor and antimicrobial peptide expression in the bovine endometrium. Reproductive Biology and Endocrinology: RB&E. 2008, 6: 53- 10.1186/1477-7827-6-53

    Article  Google Scholar 

  11. 11.

    Ushizawa K, Takahashi T, Hosoe M, Ishiwata H, Kaneyama K, Kizaki K, Hashizume K: Global gene expression analysis and regulation of the principal genes expressed in bovine placenta in relation to the transcription factor AP-2 family. Reproductive Biology and Endocrinology. 2007, 5: 17-10.1186/1477-7827-5-17. 10.1186/1477-7827-5-17

    PubMed Central  Article  PubMed  Google Scholar 

  12. 12.

    Lee H-Y, Acosta TJ, Tanikawa M, Sakumoto R, Komiyama J, Tasaki Y, Piskula M, Skarzynski DJ, Tetsuka M, Okuda K: The role of glucocorticoid in the regulation of prostaglandin biosynthesis in non-pregnant bovine endometrium. J Endocrinol. 2007, 193: 127-135. 10.1677/joe.1.06975

    Article  CAS  PubMed  Google Scholar 

  13. 13.

    Olsvik P, Lie K, Jordal A-E, Nilsen T, Hordvik I: Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Molecular Biology. 2005, 6: 21- 10.1186/1471-2199-6-21

    PubMed Central  Article  PubMed  Google Scholar 

  14. 14.

    Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Garnsworthy PC, Sinclair KD, Webb R: Integration of Physiological Mechanisms That Influence Fertility in Dairy Cows. animal. 2008, 2: 1144-1152. 10.1017/S1751731108002358..

    Article  CAS  PubMed  Google Scholar 

  16. 16.

    Diskin MG, Murphy JJ, Sreenan JM: Embryo survival in dairy cows managed under pastoral conditions. Animal Reproduction Science. 2006, 96: 297-311. 10.1016/j.anireprosci.2006.08.008

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Takemura Y, Osuga Y, Yamauchi T, Kobayashi M, Harada M, Hirata T, Morimoto C, Hirota Y, Yoshino O, Koga K, et al: Expression of Adiponectin Receptors and Its Possible Implication in the Human Endometrium. Endocrinology. 2006, 147: 3203-3210. 10.1210/en.2005-1510

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Ushizawa K, Takahashi T, Hosoe M, Kizaki K, Abe Y, Sasada H, Sato E, Hashizume K: Gene expression profiles of novel caprine placental prolactin-related proteins similar to bovine placental prolactin-related proteins. BMC Developmental Biology. 2007, 7: 16- 10.1186/1471-213X-7-16

    PubMed Central  Article  PubMed  Google Scholar 

  19. 19.

    Bionaz M, Loor JJ: Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. Physiol Genomics. 2007, 29: 312-319. 10.1152/physiolgenomics.00223.2006

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Andersen CL, Jensen JL, Ørntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Research. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496

    Article  CAS  PubMed  Google Scholar 

  21. 21.

    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 Biology. 2002, 3: research0034.0031-research0034.0011-research0034.0031-research0034.0011

    Google Scholar 

  22. 22.

    Atkinson BA, King GJ, Amoroso EC: Development of the caruncular and intercaruncular regions in the bovine endometrium. Biol Reprod. 1984, 30: 763-774. 10.1095/biolreprod30.3.763

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    Dhaliwal GS, Murray RD, Rees EM, Howard CV, Beech DJ: Quantitative unbiased estimates of endometrial gland surface area and volume in cycling cows and heifers. Research in Veterinary Science. 2002, 73: 259-265. 10.1016/S0034-5288(02)00098-X

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Asselin E, Drolet P, Fortier MA: In Vitro Response to Oxytocin and Interferon-Tau in Bovine Endometrial Cells from Caruncular and Inter-Caruncular Areas. Biol Reprod. 1998, 59: 241-247. 10.1095/biolreprod59.2.241

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Klein C, Bauersachs S, Ulbrich SE, Einspanier R, Meyer HHD, Schmidt SEM, Reichenbach H-D, Vermehren M, Sinowatz F, Blum H, Wolf E: Monozygotic Twin Model Reveals Novel Embryo-Induced Transcriptome Changes of Bovine Endometrium in the Preattachment Period. Biol Reprod. 2006, 74: 253-264. 10.1095/biolreprod.105.046748

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    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: 225-240. 10.1530/REP-07-0415

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Balaguer SA, Pershing RA, Rodriguez-Sallaberry C, Thatcher WW, Badinga L: Effects of Bovine Somatotropin on Uterine Genes Related to the Prostaglandin Cascade in Lactating Dairy Cows. J Dairy Sci. 2005, 88: 543-552.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Gimpl G, Fahrenholz F: The Oxytocin Receptor System: Structure, Function, and Regulation. Physiol Rev. 2001, 81: 629-683.

    CAS  PubMed  Google Scholar 

  29. 29.

    Spencer TE, Sandra O, Wolf E: Genes involved in conceptus-endometrial interactions in ruminants: insights from reductionism and thoughts on holistic approaches. Reproduction. 2008, 135: 165-179. 10.1530/REP-07-0327

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Bauersachs S, Ulbrich SE, Gross K, Schmidt SEM, Meyer HHD, Wenigerkind H, Vermehren M, Sinowatz F, Blum H, Wolf E: Embryo-induced transcriptome changes in bovine endometrium reveal species-specific and common molecular markers of uterine receptivity. Reproduction. 2006, 132: 319-331. 10.1530/rep.1.00996

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Horan B, Mee JF, O'Connor P, Rath M, Dillon P: The effect of strain of Holstein-Friesian cow and feeding system on postpartum ovarian function, animal production and conception rate to first service. Theriogenology. 2005, 63: 950-971. 10.1016/j.theriogenology.2004.05.014

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Meier S, Peterson AJ, Mitchell MD, Littlejohn M, Walker CG, Roche JR: Genetic strain and reproductive status affect endometrial fatty acid concentrations. J Dairy Sci. 2009, 92: 3723-3730. 10.3168/jds.2009-2047

    Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by funding from New Zealand Dairy Farmers through DairyNZ Inc. (New Zealand; Project AN708), and the Foundation for Research, Science, and Technology (New Zealand, DRCX 0301 and DRCX0202). The authors acknowledge the statistical advice and assistance of Barbara Dow

Author information



Corresponding author

Correspondence to Caroline G Walker.

Additional information

Authors' contributions

CGW was involved in experimental design, performed the experimental work and statistical analysis, and drafted the manuscript. SM designed the animal trial, carried out sample collection, project manager and funding of the project. JRR was involved in critical analysis of the manuscript and general supervision of the project. MDM was involved in critical analysis of the manuscript and general supervision of the project. MDL was involved in experimental design and critical analysis of the manuscript, and general supervision of the project.

All Authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Walker, C.G., Meier, S., Mitchell, M.D. et al. Evaluation of real-time PCR endogenous control genes for analysis of gene expression in bovine endometrium. BMC Molecular Biol 10, 100 (2009).

Download citation


  • Endometrial Tissue
  • Endometrial Sample
  • Endogenous Control Gene
  • North American
  • Universal Probe Library