Evaluation of real-time PCR endogenous control genes for analysis of gene expression in bovine endometrium
© Walker et al; licensee BioMed Central Ltd. 2009
Received: 27 May 2009
Accepted: 01 November 2009
Published: 01 November 2009
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 [3–5]. 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 . The most widely used endogenous control gene in studies of endometrial gene expression is GAPDH [9–12]. 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 , making it a questionable choice for RT-PCR normalisation .
Over the past three decades, genetic selection for milk production has resulted in a significant decline in dairy cattle fertility . 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 , 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 [17–19], 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  and Normfinder  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 http://blast.ncbi.nlm.nih.gov/Blast.cgi.
Quantification cycle (Cq) values and statistics for 15 candidate endogenous control genes assayed across 44 bovine endometrial samples.
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 [22–24]. 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[25–27]. 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  and Normfinder  algorithms.
Stability ranking of candidate endogenous control genes from Normfinder and GeNorm analyses.
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.
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 ) 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).
Tissue categories: Two different tissue types each represented by four different groups of Holstein Friesians.
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).
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 and Normfinder  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
Characteristics of gene specific real-time PCR assays.
Experimentally derived genes
Amplicon size (bp)
DnaJ (Hsp40) homolog, subfamily C, member 17
histone deacetylase 1
RAN binding protein 10
CCR4-NOT transcription complex, subunit 7
ARP1 actin-related protein 1 homolog A, centractin alpha (yeast)
Bardet-Biedl syndrome 2
suppressor of zeste 12 homolog (Drosophila)
zinc finger protein 131
solute carrier family 30 (zinc transporter), member
LOC506268 similar to Uncharacterized protein C2orf29
Literature derived genes
ribosomal protein S9
peptidylprolyl isomerase A (cyclophilin A)
ribosomal protein S15a
Gene of interest
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.
Statistical analyses were performed using the Excel applets GeNorm (version 3.5)and Normfinder  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.
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
- 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..View ArticlePubMedGoogle Scholar
- Valasek MA, Repa JJ: The power of real-time PCR. Advan Physiol Edu. 2005, 29: 151-159. 10.1152/advan.00019.2005..View ArticleGoogle Scholar
- 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-59PubMed CentralView ArticlePubMedGoogle Scholar
- Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nature Protocols. 2006, 1: 1559-1582. 10.1038/nprot.2006.236View ArticlePubMedGoogle Scholar
- 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/bpo126PubMed CentralView ArticlePubMedGoogle Scholar
- 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.01755View ArticlePubMedGoogle Scholar
- 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.0290023View ArticlePubMedGoogle Scholar
- 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-119PubMedGoogle Scholar
- 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-66View ArticleGoogle Scholar
- 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-53View ArticleGoogle Scholar
- 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-17PubMed CentralView ArticlePubMedGoogle Scholar
- 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.06975View ArticlePubMedGoogle Scholar
- 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-21PubMed CentralView ArticlePubMedGoogle Scholar
- 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.6364190View ArticlePubMedGoogle Scholar
- Garnsworthy PC, Sinclair KD, Webb R: Integration of Physiological Mechanisms That Influence Fertility in Dairy Cows. animal. 2008, 2: 1144-1152. 10.1017/S1751731108002358..View ArticlePubMedGoogle Scholar
- 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.008View ArticlePubMedGoogle Scholar
- 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-1510View ArticlePubMedGoogle Scholar
- 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-16PubMed CentralView ArticlePubMedGoogle Scholar
- 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.2006View ArticlePubMedGoogle Scholar
- 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-0496View ArticlePubMedGoogle Scholar
- 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.0011Google Scholar
- 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.763View ArticlePubMedGoogle Scholar
- 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-XView ArticlePubMedGoogle Scholar
- 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.241View ArticlePubMedGoogle Scholar
- 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.046748View ArticlePubMedGoogle Scholar
- 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-0415View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gimpl G, Fahrenholz F: The Oxytocin Receptor System: Structure, Function, and Regulation. Physiol Rev. 2001, 81: 629-683.PubMedGoogle Scholar
- 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-0327View ArticlePubMedGoogle Scholar
- 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.00996View ArticlePubMedGoogle Scholar
- 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.014View ArticlePubMedGoogle Scholar
- 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-2047View ArticlePubMedGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.