Selection of a set of reliable reference genes for quantitative real-time PCR in normal equine skin and in equine sarcoids
© Bogaert et al; licensee BioMed Central Ltd. 2006
Received: 19 December 2005
Accepted: 27 April 2006
Published: 27 April 2006
Real-time quantitative PCR can be a very powerful and accurate technique to examine gene transcription patterns in different biological conditions. One of the critical steps in comparing transcription profiles is accurate normalisation. In most of the studies published on real-time PCR in horses, normalisation occurred against only one reference gene, usually GAPDH or ACTB, without validation of its expression stability. This might result in unreliable conclusions, because it has been demonstrated that the expression levels of so called "housekeeping genes" may vary considerably in different tissues, cell types or disease stages, particularly in clinical samples associated with malignant disease. The goal of this study was to establish a reliable set of reference genes for studies concerning normal equine skin and equine sarcoids, which are the most common skin tumour in horses.
In the present study the gene transcription levels of 6 commonly used reference genes (ACTB, B2M, HPRT1, UBB, TUBA1 and RPL32) were determined in normal equine skin and in equine sarcoids. After applying the geNorm applet to this set of genes, TUBA1, ACTB and UBB were found to be most stable in normal skin and B2M, ACTB and UBB in equine sarcoids.
Based on these results, TUBA1, ACTB and UBB, respectively B2M, ACTB and UBB can be proposed as reference gene panels for accurate normalisation of quantitative data for normal equine skin, respectively equine sarcoids. When normal skin and equine sarcoids are compared, the use of the geometric mean of UBB, ACTB and B2M can be recommended as a reliable and accurate normalisation factor.
Gene expression analysis has become increasingly important in biological research where e.g. gene expression profiles from tissues associated with diseases and disorders have to be compared with each other and with those from normal tissues. One of the most powerful tools in this area is real-time quantitative reverse transcription PCR (qRT-PCR). To account for differences in starting material, RNA preparation, RNA quality and cDNA synthesis, adequate normalisation is frequently performed by comparing expression profiles of the genes of interest with those of constitutively expressed genes (= reference genes). Housekeeping genes are most widely used as reference genes, based on the assumption that they are constitutively expressed in most tissues and under certain circumstances, and that they are more or less resistant to cell cycle fluctuations [1, 2]. However, it has been demonstrated that the expression levels of these genes may vary considerably in different tissues, different cell types and different disease stages, particularly in clinical samples associated with malignant disease [3, 4]. Housekeeping genes are not only involved in the basal cell metabolism, but appear to participate in other functions too, and therefore are prone to regulation [5–7]. Especially in tumours, the metabolism is generally elevated because of permanent proliferation and expansion. Moreover, some housekeeping genes may have a specific function essential for the tumour metabolism and therefore be up or down regulated . Because of these findings, Vandesompele et al  proposed to identify a set of stable housekeeping genes in the tissue of interest and use them as internal reference genes for accurate normalisation.
Up till now, only a few gene expression studies using real-time qRT-PCR have been performed in horses. GAPDH or ACTB were commonly used as a single non-validated reference gene [10–15], since at that time not much information was available concerning this issue. Recently, Waguespack et al  compared 4 housekeeping genes (ACTB, B2M, GAPDH and TBP) in the lamellae of the hoof in horses.
In this study, 6 commonly used reference genes in both human and animal studies were investigated, both in normal skin and in equine sarcoids of horses. Equine sarcoids are fibroblastic skin tumours and are the most common tumours in horses. The disease not only induces esthetical defects, but also diminishes the economical value of affected horses . Moreover affected horses show a genetic predisposition for the development of equine sarcoids, through which the breeding value of an animal with sarcoids sharply declines [18–20]. The bovine papillomavirus (BPV) plays an important role in the aetiology of equine sarcoids [21–23]. Several clinical types exist, ranging from small, stable patches to large, aggressive and fast growing tumours . To be able to examine the gene expression profile of BPV in these different clinical types and to compare equine sarcoids with normal skin asymptomatically infected with BPV, a well suited internal control should first of all be established.
Results and discussion
Transcription profiling of the candidate genes
cDNA was synthesised from DNA-free RNA (checked with minus RT control) isolated from 8 normal equine skin and 8 equine sarcoid samples. A real-time PCR assay, based on SYBR® Green detection, was designed for the transcription profiling of six frequently used reference genes (ACTB, B2M, HPRT1, UBB, TUBA1 and RPL32 )in these cDNA samples. During optimisation of the protocol, real-time PCR products were visualised by gelelectrophoresis and sequenced for verification. For every assay, a single amplicon with the expected size was generated without primer dimer formation. Indeed, the formation of primer dimers and unspecific amplification, which can falsely increase the gene expression levels, is a major point of attention, particularly when using intercalating dyes such as SYBR® Green.
Amplicon sequences of ACTB, B2M, HPRT1, and TUBA1 were 100% identical with the described sequences on which primer design was based. The sequencing of UBB revealed 1 gap and 3 SNPs compared to the original sequence (98% identity). The gap is probably due to sequencing errors; the SNPs did not result in an amino acid variation. When translated to amino acid sequence, the RPL32 sequence was 100% identical to publicly accessible horse ESTs and human RPL32 amino acid sequences.
After optimisation, gene-specific amplification was confirmed by a single peak in melt-curve analysis. For each assay, a standard curve was generated by using 10-fold serial dilutions of pooled cDNA, generated of both normal skin and equine sarcoid tissue, characterised by a linear correlation coefficient (R2) varying from 0.991 to 0.998 and a PCR efficiency between 88.1 and 104.6%. These findings showed that these assays are suitable for quantitative purposes.
In order to select a reliable set of reference genes, each assay was performed in duplicate and included the appropriate control samples. To compare the transcription level of the selected genes across the different samples, the Ct values, ranging from 16.4 to 31.4 were converted into raw data based on the PCR efficiency, gathered by standard curve analyses.
Ranking of the reference genes. The reference genes are ranked in order of their expression stability in normal equine skin, in equine sarcoids and in the combination of both kinds of samples, decreasing from top to bottom. The reference genes chosen to calculate the normalisation factor, used for comparing equine sarcoids and normal skin, are printed in bold.
When calculating a normalisation factor (NF), a careful choice of the number of reference genes should be made. The more genes included, the more accurate the NF is. However, including too many genes may increase the risk of using unsuitable genes, and is also impractical. On the other hand, if the cut off is made too stringent, stably expressed reference genes may be excluded and accuracy might drop. In order to determine how many reference genes should be included, normalisation factors (NFn), based on the geometric mean of the expression levels of the n best reference genes, were calculated by inclusion of an extra, less stable, reference gene according to Vandesompele et al . Figures 1A(b), 1B(b) and 1C(b) show the pairwise variation Vn/Vn+1 between 2 sequential normalisation factors NFn and NFn+1 for normal skin, equine sarcoids and the combination of both kinds of samples. In all 3 cases, the inclusion of a 4th gene had no significant contribution (low V3/4 value) to the NF. The 3 member sets as described above are a good choice for the calculation of the NF.
Implementation of results in clinical research
For the reasons discussed above, we have confidence that our gene expression results are accurate and reliable. The described set of reference genes can be used in gene expression studies both in normal skin, equine sarcoids and the combination of both. One of the points of interest in veterinary medicine is the expression level of BPV in different kinds of equine sarcoids. With the normalisation technique described in this study, reliable results can be obtained. Another research topic in this domain is the study of BPV expression in normal skin, showing latent infection with BPV. Also, the expression level of specific horse genes with a putative role in tumourigenesis, can be investigated.
In conclusion, a method for genomic DNA-free RNA extraction from normal equine skin and equine sarcoids was optimised and a reference gene assay for reliable normalisation of real-time qPCR data, obtained from normal skin, equine sarcoids and the combination of both, was designed. The profiling of the gene expression pattern of 6 putative reference genes showed that 3 reference genes should be used. ACTB, TUBA1 and UBB can be used in normal skin, while ACTB, B2M and UBB are the best choice in equine sarcoids. If normal skin and equine sarcoids have to be compared, the same member set as proposed for equine sarcoids can be used.
Eight equine sarcoid samples were obtained from surgically treated horses. A whole range of tumours were sampled (one occult type, one verrucous type, one nodular type and five fibroblastic types). Three sarcoids were located on the medial part of the thigh, two in the axilla and two on the ventral side of the abdomen. Care was taken to obtain only tumoural tissue, without underlying normal stroma. Normal skin samples were obtained from healthy horses undergoing elective surgery or euthanasia (umbilical hernia, castration, osteochondrosis dissecans, dorsal displacement of the soft palatinum). Six of the samples were collected from the ventral abdomen, one from the throat region and one from the shoulder. All samples were freshly collected and stored immediately in RNAlater (Ambion). After overnight incubation at +4°C, the samples were frozen at -18°C until RNA extraction.
RNA extraction and cDNA synthesis
Total RNA was isolated from the samples using TRIR (ABgene) according to the manufacturer's instructions. Subsequently, approximately 4 μg of the total RNA solution, measured with the BioPhotometer (Eppendorf), were treated with 3 units DNase I (Ambion) to remove genomic DNA. This was followed by a spin-column purification (Microcon YM-10, Millipore). A minus RT control with primers for GAPDH was performed to check for successful removal of all the contaminating DNA. These primers were designed by the Primer 3 software  using publicly available sequences from the Nucleotide Sequence Database from the National Center for Biotechnology Information (NCBI) . The initial denaturation was performed at 94°C for 10 minutes. Thirty-five cycles of amplification were performed. Each cycle involved a denaturation step of 30 seconds at 94°C, followed by 30 seconds primer annealing at 61°C and 60 seconds primer extension at 72°C. After the last cycle, the PCR-mix was heated during 10 minutes at 72°C to become extension of the partially elongated primers.
First strand cDNA synthesis was carried out on approximately 2.6 μg of the total RNA solution with Superscript™ II Reverse Transcriptase (Invitrogen), an engineered version of M-MLV RT, and a combination of random primers (Invitrogen) and oligo(dT)20 primers (Invitrogen) in a total volume of 20 μl, following the manufacturer's instructions. After this step, a PCR was performed with primers for GAPDH to check for the presence of cDNA. The reaction conditions used in this PCR were identical to the PCR described for the minus RT control.
Reference gene selection and primer design
Name and function of genes mentioned in the text
Cytoskeletal structural protein
Beta-chain of major histocompatibility complex class I molecules
hypoxanthine phosphoribosyltransferase 1
Purine synthesis in salvage pathway
tubulin, alpha 1 (testis specific)
ribosomal protein L32
Member of the 80 different ribosome proteins
TATA box binding protein
General RNA polymerase II transcription factor
Information on the primers used in this study
Genbank accession number
Product size (bp)
Real-time quantitative PCR
Eight normal skin samples and eight equine sarcoid samples were used for quantification of reference genes.
PCR reactions were performed in a 15 μl reaction volume on the iCycler iQ Real-Time PCR Detection System (Bio-Rad) using the Platinum® SYBR® Green qPCR SuperMix UDG (Invitrogen) supplemented with 0.02 μM fluorescein and 20 ng of cDNA. Primer concentration varied according to the primers used. A blank was incorporated in each assay.
First, an UDG-treatment was done at 50°C to prevent cross contamination. The initial denaturation was performed at 95°C for 2 minutes to activate the Taq DNA polymerase, followed by 40 cycles of denaturation at 95°C for 20 seconds and a combined primer annealing/extension at the specific annealing temperature (Table 3) for 40 seconds during which fluorescence was measured. A melt curve was generated to confirm a single gene-specific peak and to detect primer dimer formation by heating the samples from 70 to 95°C in 0.5°C increments with a dwell time at each temperature of 10 seconds while continuously monitoring the fluorescence.
PCR efficiencies were calculated using a relative standard curve derived from a pooled cDNA mixture (a ten-fold dilution series with four measuring points). The pooled cDNA was obtained from normal equine skin and equine sarcoids, using the same RNA isolation and cDNA synthesis protocols as described above.
Each reaction was run in duplicate, whereby a no-template control was included.
During optimisation of the protocol, the PCR products were loaded on a 3% agarose gel after each run to confirm specific gene amplification and the absence of primer dimer formation.
Determination of reference gene expression stability
To determine the stability of the selected reference genes, the geNorm Visual Basic application for Microsoft Excel was used as described by Vandesompele et al . This approach relies on the principle that the expression ratio of two perfect reference genes should be identical in all samples, independent of the experimental condition or cell type.
This research was funded by the Special Research Fund, Ghent University, grant no. 01D01103.
- Radonic A, Thulke S, Bae HG, Müller MA, Siegert W, Nitsche A: Reference gene selection for quantitative real-time PCR analysis in virus infected cells: SARS corona virus, Yellow fever virus, Human Herpesvirus-6 and Cytomegalovirus infections. Virol J. 2005, 2: 7-10.1186/1743-422X-2-7.View ArticleGoogle Scholar
- Zhang X, Ding L, Sandford AJ: Selection of reference genes for gene expression studies in human neutrophils by real-time PCR. BMC Mol Biol. 2005, 6: 4-10.1186/1471-2199-6-4.View ArticleGoogle Scholar
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Gristar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7.View ArticleGoogle Scholar
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29: 332-337.Google Scholar
- Petersen BH, Rapaport F, Henry DP, Huseman C, Moore DP: Effect of treatment with biosynthetic human growth hormone (GH) on peripheral blood lymphocyte populations and function in growth hormone-deficient children. J Clin Endocrinol Metab. 1990, 70: 1756-1760.View ArticleGoogle Scholar
- Singh R, Green MR: Sequence specific binding of transfer RNA by glyceraldehydes-3-phosphate dehydrogenase. Science. 1993, 259: 365-368.View ArticleGoogle Scholar
- Ishitani R, Sunaga K, Hirano A, Saunders P, Katsube N, Chang DM: Evidence that glyceraldehydes-3-phosphate dehydrogenase is involved in age-induced apoptosis in mature cerebellar neurons in culture. J Neurochem. 1996, 66: 928-935.View ArticleGoogle Scholar
- Neuvians TP, Gashaw I, Sauer CG, von Ostau C, Kliesch S, Bergmann M, Häcker A, Grobholz R: Standardization strategy for quantitative PCR in human seminoma and normal testis. J Biotechnol. 2005, 117: 163-171. 10.1016/j.jbiotec.2005.01.011.View 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 by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: 34-10.1186/gb-2002-3-7-research0034. [http://medgen.ugent.be/~jvdesomp/genorm/]View ArticleGoogle Scholar
- Giguère S, Prescott JF: Quantitation of equine cytokine mRNA expression by reverse transcription-competitive polymerase chain reaction. Vet Immunol Immunop. 1999, 67: 1-15. 10.1016/S0165-2427(98)00212-8.View ArticleGoogle Scholar
- Leutenegger CM, von Rechenberg B, Huder JB, Zlinsky K, Mislin C, Akens MK, Auer J, Lutz H: Quantitative real-time PCR for equine cytokine mRNA in nondecalcified bone tissue embedded in methyl methacrylate. Calcif Tissue Int. 1999, 65: 378-383. 10.1007/s002239900717.View ArticleGoogle Scholar
- Swiderski C, Klei TR, Horohov DW: Quantitative measurement of equine cytokine mRNA expression by polymerase chain reaction using target-specific standard curves. J Immunol Methods. 1999, 222: 155-169. 10.1016/S0022-1759(98)00193-8.View ArticleGoogle Scholar
- von Rechenberg B, Leutenegger C, Zlinsky K, McIlwraith CW, Akens MK, Auer JA: Upregulation of mRNA of interleukin-1 and -6 in subchondral cystic lesions of four horses. Equine Vet J. 2001, 33: 143-149.View ArticleGoogle Scholar
- Ainsworth DM, Appleton JA, Eicker SW, Luce R, Flaminio MJ, Antzack DF: The effect of strenuous exercise on mRNA concentrations of interleukin-12, interferon-gamma and interleukin-4 in equine pulmonary and peripheral blood mononuclear cells. Vet Immunol Immunop. 2003, 91: 61-71. 10.1016/S0165-2427(02)00274-X.View ArticleGoogle Scholar
- Lim WS, Edwards JF, Boyd NK, Payne SL, Ball JM: Simultaneous quantitation of equine cytokine mRNAs using a multi-probe ribonuclease protection assay. Vet Immunol Immunop. 2003, 91: 45-51. 10.1016/S0165-2427(02)00263-5.View ArticleGoogle Scholar
- Waguespack RW, Kemppainen RJ, Cochran A, Lin HC, Belknap JK: Increased expression of MAIL, a cytokine-associated nuclear protein, in the prodromal stage of black walnut-induced laminitis. Equine Vet J. 2004, 36: 285-291.View ArticleGoogle Scholar
- Gerber H: The genetic basis of some equine diseases. Equine Vet J. 1989, 21: 244-248.View ArticleGoogle Scholar
- Lazary S, Gerber H, Glatt PA, Straub R: Equine leucocyte antigens in sarcoid-affected horses. Equine Vet J. 1985, 17: 283-286.View ArticleGoogle Scholar
- Meredith D, Elser AH, Wolf B, Soma LR, Donawick WJ, Lazary F: Equine Leucocyte Antigens: relationships with sarcoid tumors and laminitis in two pure breeds. Immunogenetics. 1986, 23: 221-225. 10.1007/BF00373016.View ArticleGoogle Scholar
- Gerber H, Dubath M-L, Lazary S: Association between predisposition to equine sarcoid and MHC in multiple-case families. Proceedings of the 5th International Conference on Equine Infectious Diseases. Edited by: Powell DG. 1988, The University Press of Kentucky, 272-277.Google Scholar
- Otten N, von Tscharner C, Lazary S, Antczak DF, Gerber H: DNA of bovine papillomavirus type 1 and 2 in equine sarcoids: PCR detection and direct sequencing. Arch Virol. 1993, 132: 121-131. 10.1007/BF01309847.View ArticleGoogle Scholar
- Teifke JP, Hardt M, Weiss E: Detection of bovine papillomavirus DNA in formalin-fixed and paraffin-embedded equine sarcoids by polymerase chain reaction and non-radioactive in situ hybridization. European Journal of Veterinary Pathology. 1994, 1: 5-10.Google Scholar
- Nasir L, Reid SWJ: Bovine papillomaviral gene expression in equine sarcoid tumours. Virus Res. 1999, 61: 171-175. 10.1016/S0168-1702(99)00022-2.View ArticleGoogle Scholar
- Pascoe RR, Knottenbelt DC: Neoplastic conditions. Manual of equine dermatology. Edited by: Pascoe RR, Knottenbelt DC. 1999, London: Saunders, 244-252.Google Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386. [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi]Google Scholar
- National Center for Biotechnology Information. [http://www3.ncbi.nlm.nih.gov/]
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595. [http://www.bioinfo.rpi.edu/applications/mfold/old/dna/form1.cgi]View ArticleGoogle Scholar
- Wain HM, Bruford EA, Lovering RC, Lush MJ, Wright MW, Povey S: Guidelines for Human Gene Nomenclature. Genomics. 2002, 79: 464-470. 10.1006/geno.2002.6748. [http://www.gene.ucl.ac.uk/cgi-bin/nomenclature/searchgenes.pl]View 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.