- Research article
- Open Access
Evaluation of glyceraldehyde-3-phosphate, prolylpeptidyl isomerase A, and a set of stably expressed genes as reference mRNAs in urate crystal inflammation
© Pessler et al; licensee BioMed Central Ltd. 2011
- Received: 15 June 2011
- Accepted: 25 October 2011
- Published: 25 October 2011
The murine air pouch membrane represents an easily accessible tissue for studies on gene regulation in acute inflammation. Considering that acute inflammation may affect expression of molecular reference genes, we evaluated the expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and prolylpeptidyl isomerase A (PPIA) in the air pouch membrane during a complete time course of urate crystal inflammation and correlated the results with expression of interleukin (IL)-1β and hypoxia inducible factor (HIF)-1α. In addition, we aimed to identify alternate potential reference genes.
Using custom microfluidic real-time PCR arrays, the expression of 96 genes including GAPDH, PPIA, IL-1β, and HIF-1α was determined in dissected air pouch membranes 1, 4, 9, 18, 27, and 50 hours (h) after injecting monosodium urate (MSU) crystals into the pouch. One-way ANOVA was used to detect differential gene expression throughout the time course. Using the genes on these arrays as a convenience sample, alternate candidate reference genes were sought (1) with a biostatistical approach and (2) using the geNorm software tool.
Pouch leukocytes peaked at t = 9h and declined toward t = 50h. PPIA expression was not differentially regulated (p = 0.52, ANOVA). In contrast, GAPDH mRNA increased steadily after crystal injection, reaching a maximal 2.8-fold increase at t = 18h (p = 0.0006, t test), which followed a marked induction of IL-1β (max., 208-fold at t = 4h, p = 8.4 × 10-5, t test) and HIF-1α (max., 6.6-fold at t = 4h, p = 0.00025, t test). Fifteen genes were artifactually identified as "significantly regulated" when Ct values were normalized against GAPDH expression. The biostatistical approach and the geNorm analysis identified overlapping sets of candidate reference genes. Both ranked PPIA as the best candidate, followed by defender against cell death 1 (DAD1) and high-mobility group B1 (HMGB1).
GAPDH mRNA expression is up-regulated in urate crystal inflammation, possibly due to inflammation-associated hypoxia. Using GAPDH mRNA for molecular normalization resulted in significant artifacts in the calculated expression of the target mRNAs. PPIA and other stably expressed genes promise to be more appropriate reference genes in this model.
- Reference Gene
- GAPDH mRNA
- Candidate Reference Gene
- Median Absolute Deviation
- GAPDH Expression
Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and prolylpeptidyl isomerase A (PPIA) as potential reference genes
GAPDH is often used for molecular normalization of gene expression data from microarrays or real-time reverse transcriptase polymerase chain reactions (qPCR). This is based on the assumption that expression of this "housekeeping gene" does not change much during the life cycle of most cells and can thus be used as a relatively constant reference signal. However, while this notion has been substantiated in some scenarios, there are clear examples that GAPDH mRNA expression can vary (e.g., [1–3]). Notably, hypoxia can induce GAPDH mRNA levels, likely because binding of a complex of the inducible α subunit and the constitutively expressed β subunit of hypoxia inducible factor (HIF)-1 to a hypoxia response element (HRE) in the GAPDH promoter region can increase transcription of this gene [4–6]. Considering that hypoxia is a well documented feature of inflammatory cells and inflamed tissues [7, 8], including the synovial membrane , that HIF-1 can be activated in inflammation due to toll-like receptor (TLR) signaling [8, 10], and that HIF-1α is expressed widely in inflamed synovial membranes , GAPDH may not be a suitable reference gene for molecular normalization in gene expression studies of acute synovitis, including crystal inflammation.
PPIA (also known as cyclophilin A) is a ubiquitously expressed intermediate factor of calcium/calmodulin signaling. Its activity is regulated predominantly at the post-transcriptional level. While it has been validated as a useful reference gene for qPCR in specific scenarios [12, 13], its transcriptional regulation has been demonstrated in hypoxic cells  and in at least one example of a chronically inflamed tissue . Thus, there is reason to suspect that it, too, may be a suboptimal reference gene in studies on inflammation.
The murine air pouch model of inflammation
The murine air pouch is a bursa-like structure that is lined with a membrane resembling the lining of human joints, both histologically and biochemically . The pouch lumen is an easily accessible space, and different inflammatory processes can be elicited easily by injecting the respective pro-inflammatory agent. The air pouch membrane can be removed from the mouse nearly quantitatively by blunt dissection and thus provides an attractive system for studying inflammation-related gene expression changes in a synovium-like tissue, while minimizing transcriptional noise hailing from the adjacent structures . However, the usefulness of common reference genes for real-time PCR analysis of this tissue has not been examined. Considering the fulminant inflammatory reaction that ensues after injecting the crystals, changes in basic elements of cellular metabolism affecting expression of otherwise stably expressed genes appear likely. Using a time course experiment spanning initiation, peak and resolution of inflammation in the air pouch, we have therefore evaluated GAPDH and PPIA as reference genes for molecular normalization in this model. We find relatively stable expression of PPIA, but a steady rise of GAPDH mRNA levels that peak after maximal leukocyte accumulation in the pouch lumen, and we speculate that inflammation-associated hypoxia and/or oxidative stress are causative factors.
Murine air pouch model
Quantitative PCR analysis
Membranes were homogenized, total RNA extracted with RNeasy spin columns (Qiagen, Valencia, CA), and cDNA synthesized with Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA), following the recommendations of the manufacturers. Expression of 96 target genes was quantified on an HT7900 Sequence Analyzer (Applied Biosystems [ABI], Foster City, CA) using custom microfluidic Taqman PCR arrays (ABI) following the manufacturer's protocols (http://www.appliedbiosystems.com). The arrays had been custom designed to contain three reference genes (18S rRNA, PPIA, and GAPDH) and 93 target genes relating to several functional pathways and processes (cell cycle, apoptosis, prostaglandin metabolism, inflammasome components, cytokines, defensins, matrix metalloproteinases, markers for specific cell types, and several miscellaneous targets of interest) to study their regulation, if any, in urate crystal inflammation. While the present manuscript focuses on reference gene selection, the analysis of expression of the functional pathways will be the subject of a separate publication.
Pre-analytic data processing
Ten of the 96 genes on the arrays were removed from the analysis because two or more values were missing and this lack of detection could not be explained by a biologically plausible absence of expression at beginning or end of the time course. A Ct of 40 was assigned for missing values when absence of detection was likely due to physiological lack of expression at beginning or end of the time course; the mean Ct value of the time point was used when missing values were likely due to technical failure. All Ct values were subtracted from 40 such that higher Ct values corresponded to higher gene expression. To avoid the value 0 for samples with Ct = 40, the Laplace correction (adding a pseudocount of 1 to all values of the data set)  was used, thus obtaining (41-Ct). 18S rRNA, which is routinely included on these custom PCR arrays by the manufacturer, had to be excluded from the analysis because a batch effect was noted in principal component analysis (results not shown), which was due to low detection levels for this gene on three of the arrays. The Ct values of 18S rRNA did not correlate with mean Ct values across the arrays. Therefore, they did not reflect low overall gene expression in the affected samples, but were likely due to a manufacturing error. Eighty-five of the original 96 genes thus passed the preanalytic quality screen and were included in the analyses.
Differential regulation in the time course of target genes of interest was determined with one-way ANOVA and the Kruskal-Wallis test, comparing t = 0h against the subsequent time points. P values were adjusted for multiple testing using the Benjamini-Hochberg (BH) correction method. Similar results were obtained with the two methods, but only p values obtained with ANOVA are shown for the sake of simplicity. For each target, a t test was performed between t = 0h and the time point associated with the greatest change in expression. Variation in gene expression was also assessed by median absolute deviation from the median (MAD) . MAD is a more robust measure of data variability than standard deviation, which is obtained by multiplying MAD by a constant scaling factor. All statistical analyses were performed using code written in the R programming language (R project for statistical computing; http://www.r-project.org). The software tool geNorm  (available at http://www.biogazelle.com/genormplus/genorm) was also used to identify potential reference genes. It ranks a group of candidate reference genes according to the M value, i.e. the mean of all the standard deviation values between the Ct values of the gene in question and the Ct values of each of the other potential reference genes under consideration. A lower M value denotes a higher rank among the candidates.
Up-regulation of GAPDH, but not PPIA, mRNA in urate crystal inflammation
To exclude that the rise in GAPDH mRNA merely reflected a general induction of transcription at these time points or systematic differences in RNA or cDNA amount or quality, the mean Ct value across all 85 targets and samples was calculated for each time point. This overall mean transcription peaked at t = 4h and thus clearly followed different kinetics than regulation of GAPDH.
Artifacts resulting from normalization against GAPDH expression
Identification of alternate candidate reference mRNAs
In an independent approach, the same set of 23 non-regulated genes was analyzed with the software program geNorm. Thirteen genes had M values < 0.5, a generally accepted cut-off for the identification of potentially useful reference genes  (Figure 4B). Of note, this unbiased screen identified the same top three genes as the MAD/ANOVA approach, PPIA (M = 0.23) being the best candidate, followed by DAD1 (M = 0.25) and HMGB1 (M = 0.26).
Inutility of GAPDH as reference gene in MSU crystal inflammation
We detected significant up-regulation of GAPDH expression in the anticlimactic phase of acute urate crystal inflammation in this simple tissue resembling the human synovial membrane. These results suggest strongly that GAPDH should not be used as a reference gene in urate crystal inflammation and, likely, similar forms of acute inflammation driven by innate immunity. It remains to be determined whether this also applies to chronic or less fulminant inflammatory processes, but caution is indicated. In any case, these results underscore how important it is to test for any systematic regulation of GAPDH mRNA expression before considering it for molecular normalization in a new experimental system.
Potential mechanisms of GAPDH induction in MSU crystal inflammation
What might be the underlying mechanism that increased transcription of this enzyme? The presented data are the first report of HIF-1α induction in urate crystal inflammation, but the available evidence from related inflammatory systems allows us to speculate that the increased GAPDH expression is due to increased HIF-1α activity. This might result from at least two synergizing processes, (1) higher HIF-1α protein levels due to increased stability resulting from proline hydroxylation (which could not be assessed in our mRNA-based analysis), and (2) increased transcription of the HIF-1α gene within the first hours of crystal injection (Figure 2). These changes may arise from activation of more than one molecular pathway. Urate crystals induce fulminant inflammation in the air pouch that is characterized by massive induction of innate immunity. The crystals are sensed as danger-associated molecular patterns (DAMPS) by the NALP3 inflammasome (leading to the proteolytic processing and release of IL-1β), and by other receptors including TLR 2 and 4 . Notably, TLR 2 and 4 (as well as TLR 7 and 8, but these are not known to be activated by urate crystals) have been shown to activate HIF-1α activity [8, 10]. This can be seen as part of the intracellular adaptation to hypoxic and oxidative stress that infiltrating inflammatory cells are invariably exposed to. Injecting urate crystals into the air pouch causes a major influx of innate immune cells, consisting of neutrophilic granulocytes, monocytes (which express TLR 2 and 4), and mast cells  (whose activation also leads to increased HIF-1α activity ), all of which ultimately results in the release of major proinflammatory cytokines and reactive oxygen species (ROS). Moreover, focal hypoxia is a common feature of inflamed tissues , including the synovial membrane , likely due to increased oxygen demand in the tissue in the presence of vascular congestion. A recent report in this journal on expression of GAPDH in tumor and non-tumor cell lines clearly demonstrated that GAPDH expression increases under hypoxic conditions in a cell-type specific manner , likely due to HIF-1α induction. GAPDH is a key enzyme of glycolysis. A HIF-1α-induced increase in GAPDH activity would activate this pathway, thus improving the intracellular energy balance by providing additional ATP molecules. In addition, intracellular ROS can induce HIF-1α activity . Besides these intracellular ROS, which serve as signaling molecules in innate immune cells during activation, neutrophil and mast cell degranulation leads to additional oxidative stress in the tissue, and increased GAPDH activity may increase survival of affected cells by increasing the availability of ATP and reducing equivalents. Thus, the observed induction of GAPDH in the presented study may be explained biochemically in two, not mutually exclusive, ways: (1) as an inevitable result of TLR signal transduction triggered by urate crystals and (2) as a central feature of the tissue adaptation to hypoxia and oxidative stress. Future studies should be directed at testing these models experimentally.
PPIA, DAD1 and HMGB1 as potential reference genes for urate crystal inflammation
The MAD and geNorm analyses highlighted several stably expressed genes, and both of these unbiased approaches identified PPIA as the best candidate reference gene. PPIA has been validated for molecular normalization in qPCR analysis of, for instance, colorectal and breast cancer tissue [12, 13]. It now promises to be a useful reference gene also in urate crystal inflammation and, possibly, similar states of inflammation. DAD1 is a subunit of the enzyme oligosaccharyltransferase and can counteract apoptotic cell death . Interestingly, MSU crystals inhibit neutrophil apoptosis . While it is not known whether this process involves DAD1 activity, our results indicate that the crystals at least do not affect transcription of the DAD1 gene. HMGB1 is a multifunctional protein that usually resides in the nucleus and functions in DNA repair . But it also plays important extracellular roles in innate immunity as an "alarmin", where it is released from cells through activation of the very NALP3 inflammasome that is activated by MSU crystals . Although this has not been addressed experimentally, induction of HMGB1 activity in response to MSU crystal is highly likely. Its identification as a stably expressed gene in MSU crystal inflammation was, therefore, unexpected. A likely explanation is that its activity is regulated predominantly, if not exclusively, at the post-transcriptional level in this type of acute but self-limiting process.
Using the geometric mean of a combination of two or more reference genes has become a commonly recommended strategy for molecular normalization of RT-PCR data sets . The geNorm analysis suggested that the combination of PPIA, DAD1 and HMGB1 would constitute the optimal normalization strategy in the current data set. Using this combination affected the number of genes identified as "significantly regulated" in a relatively minor way, which might be viewed as "fine tuning" of normalization against PPIA alone. Additional changes due to using this combination may be more subtle and, for instance, affect the measured kinetics of some of the regulated genes.
This study is limited by the fact that reference genes that are commonly used in inflammation-related gene expression studies, notably β-actin and 18S and 28S rRNA, were not included in the analysis and that no attempt was made to elucidate the mechanisms underlying GAPDH induction or to verify expression of the protein products of the studied mRNAs. Moreover, DAD1 and HMGB1 and the other stably expressed genes other than PPIA were included on the arrays a priori not as potential reference genes but in order to test whether their expression might be functionally regulated in urate crystal inflammation. However, our results do suggest that a post hoc data mining approach of PCR arrays may be useful to screen for additional candidate reference genes. Clearly, future studies of urate crystal inflammation should be designed to (1) further validate PPIA, DAD1 and HMGB1 as reference genes and to compare their performance directly with that of commonly used reference genes and (2) to provide more experimental evidence for the hypothetical model of HIF-1α-associated GAPDH induction outlined at the beginning of this section.
GAPDH mRNA expression is up-regulated in urate crystal inflammation, possibly due to inflammation-associated hypoxia and/or oxidative stress. Using GAPDH mRNA for molecular normalization resulted in significant artifacts in the calculated expression of the target mRNAs, and it should therefore not be used as a reference mRNA in studies on MSU crystal inflammation. PPIA and other stably expressed genes, including DAD1 and HMGB1, promise to be more appropriate reference genes in this model.
We thank Nicole Chartrain and Marlin Camps-Ramirez (Merck Research Laboratories, Rahway, NJ) for expert technical assistance, and Manas Akmatov (Helmholtz Centre for Infection Research) for helpful discussion and a critical reading of the manuscript. We also thank Christian T. Mayer (Twincore GmbH, Hannover, Germany) and John Tobias (University of Pennsylvania, Philadelphia, PA) for their help in the early phases of the project. The study was supported by funds from the Philadelphia VA Medical Center Research Foundation, by internal funds from the Helmholtz Association (Program Infection and Immunity) to FP, and by European Regional Development Funds from NBank (Niedersächsische Investitions- und Förderbank) to CDB.
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29: 332-337.PubMedGoogle Scholar
- Said HM, Polat B, Hagemann C, Anacker J, Flentje M, Vordermark D: Absence of GAPDH regulation in tumor-cells of different origin under hypoxic conditions in - vitro. BMC Res Notes. 2009, 2: 8-10.1186/1756-0500-2-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Fink T, Lund P, Pilgaard L, Rasmussen JG, Duroux M, Zachar V: Instability of standard PCR reference genes in adipose-derived stem cells during propagation, differentiation and hypoxic exposure. BMC Mol Biol. 2008, 9: 98-10.1186/1471-2199-9-98.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhong H, Simons JW: Direct comparison of GAPDH, beta-actin, cyclophilin, and 28S rRNA as internal standards for quantifying RNA levels under hypoxia. Biochem Biophys Res Commun. 1999, 259: 523-526. 10.1006/bbrc.1999.0815.PubMedView ArticleGoogle Scholar
- Graven KK, Yu Q, Pan D, Roncarati JS, Farber HW: Identification of an oxygen responsive enhancer element in the glyceraldehyde-3-phosphate dehydrogenase gene. Biochim Biophys Acta. 1999, 1447: 208-218.PubMedView ArticleGoogle Scholar
- Higashimura Y, Nakajima Y, Yamaji R, Harada N, Shibasaki F, Nakano Y, Inui H: Up-regulation of glyceraldehyde-3-phosphate dehydrogenase gene expression by HIF-1 activity depending on Sp1 in hypoxic breast cancer cells. Arch Biochem Biophys. 2011, 509: 1-8. 10.1016/j.abb.2011.02.011.PubMedView ArticleGoogle Scholar
- Eltzschig HK, Carmeliet P: Hypoxia and inflammation. N Engl J Med. 2011, 364: 656-665. 10.1056/NEJMra0910283.PubMedPubMed CentralView ArticleGoogle Scholar
- Sumbayev VV, Nicholas SA: Hypoxia-inducible factor 1 as one of the "signaling drivers" of Toll-like receptor-dependent and allergic inflammation. Arch Immunol Ther Exp (Warsz). 2010, 58: 287-294. 10.1007/s00005-010-0083-0.View ArticleGoogle Scholar
- Ng CT, Biniecka M, Kennedy A, McCormick J, Fitzgerald O, Bresnihan B, Buggy D, Taylor CT, O'Sullivan J, Fearon U, Veale DJ: Synovial tissue hypoxia and inflammation in vivo. Ann Rheum Dis. 2010, 69: 1389-1395. 10.1136/ard.2009.119776.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicholas SA, Sumbayev VV: The involvement of hypoxia-inducible factor 1 alpha in Toll-like receptor 7/8-mediated inflammatory response. Cell Res. 2009, 19: 973-983. 10.1038/cr.2009.44.PubMedView ArticleGoogle Scholar
- Brouwer E, Gouw AS, Posthumus MD, van Leeuwen MA, Boerboom AL, Bijzet J, Bos R, Limburg PC, Kallenberg CG, Westra J: Hypoxia inducible factor-1-alpha (HIF-1alpha) is related to both angiogenesis and inflammation in rheumatoid arthritis. Clin Exp Rheumatol. 2009, 27: 945-951.PubMedGoogle Scholar
- Kheirelseid EA, Chang KH, Newell J, Kerin MJ, Miller N: Identification of endogenous control genes for normalisation of real-time quantitative PCR data in colorectal cancer. BMC Mol Biol. 2010, 11: 12-10.1186/1471-2199-11-12.PubMedPubMed CentralView ArticleGoogle Scholar
- McNeill RE, Miller N, Kerin MJ: Evaluation and validation of candidate endogenous control genes for real-time quantitative PCR studies of breast cancer. BMC Mol Biol. 2007, 8: 107-10.1186/1471-2199-8-107.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim H, Kim WJ, Jeon ST, Koh EM, Cha HS, Ahn KS, Lee WH: Cyclophilin A may contribute to the inflammatory processes in rheumatoid arthritis through induction of matrix degrading enzymes and inflammatory cytokines from macrophages. Clin Immunol. 2005, 116: 217-224. 10.1016/j.clim.2005.05.004.PubMedView ArticleGoogle Scholar
- Edwards JC, Sedgwick AD, Willoughby DA: The formation of a structure with the features of synovial lining by subcutaneous injection of air: an in vivo tissue culture system. J Pathol. 1981, 134: 147-156. 10.1002/path.1711340205.PubMedView ArticleGoogle Scholar
- Pessler F, Mayer CT, Jung SM, Behrens EM, Dai L, Menetski JP, Schumacher HR: Identification of novel monosodium urate crystal regulated mRNAs by transcript profiling of dissected murine air pouch membranes. Arthritis Res Ther. 2008, 10: R64-10.1186/ar2435.PubMedPubMed CentralView ArticleGoogle Scholar
- Jung SM, Schumacher HR, Kim H, Kim M, Lee SH, Pessler F: Reduction of urate crystal-induced inflammation by root extracts from traditional oriental medicinal plants: elevation of prostaglandin D2 levels. Arthritis Res Ther. 2007, 9: R64-10.1186/ar2222.PubMedPubMed CentralView ArticleGoogle Scholar
- McCarty DJ, Faires JS: A comparison of the duration of local anti-inflammatory effect of several adrenocorticosteroid esters--a bioassay technique. Curr Ther Res Clin Exp. 1963, 5: 284-290.PubMedGoogle Scholar
- Berthold Michael, Christian Borgelt, Frank Höppner, Klawonn F: Descriptive Statistics. Guide to Intelligent Data Analysis: How to Intelligently Make Sense of Real Data. 2010, Springer, 315-1View ArticleGoogle Scholar
- Hoaglin DC: Understanding Robust and Exploratory Data Analysis. 1983, John Wiley & SonsGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8: R19-10.1186/gb-2007-8-2-r19.PubMedPubMed CentralView ArticleGoogle Scholar
- Martinon F: Mechanisms of uric acid crystal-mediated autoinflammation. Immunol Rev. 2010, 233: 218-232. 10.1111/j.0105-2896.2009.00860.x.PubMedView ArticleGoogle Scholar
- Schiltz C, Lioté F, Prudhommeaux F, Meunier A, Champy R, Callebert J, Bardin T: Monosodium urate monohydrate crystal-induced inflammation in vivo: quantitative histomorphometric analysis of cellular events. Arthritis Rheum. 2002, 46: 1643-1650. 10.1002/art.10326.PubMedView ArticleGoogle Scholar
- Walczak-Drzewiecka A, Ratajewski M, Wagner W, Dastych J: HIF-1alpha is up-regulated in activated mast cells by a process that involves calcineurin and NFAT. J Immunol. 2008, 181: 1665-1672.PubMedView ArticleGoogle Scholar
- Kelleher DJ, Gilmore R: DAD1, the defender against apoptotic cell death, is a subunit of the mammalian oligosaccharyltransferase. Proc Natl Acad Sci USA. 1997, 94: 4994-4999. 10.1073/pnas.94.10.4994.PubMedPubMed CentralView ArticleGoogle Scholar
- Tudan C, Fong D, Duronio V, Burt HM, Jackson JK: The inhibition of spontaneous and tumor necrosis factor-alpha induced neutrophil apoptosis by crystals of calcium pyrophosphate dihydrate and monosodium urate monohydrate. J Rheumatol. 2000, 27: 2463-2472.PubMedGoogle Scholar
- Liu Y, Prasad R, Wilson SH: HMGB1: roles in base excision repair and related function. Biochim Biophys Acta. 2010, 1799: 119-130.PubMedPubMed CentralView ArticleGoogle Scholar
- Vande Walle L, Kanneganti TD, Lamkanfi M: HMGB1 release by inflammasomes. Virulence. 2011, 2: 162-165. 10.4161/viru.2.2.15480.PubMedPubMed CentralView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.