Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Genomic selection of reference genes for real-time PCR in human myocardium

BMC Medical Genomics volume 1, Article number: 64 (2008) Cite this article

Abstract

Background

Reliability of real-time PCR (RT-qPCR) data is dependent on the use of appropriate reference gene(s) for normalization. To date, no validated reference genes have been reported for normalizing gene expression in human myocardium. This study aimed to identify validated reference genes for use in gene expression studies of failed and non-failed human myocardium.

Methods

Bioinformatic analysis of published human heart gene expression arrays (195 failed hearts, 16 donor hearts) was used to identify 10 stable and abundant genes for further testing. The expression stability of these genes was investigated in 28 failed and 28 non-failed human myocardium samples by RT-qPCR using geNorm software.

Results

Signal recognition particle 14 kDa (SRP14), tumor protein, translationally-controlled 1 (TPT1) and eukaryotic elongation factor 1A1 (EEF1A1) were ranked the most stable genes. The commonly used reference gene, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was ranked the least stable of the genes tested. The normalization strategy was tested by comparing RT-qPCR data of both normalized and raw expression levels of brain natriuretic peptide precursor (NPPB), a gene known to be up-regulated in heart failure. Non-normalized levels of NPPB exhibited a marginally significant difference between failed and non-failed samples (p = 0.058). In contrast, normalized NPPB expression levels were significantly higher in heart-failed patients compared with controls (p = 0.023).

Conclusion

This study used publicly available gene array data to identify a strategy for normalization involving two reference genes in combination that may have broad application for accurate and reliable normalization of RT-qPCR data in failed and non-failed human myocardium.

Peer Review reports

Background

Analysis of gene expression levels has become increasingly important in biological research, providing insight into the complex regulatory networks that underlie health and disease [1]. Quantitative real-time PCR (RT-qPCR) is one of the most sensitive, specific, and widely-used methods for measuring the expression levels of selected genes [2, 3]. Accuracy of RT-qPCR data is reliant on appropriate normalization of individual samples, and sub-optimal normalization is a common pitfall of RT-qPCR analyses [1, 2, 4]. There are many sources of variation in gene expression measurements, including differences in cell number, transcriptional efficiency, mRNA integrity, and differences in efficiency of RNA extraction and variability of cDNA synthesis between samples. While a number of strategies have been applied to account for inter-sample variation, the use of two or more reference genes is currently the favored approach for accurate and reliable normalization of the mRNA fraction [2, 5].

Among the most commonly used reference genes in RT-qPCR are glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin [6]. These genes were introduced for normalization when mRNA expression was assessed using semi-quantitative methods, such as Northern blotting and RNase protection, because of their relatively high abundance in numerous cell-types [1, 4]. However, the increased sensitivity of RT-qPCR demands a significantly greater stringency for reference gene selection [2, 5]. Previous studies have demonstrated that the expression levels of classical reference genes can vary markedly across cells, tissues and experimental conditions, including heart cells and tissues [7, 8]. More stable alternatives exist, although these tend to be specific to the tissue under investigation [1, 6, 9, 10]. To date, no validated universal or human heart-specific reference genes have been identified.

To select optimal reference genes for RT-qPCR analysis of human myocardium, we screened human heart gene expression data from a published microarray study and investigated gene abundance and expression stability. We validated the top-ranked candidate reference genes in an independent set of human myocardium samples obtained from heart transplant patients and heart-healthy donors, and identified the minimum set required for accurate and reliable normalization of RT-qPCR data in failed and non-failed human myocardium.

Methods

Candidate gene selection

