Age-related changes in relative expression stability of commonly used housekeeping genes in selected porcine tissues
© Schellander et al; licensee BioMed Central Ltd. 2011
Received: 28 June 2011
Accepted: 24 October 2011
Published: 24 October 2011
Gene expression analysis using real-time RT-PCR (qRT-PCR) is increasingly important in biological research due to the high-throughput and accuracy of qRT-PCR. For accurate and reliable gene expression analysis, normalization of gene expression data against housekeeping genes or internal control genes is required. The stability of reference genes has a tremendous effect on the results of relative quantification of gene expression by qRT-PCR. The expression stability of reference genes could vary according to tissues, age of individuals and experimental conditions. In the pig however, very little information is available on the expression stability of reference genes. The aim of this research was therefore to develop a new set of reference genes which can be used for normalization of mRNA expression data of genes expressed in varieties of porcine tissues at different ages.
The mRNA expression stability of nine commonly used reference genes (B2M, BLM, GAPDH, HPRT1, PPIA, RPL4, SDHA, TBP and YWHAZ) was determined in varieties of tissues collected from newborn, young and adult pigs. geNorm, NormFinder and BestKeeper software were used to rank the genes according to their stability. geNorm software revealed that RPL4, PPIA and YWHAZ showed high stability in newborn and adult pigs, while B2M, YWHAZ and SDHA showed high stability in young pigs. In all cases, GAPDH showed the least stability in geNorm. NormFinder revealed that TBP was the most stable gene in newborn and young pigs, while PPIA was most stable in adult pigs. Moreover, geNorm software suggested that the geometric mean of three most stable gene would be the suitable combination for accurate normalization of gene expression study.
Although, there was discrepancy in the ranking order of reference genes obtained by different analysing software methods, the geometric mean of the RPL4, PPIA and YWHAZ seems to be the most appropriate combination of housekeeping genes for accurate normalization of gene expression data in different porcine tissues at different ages.
The pig is one of the most studied organism in research community as a food as well as a model animal, and many projects in pigs require the quantification of genes for many purposes. Real-time quantitative PCR (qRT-PCR) is the most frequently used method for gene quantification nowadays. qRT-PCR is an efficient method for quantification of mRNA transcript levels due to its high sensitivity, reproducibility and large dynamic range. Furthermore, it is fast, easy to use and provides simultaneous measurement of gene expression in many different samples for a limited number of genes [1–3]. In case of qRT-PCR, when analyzing data for relative quantification, results are normalized to a reference. The most accepted approach to mRNA quantification is normalization of the expression level of a gene of interest (target gene) to the expression level of an internal stably expressed gene (control gene) [4–6]. The control gene, often termed reference gene or housekeeping gene, is a stably expressed gene that is experimentally verified in given species and tissues under given experimental conditions [3, 7–9]. Normalizing to a reference gene is a widely used method because it is simple in theory. The normalization adjusts for differences in the quality or quantity of template RNA or starting material and differences in RNA preparation and cDNA synthesis, since the reference gene is exposed to the same preparation steps as the gene of interest. This allows the direct comparison of normalized transcript expression levels between samples. However, this approach requires the selection of at least one reference gene for validation of a corresponding qRT-PCR method. Normalization is extremely important to allow accurate comparison of the results between different samples and conditions in gene expression studies . For instance, the commonly used reference genes such as GAPDH and β-actin are unfortunately often used without prior validation of their expression stability under the specific study conditions, but a number of studies have shown that the expression of those genes is significantly altered in some experimental conditions [10–12]. It is therefore necessary to validate the expression stability of reference genes prior to their use in an experimental protocol. Recently it has been recommended that a combination of reference genes should be used to obtain a more stable reference  and the use of a single reference gene is nowadays discouraged by more and more authors [4, 6, 13]. Because, a variability or alteration in the chosen reference gene by the experiment, however, may change the obtained results entirely and could be incorrect. Therefore, the validation of potential reference genes is essential.
An ideal reference gene should be stably expressed and unaffected by experimental protocol or status . But, the recent studies showed that the housekeeping gene expressions could be changed according to the type of tissues [3, 8, 15] breeds , experimental condition (such as treatment or disease) [16–19] and age [15, 20, 21]. A set of reference genes are suggested on the basis of their stability over tissues in pigs [3, 7, 15, 22, 23] but studies for expression stability of housekeeping genes in varieties of porcine tissue collected from different age of pigs are rare. Therefore, this study was aimed to explore the expressions of nine mostly used housekeeping genes in 14 different tissues collected from three different ages of pigs (1 day old piglet, 2 months old young and 5 months old adult pigs) in order to select the suitable set of housekeeping genes that could be used as an internal control to normalize gene expression in pigs.
A total of nine clinically healthy pigs of three age group were selected: neonatal (one day old), young (2 months old) and adult (5 months old) for this experiment. Each age group consisted of three animals of Pietrain, and all the animals were male and from the same batch. All pigs were kept at the Frankenforst experimental research farm at the University of Bonn (Germany). The animals were reared and slaughtered according to the rules of German performance stations . The animals were fed same diet ad libitum during the whole experimental period. Blood was collected for peripheral blood mononuclear cells (PBMC) isolation. Lymph nodes (cervical and mesenteric), intestinal mucosa from duodenum, jejunum and ileum, tissues from stomach, liver, spleen, thymus, lung, kidney, heart and skin from ear were collected for mRNA isolation after slaughter. For mRNA isolation from tissues, samples were directly put into liquid nitrogen after washing in PBS. PBMC was isolated from whole blood using Ficoll-Histopaque (Sigma) following manufacturer's protocol. All samples were kept in -80°C till used.
RNA isolation and cDNA synthesis
Total RNA was isolated from individual samples by using Tri-Reagent (Sigma-Aldrich, Munich, Germany) according to the standard protocol. In brief, sample was first grinded in a mortar, then mixed and homogenized with 1 ml Tri-Reagent using electric homogenizer. To ensure complete dissociation of nucleoprotein complexes, the sample was allowed to stand for 5 min before adding 0.2 ml of chloroform. The mixture was shaken and left at room temperature for 10 min and centrifuged at 12,000 × g for 15 min at 4°C. The upper aqueous phase was transferred to another fresh centrifuge tube and RNA was precipitated with 0.5 ml of isopropanol. After being incubated at room temperature for 10 min, the sample was centrifuged at 12,000 × g for 10 min at 4°C to get the RNA pellet, which was subsequently washed by 75% (v/v) ethanol. Centrifugation was then performed and the RNA pellet was air-dried and resuspended in 25 μl of DEPC treated water. RNA was isolated from PBMC using Picopure RNA isolation kit (Cat.# KIT0202; Arcturus). All samples were kept at -80°C until cleanup.
In order to remove possible contaminating genomic DNA, the extracted RNA was treated with 5 μl RQ1 DNase buffer, 5 units DNase and 40 units of RNase inhibitor in a 40 μl reaction volume. The mixture was incubated at 37°C for 1h followed by purification with the RNeasy Mini Kit (Qiagen, Hilden, Germany). Concentration of clean-up RNA was determined spectrophotometrically by using the NanoDrop (ND-8000) instrument; the purity of RNA was estimated by the ratio A260/A280 with respect to contaminants that absorb in the UV. Additional examination of integrity was done by denaturing agarose gel electrophoresis and ethidium bromide staining. Finally, the purified RNA was stored at -80°C for further analysis.
Approximately 1.5 μg of total RNA for each sample was transcribed into cDNA. cDNA was synthesised using GoScript (Cat.#A5000) reverse Transcription System (Promega, Germany) combined with OligoDT15 Primers, Recombinant RNasin® Ribonuclease Inhibitor and GoScript™ Reverse Transcriptase according to the manufacturer's specification and protocol. cDNA was stored at -80°C until further use.
Selection of reference genes and primer design
Selected candidate reference genes, primers, and PCR reactions efficiencies
GeneBank accession number
Primer sequence (forward/reverse)
Amplicon length (bp)
Amplification efficiency (%)
Average Ct of cDNA
Nine-fold serial dilution of plasmids DNA were prepared and used as template for the generation of the standard curve. In each run, the 96-well microtiter plate contained each cDNA sample, plasmid standards for the standard curves and no-template control. A no-template control (NTC) was included in each run for each gene to check for contamination. Quantitative real-time RT-PCR (qRT-PCR) was set up using 2 μl first-strand cDNA template, 7.4 μl deionized H2O, 0.3 μM of upstream and downstream primers and 10 μl 1× Power SYBR Green I master mix with ROX as reference dye (Bio-Rad). The thermal cycling conditions were 3 min at 95°C followed by 15 s at 95°C (40 cycles) and 1 min at 60°C. Experiments were performed using the StepOnePlus™ Real-Time PCR System (Applied Biosystems). Based on the Ct values for all dilution points in a series, a standard curve was generated using linear regression and the slope and the PCR amplification efficiency of each primer pair is calculated from the slope of a standard curve . Melting curve analysis was constructed to verify the presence of gene-specific peak and the absence of primer dimer. Agarose gel electrophoresis was performed to test for the specificity of the amplicons. To ensure repeatability of the experiments, all the reactions were executed in triplicate and the average was used for further analysis.
Determination of reference gene expression stability
The raw qRT-PCR amplification data was exported from the StepOne® software (Applied Biosystem) to Microsoft® Excel. The averages of the Ct-values for each triplicate were used for stability comparison of candidate reference genes in the NormFinder, geNorm and BestKeeper. The Proc GLM (ver9.2; SAS, SAS Institute Inc., Cary, NC, USA) analysis was performed to detect the effect of age and organs on the expression of housekeeping genes. Differences in gene expression levels between age groups within tissues were determined using t-test in SAS. P < 0.05 was considered statistically significant.
Ct values of all samples were exported to Excel, ordered for use in geNormPlus software (15 days free trial version qBasePlus; http://www.biogazelle.com) and transformed to relative quantities using the gene-specific PCR amplification efficiency . These relative quantities were then exported to geNormPlus to analyze gene expression stability . The approach of reference gene selection implemented in geNorm relies on the principle that the expression ratio of two ideal reference genes should be identical in all samples, independent of the treatment, condition, or tissue type. Increasing variations in the expression ratio between two genes correspond to lower expression stability across samples. geNorm calculates the stability using a pairwise comparison model  and determines the level of pairwise variation for each reference gene with all other reference genes as the standard deviation of the logarithmically transformed expression ratios. In this way, the reference gene expression stability measure (M value) was calculated as the average pairwise variation of a particular gene with all other control genes included in the analysis [6, 8]. Lower M values represent higher expression stabilities, whereas the least stable gene showed the highest M value generates a ranking of genes according to their M values resulting in the identification of the genes with the most stable expression in the samples under analysis. geNorm was also used to estimate the normalization factor (NF n ) by calculating the geometric mean of the expression levels of the n best reference genes . The optimisation of the number of reference genes starts with the inclusion of the two genes with the lowest M value, and continues by sequentially adding genes with increasing values of M. Thus, geNorm calculates the pairwise variation V n /Vn+1between two sequential normalization factors NF n and NFn+1containing an increasing number of reference genes . A large variation means that the added gene has a significant effect on the normalization and should preferably be included for calculation of a reliable normalization factor. Ideally, extra reference genes are included until the variation V n /Vn+1drops below a given threshold. According to the geNorm, if Vn/Vn+1 < 0.15 the inclusion of an additional reference gene is not required and the recommended number of reference genes is given by n .
NormFinder uses an ANOVA-based model . The software calculates a stability value for all candidate reference genes tested. The stability value is based on the combined estimate of intra- and inter-group expression variations of the genes studied . For each gene, the average Ct value of each triplicate reaction was converted to relative quantity data as described for geNorm, to calculate the stability value with NormFinder program . The NormFinder reference tool was applied to rank the candidate reference gene expression stability for all samples with no subgroup determination as well as with age as subgroup. A low stability value, indicating a low combined intra- and inter-group variation, indicates high expression stability .
The average Ct value of each triplicate reaction was used (without conversion to relative quantity) to analyze the stability value of studied genes via BestKeeper  which creates a pairwise correlation coefficient between each gene and the BestKeeper index (BI). This index is the geometric mean of the Ct values of all candidate reference genes. BestKeeper also calculates standard deviation (SD) of the Ct values between the whole data set. The gene with the highest coefficient of correlation with the BI indicates the highest stability .
Purity, quantity of extracted RNA and verification of amplicons
Expression levels of candidate reference genes
Identification of optimal reference genes
Expression stability of nine candidate reference gens evaluated by BestKeeper software
Irrespective of age
SD [± Ct]
CV [% Ct]
SD [± Ct]
CV [% Ct]
SD [± Ct]
CV [% Ct]
SD [± Ct]
CV [% Ct]
Determination of the optimal number of reference genes for normalization
For an exact comparison of mRNA transcription in different samples or tissues it is crucial to choose the appropriate reference gene. The optimal reference gene should be constantly transcribed in all types of cells at any time in cell cycle and differentiation. Moreover the transcription of such a gene should not be regulated by internal or external influences, at least not more than the general variation in RNA synthesis . The reference gene used for normalization of gene expression in qRT-PCR studies should also pass through the same steps of analysis as the gene to be quantified. However, such a perfect reference gene does probably not exist. Recent research has demonstrated that the expression of housekeeping genes may be altered due to differences in tissue types [3, 15, 22], breeds , ages [21, 23] and experimental condition or treatment [6, 16–19]. Therefore, it is critical to elucidate differences that may exist in housekeeping genes between younger and older adults. As an increasing volume of data continues to be published exploring mRNA expression in cases of age-depended disease, there has been a greater interest in evaluating the commonly used, widely expressed housekeeping genes for comparisons between ages. Without this information, age-dependent comparisons are very difficult to make. Therefore, it is necessary to investigate the validity and reliability of measuring the expression of various housekeeping genes in porcine tissues at different ages using qRT-PCR. To the author's knowledge, this study is the first to report that aging can influence the expression of certain housekeeping genes in pigs.
Numerous studies have been carried out in order to evaluate reference genes in specific tissues in several species. The majority of these studies are directed towards specific tissues in pigs [3, 7, 29, 35, 36]. Taken together, it is very difficult to find a 'universal' reference gene having stable expression in all cell types and tissues, and in particular to find reference genes that remain stable between samples taken at different ages under different experimental conditions. According to the NCBI-PubMed statistics , GAPDH and ACTB are the two mostly used porcine housekeeping genes. But they have been shown to vary considerably and are consequently unsuitable as reference genes for normalization of gene expression analysis in some cases [10–12]. Also the low expressed reference gene TBP is highly regulated in pigs . The first priority, however, is to identify genes with stable expression preferably across cell types since many qRT-PCR studies are performed on cDNA isolated from tissues with a mixed cell population. Presently, only few major publications describe the stability of housekeeping genes in pig and are based on limited samples of specific categories [3, 7, 29, 35, 36]. Our comprehensive set of representative tissue samples and selected housekeeping genes provide valuable recommendations for the choice of endogenous control genes for the study of gene expression patterns in normal tissues. Notably, our results coincided with the finding of Gu et al  reported that YWHAZ is one of the most stably expressed reference genes across tissues in healthy pigs. Nygard et al.  reported that RPL4, TBP and YWHAZ have the highest stability across tissues collected from healthy pigs which are in good agreement with our findings. In this study, geNorm showed that PPIA, YWHAZ and RPL4 are the most stable housekeeping genes across tissues in case of newborn piglets, adults and in irrespective of ages. Additionally, PPIA, RPL4 and YWHAZ are detected to be the most stably expressed genes across the tissues by NormFinder.
geNorm finding is contradictory to the findings of Erkens et al.  who reported that TBP is one of the most stable housekeeping gene in porcine backfat and muscle (longissimus dorsi) while SDHA is reported as an unstable gene. Kuijk et al.  reported that GAPDH and B2M are the most and least stably genes, respectively in porcine oocytes and perimplantation embryos. Although tissues are different, the finding of this study is in good agreement with Piorkowska et al  who recently reported that GAPDH is the least stable reference candidate gene in porcine adipose tissues collected from different pig breeds. The findings of this study that commonly used housekeeping genes studied are expressed differentially across porcine tissues is supported by Svobodova et al  in pigs. Moreover, Svobodova et al  found that GAPDH expression was unstable across porcine tissues which is in good agreement with our result. Svobodova et al  reported that HPRT1 has the highest stability, whereas this study found that according to the geNorm HPRT1 is moderately stable across the porcine tissues (figure 3a-3d) but unstable according to the NormFinder (figure 3e-3h). Pierzchala et al.  recently reported that HPRT1 and TBP are the most stable housekeeping genes in porcine liver and in three different muscle tissues which is partially supporting the NormFinder result as well as conflicting to the geNorm result. In this study we found that HPRT1 and TBP are moderately stable genes in geNorm analysis (figure 3a-3d), whereas TBP is a stable gene but HPRT1 is an unstable gene in NormFinder analysis (figure 3e-3h). RPL4, HPRT1 and B2M are reported as stably expressed and suitable candidate genes in intestinal tissues collected from healthy pig and from pigs with enteritis . Reportedly, GAPDH is the least stable gene while RPL27 is most stable housekeeping gene in porcine stomach tissue . However, different housekeeping genes are identified between the previous studies and our study, as the samples varied in their cell, tissue, sex and developmental stage specificities, and different catalogues of selected housekeeping genes are chosen.
According to the BestKeeper analysis software, all the studied reference candidate genes, except YWHAZ at 2 months old young pigs tissues, are less suitable. Several studies previously reported similar findings for BestKeeper [8, 29, 34] and few studies followed the BestKeeper analysis method compared to geNorm and NormFinder. It is important to note that very similar discrepancies between the different algorithms have been observed in previous studies comparing statistical analysis methods [8, 16, 29, 34, 38, 39]. However, we found that the first three most stable reference genes in most cases were consistently the same when using geNorm and NormFinder, even if they were not in the exact same ranking order. Similar findings are reported by previous studies in horse, human and plants [8, 16, 38, 40]. Such discrepancy could be explained by genes' coregulation. Indeed, coregulated genes may become highly ranked independently of their expression stabilities with geNorm software . Moreover, NormFinder takes into account variation across subgroups, thus avoiding artificial selection of coregulated genes by analyzing the expression stability of candidate genes independently from each other . However, no studies dealing with porcine reference genes stability used different analysis methods except geNorm [3, 7, 22, 23, 29].
As described above, geNorm also provides a measure for the best number of reference genes that should be used for optimal normalization. In agreement with several previous studies, we postulate that the use of more than one reference gene allows for a more accurate normalization than the use of only one reference gene [4, 6, 16, 22, 32]. Based on a cut-off point for the V value, as described by Vandesompele et al , a combination of the six most stable reference genes was calculated as being optimal for gene expression studies in different porcine tissues over ages (figure 4). However, as we described above and other studies [6, 22] recommended, the combination of the most three stable genes seems to be appropriate for accurate normalization.
This investigation found evidence that there can be variafition in the expression of commonly used housekeeping genes with populations of different ages. Due to the new influx of data suggesting alterations in mRNA expression according to ages, we feel that beside therapy uses or experimental condition, the selection of housekeeping genes based upon the age of populations used should be taken into consideration. This shows again that the choice of reference genes cannot be transposed from one study to the other without validation for the specifics of each experimental protocol. In general, we recommend using the geometric mean of RPL4, PPIA and YWHAZ to guarantee suitable normalization across the porcine tissues obtained from pigs of different ages.
List of abbreviations
quantitative real-time reverse transcriptase polymerase chain reaction
- B2M :
- BLM :
Bloom syndrome, RecQ helicase-like
- GAPDH :
glyceraldehyde 3-phosphate dehydrogenase
- HPRT1 :
hypoxanthine phosphoribosyltransferase 1
- PPIA :
peptidylprolyl isomerase A (cyclophilin A)
- RPL4 :
ribosomal protein L4
- SDHA :
succinate dehydrogenase complex subunit A flavoprotein
- TBP :
TATA box binding protein
- YWHAZ :
tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide
This project was supported by the Gene Dialog project, FUGATO Plus, BMBF, grant no: 0315130C, Germany. The authors are indebted to Prof. Dr. Florian M. W. Grundler, Molecular Phytomedicine, University of Bonn, Germany for providing with StepOnePlus Real-time PCR system (Applied Biosystem) during experiment.
- Arya M, Shergill IS, Williamson M, Gommersall L, Arya N, Patel HR: Basic principles of real-time quantitative PCR. Expert Rev Mol Diagn. 2005, 5 (2): 209-219. 10.1586/14737220.127.116.11.PubMedView ArticleGoogle Scholar
- Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nat Protoc. 2006, 1 (3): 1559-1582. 10.1038/nprot.2006.236.PubMedView ArticleGoogle Scholar
- Nygard AB, Jorgensen CB, Cirera S, Fredholm M: Selection of reference genes for gene expression studies in pig tissues using SYBR green qPCR. BMC Mol Biol. 2007, 8: 67-10.1186/1471-2199-8-67.PubMedPubMed CentralView ArticleGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6 (4): 279-284. 10.1038/sj.gene.6364190.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle 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 Biol. 2002, 3 (7): RESEARCH0034-PubMedPubMed CentralView ArticleGoogle Scholar
- Erkens T, Van Poucke M, Vandesompele J, Goossens K, Van Zeveren A, Peelman LJ: Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A. BMC Biotechnol. 2006, 6: 41-10.1186/1472-6750-6-41.PubMedPubMed CentralView ArticleGoogle Scholar
- Maroufi A, Van Bockstaele E, De Loose M: Validation of reference genes for gene expression analysis in chicory (Cichorium intybus) using quantitative real-time PCR. BMC Mol Biol. 2010, 11: 15-10.1186/1471-2199-11-15.PubMedPubMed CentralView ArticleGoogle Scholar
- Lovdal T, Lillo C: Reference gene selection for quantitative real-time PCR normalization in tomato subjected to nitrogen, cold, and light stress. Anal Biochem. 2009, 387 (2): 238-242. 10.1016/j.ab.2009.01.024.PubMedView ArticleGoogle Scholar
- Barber RD, Harmer DW, Coleman RA, Clark BJ: GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiol Genomics. 2005, 21 (3): 389-395. 10.1152/physiolgenomics.00025.2005.PubMedView ArticleGoogle Scholar
- Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, Kristiansen G, Jung K: In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC Mol Biol. 2007, 8: 47-10.1186/1471-2199-8-47.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Anal Biochem. 2002, 309 (2): 293-300. 10.1016/S0003-2697(02)00311-1.PubMedView ArticleGoogle Scholar
- Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods. 2000, 46 (1-2): 69-81. 10.1016/S0165-022X(00)00129-9.PubMedView ArticleGoogle Scholar
- Pierzchala M, Lisowski P, Urbanski P, Pareek CS, Cooper RG, Kuryl J: Evaluation based selection of housekeeping genes for studies of gene expression in the porcine muscle and liver tissues. Journal of Animal and Veterinary Advances. 2011, 10 (4): 401-405.View ArticleGoogle Scholar
- Beekman L, Tohver T, Dardari R, Leguillette R: Evaluation of suitable reference genes for gene expression studies in bronchoalveolar lavage cells from horses with inflammatory airway disease. BMC Mol Biol. 2011, 12: 5-10.1186/1471-2199-12-5.PubMedPubMed CentralView ArticleGoogle Scholar
- De Boever S, Vangestel C, De Backer P, Croubels S, Sys SU: Identification and validation of housekeeping genes as internal control for gene expression in an intravenous LPS inflammation model in chickens. Vet Immunol Immunopathol. 2008, 122 (3-4): 312-317. 10.1016/j.vetimm.2007.12.002.PubMedView ArticleGoogle Scholar
- Maccoux LJ, Clements DN, Salway F, Day PJ: Identification of new reference genes for the normalisation of canine osteoarthritic joint tissue transcripts from microarray data. BMC Mol Biol. 2007, 8: 62-10.1186/1471-2199-8-62.PubMedPubMed CentralView ArticleGoogle Scholar
- Penning LC, Vrieling HE, Brinkhof B, Riemers FM, Rothuizen J, Rutteman GR, Hazewinkel HA: A validation of 10 feline reference genes for gene expression measurements in snap-frozen tissues. Vet Immunol Immunopathol. 2007, 120 (3-4): 212-222. 10.1016/j.vetimm.2007.08.006.PubMedView ArticleGoogle Scholar
- Al-Bader MD, Al-Sarraf HA: Housekeeping gene expression during fetal brain development in the rat-validation by semi-quantitative RT-PCR. Brain Res Dev Brain Res. 2005, 156 (1): 38-45.PubMedView ArticleGoogle Scholar
- Touchberry CD, Wacker MJ, Richmond SR, Whitman SA, Godard MP: Age-related changes in relative expression of real-time PCR housekeeping genes in human skeletal muscle. J Biomol Tech. 2006, 17 (2): 157-162.PubMedPubMed CentralGoogle Scholar
- Gu Y, Li M, Zhang K, Chen L, Jiang A, Wang J, Li X: Evaluation of endogenous control genes for gene expression studies across multiple tissues and in the specific sets of fat- and muscle-type samples of the pig. Journal of Animal Breeding and Genetics. 2011, 128:Google Scholar
- Piorkowska K, Oczkowicz M, Rozycki M, Ropka-Molik K, Piestrzynska-Kajtoch A: Novel porcine housekeeping genes for real-time RT-PCR experiments normalization in adipose tissue: assessment of leptin mRNA quantity in different pig breeds. Meat Sci. 87 (3): 191-195.
- Zentralverband der Deutschen Schweineproduktion (ZDS): Richtlinie für die Stationsprüfung auf Mastleistung, Schlachtkörperwert und Fleischbeschaffenheit beim Schwein, 10.12.2003.
- Kaewmala K, Uddin MJ, Cinar MU, Grosse-Brinkhaus C, Jonas E, Tesfaye D, Phatsara C, Tholen E, Looft C, Schellander K: Association study and expression analysis of CD9 as candidate gene for boar sperm quality and fertility traits. Anim Reprod Sci. 2011, 125 (1-4): 170-179. 10.1016/j.anireprosci.2011.02.017.PubMedView ArticleGoogle Scholar
- Kayan A, Cinar MU, Uddin MJ, Phatsara C, Wimmers K, Ponsuksili S, Tesfaye D, Looft C, Juengst H, Tholen E, et al: Polymorphism and expression of the porcine Tenascin C gene associated with meat and carcass quality. Meat Sci. 2011, 89 (1): 76-83. 10.1016/j.meatsci.2011.04.001.PubMedView ArticleGoogle Scholar
- Kayan A, Uddin MJ, Cinar MU, Grosse-Brinkhaus C, Phatsara C, Wimmers K, Ponsuksili S, Tesfaye D, Looft C, Juengst H: Investigation on interferon alpha-inducible protein 6 (IFI6) gene as a candidate for meat and carcass quality in pig. Meat Sci. 88 (4): 755-760.
- Laenoi W, Uddin MJ, Cinar MU, Phatsara C, Tesfaye D, Scholz AM, Tholen E, Looft C, Mielenz M, Sauerwein H: Molecular characterization and methylation study of matrix gla protein in articular cartilage from pig with osteochondrosis. Gene. 459 (1-2): 24-31.
- Oczkowicz M, Różycki M, Piórkowska K, Piestrzyńska-Kajtoch A, Rejduch B: A new set of endogenous reference genes for gene expression studies of porcine stomach. Journal of Animal and Feed Sciences. 2010, 19: 570-576.Google Scholar
- Rozen S, Skaletsky H: Primer3 on the. Methods Mol Biol. 2000, 132: 365-386. [http://www.for/general/users/and/for/biologist/programmers]PubMedGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8 (2): R19-10.1186/gb-2007-8-2-r19.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft 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 Res. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26 (6): 509-515.PubMedView ArticleGoogle Scholar
- Stern-Straeter J, Bonaterra GA, Hormann K, Kinscherf R, Goessler UR: Identification of valid reference genes during the differentiation of human myoblasts. BMC Mol Biol. 2009, 10: 66-10.1186/1471-2199-10-66.PubMedPubMed CentralView ArticleGoogle Scholar
- Svobodova K, Bilek K, Knoll A: Verification of reference genes for relative quantification of gene expression by real-time reverse transcription PCR in the pig. J Appl Genet. 2008, 49 (3): 263-265. 10.1007/BF03195623.PubMedView ArticleGoogle Scholar
- Kuijk EW, du Puy L, van Tol HT, Haagsman HP, Colenbrander B, Roelen BA: Validation of reference genes for quantitative RT-PCR studies in porcine oocytes and preimplantation embryos. BMC Dev Biol. 2007, 7: 58-10.1186/1471-213X-7-58.PubMedPubMed CentralView ArticleGoogle Scholar
- Schroyen Martine., Verhelst R, Niewold T, Buys N: 10th Gut Day Symposium (The 10th Darmendag): 2008; Friday November 14th 2008 in the Boothzaal of the Utrecht University Library, Utrecht. 2008Google Scholar
- Cappelli K, Felicetti M, Capomaccio S, Spinsanti G, Silvestrelli M, Supplizi AV: Exercise induced stress in horses: selection of the most stable reference genes for quantitative RT-PCR normalization. BMC Mol Biol. 2008, 9: 49-10.1186/1471-2199-9-49.PubMedPubMed CentralView ArticleGoogle Scholar
- Hosseini A, Sauerwein H, Mielenz M: Putative reference genes for gene expression studies in propionate and beta-hydroxybutyrate treated bovine adipose tissue explants. J Anim Physiol Anim Nutr (Berl). 94 (5): e178-184.
- Kriegova E, Arakelyan A, Fillerova R, Zatloukal J, Mrazek F, Navratilova Z, Kolek V, du Bois RM, Petrek M: PSMB2 and RPL32 are suitable denominators to normalize gene expression profiles in bronchoalveolar cells. BMC Mol Biol. 2008, 9: 69-10.1186/1471-2199-9-69.PubMedPubMed CentralView ArticleGoogle 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.