|Home | About | Journals | Submit | Contact Us | Français|
Reliable reference genes are a vital prerequisite for any functional study employing quantitative real-time RT-PCR (RT-qPCR) for analyzing gene expression. Yet a proper selection and assessment of the chosen reference genes is only rarely included into a study. To date, no reference genes have been validated for differentiation of THP-1 monocytes. Here we report on the selection of validated reference genes during differentiation of THP-1 monocytes into macrophages induced by phorbol 12-myristate 13-acetate (PMA).
The mRNA expression of 21 preselected potential reference genes was measured by RT-qPCR at several time-points over six days of PMA-induced THP-1 monocyte-to-macrophage differentiation. A ranking according to expression stability was calculated. Calculations were performed using Microsoft Excel-based applets GeNorm, NormFinder and BestKeeper. Our results indicated ACTB (β-actin) (Cq ± SD, 14.1 ± 0.3) and RPL37A (ribosomal protein L37a) (14.5 ± 0.3) as the most stable genes. While other frequently used reference genes such as GAPDH (glycereraldehyde-3-phosphate dehydrogenase) (20.8 ± 0.8) or G6PD (glucose-6-phophate dehydrogenase) (16.1 ± 1.0) were found to be not as reliable and were therefore unsuited for use as reference genes. These findings were validated by investigating mRNA expression of macrophage scavenger receptor CD36, known to be regulated during monocyte-to-macrophage differentiation. Using ACTB and RPL37A as reference genes a profound and significant regulation of CD36 could be demonstrated, while use of G6PD resulted in a much less pronounced apparent regulation of CD36.
Consequently, it is recommended to normalize any real-time PCR-based expression data obtained during THP-1 monocyte differentiation using ACTB and RPL37A.
Macrophages are the phagocytic cells of the immune system which play a pivotal role in many disease processes . Upon a local stimulus circulating blood monocytes immigrate into the respective tissue where they differentiate into mature macrophages. Due to the very complex interactions of macrophages with the cells of the surrounding tissue and their manifold activities, which include removal of necrotic and apoptotic tissue or invaded microorganisms, contribution to wound healing, and presentation of antigens [1,2], intensive effort in basic and clinical research has been spend on unraveling the biology of macrophages and their behavior.
Although pure human primary monocytes can be obtained either by dextran sedimentation followed by Ficoll density centrifugation , or leukapharesis combined with counter current elutriation , the number of primary monocytes for functional studies is limited due to insufficient proliferation . Therefore immortalized proliferating cell lines are often used instead, such as the human THP-1 monocytic leukemia cell line. This cell line is a well-established model, for example in toxicology, immunology and atherosclerosis research, regarding monocyte and macrophage function and biology [6-8].
Phorbol esters, such as phorbol 12-myristate 13-acetate (PMA), are frequently used to elicit the differentiation of THP-1 monocytes into macrophage-like cells which mimic many characteristic features of human primary macrophages [6,9]. The process of differentiation is accompanied by profound changes within the cells, as the cells become adherent and adjust their morphology and physiology [10-12].
Quantitative real-time RT-PCR (RT-qPCR) is a powerful tool for quantifying RNA expression and determining differences in expression levels. Yet, in order to correctly assess the results obtained, a reliable reference is strictly required. As RT-qPCR expression data depend on a multitude of factors, such as amount and quality of isolated RNA, efficiencies of enzymes (reverse transcriptase and DNA-dependent DNA polymerase), and overall variability in transcriptional activity between samples [13-15], it is common procedure to normalize expression of a gene of interest using an internal standard, i.e. a reference gene. Usually genes known to be stably expressed in general, such as GAPDH or ACTB, are chosen as internal standard. Yet the validity of this assumption is rarely verified, although over the past years evidence has emerged that genes previously thought to be stably expressed might actually be regulated under certain conditions [16,17]. Therefore the recently published MIQE (Minimum Information for publication of Quantitative real-time PCR Experiments) guidelines propose reference gene validation for all RT-qPCR experiments . In this study we aimed at identifying suitable reference genes during differentiation of human THP-1 monocytes to macrophages. To the best of our knowledge, this is the first study establishing reference genes for the widely used THP-1 cell culture model of monocyte maturation to macrophages.
Selection of reference genes and proper assessment of the stability of selected genes is a circular problem since normalization is required in order to verify the stability of expression of a particular gene. This problem can be circumvented by using several reference genes in combination [13-15]. Based on this approach different Microsoft Excel-based tools have been developed, which allow selection of the most stably expressed genes from a set of expression data of several potential reference genes. For our study we used the GeNorm applet published by Vandesompele et al. , the BestKeeper applet developed by Pfaffl et al. , and the NormFinder applet published by Lindbjerg et al. .
The aim of our study was to identify and validate reliable reference genes for the THP-1 cell culture model of monocyte maturation. For this we have preselected 21 genes from a list of reference genes, which had shown promising little variation in previous experiments , as well as from a list of genes often used in studies involving human macrophages [20-23], and thus were considered potential reference genes. Background information on the preselected potential 21 reference genes is provided as Additional file 1, Table S1. The observed range of Cq values was distributed over a fairly large range representing highly expressed genes (e.g. B2M [Cq ± SD, 13.1 ± 0.7] and ACTB [14.1 ± 0.3]), as well as less abundant mRNAs (e.g. UBE2D2 [21.3 ± 0.6] and GAPDH [20.8 ± 0.8]) (Figure (Figure1).1). Diagrams representing the variance of Cq values at each point of measurement for every single gene are available as Additional file 2, Figure S1.
Obtained raw Cq values were manually transformed into GeNorm data input format and afterwards analyzed by GeNorm. First level analysis performed by GeNorm creates a ranking of all reference genes based on their pairwise variances (Figure (Figure2).2). This identified the pair of ACTB and RPL37A as the most stable genes of the entire set during differentiation of THP-1 monocytes into macrophages; PPARD and PPARG were the least stable genes (Figure (Figure2),2), which also show highest variation of Cq values during maturation of THP-1 cells (Figure (Figure11).
On the second level of analysis the number of reference genes, which need to be included into the assessment in order to provide a reliable normalization, was determined. For this purpose pairwise variances of the normalization factors were calculated when an additional reference gene is included or not. Vandesompele et al. suggested a cut-off at 0.15 , i.e. if the inclusion of a further reference gene into the calculation of the normalization factors causes a pairwise variance of less than 0.15 as opposed to the exclusion of the respective reference gene, than this gene has no significant contribution to the normalization factor and does not need to be included into the calculation. On the other hand if the variance exceeds 0.15 there is a significant gain in stability and reliability and the respective reference gene should be considered for the calculations. A graphical representation of the calculated pairwise variances is therefore given as Additional file 3, Figure S2. According to these pairwise variances we recommend to include the first two reference genes (ACTB and RPL37A) only into the normalization factor, since none of the pairwise variances actually exceeds the threshold of 0.15.
In order to demonstrate the importance of combining several reliable reference genes into a set of normalization factors for normalization of measured Cq values, we show in Figure Figure33 the results of two different normalizations for a single set of measurements. CD36, the gene encoding the scavenger receptor CD36, is known to be regulated during THP-1 monocyte differentiation . The CD36 protein contributes to the uptake of oxidized LDL particles and fatty acids by macrophages [25,26]. The obtained Cq values were normalized using the two most reliable reference genes identified as described above (ACTB and RPL37A; Figure Figure3A),3A), and the widely used reference gene G6PD (Figure (Figure3B),3B), respectively. As shown in Figure Figure3,3, normalization with G6PD in comparison to ACTB and RPL37A indicates a much less pronounced regulatory effect.
The stability ranking of the reference genes determined by GeNorm was assessed by comparison with two alternative Microsoft Excel-based applets, termed NormFinder , and BestKeeper . In Table Table11 the rankings produced by these applets are listed next to the GeNorm ranking obtained from our analyses. Interestingly, these rankings differ from each other in the midfield positions yet the top and bottom ranked genes are fairly constant. Especially the genes ACTB and RPL37A top ranked by GeNorm have been reproducibly listed among the four most reliable genes. For an overview of the ranking calculated by NormFinder see Additional file 4, Figure S3.
Recent studies provide clear evidence that traditional non-validated "housekeeping" genes such as GAPDH and ACTB are not stably expressed in many experimental setups [27-30], and are therefore not suitable for normalization in a broad range of cell models. Despite this knowledge, the majority of published RT-qPCR studies still lack proper validation of the reference genes. To reduce the risk of erroneous results due to instable reference gene expression we aimed at validating reference genes for the widely used THP-1 cell line which allow reliable RT-qPCR analyses of mRNA expression during maturation of monocytes into macrophages.
Our analyses of 21 preselected potential reference genes stably expressed in a broad range of tissues revealed that during differentiation of THP-1 monocytes to macrophages the normalization using a GeNorm normalization factor calculated from ACTB and RPL37A mRNA expression data is the most stable and reliable option and is sufficient for an accurate assessment of relative changes in gene expression. In contrast to this, other very frequently used reference genes, such as GAPDH and G6PD, are much less reliable. Additionally to the lack of reliability a normalization using G6PD, as shown in Figure Figure3,3, yields considerably different results when compared to normalization using the reference genes ACTB and RPL37A as recommended by GeNorm. Thus, the frequently used reference gene G6PD is actually not suitable for normalization and indeed may cause misinterpretations of results.
A comparison with recently published studies involving monocytes and macrophages further highlights the importance of individual validation of reference genes even if comparable cell models are used. For example, Piehler et al. validated reference genes in lipopolysaccharide-stimulated primary human monocytes isolated from peripheral blood . In contrast to our results, Piehler and colleagues identified PPIB and B2M as most stably expressed genes in primary human monocytes, while ACTB was found to be inapplicable as reference gene. In our THP-1 monocyte maturation model, ACTB is the most stably expressed gene whereas PPIB and B2M are inappropriate reference genes. Furthermore, ABL1 und GUSB, two reference genes often used in studies involving leukocytes and monocytes [20-22], failed both as reliable reference genes in our hands. Our findings emphasize the need of appropriate and careful validation of reference genes for cell culture models in general and the THP-1 model system in particular.
We also asked the question whether further promising reference gene candidates stably expressed during THP-1 maturation may exist. To answer this question we compared our findings with microarray data on mRNA expression in monocytes and macrophages of different origin that were available in the NCBI GEO database at the time of writing (20th September 2010). Since distinct model systems and experimental settings require individual validation of reference genes, we considered studies only closely resembling our own experimental setup. Three studies using a monocyte-to-macrophage maturation model involving either THP-1 or primary human cells were available (GEO entries GDS3554, GDS3203 and GDS2430) [32-34]. Unfortunately, these microarray raw data online provide average signal intensities for more than 20.000 probes but lack information on the specificity of the signal and whether the respective probe of the corresponding Affymetrix microarray was called present or absent. Thus, a comparative assessment of our selection of reference genes using these microarray data has severe limitations, because many of the low abundant genes detected with the Affymetrix microarrays may represent non-specific background. Nevertheless, we found that ACTB was ranked among the genes showing least variances of average signal intensity across all samples, while RPL37A was ranked in a midfield position and G6PD was among the probes with highest variances of average signal intensity. With the limitations in mind, these findings indicate that our selection of reference genes holds true also on a broader scale. However, further stable reference genes may be available for THP-1 cells in global gene expression data in public databases but further experimental validation is required prior using them for normalization of RT-qPCR data.
Apart from GeNorm other Microsoft Excel-based applications for determination of appropriate reference genes have been published; NormFinder  and Bestkeeper  were therefore used to assess the quality of the ranking obtained by GeNorm. Comparing the results of the three different applets it is remarkable that the top and the bottom ranked genes more or less retain their ranking independent of the applet used, while some of the genes ranked in the middle positions have changed their ranking quite considerably. The reproducibility of the top rankings gives high confidence regarding the actual stability of the selected reference genes. The lack of confidence within the middle positions might be due to only very slight differences in expression stability within that range (as seen in Figure Figure2).2). We therefore conclude that minimal differences within the calculation algorithms of each applet may cause the observed differences in ranking positions.
To sum up, our study provides clear evidence for the necessity to carefully validate reference genes for normalization of gene expression data obtained by RT-qPCR. As demonstrated by normalizing CD36 expression data using validated reference genes (ACTB and RPL37A) as well as the non-validated reference gene G6PD, application of appropriate reference genes may have a significant impact on experimental results. Our study has clearly shown that a GeNorm normalization factor calculated from ACTB and RPL37A mRNA expression is the most stable and reliable option for any RT-qPCR-based expression data obtained during THP-1 monocyte differentiation.
THP-1 monocytes were obtained from ATCC (Manassas, Virginia, USA) and cultured according to the supplier's recommendations. Differentiation of THP-1 monocytes into macrophages was initiated by adding 100 ng/ml phorbol-12-myristate-13-acetate (PMA, Sigma-Aldrich, Seelze, Germany) and 50 μM β-mercaptoethanol (Sigma-Aldrich) to RPMI 1640 medium (PAA, Cölbe, Germany) according to standard protocols . Cells were cultured at 37°C in a humidified 5% (v/v) CO2 air atmosphere. Cells were allowed to differentiate for 1, 2, 4 or 6 days, before they were harvested for RNA isolation. A sample of undifferentiated monocytes was included as well.
Total RNA was prepared from cell lysates using Qiagen RNeasy Mini kit (Hilden, Germany) according to manufacturer's instructions including an on-column DNase I (Qiagen) digestion as previously reported . Adequate RNA quality was assessed both by agarose gel electrophoresis and photometrically (Additional file 5, Table S2) and was in each case found to be appropriate; all 260/280 ratios were above 2.0. 260/230 ratios were between 1.0 and 2.0. In principle 260/230 ratios should be larger than 2.0 for pure RNA because lower values indicate the presence of contaminations such as guanidine thiocyanate and phenol-chloroform . Since we did not apply phenol extraction absorbance at 230 nm is caused by guanidine thiocyanate only. Previous studies have shown that even very small quantities of guanidine thiocyanate already significantly influence 260/230 ratios, albeit have no measurable effect on downstream applications such as RT-qPCR until concentrations of more than 100 mM . According to Qiagen's technical documentation our samples are suitable for reliable RT-qPCR analyses .
Synthesis of cDNA was performed using Revert Aid First strand cDNA synthesis kit (Fermentas, St. Leon-Rot, Germany). For each cDNA synthesis 5 μg of total RNA and 0.5 μg oligo(dT) primers were used as previously described . Complementary DNA was diluted tenfold and stored at -30°C prior to PCR analyses.
Real-time PCR analyses were performed on a LightCycler 480 II instrument (Roche Diagnostics, Mannheim, Germany) using the Qiagen QuantiTect SYBR Green PCR kit as previously described . Primers were designed using PrimerExpress software version 2.0.0 (Applied Biosystems, Weiterstadt, Germany) and are listed in Table Table2.2. All primer pairs were designed to have melting temperatures of about 60°C and both primers of a primer pair are located in different exons. Primers were purchased from Invitrogen (Karlsruhe, Germany). PCR runs included a 15 min pre-incubation at 95°C to allow heat activation of polymerase, followed by 40 cycles of a two-step PCR consisting of a denaturing phase at 94°C for 15 s and a combined annealing and extension phase at 60°C for 60 s. After completion of PCR a melting curve was recorded. PCR results were analyzed using the LightCycler software release version 18.104.22.168 (Roche Diagnostics). Quality of PCR amplicons was assessed by melting point analysis in order to exclude the formation of primer dimers or other by-products. Size of PCR products was confirmed by agarose gel electrophoresis. The fit point algorithm of the LightCycler software was used to calculate Cq values. In order to determine expression levels and their changes during maturation of THP-1 cells relative quantifications were calculated using GeNorm normalization factor. In order to calculate the normalization factor raw Cq values were converted into fold changes, according to the GeNorm manual the fold change of the sample with the highest expression was set to 1. From these values the normalization factor was calculated as the geometric mean of ACTB and RPL37A expression data or G6PD only . After normalization of CD36 data using the calculated GeNorm normalization factor these values were converted into conventional fold changes, by setting the value of the monocyte sample to 1 and scaling all other values proportionally. Samples were prepared in biological triplicates and always measured in technical duplicates.
Statistical analyses were performed using Microsoft Excel 2007 and three different Microsoft Excel-based applets: GeNorm, BestKeeper and NormFinder. Applets were used according to the instructions provided by the respective suppliers; this included a transformation of the raw Cq values into required data input formats for the applets GeNorm , and NormFinder . BestKeeper  analyses were based on raw Cq values without transformation. In order to test for statistical significance a one-factor ANOVA with a post-hoc Dunnett test was performed; this test is known to be rather robust against a violation of the normality condition. Statistical significances indicated by the Dunnett test were confirmed by a non-parametric permutation assay.
MBM and SS designed and performed the experiments, analyzed the data and drafted the manuscript. SL designed the study, was involved with the analysis and interpretation of data and critically helped to draft the manuscript. All authors read and approved the final manuscript.
Table S1. Gene product function and GO annotation of the genes used as potential reference genes and CD36. Information was obtained from the NCBI resource http://www.ncbi.nlm.nih.gov/ and the Gene Ontology website http://www.geneontology.org/.
Figure S1. Summary of observed variances of Cq values of the 21 preselected potential reference genes. For each potential reference gene the observed variances of Cq values at each day of differentiation of THP-1 monocytes to macrophages are shown. Squares indicate mean values. Bars represent standard deviations.
Figure S2. Determination of the number of genes required for calculating GeNorm normalization factor. Variances of pairwise combined normalization factors were calculated in order to determine which genes had to be considered for inclusion into the GeNorm normalization factor. Each bar represents the variance of the normalization factors when an additional gene is included into the calculation; the starting set of normalization factors is calculated from the two most stable genes (ACTB and RPL37A). Further genes are included according to the stability ranking calculated previously; according to Vandesompele et al. further genes are recommended to be included until the variance is below 0.15 . Reference  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:research0034.1-0034.11.
Figure S3. NormFinder analysis showing logarithmic intergroup and intragroup variances of the 21 preselected reference genes and CD36. NormFinder application was used to calculate inter- and intragroup variances as an estimate of gene stability . Squares indicate intergroup variance. Bars represent intragroup variance. Two distinct groups were defined: Group 1 is constituted of expression data measured for undifferentiated THP-1 monocytes; group 2 combines all expression data of differentiating and differentiated THP-1 macrophages. A gene's stability is represented by the distance of the respective square from the horizontal line at 0. Reference  Lindbjerg CA, Jensen JL, Ørntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 2004, 64:5245-5250.
Table S2. 260/280 and 260/230 ratios for assessment of RNA quality. RNA quality was assessed photometrically using an Eppendorf BioPhotometer plus.
This work was supported by grants from the Deutsche Infarktforschungshilfe, the Ernest-Solvay-Stiftung, the Thüringer Ministerium für Bildung, Wissenschaft und Kultur and the Wilhelm-Vaillant-Stiftung to SL. We thank Maria Braun and Thomas Böking for their excellent technical assistance.