To select candidate reference genes for subsequent testing by RT-qPCR, publicly available microarray data from explanted human myocardium was screened to identify transcripts with stable and abundant expression. Abundant transcripts were selected to ensure robust measurement by RT-qPCR in all samples. Affymetrix gene expression profiles generated from left ventricle myocardium from 195 heart transplant recipients with advanced ischemic or idiopathic cardiomyopathy and 16 unmatched heart donors [11] were obtained from the NCBI GEO database (GEO accession GSE5406, http://www.ncbi.nlm.nih.gov). These data were derived from Affymetrix HG-U133A microarrays, each consisting of 22,283 probe sets representing approximately 13,000 genes. MAS 5.0 [12] normalized data were used to calculate the mean and standard deviation of the expression intensity for each probe set across all samples, using the R software package [13]. The 5% of probe sets with the smallest standard deviations (i.e., the least variable) were then identified, and of those, any probe set with a mean expression intensity of greater than 10 on the log2 scale was selected for assessment as a candidate reference for RT-qPCR analysis. For comparison, the expression stability and abundance of 50 genes previously used or proposed by prior studies as reference genes [6, 10, 14] were also screened for abundance and stability across all probes for identification of candidate reference genes.

Patient samples

Tissue from the left ventricular free wall of the myocardium was obtained from explanted failed hearts from heart transplant recipients at the Cleveland Clinic Foundation (n = 28). Non-failed heart tissue was obtained from unmatched organ donors through Lifebanc of Northeast Ohio (n = 28). Tissue handling was as previously described [15, 16]. The investigation was approved by the Cleveland Clinic Internal Review Board (ethics approval IRB 2378) and all patients provided informed consent. The study adhered to the principles outlined in the Declaration of Helsinki and Title 45, US Code of Federal Regulation, Part 46, Protection of Human Subjects, revised November 13, 2001, effective December 13, 2001. All procedures followed were in accordance with institutional guidelines.

Sample preparation

Frozen tissue blocks (previously stored at -80°C, mean weight ± standard deviation = 181 ± 81 mg) were placed in pre-chilled tubes containing TRIzol® (Invitrogen, Carlsbad, CA) and immediately subjected to automated grinding in a Mixer Mill MM301 (Retsch, Haan, Germany). Total RNA was isolated with chloroform and purified using RNeasy Midi columns in a total volume of 300 μL (Qiagen, Valencia, CA) according to the manufacturer's instructions. RNA yield and purity was determined using a Nanodrop spectrophotometer (Nanodrop Technologies, Montanin, DE). The mean RNA concentration and 260:280 ratio (± standard deviation) was 749 ± 271 ng/μL and 2.05 ± 0.03, respectively. The integrity of the RNA was assessed with gel electrophoresis. All samples were digested with DNase I (Invitrogen, Carlsbad, CA). First strand cDNA synthesis was performed from 2 μg of total RNA with oligodT primers and Superscript III, followed by RNase H digestion (Invitrogen, Carlsbad, CA) as previously described [17].

Real-time PCR

Primers for all 10 candidate reference genes were designed using Primer3 software [18] within regions spanned by the Affymetrix probe set target sequences (Table 1). At least one primer in every primer pair was designed across an exon/exon boundary to prevent amplification of any dsDNA transcripts not removed by digestion with DNase I. Affymetrix probe sequences were obtained from the Affymetrix NETAFFX Analysis Centre http://www.affymetrix.com/analysis/index.affx. Alignments of Affymetrix probe set target sequences with mRNA and genomic DNA sequences for each gene were performed with Genious version 2.0.01 software. RT-qPCR reactions were performed in duplicate in a Rotor-Gene 3000 (Corbett Research, Sydney, Australia) using SYBR Green I detection of dsDNA synthesis. Reactions (20 μL) contained 1 μL template cDNA, 1× PCR buffer, 1.5 mM MgCl2, 0.2 mM dNTPs (Fermentas, Glen Burnie, MD), 0.5 mM forward and reverse primers, 5× SYBR Green I (Roche Diagnostics, Mannheim, Germany) and 1 U Taq-Ti DNA polymerase (Fisher Biotec, West Perth, Australia). The cycling conditions comprised 2 mins polymerase activation at 94°C followed by 30 cycles of 94°C for 10s, 56°C for 20s, and 68°C for 10s. On completion, amplimers were exposed to a temperature gradient from 79–95°C (melt curve) to confirm that a single product had been amplified. For each sample the Ct value (the fluorescent point at which the reactions are compared) was fitted to a standard curve consisting of five serial dilution points (in triplicate) of purified DNA template (amplicon derived, copy number ranging from ~109 to ~105 copies) and a no-template control. The reaction efficiency for each primer set is detailed in Table 1. The mean standard deviation of Ct values for duplicated samples was 0.10. Quantification of each sample was performed using Rotor-Gene software version 6.1.

Table 1 Real-time PCR primers
Full size table

NPPB expression levels were determined using a Taqman gene expression assay with inventoried probes (assay id # Hs00173590_m1, Applied Biosystems, Foster City, CA). Reactions (20 μL) were performed in triplicate on a 7500 Fast real-time PCR system (Applied Biosystems) in standard mode according to manufacturer's instructions. Samples were quantified using a standard curve consisting of five serial dilution points (in triplicate) of purified DNA template (amplicon derived, copy number ranging from 1.13 × 108 – 1.13 × 104 copies) and a no-template control. Reaction efficiency was 101.2%. The mean standard deviation of Ct values for duplicate samples was 0.16. Quantification was performed with 7500 Software version 2.0 (Applied Biosystems).

Evaluation of candidate reference genes using GeNorm

The Visual Basic Application geNorm version 3.4 [5, 14] was used to evaluate the expression stability of each gene, determine the ideal number of genes required for normalization, and calculate individual normalization factors based on the expression levels of the best reference genes. In geNorm analysis, the calculated gene stability measure (M) relies on the principle that the ratio of two ideal reference genes will be identical in all samples. M is calculated as the average pairwise variation between a particular gene and all other genes in the analysis. Decreasing expression stability occurs when one or more of the genes is not constantly expressed. At each step in the analysis the least stable gene (highest M-value) is excluded and M is recalculated. This process continues in a step-wise manner until the two most stable genes remain. To determine the optimum number of genes for normalization, the two most stable genes are used as a starting point. Genes are added sequentially (from most stable to least stable) until the pairwise variation between two sequential normalization factors drops below a set threshold, indicating that the prior gene set is sufficient for accurate normalization [5]. A large variation indicates that the added gene has a significant effect and should preferably be included, provided it is stably expressed. However, genes that are less stably expressed also give high pairwise variation because they have a different pattern of expression compared with the more stably expressed genes. Therefore, to achieve accurate normalization, less stable genes should not be included in calculation of the normalization factor, despite their high pairwise variation values. Based on their findings Vandesompele et. al. recommended setting the pairwise variation threshold at 0.15 [5].

Using the geNorm algorithm, normalization factors were calculated by transforming Ct values into quantities (using standard curves or the comparative Ct method), setting the highest relative quantity for each gene to 1, and scaling the expression value for all other samples to a proportion relative to this highest value. Normalization factors using multiple genes were determined by calculating the geometric mean of the scaled reference genes.

Statistical analyses

Univariate analyses were performed using analysis of variance to compare gene expression levels in heart-failed and non heart-failed groups. Gene expression data displayed consistently skewed distributions and hence these data were log-transformed and geometric means with 95% confidence intervals have been reported. Correlations between expression levels of ribosomal genes and between normalization factors generated with two, three, four or five genes were calculated using Pearson correlation coefficients. All statistical analyses were performed with SPSS version 16 (SPSS Inc., Chicago, IL). A p-value < 0.05 was taken to indicate statistical significance.

Results

Twenty genes were identified as being represented by at least one probe set exhibiting low variation and high abundance among all samples (Figure 1). From this analysis 7 genes were selected as candidate reference genes for further testing with RT-qPCR based on whether other probe sets for that gene also had small variance and high abundance across samples. The genes selected were GAPDH, RPL22, RPS4X, TPT1, RPL23A, RPL41 and SRP14 (see Table 2 for full gene name, chromosomal location, function and expression levels in heart-failed and control tissue). In addition, analysis of 50 genes previously used or proposed by prior studies as reference genes [6, 10, 14] identified an additional three genes with small variance across probes: RPL13A, EEF1A1 and RNPS1 (Table 2, see Additional file 1 for Affymetrix probe intensities for all 50 genes). Thus a total of 10 genes were selected as candidates for testing, of which 4 encoded ribosomal proteins (RPL22, RPL23A, RPL41 and RPL13A) and are involved in protein biosynthesis. Because the expression levels of the ribosomal genes were tightly correlated (Table 3) and co-regulation of genes may confound geNorm analyses, only one of the ribosomal proteins could be included in the analysis. A preliminary analysis identified RPL13A as the most stable ribosomal gene for inclusion in the geNorm analysis (data not shown). The remaining ribosomal genes were excluded. There was no relationship between the functional class of any of the other genes (Table 2).

Figure 1
figure 1

Selection of candidate reference genes. Expression levels of probe sets for genes with one or more probe sets among the top 20 most stable and abundantly expressed across all samples, identified from screening publicly available Affymetrix gene expression profiles of left ventricle myocardium from 195 heart transplant recipients and 16 unmatched heart donors [11] (GEO accession GSE5406, http://www.ncbi.nlm.nih.gov). GAPDH, RPL22, RPS4X, TPT1, RPL23A, RPL41 and SRP14 genes were selected as candidate reference genes (shown in red), based on their high abundance and consistent expression across the majority of probe sets. Affymetrix probe set nomenclature is preceded by the gene symbol. Boxes indicate median and interquartile range, whisker length is 150% of the interquartile range. Observations beyond the whiskers are denoted by open circles.

Full size image
Table 2 Candidate cardiac reference genes ranked in order of abundance
Full size table
Table 3 Pearson correlations of expression levels of ribosomal genes
Full size table

Genes were ranked from least stable to most stable by geNorm analysis: GAPDH, RPS4X, RPL13A, RNPS1, EEF1A1, TPT1/SRP14 (Figure 2). All genes had high expression stability (M-values <1), below the default limit of M = 1.5, except for GAPDH. GAPDH was considerably more variably expressed (M-value = 1.94). GeNorm analysis indicated that the top five most stable genes would be needed for reliable normalization of RT-qPCR data (geNorm recommended threshold = 0.15, Figure 3).

Figure 2
figure 2

Average expression variability of reference genes (M), during stepwise exclusion of the least stable gene. At each step the most variable gene (highest M-value) is excluded and M is recalculated. This process continues in a step-wise manner until the two most stable genes remain. The figure indicates that GAPDH is least stably expressed relative to the other genes, and that SRP14 and TPT1 have the most stable expression.

Full size image
Figure 3
figure 3

Determination of the ideal number of reference genes for normalization. Using the two most stable genes as a starting point (SRP14 and TPT1), genes are added sequentially (from most stable to least stable) until the pairwise variation between two sequential normalization factors drops below the recommended threshold of 0.15, indicating that the prior gene set is sufficient for accurate normalization. Variation greater than 0.15 indicates that the added gene has a significant effect and should preferably be included, provided it is stably expressed. This figure shows that the sixth most stable gene (RPS4X) is not required for reliable normalization of RT-PCR data.

Full size image

Univariate analysis showed that the expression levels of all candidate reference genes except for TPT1 (p = 0.080), EEF1A1 (p = 0.116) and RNPS1 (p = 0.061) were significantly lower in the myocardium of end-stage heart failure patients compared with non heart-failed donors (Table 2). Thus, these genes have utility as a reference within heart-failed myocardium samples or non-failed donor myocardium samples separately, but not for comparing gene expression between these groups. Notably expression levels of the commonly used reference gene, GAPDH, differed significantly between failed and non-failed myocardium (p = 0.015, Table 2).

Within heart-failed and non-failed myocardium samples separately, we explored whether fewer than five reference genes could be used for normalization. A series of normalization factors were generated for each sample using geometric means of the top two, three or four most stable genes and correlated with the original set of normalization factors generated using all five top-ranked genes. Very tight correlation was observed between normalization factors for all gene sets (Pearson correlation coefficients, heart-failed samples: 0.976–0.995, p < 0.001; non-failed controls: 0.973–0.992, p < 0.001). Figure 4 illustrates the tight correlation between normalization factors generated with the top two (SRP14 and TPT1) and the top five most stable genes, indicating that as few as two genes may be sufficient for reliable normalization within failed and non-failed myocardium. Normalization factors generated using SRP14 and TPT1 were not correlated with GAPDH levels, in failed myocardium (Pearson correlation coefficient = -0.153, p = 0.436) or non-failed myocardium (Pearson correlation coefficient = -0.181, p = 0.356), suggesting that normalization of gene expression data with GAPDH in these samples would be inappropriate.

Figure 4
figure 4

Correlation of individual normalization factors generated using the top five most stable genes and the top two most stable genes. Tight correlation between normalization factors indicates that as few as two genes may be sufficient for reliable normalization of RT-qPCR data for (A) heart-failed myocardium and (B) non-failed donor myocardium. PCC = Pearson correlation coefficient.

Full size image

The expression levels of TPT1, EEF1A1 and RNPS1, did not differ significantly between failed and non-failed human myocardium samples (Table 2). These genes ranked among the top four most stable genes of those tested (Figure 2). Normalization of gene expression in failed and non-failed myocardium using these genes was tested by comparing normalized and raw expression levels of NPPB, a gene known to be up-regulated in heart-failed heart tissue. In the raw data NPPB levels were higher in heart-failed patients compared with controls, but the difference did not reach statistical significance (p = 0.058, Table 4). In normalized data, NPPB levels were significantly higher in heart-failed patients compared with controls (p = 0.023, Table 4). Removing the least stable gene, RNPS1, from the pool of reference genes only marginally affected the results (p = 0.024, Table 4). Mean NPPB levels were almost identical to data normalized with RNPS1 included, indicating that as few as two genes (TPT1 and EEF1A1) may be sufficient for reliable normalization. Normalization factors generated using TPT1, EEF1A1 and RNPS1 were not correlated with GAPDH levels (Pearson correlation coefficient = -0.085, p = 0.535). Normalizing NPPB to GAPDH alone gave a larger difference between failed and non-failed myocardium (4.2-fold, p = 0.012, Table 4) compared with data normalized with TPT1, EEF1A1 and RNPS1 combined (3.4-fold, p = 0.023, Table 4).

Table 4 Normalization of NPPB using fewer than three reference genes
Full size table

Discussion

Validation of tissue-specific reference genes is a fundamental first-step in RT-qPCR analysis. We have identified two genes, SRP14 and TPT1, that in combination can be used for accurate and reliable normalization of RT-qPCR gene expression data within failed human myocardium and non-failed human myocardium separately, and two genes, TPT1 and EEF1A1 that are more stably expressed than GAPDH and may be used for normalization of RT-qPCR data when comparing gene expression levels between failed and non-failed human myocardium samples.

The physiological functions of SRP14, TPT1 and EEF1A1 relate to important independent cellular processes, emphasizing their utility as reference genes. SRP14 is a component of the signal recognition particle, a universally conserved ribonucleoprotein complex that mediates targeting of membrane and secretory proteins to the endoplasmic reticulum [19]. TPT1 encodes the abundant and highly conserved translationally controlled tumor protein (TCTP), which is reported to play a role in cell growth and cell cycle progression [20]. EEF1A1 is a component of the alpha subunit of the elongation factor-1 complex, which facilitates the enzymatic delivery of aminoacyl tRNAs to the ribosome [21]. Notably, TPT1 is expressed in a tissue- and development-specific manner, suggesting that its utility as a reference gene in human myocardium may be limited to analyses in adult hearts [22].

To date, reference genes have been identified for a wide range of species, tissue types, experimental conditions and disease states. The typical methodological approach used in these studies has been to test a selection of commonly used reference genes such as GAPDH, β-actin, α-tubulin, hypoxanthine phosphoribosyl-transferase 1 and 18S RNA to identify the most stable combination specific to the sample and experimental conditions under investigation [23–29]. A small number of studies have reported that none of the traditional reference genes tested were appropriate for use in the particular setting investigated [30–32]. More recently, a genome-wide approach has been employed, where candidate reference genes selected from publicly available microarray data has lead to identification of novel transcripts that are more stably expressed than commonly used reference genes. These studies can be broadly classified into two groups: those that screen expression data from a diverse range of tissue types and species to identify transcripts that are stably and ubiquitously expressed and may therefore have board utility as reference genes [6, 10, 33–35], and those that screen array data from a specific tissue or disease setting to identify a set of candidate genes for validation in an independent sample set, as in the present study. To date, studies in the latter group have validated panels of novel reference genes in human normal and cancerous lung tissues [36, 37], human normal and tumor gastrointestinal tissues [38], human embryonic stem cells [39], porcine mammary tissue [40], canine osteoarthritic joint tissue [41], barley [42] and lymphoblastoid cell lines and fibroblasts [43]. Notably, Peltier et. al. have identified a set of microRNA references for use in RT-qPCR analysis of microRNAs in normal and cancerous human solid tumors [44]. However, no validated universal reference genes have been identified.

In the present study, screening of Affymetrix cardiac expression data enabled us to identify four novel genes and six genes commonly used, or recommended for use, as references in human tissue, as potentially suitable for use in heart tissue RT-qPCR analyses. Interestingly, there were striking similarities between the genes we identified as stably expressed in cardiac microarrays and those identified in a recent meta-analysis of 13,629 human gene arrays of a diverse range of cell types and experimental conditions [6]. Both studies identified SRP14, RPL13A and RPL22 as being highly stably expressed, and found an over-representation of ribosomal genes among the top ranked transcripts. However, we also identified novel candidate genes specific to heart tissue, including TPT1, one of the most stably expressed genes in our analysis. Our finding that all genes selected were more stably expressed than GAPDH is consistent with an analogous study in gastrointestinal normal and tumor tissues [38], where GAPDH was the less stably expressed genes in tissue from stomach, small intestine, liver, and lymph nodes compared with expression of 8 candidate reference genes selected from analysis of publicly available Affymetrix array data.

GeNorm [5] is one of several statistical algorithms that evaluate the relative expression stability of genes on the basis of non-normalized expression levels. Other programs for assessing gene stability include Best keeper [45], NormFinder [46], Global Pattern Recognition [47] and equivalence tests [48]. These programs have provided a rational basis for selection of reference genes for normalization of RT-qPCR data. While the relative gene stability rankings may vary subtly depending on which algorithm is applied [41, 49, 50], use of multiple reference genes provides a considerably more robust result compared with relying on a single RNA transcript [5]. This reflects the heterogeneity inherent in human samples related to age, gender, ethnicity, and gene-environment interactions. In the present study all donors were on life-support as a result of head trauma or cerebral vascular accident and thus the gene expression profile of the donor heart samples may have been affected by the traumatic events and acute drug treatments that preceded the donation of tissue. Similarly, gene expression in heart-failed tissue will vary depending on the etiology and duration of heart failure and chronic and acute drug treatments.

GeNorm analysis enables the ideal number of reference genes required for reliable normalization to be determined. Our analysis indicated that the top five genes would be required for accurate normalization of RT-qPCR data within failed and non-failed human myocardium, however, the number of genes used needs to be a balance between accuracy and practical considerations. Realistically, if two relatively stable genes have been identified, it is unnecessary to add more genes if the normalization factor does not change markedly when they are included. We were able to demonstrate that as few as two genes may be sufficient for normalization in human myocardium (SRP14 and TPT1 for analysis within failed or non-failed tissue and TPT1 and EEF1A1 for comparisons of gene expression between these groups), as the variation between normalization factors generated using two, three, four or five genes was very small. Furthermore, when comparing NPPB levels between failed and non-failed myocardium, removal of the least stable gene, RNPS1, impacted minimally on normalized NPPB levels, suggesting that TPT1 and EEF1A1 may be sufficient for normalization. This reflects the uniformity in expression stability across the top-ranked genes (Figure 2) and is equivalent to using a pairwise variation threshold of 0.23.

Several studies have demonstrated how a single, unvalidated reference gene can generate biased results if it is itself altered by the experimental conditions [5, 51–53], In contrast to GAPDH, the expression levels of TPT1 and EEF1A1 did not differ significantly between failed and non-failed human myocardium samples, and they were ranked as being considerably more stably expressed. However, because expression of both genes tended to be lower in failed myocardium compared with non-failed myocardium, TPT1 and EEF1A1 require further validation in a larger sample to confirm that their expression levels do not differ between these groups. Normalizing NPPB expression to GAPDH alone gave a greater difference in NPPB levels between failed and non-failed myocardium than the data normalized to TPT1 and EEF1A1. However, the expression levels of GAPDH were lower in failed heart samples compared with healthy hearts, and this will have artificially skewed the data. Normalizing NPPB expression to TPT1 and EEF1A1 gave a modestly significant difference in NPPB levels between failed and non-failed myocardium that is more likely to reflect the true level of expression of NPPB in these two groups. In particular our finding that GAPDH is differentially expressed in failing and non-failing myocardium has implications for previous RT-qPCR studies of human myocardium that have used GAPDH for normalization without validation.

We have identified reference genes for both failing and non-failing human myocardium that have significantly improved stability compared with the commonly used reference gene, GAPDH (SRP14 and TPT1 for analysis within failed or non-failed tissue and TPT1 and EEF1A1 for comparisons of gene expression between these groups). Our findings suggest that routine normalization of RT-qPCR data in human myocardium with GAPDH should be avoided as it was the least stably expressed of all genes tested. This highlights the importance of validating reference genes for normalization of RT-qPCR data. The reference genes identified in this study will enable more reliable interpretation of RT-qPCR results in these tissues. Whilst TPT1 and EEF1A1 require further validation in a larger study, these findings serve as a basic guideline for reference gene selection in human myocardium.

Conclusion

This study has identified a strategy for normalization of cardiac gene expression involving two reference genes, SRP14 and TPT1, that in combination may have broad application for accurate and reliable normalization of RT-qPCR data within failed human myocardium and non-failed human myocardium separately. In addition, for comparisons of gene expression between failed and non-failed human myocardium TPT1 and EEF1A1, in combination, may provide a more reliable reference than GAPDH for normalization of RT-qPCR data for these analyses.

Abbreviations

cDNA:

complementary deoxyribonucleic acid

dsDNA:

double stranded deoxyribonucleic acid

EEF1A1:

eukaryotic elongation factor 1A1

GAPDH:

glyceraldehyde-3-phosphate dehydrogenase

GEO:

Gene Expression Omnibus

mRNA:

messenger ribonucleic acid

NCBI:

National Center for Biotechnology Information

NPPB:

brain natriuretic peptide precursor

PCR:

polymerase chain reaction

RNA:

ribonucleic acid

RNPS1:

ribonucleic acid binding protein S1

RPL13A:

ribosomal protein L13a

RPL22:

ribosomal protein L22

RPL23A:

ribosomal protein L23a

RPL41:

ribosomal protein L41

RPS4X:

ribosomal protein S4: X-linked

RT-qPCR:

real-time polymerase chain reaction

SRP14:

signal recognition particle 14

TCTP:

translationally controlled tumor protein

TPT1:

tumor protein: translationally-controlled 1.

References

  1. Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of molecular endocrinology. 2000, 25: 169-193. 10.1677/jme.0.0250169.

    Article  CAS  PubMed  Google Scholar 

  2. Hendriks-Balk MC, Michel MC, Alewijnse AE: Pitfalls in the normalization of real-time polymerase chain reaction data. Basic research in cardiology. 2007, 102: 195-197. 10.1007/s00395-007-0649-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Wang T, Brown MJ: mRNA quantification by real time TaqMan polymerase chain reaction: validation and comparison with RNase protection. Analytical biochemistry. 1999, 269: 198-201. 10.1006/abio.1999.4022.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  5. 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-10.1186/gb-2002-3-7-research0034.

    Article  PubMed  PubMed Central  Google Scholar 

  6. de Jonge HJ, Fehrmann RS, de Bont ES, Hofstra RM, Gerbens F, Kamps WA, de Vries EG, Zee van der AG, te Meerman GJ, ter Elst A: Evidence based selection of housekeeping genes. PLoS ONE. 2007, 2: e898-10.1371/journal.pone.0000898.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Brattelid T, Tveit K, Birkeland JA, Sjaastad I, Qvigstad E, Krobert KA, Hussain RI, Skomedal T, Osnes JB, Levy FO: Expression of mRNA encoding G protein-coupled receptors involved in congestive heart failure – a quantitative RT-PCR study and the question of normalisation. Basic research in cardiology. 2007, 102: 198-208. 10.1007/s00395-007-0648-1.

    Article  CAS  PubMed  Google Scholar 

  8. Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. Journal of biochemical and biophysical methods. 2000, 46: 69-81. 10.1016/S0165-022X(00)00129-9.

    Article  CAS  PubMed  Google Scholar 

  9. Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. BioTechniques. 2000, 29: 332-337.

    CAS  PubMed  Google Scholar 

  10. Warrington JA, Nair A, Mahadevappa M, Tsyganskaya M: Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiological genomics. 2000, 2: 143-147.

    CAS  PubMed  Google Scholar 

  11. Hannenhalli S, Putt ME, Gilmore JM, Wang J, Parmacek MS, Epstein JA, Morrisey EE, Margulies KB, Cappola TP: Transcriptional genomics associates FOX transcription factors with human heart failure. Circulation. 2006, 114: 1269-1276. 10.1161/CIRCULATIONAHA.106.632430.

    Article  CAS  PubMed  Google Scholar 

  12. Hubbell E, Liu WM, Mei R: Robust estimators for expression analysis. Bioinformatics (Oxford, England). 2002, 18: 1585-1592. 10.1093/bioinformatics/18.12.1585.

    Article  CAS  Google Scholar 

  13. A Language and Environment for Statistical Computing. [http://www.R-project.org]

  14. geNorm. [http://medgen.ugent.be/~jvdesomp/genorm/]

  15. Yang J, Moravec CS, Sussman MA, DiPaola NR, Fu D, Hawthorn L, Mitchell CA, Young JB, Francis GS, McCarthy PM, Bond M: Decreased SLIM1 expression and increased gelsolin expression in failing human hearts measured by high-density oligonucleotide arrays. Circulation. 2000, 102: 3046-3052.

    Article  CAS  PubMed  Google Scholar 

  16. Zakhary DR, Moravec CS, Stewart RW, Bond M: Protein kinase A (PKA)-dependent troponin-I phosphorylation and PKA regulatory subunits are decreased in human dilated cardiomyopathy. Circulation. 1999, 99: 505-510.

    Article  CAS  PubMed  Google Scholar 

  17. Ellmers LJ, Scott NJ, Piuhola J, Maeda N, Smithies O, Frampton CM, Richards AM, Cameron VA: Npr1-regulated gene pathways contributing to cardiac hypertrophy and fibrosis. Journal of molecular endocrinology. 2007, 38: 245-257. 10.1677/jme.1.02138.

    Article  CAS  PubMed  Google Scholar 

  18. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  19. Lakkaraju AK, Mary C, Scherrer A, Johnson AE, Strub K: SRP keeps polypeptides translocation-competent by slowing translation to match limiting ER-targeting sites. Cell. 2008, 133: 440-451. 10.1016/j.cell.2008.02.049.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bommer UA, Thiele BJ: The translationally controlled tumour protein (TCTP). The international journal of biochemistry & cell biology. 2004, 36: 379-385. 10.1016/S1357-2725(03)00213-9.

    Article  CAS  Google Scholar 

  21. Brands JH, Maassen JA, van Hemert FJ, Amons R, Moller W: The primary structure of the alpha subunit of human elongation factor 1. Structural aspects of guanine-nucleotide-binding sites. European journal of biochemistry/FEBS. 1986, 155: 167-171. 10.1111/j.1432-1033.1986.tb09472.x.

    Article  CAS  PubMed  Google Scholar 

  22. Thiele H, Berger M, Skalweit A, Thiele BJ: Expression of the gene and processed pseudogenes encoding the human and rabbit translationally controlled tumour protein (TCTP). European journal of biochemistry/FEBS. 2000, 267: 5473-5481. 10.1046/j.1432-1327.2000.01609.x.

    Article  CAS  PubMed  Google Scholar 

  23. Ahn K, Huh JW, Park SJ, Kim DS, Ha HS, Kim YJ, Lee JR, Chang KT, Kim HS: Selection of internal reference genes for SYBR green qRT-PCR studies of rhesus monkey (Macaca mulatta) tissues. BMC molecular biology. 2008, 9: 78-10.1186/1471-2199-9-78.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Daud II, Scott ME: Validation of reference genes in cervical cell samples from human papillomavirus-infected and -uninfected women for quantitative reverse transcription-PCR assays. Clin Vaccine Immunol. 2008, 15: 1369-1373. 10.1128/CVI.00074-08.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Langnaese K, John R, Schweizer H, Ebmeyer U, Keilhoff G: Selection of reference genes for quantitative real-time PCR in a rat asphyxial cardiac arrest model. BMC molecular biology. 2008, 9: 53-10.1186/1471-2199-9-53.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Olbrich M, Gerstner E, Welzl G, Fleischmann F, Osswald W, Bahnweg G, Ernst D: Quantification of mRNAs and housekeeping gene selection for quantitative real-time RT-PCR normalization in European beech (Fagus sylvatica L.) during abiotic and biotic stress. Zeitschrift fur Naturforschung. 2008, 63: 574-582.

    CAS  PubMed  Google Scholar 

  27. Perez R, Tupac-Yupanqui I, Dunner S: Evaluation of suitable reference genes for gene expression studies in bovine muscular tissue. BMC molecular biology. 2008, 9: 79-10.1186/1471-2199-9-79.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Waxman S, Wurmbach E: De-regulation of common housekeeping genes in hepatocellular carcinoma. BMC genomics. 2007, 8: 243-10.1186/1471-2164-8-243.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Yoo WG, Kim TI, Li S, Kwon OS, Cho PY, Kim TS, Kim K, Hong SJ: Reference genes for quantitative analysis on Clonorchis sinensis gene expression by real-time PCR. Parasitology research. 2008

    Google Scholar 

  30. Giricz O, Lauer-Fields JL, Fields GB: The normalization of gene expression data in melanoma: investigating the use of glyceraldehyde 3-phosphate dehydrogenase and 18S ribosomal RNA as internal reference genes for quantitative real-time PCR. Analytical biochemistry. 2008, 380: 137-139. 10.1016/j.ab.2008.05.024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Gutierrez L, Mauriat M, Guenin S, Pelloux J, Lefebvre JF, Louvet R, Rusterucci C, Moritz T, Guerineau F, Bellini C, van Wuytswinkel O: The lack of a systematic validation of reference genes: a serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant biotechnology journal. 2008, 6: 609-618. 10.1111/j.1467-7652.2008.00346.x.

    Article  CAS  PubMed  Google Scholar 

  32. Sellars MJ, Vuocolo T, Leeton LA, Coman GJ, Degnan BM, Preston NP: Real-time RT-PCR quantification of Kuruma shrimp transcripts: a comparison of relative and absolute quantification procedures. Journal of biotechnology. 2007, 129: 391-399. 10.1016/j.jbiotec.2007.01.029.

    Article  CAS  PubMed  Google Scholar 

  33. Frericks M, Esser C: A toolbox of novel murine house-keeping genes identified by meta-analysis of large scale gene expression profiles. Biochimica et biophysica acta. 2008

    Google Scholar 

  34. Lee S, Jo M, Lee J, Koh SS, Kim S: Identification of novel universal housekeeping genes by statistical analysis of microarray data. J Biochem Mol Biol. 2007, 40 (2): 226-231.

    Article  CAS  PubMed  Google Scholar 

  35. Zhu J, He F, Song S, Wang J, Yu J: How many human genes can be defined as housekeeping with current expression data?. BMC genomics. 2008, 9: 172-10.1186/1471-2164-9-172.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Saviozzi S, Cordero F, Lo Iacono M, Novello S, Scagliotti GV, Calogero RA: Selection of suitable reference genes for accurate normalization of gene expression profile studies in non-small cell lung cancer. BMC cancer. 2006, 6: 200-10.1186/1471-2407-6-200.

    Article  PubMed  Google Scholar 

  37. Su LJ, Chang CW, Wu YC, Chen KC, Lin CJ, Liang SC, Lin CH, Whang-Peng J, Hsu SL, Chen CH, Huang CY: Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC genomics. 2007, 8: 140-10.1186/1471-2164-8-140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kidd M, Nadler B, Mane S, Eick G, Malfertheiner M, Champaneria M, Pfragner R, Modlin I: GeneChip, geNorm, and gastrointestinal tumors: novel reference genes for real-time PCR. Physiological genomics. 2007, 30: 363-370. 10.1152/physiolgenomics.00251.2006.

    Article  CAS  PubMed  Google Scholar 

  39. Synnergren J, Giesler TL, Adak S, Tandon R, Noaksson K, Lindahl A, Nilsson P, Nelson D, Olsson B, Englund MC, Abbot S, Sartipy P: Differentiating human embryonic stem cells express a unique housekeeping gene signature. Stem cells (Dayton, Ohio). 2007, 25: 473-480. 10.1634/stemcells.2006-0247.

    Article  CAS  Google Scholar 

  40. Tramontana S, Bionaz M, Sharma A, Graugnard DE, Cutler EA, Ajmone-Marsan P, Hurley WL, Loor JJ: Internal controls for quantitative polymerase chain reaction of swine mammary glands during pregnancy and lactation. Journal of dairy science. 2008, 91: 3057-3066. 10.3168/jds.2008-1164.

    Article  CAS  PubMed  Google Scholar 

  41. 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 molecular biology. 2007, 8: 62-10.1186/1471-2199-8-62.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Faccioli P, Ciceri GP, Provero P, Stanca AM, Morcia C, Terzi V: A combined strategy of "in silico" transcriptome analysis and web search engine optimization allows an agile identification of reference genes suitable for normalization in gene expression studies. Plant molecular biology. 2007, 63: 679-688. 10.1007/s11103-006-9116-9.

    Article  CAS  PubMed  Google Scholar 

  43. de Brouwer AP, van Bokhoven H, Kremer H: Comparison of 12 reference genes for normalization of gene expression levels in Epstein-Barr virus-transformed lymphoblastoid cell lines and fibroblasts. Molecular diagnosis & therapy. 2006, 10: 197-204.

    Article  CAS  Google Scholar 

  44. Peltier HJ, Latham GJ: Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA. 2008, 14 (5): 844-852. 10.1261/rna.939908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. 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. Biotechnology letters. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47.

    Article  CAS  PubMed  Google Scholar 

  46. 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 research. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.

    Article  CAS  PubMed  Google Scholar 

  47. Akilesh S, Shaffer DJ, Roopenian D: Customized molecular phenotyping by quantitative gene expression and pattern recognition analysis. Genome research. 2003, 13: 1719-1727. 10.1101/gr.533003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Haller F, Kulle B, Schwager S, Gunawan B, von Heydebreck A, Sultmann H, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Analytical biochemistry. 2004, 335: 1-9. 10.1016/j.ab.2004.08.024.

    Article  CAS  PubMed  Google Scholar 

  49. Robinson TL, Sutherland IA, Sutherland J: Validation of candidate bovine reference genes for use with real-time PCR. Veterinary immunology and immunopathology. 2007, 115: 160-165. 10.1016/j.vetimm.2006.09.012.

    Article  CAS  PubMed  Google Scholar 

  50. Spinsanti G, Panti C, Lazzeri E, Marsili L, Casini S, Frati F, Fossi CM: Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies. BMC molecular biology. 2006, 7: 32-10.1186/1471-2199-7-32.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Bas A, Forsberg G, Hammarstrom S, Hammarstrom ML: Utility of the housekeeping genes 18S rRNA, beta-actin and glyceraldehyde-3-phosphate-dehydrogenase for normalization in real-time quantitative reverse transcriptase-polymerase chain reaction analysis of gene expression in human T lymphocytes. Scandinavian journal of immunology. 2004, 59: 566-573. 10.1111/j.0300-9475.2004.01440.x.

    Article  CAS  PubMed  Google Scholar 

  52. Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GA, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Analytical biochemistry. 2005, 344: 141-143. 10.1016/j.ab.2005.05.022.

    Article  CAS  PubMed  Google Scholar 

  53. 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. Analytical biochemistry. 2002, 309: 293-300. 10.1016/S0003-2697(02)00311-1.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We gratefully acknowledge the donation of human myocardium by heart transplant patients and the families of heart donors for research purposes. This work was funded by the National Heart Foundation of New Zealand, the New Zealand Foundation of Research, Science and Technology and the Health Research Council of New Zealand.

Author information

Authors and Affiliations

  1. Christchurch Cardioendocrine Research Group, Department of Medicine, University of Otago-Christchurch, PO Box 4345, Christchurch, 8014, New Zealand

    Anna P Pilbrow, Leigh J Ellmers, Richard W Troughton, A Mark Richards, Chris M Frampton & Vicky A Cameron

  2. Department of Biochemistry, University of Otago, PO Box 56, Dunedin, 9054, New Zealand

    Michael A Black

  3. Kaufman Center for Heart Failure, Department of Cardiovascular Medicine, Cleveland Clinic, 9500 Euclid Avenue, Cleveland, Ohio, 44195, USA

    Christine S Moravec & Wendy E Sweet

Authors
  1. Anna P Pilbrow
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Leigh J Ellmers
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Michael A Black
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Christine S Moravec
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Wendy E Sweet
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Richard W Troughton
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. A Mark Richards
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Chris M Frampton
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Vicky A Cameron
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Vicky A Cameron.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

APP was involved in the conception and design of the study, performed the gene expression studies and statistical analysis, and wrote the first draft of the manuscript. LJE participated in the data acquisition and provided critical revision of the manuscript. MAB performed the bioinformatics analysis and provided critical revision of the manuscript. CSM was involved in the conception, design and coordination of the study and management of the Cleveland heart tissue bank. WES oversaw tissue acquisition and data coordination. RWT was involved in the conception and coordination of the study and acquisition of funding. AMR was involved with acquisition of funding and provided critical revision of the manuscript. CMF contributed to the study design and performed statistical analysis. VAC was involved in the conception, design and coordination of the study, acquisition of funding, interpretation of the data and provided critical revision of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12920_2008_64_MOESM1_ESM.pdf

Additional file 1: Selection of candidate reference genes. Expression levels of probe sets for 50 genes previously used or recommended for use as reference genes from publically available Affymetrix gene expression profiles of left ventricle myocardium from 195 heart transplant recipients and 16 unmatched heart donors [11] (GEO accession GSE5406, http://www.ncbi.nlm.nih.gov). RPL13A, EEF1A1 and RNPS1 (shown in red) were selected as candidate reference genes from this analysis, based on their high abundance and consistent expression across the majority of probes. Affymetrix probe set nomenclature is preceded by the gene symbol. Boxes indicate median and interquartile range, whisker length is 150% of the interquartile range. Observations beyond the whiskers are denoted by open circles. (PDF 7 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), 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

Pilbrow, A.P., Ellmers, L.J., Black, M.A. et al. Genomic selection of reference genes for real-time PCR in human myocardium. BMC Med Genomics 1, 64 (2008). https://doi.org/10.1186/1755-8794-1-64

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1755-8794-1-64

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us