Selection of reliable reference genes during THP-1 monocyte differentiation into macrophages
© Maeß et al; licensee BioMed Central Ltd. 2010
Received: 25 May 2010
Accepted: 1 December 2010
Published: 1 December 2010
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 (M inimum I nformation for publication of Q uantitative real-time PCR E xperiments) 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. .
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.
Comparison of stability ranking of the 21 preselected potential reference genes by BestKeeper, GeNorm and NormFinder.
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 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 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.
PCR primers used in this study.
GenBank accession no.
Amplicon Size [bp]
PCR efficiency [%]
c-abl oncogene 1
NM_001001548.2, NM_001001547.2, NM_000072.3, NM_001127443.1, NM_001127444.1
Eukaryotic translation initiation factor 2B2
H3 histone, family 3A
Polymerase (RNA) II (DNA-directed) polypeptide K
Peroxisome proliferator-activated receptor α
Peroxisome proliferator-activated receptor δ
Peroxisome proliferator-activated receptor γ
NM_138712.3, NM_015869.4, NM_138711.3, NM_005037.5
Peptidylprolyl isomerase B (cyclophilin B)
Proteasome subunit, β type 2
Proteasome subunit, β type 6
Ribosomal protein L37a
Signal recognition particle 14 kDa
NM_006009.2, NM_006082.2, NM_032704.3
Thioredoxin reductase 1
NM_003330.2, NM_182742.1, NM_182729.1, NM_182743.1, NM_001093771.1
Ubiquitin-conjugating enzyme E2D 2
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.
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.
- Lewis CE, McGee JO: The Macrophage. 2003, Oxford: Oxford University Press
- Ricardo SD, van Goor H, Eddy AA: Macrophage diversity in renal injury and repair. J Clin Invest. 2008, 118: 3522-3530. 10.1172/JCI36150View ArticlePubMedPubMed Central
- Lee WY, Chin AC, Voss S, Parkos CA: In vitro neutrophil transepithelial migration. Methods Mol Biol. 2006, 341: 205-215.PubMed
- Schnoor M, Cullen P, Lorkowski J, Stolle K, Robenek H, Troyer D, Rauterberg J, Lorkowski S: Production of type VI collagen by human macrophages: a new dimension in macrophage functional heterogeneity. J Immunol. 2008, 180: 5707-5719.View ArticlePubMed
- Lorkowski S, Kratz M, Wenner C, Schmidt R, Weitkamp B, Fobker M, Reinhardt J, Rauterberg J, Galinski EA, Cullen P: Expression of the ATP-binding cassette transporter gene ABCG1 (ABC8) in Tangier disease. Biochem Biophys Res Commun. 2001, 283: 821-830. 10.1006/bbrc.2001.4863View ArticlePubMed
- Auwerx J: The human leukemia cell line, THP-1: a multifaceted model for the study of monocyte-macrophage differentiation. Experientia. 1991, 47: 22-31. 10.1007/BF02041244View ArticlePubMed
- Cassol E, Alfano M, Biswas P, Poli G: Monocyte-derived macrophages and myeloid cell lines as targets of HIV-1 replication and persistence. J Leukoc Biol. 2006, 80: 1018-1030. 10.1189/jlb.0306150View ArticlePubMed
- Preiss S, Namgaladze D, Brune B: Critical role for classical PKC in activating Akt by phospholipase A2-modified LDL in monocytic cells. Cardiovasc Res. 2007, 73: 833-840. 10.1016/j.cardiores.2006.12.019View ArticlePubMed
- Tsuchiya S, Kobayashi Y, Goto Y, Okumura H, Nakae S, Konno T, Tada K: Induction of maturation in cultured human monocytic leukemia cells by a phorbol diester. Cancer Res. 1982, 42: 1530-1536.PubMed
- Sokol RJ, Hudson G, James NT, Frost IJ, Wales J: Human macrophage development: a morphometric study. J Anat. 1987, 151: 27-35.PubMedPubMed Central
- Sokol RJ, Hudson G, Wales JM, Goldstein DJ, James NT: Quantitative enzyme cytochemistry during human macrophages development. J Anat. 1993, 183: 97-101.PubMedPubMed Central
- Chang ZL: Recent development of the mononuclear phagocyte system: In memory of Metchnikoff and Ehrlich on the 100th Anniversary of the 1908 Nobel Prize in Physiology or Medicine. Biol Cell. 2009, 101: 709-21. 10.1042/BC20080227View ArticlePubMed
- 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. 10.1158/0008-5472.CAN-04-0496View Article
- Pfaffl MW, Tichopad A, Neuvians P, Neuvians TP: Determination of stable housekeeping genes, differentially regulated genes and sample integrity: BestKeeper - Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47View ArticlePubMed
- 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. 10.1186/gb-2002-3-7-research0034. 10.1186/gb-2002-3-7-research0034View Article
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrin. 2000, 25: 169-193. 10.1677/jme.0.0250169. 10.1677/jme.0.0250169View Article
- Thellin O, Zorzi W, Lakaye B, De Bormann B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotech. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7. 10.1016/S0168-1656(99)00163-7View Article
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE Guidelines: M inimum I nformation for Publication of Q uantitative Real-Time PCR E xperiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797View ArticlePubMed
- Warrington JA, Nair A, Mahadevappa M, Tsyganskaya M: Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiol Genomics. 2000, 2: 143-147.PubMed
- Zampieri M, Ciccarone F, Gustafierro T, Bacalini MG, Calabrese R, Moreno-Villanueva M, Reale A, Chevanne M, Bürkle A, Caiafa : Validation of suitable internal control genes for expression studies in aging. Mech Ageing Dev. 2010, 131: 89-95. 10.1016/j.mad.2009.12.005View ArticlePubMed
- Yin R, Tian F, Frankenberger B, de Angelis MH, Stoeger T: Selection and evaluation of stable housekeeping genes for gene expression normalization in carbon nanoparticle-induced mice acute pulmonary inflammation. Biochem Biophys Res Commun. 2010, 399: 531-536. 10.1016/j.bbrc.2010.07.104View ArticlePubMed
- Labouyrie E, Dubus P, Groppi A, Mahon FX, Ferrer J, Parrens M, Reiffers J, de Mascarel A, Merlio JP: Expression of Neurotrophins and their receptors in human bone marrow. Am J Pathol. 1999, 154: 405-415.View ArticlePubMedPubMed Central
- Chinetti G, Fruchart J-C, Staels B: Peroxisome proliferator-activated receptors: new targets for the pharmacological modulation of macrophage gene expression and function. Curr Opin Lipidol. 2003, 14: 459-468. 10.1097/00041433-200310000-00006View ArticlePubMed
- Hayden JM, Brachova L, Higgins K, Obermiller L, Sevanian A, Khandrika S, Reaven PD: Induction of monocyte differentiation and foam cell formation in vitro by 7-ketocholesterol. J Lipid Res. 2002, 43: 26-35.PubMed
- Collins RF, Touret N, Kuwata H, Tandon NN, Grinstein S, Trimble WS: Uptake of oxidized low density lipoprotein by CD36 occurs by an Actin-dependent pathway distinct from macropinocytosis. J Biol Chem. 2009, 284: 30288-30297. 10.1074/jbc.M109.045104View ArticlePubMedPubMed Central
- Zhou J, Febbraio M, Wada T, Zhai Y, Kuruba R, He J, Lee JH, Khadem S, Ren S, Li S, Silverstein RL, Xie W: Hepatic fatty acid transporter Cd36 is a common target of LXR, PXR and PPARγ in promoting steatosis. Gastroenterology. 2008, 134: 556-567. 10.1053/j.gastro.2007.11.037View ArticlePubMed
- Pilbrow AP, Ellmers LJ, Black MA, Moravec CS, Sweet WE, Troughton RW, Richards AM, Frampton CM, Cameron VA: Genomic selection of reference genes for real-time PCR in human myocardium. BMC Med Genomics. 2008, 1: 64- 10.1186/1755-8794-1-64View ArticlePubMedPubMed Central
- Stern-Straeter J, Bonaterra GA, Hörmann K, Kinscherf R, Goessler UR: Identification of valid reference genes during differentiation of human myoblasts. BMC Mol Biol. 2009, 10: 66- 10.1186/1471-2199-10-66View ArticlePubMedPubMed Central
- Fu L-Y, Jia H-L, Dong Q-Z, Wu J-C, Zhao Y, Zhou H-J, Ren N, Ye Q-H, Qin L-X: Suitable reference genes for real-time PCR in human HBV-related hepatocellular carcinoma with different clinical prognoses. BMC Cancer. 2009, 9: 49- 10.1186/1471-2407-9-49View ArticlePubMedPubMed Central
- Tong Z, Gao Z, Wang F, Zhou J, Zhang Z: Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol Biol. 2009, 10: 71- 10.1186/1471-2199-10-71View ArticlePubMedPubMed Central
- Piehler AP, Grimholt RM, Øvstebø R, Berg JP: Gene expression results in lipopolysaccharide-stimulated monocytes depend significantly on the choice of reference genes. BMC Immunol. 2010, 11: 21- 10.1186/1471-2172-11-21View ArticlePubMedPubMed Central
- Liu H, Shi B, Huang C-C, Eksarko P, Pope RM: Transcriptional diversity during monocyte to macrophage differentiation. Immunol Lett. 2008, 117: 70-80. 10.1016/j.imlet.2007.12.012View ArticlePubMedPubMed Central
- Mauche S, Poirier O, Godefroy T, Olaso R, Gut I, Collet J-P, Montalescot G, Cambien F: Performance comparison of two microarray platforms to assess differential gene expression in human monocyte and macrophage cells. BMC Genomics. 2008, 9: 302- 10.1186/1471-2164-9-302View Article
- Martinez FO, Gordon S, Locati M, Mantovani A: Transcriptional profiling of the human monocyte-to-macrophage differentiation and polarization: New molecules and patterns of gene expression. J Immunol. 2006, 177: 7303-7311.View ArticlePubMed
- Schnoor M, Buers I, Sietmann A, Brodde MF, Hofnagel O, Robenek H, Lorkowski S: Efficient non-viral transfection of THP-1 cells. J Immunol Methods. 2009, 344: 109-115. 10.1016/j.jim.2009.03.014View ArticlePubMed
- Stolle K, Schnoor M, Fuellen G, Spitzer M, Cullen P, Lorkowski S: Cloning, genomic organization, and tissue-specific expression of the RASL11B gene. Biochim Biophys Acta. 2007, 1769: 514-524.View ArticlePubMed
- Von Ahlfen S, Schlumpberger M: Effects of low A 260 /A 230 ratios in RNA preparations on downstream applications. QIAGEN Gene Expression Newsletter. 2010, 15: 6-7.
- Stolle K, Schnoor M, Fuellen G, Spitzer M, Engel T, Spener F, Cullen P, Lorkowski S: Cloning, cellular localization, genomic organization, and tissue-specific expression of the TGFbeta1-inducible SMAP-5 gene. Gene. 2005, 351: 119-130. 10.1016/j.gene.2005.03.012View ArticlePubMed
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.