Multiplex transcriptional analysis of paraffin-embedded liver needle biopsy from patients with liver fibrosis
© Staten et al.; licensee BioMed Central Ltd. 2012
Received: 18 October 2012
Accepted: 11 December 2012
Published: 27 December 2012
The possibility of extracting RNA and measuring RNA expression from paraffin sections can allow extensive investigations on stored paraffin samples obtained from diseased livers and could help with studies of the natural history of liver fibrosis and inflammation, and in particular, correlate basic mechanisms to clinical outcomes.
To address this issue, a pilot study of multiplex gene expression using branched-chain DNA technology was conducted to directly measure mRNA expression in formalin-fixed paraffin-embedded needle biopsy samples of human liver. Twenty-five genes were selected for evaluation based on evidence obtained from human fibrotic liver, a rat BDL model and in vitro cultures of immortalized human hepatic stellate cells. The expression levels of these 25 genes were then correlated with liver fibrosis and inflammation activity scores. Statistical analysis revealed that three genes (COL3A1, KRT18, and TUBB) could separate fibrotic from non-fibrotic samples and that the expression of ten genes (ANXA2, TIMP1, CTGF, COL4A1, KRT18, COL1A1, COL3A1, ACTA2, TGFB1, LOXL2) were positively correlated with the level of liver inflammation activity.
This is the first report describing this multiplex technique for liver fibrosis and has provided the proof of concept of the suitability of RNA extracted from paraffin sections for investigating the modulation of a panel of proinflammatory and profibrogenic genes. This pilot study suggests that this technique will allow extensive investigations on paraffin samples from diseased livers and possibly from any other tissue. Using identical or other genes, this multiplex expression technique could be applied to samples obtained from extensive patient cohorts with stored paraffin samples in order to correlate gene expression with valuable clinically relevant information. This method could be used to provide a better understanding of the mechanisms of liver fibrosis and inflammation, its progression, and help development of new therapeutic approaches for this indication.
KeywordsLiver fibrosis Liver inflammation Multiplex gene expression
Conducting gene expression analysis at the level of the organ of interest and correlating specific expression patterns to physiopathology, prognosis or early response to therapy remains a topic under intense investigation in liver fibrosis . The ability to analyze gene expression from formalin-fixed paraffin-embedded (FFPE) tissue would be logistically useful because of access to banked samples which are usually paraffin-embedded specimens. For example, the description of this approach for hepatocellular carcinoma on paraffin-fixed tissue  has underscored the importance of the technology.
Here, we describe a multiplex gene expression pilot study using the bDNA technology QuantiGene™ Reagent System (Affymetrix, Santa Clara, CA, USA) for quantifying mRNA expression in FFPE liver samples . Using this single multiplex assay, the expression levels of 25 genes were correlated with the extent of disease activity in fibrogenic chronic liver disorders. Included in the panel were genes known to be associated with liver fibrosis and hepatic stellate cell (HSC) activation, which have both been shown to correlate with inflammation . In addition, several putative HSC housekeeping genes were chosen to assess for a correlation of HSC activity to stage of fibrosis or inflammation. Proof of concept was established showing that RNA extracted from paraffin sections could be used to investigate the modulation of a panel of proinflammatory and profibrogenic genes with this multiplex assay.
Identification of gene expression panel
A 30-plex gene expression panel, including five human gene normalization controls and the two human ribosomal (DNA and RNA) controls used in singleplex
18S ribosomal DNA
Tissue Quality & Quantity Assessment
28S ribosomal RNA
Peptidylprolyl isomerase A (cyclophilin A)
Ribosomal protein L13a
TATA box binding protein
Actin, alpha 2, smooth muscle, aorta
Fibrosis and/or HSC Activation Marker
Aryl-hydrocarbon receptor nuclear translocator 2
Collagen, type I, alpha 1
Fibrosis and/or HSC Activation Marker
Collagen, type III, alpha 1
Collagen, type IV, alpha 1
Cytoplasmic polyadenylation element binding protein 1
Connective tissue growth factor
Fibrosis and/or HSC Activation Marker
Collagen triple helix repeat containing 1
Chemokine (C-X-C motif) ligand 2 ; CINC-2a, GROb, Gro2, MIP-2, MIP-2a
Forkhead box F1
Hepcidin antimicrobial peptide
Fibrosis and/or HSC Activation Marker
Heat shock 70 kDa protein 1A
Insulin-like growth factor binding protein 6
Keratin 18; cytokeratin 18, cell proliferation inducing protein 46
Fibrosis and/or HSC Activation Marker
Lysyl oxidase-like 2
MAS-related GPR, member F
Neuronal protein 3.1, C5orf13
Fibrosis and/or HSC Activation Marker
Transforming growth factor, beta 1
Transforming growth factor, beta 2
TIMP metallopeptidase inhibitor 1
Tumor necrosis factor (TNF superfamily, member 2)
HSC Housekeeping & Fibrosis Marker
HSC, hepatic stellate cells
Performance of the multiplex assay with FFPE liver samples
The technique was first standardized on cadaver liver tissue blocks. Five sections (60 mm2) of 8-μm thick sections from cadaver liver tissue blocks resulted in the most reproducible results. The 30-plex assay was then used to assess gene expression in all 25 needle biopsies. The optimum quantity of liver tissue for best sensitivity and reproducibility was equivalent to two to three needle biopsy sections (approximately 20 to 30 mm2) 8-μm thick for detection of all 30 genes tested. QuantiGene™ 2.0 multiplex intraassay coefficient of variations (CVs) were less than 10% for duplicate sections and less than 15% for triplicate sections. The interassay CVs were less than 20% for technical replicates. Data for four genes from the original set of 25 were discarded due to weak signal intensities across the majority of samples (CXCL2 and SYP) or in more than 50% of samples (ARNT2 and CPEB1). In addition, one of the QuantiGene™ assay normalization genes, UBC, was not included for further analysis due to inconsistent behavior across multiple strip testing.
Six genes were identified with significantly upregulated expression in fibrotic liver: LOXL2, TUBB, COL1A1, TBP (all P <0.0001), COL3A1 (P = 0.0002), and P311 (P = 0.0008). One gene showed significantly downregulated expression in fibrotic liver: KRT18 (P = 0.0002). Among these seven differentially expressed genes, TBP was originally selected as a gene for assay normalization.
Development of a logistic model for assessing liver fibrosis progression
A logistic model was also established to determine those informative genes that could provide insight into liver fibrosis progression . The disease status as a binary response was modeled on the covariates, that is, the gene expression measurements. First, each gene was evaluated separately by fitting the simple logistic model. From this preliminary one-gene-at-a-time modeling step, nine out of 25 genes were selected, including the seven genes listed above along with two additional genes: ACTB and PPIA. A logistical model was then built by utilizing the forward selection method based on the nine genes from the previous step using the SAS PROC LOGISTIC procedure . Comparison of the modeling results with the t-test, principle component analysis (PCA), and partial least square (PLS) analyses resulted in the selection of three genes (COL3A1, KRT18, and TUBB,) as potential informative covariates. However, this logistical model needs to be tested with additional data in order to validate this model and establish it as a useful prediction model for liver fibrosis status.
Correlation between gene expression and disease activity
This report describes a pilot study evaluating a multiplex gene analysis assay, relevant for complex diseases with limited tissue availability, which proved reliable for the study of tissue samples of fibrogenic chronic liver diseases using FFPE needle biopsy samples. The optimum quantity of liver tissue for best sensitivity and reproducibility was equivalent to two to three needle biopsy sections 8-μm thick for detection of all genes tested. Using this approach, the expression of ten genes, known to be associated with liver fibrosis and/or HSCs activation, were identified as positively correlating with necroinflammatory activity, whereas two HSC housekeeping genes were identified as negatively correlating with activity, though we could not find a strong correlation between gene expression and fibrosis stage.
We used the bDNA technology QuantiGene™ Reagent System for quantifying mRNA expression. The technical advantage of the sandwich nucleic acid hybridization assay  lies in amplification of reporter signals without enzymatic amplification of the target mRNA, which allows direct measurement of mRNA in tissue homogenates, overcoming issues associated with mRNA extraction. This technology was chosen over a qPCR-based assay method for two primary reasons: 1) QuantiGene™ detection can accurately measure mRNA in FFPE samples whereas qPCR cannot reliably measure mRNA, which has an average size of 100 bases; and 2) QuantiGene™ can be multiplexed more easily than qPCR which is more suited for expression analysis of mRNA extracted from fresh or frozen samples. Prior attempts to employ qPCR on a subset of selected genes in paraffin resulted in no detectable gene expression in paraffin liver RNA and resulted in no detectable gene expression or gene expression beyond the limit of detection (Staten, data not published). Thus, we were not able to use qPCR to compare gene expression between samples or correlate between qPCR and bDNA technologies.
The selection of the genes to be evaluated was made arbitrarily on the basis of a mixed combination of data obtained from human fibrotic liver, a rat BDL model and in vitro culture of immortalized hepatic stellate cells, namely L190 and LX2, and therefore the selected genes are likely related to established inflammation and fibrosis. Other genes could be used in this multiplex assay to address complementary, but different, questions. For example we did not select all of the candidate genes strictly on the basis of a subtraction analysis done by comparing gene expression in fibrotic liver (or activated HSC) and normal liver (or quiescent stellate cells), as described by others previously [9, 10]. This alternative approach may be more appropriate to study the difference of gene expression between disease state and normal liver. However, our approach and the choice of genes used in the present study are more appropriate for investigating the modulation of inflammation and fibrosis between different samples of diseased liver, and could also be applied to verify the modulation of these genes following a period of treatment (for example, predictor of antifibrotic and anti-inflammatory response).
Prognostic factors of rapid progression of liver fibrosis to cirrhosis based on liver biopsy, but also on somatic genetic markers, have been described [11, 12]. For example, Asselah and colleagues  identified 11 genes that differentiate mild from moderate fibrosis based on their respective expression levels. But only two of these genes (COL1A1, TIMP1) overlapped with our study. For identification of potential therapeutic targets, this multiplex assay can be customized with other genes of interest. Liver fibrosis is a disease characterized by collagen deposition  and we included in the multiplex assay multiple genes coding for proteins involved in extra cellular matrix formation, deposition and turnover (COL1A1, COL3A1, COL4A1, ACTA2 and TIMP1); and one gene encoding a component of the cytoskeleton (KRT18). However, in our report, the highest correlation with inflammatory activity stage was observed for annexin 2 (ANXA2). ANXA2 has many described functions but its role in the hepato-biliary system is not fully understood. It is expressed in many cells including hepatocytes, HSCs  and cholangiocytes . ANXA2 is involved in the plasminogen/plasmin system with fibrinolytic effects promoting clot dissolution in alcoholic liver cirrhosis . In inflammatory cholangiopathy ANXA2 can mediate a compensatory effect for the impaired bicarbonate secretion of cholangiocytes through modulation of anion exchanger . ANXA2 has been shown to be expressed as a vitamin D binding protein on hepatocytes and HSCs. Some of these interactions could be important for liver fibrogenesis and inflammation through a TGFβ-independent mechanism .
In conclusion, the results of this pilot study provide feasibility evidence (proof of concept) for the use of multiplex gene expression methods on FFPE samples from chronic fibrotic liver tissues. However, because of possible limitations of our study (choice of genes, limited samples used for analysis, scoring and grading done only with METAVIR), these results are not yet conclusive for the definition of the diagnostic potential of the panel of genes employed. Indeed, we were unable to source FFPE fine needle biopsies from nonfibrotic control livers. Furthermore, our analysis includes two distinct tissue sources, that is, cadaveric samples (fibrotic and nonfibrotic livers) and needle biopsies from fibrotic livers. Thus, the time to process samples from cadavers between time of death and time of processing (4 to 10 hours post mortem) may impact the mRNA levels of specific hepatic genes. Therefore, our results should be corroborated by further studies including an adequate number of cases thoroughly characterized in terms of etiology, grading and staging, and with adequate sampling conditions. Following this demonstration, the same technical multiplex scheme could be applied to samples obtained from extensive cohorts with stored paraffin samples to obtain more precise and clinically relevant information. In particular, it will be possible to correlate gene expression (employing a customized multiplex gene panel), with histopathology (stage and activity), clinical and biochemical data, and prognostic indexes (time to cirrhosis, liver related events, survival, etcetera) in retrospective or prospective analysis.
Fixed human liver samples
Human FFPE block samples were acquired from various external academic sources in compliance with Pfizer’s human tissue policy and approval by the Institutional Review Board of the source institution. Samples were de-identified and patient identities were not traceable by Pfizer at any point in time. Needle biopsies from fibrotic livers were collected from 18 HCV and seven non-alcoholic steatohepatitis (NASH) patients (14 male and 11 females; 23 to 65 years of age). To the best of our knowledge, these biopsies were fixed within 0.5 hours of procurement and processed through surgical pathology at the home institutions. In addition, we obtained 11 nonfibrotic and five (HCV-negative) fibrotic liver samples collected from cadavers fixed with a lag-time of 4 to 10 hours postmortem. These FFPE blocks were sectioned (8 μm thick), stained with hematoxylin and eosin, and transferred into 1.5 mL microcentrifuge tubes for further processing .
Fibrosis and activity scoring
All of the hematoxylin and eosin-stained liver samples were assessed by a single pathologist at Pfizer for inflammatory activity and changes in tissue architecture using scoring systems adapted from the Knodell system for inflammation  and the METAVIR scoring system for fibrosis . Inflammation grading was indicated as: 0 = no inflammation, 1 = minimal inflammation, 2 = mild inflammation, 3 = moderate inflammation, and 4 = marked/severe inflammation, with half point scores used at the pathologist’s discretion. The presence of necrotic hepatocytes was accounted for in the overall grading scheme. Fibrosis grading was indicated as: 0 = normal architecture, 1 = minimally expanded portal tracts, 2 = periportal fibrosis, 3 = bridging fibrosis, and 4 = cirrhosis.
QuantiGene™ multiplex bDNA assay for the 30-plex gene expression panel
FFPE tissue homogenates were prepared as previously described . Briefly, FFPE sections were homogenized in tissue homogenizing solution (THS) (Affymetrix Inc. Santa Clara, CA) supplemented with 500 μg/mL Proteinase K and incubated for 6 hours at 65°C. The tissue homogenate was processed immediately or stored at −80°C until use.
For the cadaver liver tissue blocks, five sections of 60 mm2 (each 8 μm thick) resulted in the most reproducible results (data not shown). These sections were homogenized individually in 300 μL of THS and digested with 500 μg/mL Proteinase K for 6 hours at 65°C. For needle biopsy samples, two sections from each of the liver needle biopsies ranging from approximated 20 to 30 mm2 (8 μm thick) were homogenized in 150 μL THS using 500 μg/mL proteinase K for 6 hours at 65°C. The mRNA quality and quantity were assessed using the 18S rRNA probe set (1 μL homogenate) and 28S rRNA probe set (0.001 μL homogenate) in duplicate on all homogenates . Finally, 40 μL of homogenate of each sample were tested in triplicate using the 30-plex gene panel.
The multiplex bDNA assays were run as recommended in the QuantiGene™ Plex 2.0 Reagent System User Manual. Probe design software  was used to generate probe sets for singleplex or multiplex bDNA assays . Probe sets were designed and synthesized by Affymetrix, Inc. (Santa Clara, CA, USA). Samples containing 40 μL of homogenate were mixed with pooled multiplex probe sets and capture beads (2,000 beads per gene per assay) prior to hybridization in 100 μL total volume, for 20 hours at 55°C. Hybridization reactions were transferred to a 0.45 μm filter plate (Millipore, Billerica, MA, USA), followed by sequential hybridization at 55°C with the bDNA amplifier and 5′-dT(biotin)-conjugated label probe. Unbound materials were washed from beads complexed with probe set and mRNA by alternating filtration and the addition of wash buffer (0.1X SSC, 0.03% lithium lauryl sulfate). Two washes were performed after each step. After the final wash, streptavidin-conjugated R-phycoerythrin (SAPE) was added and incubated at room temperature for 30 minutes. Analysis was performed after an additional wash to remove unbound SAPE.
Background signals (relative light units) were determined in the absence of RNA and subtracted from signals obtained with RNA, setting any negative signals to zero. Background subtracted signals were used to calculate the coefficient of variation (CV) and sensitivity of all assays, as well as for normalization of the RNA expression of each gene. Sensitivity was evaluated by determining the limit of detection, defined as target RNA concentration at which the signal is three standard deviations above background. Data were normalized by four assay housekeeping genes (TBP, ACTB, PPIA, and RPL13A) using the geometric mean  for each sample. The UBC assay housekeeping gene was omitted due to inconsistent behavior across multiple strip assays.
In order to identify significantly differentially expressed liver fibrotic genes, a screening approach was first followed to compare fibrotic (sourced from needle biopsy and cadavers) and nonfibrotic liver tissue (sourced only from cadavers) using the two-sample t-test on the gene expression data of 26 genes (excluding four genes of the 30 genes for analysis: SYP, CXCL2, ARNT2, and CBEP1). Considering multiplicity testing issues, comparisons were controlled using the false discovery rate at the 5% test level . In order to correlate liver fibrosis and inflammation stages with gene expression profiles, principle component analysis (PCA) and partial least-squares analysis (PLS) (Evince software package, Umbio, Umeå, Sweden) were performed on log transformed intensities, scaled and centered to unit variance.
RA’s current address: GlaxoSmithKline (China) R&D Co., Ltd. Building 3, 898 Halei Road. Zhangjiang Hi-Tech Park, Pudony. Shangai 201203, China.
DD’s current address: Pfizer Global Research & Development, 1 Burtt Rd, Andover, MA 01810, USA.
coefficient of variation
extra cellular matrix
hepatitis C virus
hepatic stellate cells
hematoxylin and eosin
principle component analysis
partial least square
tissue homogenizing solution
tissue homogenizing solution.
We would like to thank Professor Massimo Pinzani for his critical review of our manuscript.
This study was supported by Pfizer.
- Hannivoort RA, Hernandez-gea V, Friedman SL: Genomics and proteomics in liver fibrosis and cirrhosis. Fibrogenesis & Tissue Repair. 2012, 5: 1-10.1186/1755-1536-5-1.View Article
- Hoshida Y, Villanueva A, Kobayashi M, Peix J, Chiang DY, Camargo A, Gupta S, Moore J, Wrobel MJ, Lerner J, Reich M, Chan JA, Glickman JN, Ikeda K, Hashimoto M, Watanabe G, Daidone MG, Roayaie S, Schwartz M, Thung S, Salvesen HB, Gabriel S, Mazzaferro V, Bruix J, Friedman SL, Kumada H, Llovet JM, Golub TR: Gene expression in fixed tissues and outcome in hepatocellular carcinoma. N Engl J Med. 2008, 359: 1995-2004. 10.1056/NEJMoa0804525.PubMed CentralView ArticlePubMed
- Flagella M, Bui S, Zheng Z, Nguyen CT, Zhang A, Pastor L, Ma Y, Yang W, Crawford KL, McMaster GK, Witney F, Luo Y: A multiplex branched DNA assay for parallel quantitative gene expression profiling. Anal Biochem. 2006, 352: 50-60. 10.1016/j.ab.2006.02.013.View ArticlePubMed
- Cheng K, Mahato RI: Gene modulation for treating liver fibrosis. Crit Rev Ther Drug Carrier Syst. 2007, 24: 93-146. 10.1615/CritRevTherDrugCarrierSyst.v24.i2.10.PubMed CentralView ArticlePubMed
- Shimada H, Staten NR, Rajagopalan LE: TGF-β1 mediated activation of Rho kinase induces TGF-β2 and endothelin-1 expression in human hepatic stellate cells. J Hepatol. 2011, 54: 521-528. 10.1016/j.jhep.2010.07.026.View ArticlePubMed
- Xu L, Hui AY, Albanis E, Arthur MJ, O’Byrne SM, Blaner WS, Mukherjee P, Friedman SL, Eng FJ: Human hepatic stellate cell lines, LX-1 and LX-2: new tools for analysis of hepatic fibrosis. Gut. 2005, 54: 142-151. 10.1136/gut.2004.042127.PubMed CentralView ArticlePubMed
- Agresti A: Categorical Data Analysis. 2002, Wiley, Hoboken, NJ, USA, 2View Article
- SAS Institute, Inc: The GENMOD Procedure. 1999, SAS Institute, Inc, Cary, NC
- Boers W, Aarrass S, Linthorst C, Pinzani M, Elferink RO, Bosma P: Transcriptional profiling reveals novel markers of liver fibrogenesis: gremlin and insulin-like growth factor binding proteins. J Biol Chem. 2006, 281: 16289-16295. 10.1074/jbc.M600711200.View ArticlePubMed
- De Minicis S, Seki E, Uchinami H, Klume J, Zhang Y, Brenner DA, Schwabe RF: Gene expression profiles during hepatic stellate cell activation in culture and in vivo. Gastroenterology. 2007, 132: 1937-1946. 10.1053/j.gastro.2007.02.033.View ArticlePubMed
- Missiha SB, Ostrowski M, Heathcote EJ: Disease progression in chronic hepatitis C: modifiable and non-modifiable factors. Gastroenterology. 2008, 134: 1699-1714. 10.1053/j.gastro.2008.02.069.View ArticlePubMed
- Trépo E, Potthoff A, Pradat P, Bakshi R, Young B, Lagier R, Moreno C, Verset L, Cross R, Degré D, Lemmers A, Gustot T, Berthillon P, Rosenberg W, Trépo C, Sninsky J, Adler M, Wedemeyer H: Role of cirrhosis risk score for the early prediction of fibrosis progression in hepatitis C patients with minimal liver disease. J Hepatol. 2011, 55: 38-44. 10.1016/j.jhep.2010.10.018.View ArticlePubMed
- Asselah T, Bieche I, Laurendeau I, Paradis V, Vidaud D, Degott C, Martinot M, Bedossa P, Valla D, Vidaud M, Marcellin P: Liver gene expression signature of mild fibrosis in patients with chronic hepatitis C. Gastroenterology. 2005, 129: 2064-2075. 10.1053/j.gastro.2005.09.010.View ArticlePubMed
- Chen CZC, Raghunath M: Focus on collagen: in vitro systems to study fibrogenesis and antifibrosis – state of the art. Fibrogenesis & Tissue Repair. 2009, 2: 7-10.1186/1755-1536-2-7.View Article
- Seth D, Hogg PJ, Gorrell MD, McCaughan GW, Haber PS: Direct effects of alcohol on hepatic fibrinolytic balance: implications for alcoholic liver disease. J Hepatol. 2008, 48: 614-627. 10.1016/j.jhep.2007.12.015.View ArticlePubMed
- Kido O, Fukushima K, Ueno Y, Inoue J, Jefferson DM, Shimosegawa T: Compensatory role of inducible annexin A2 for impaired biliary epithelial anion-exchange activity of imflammatory cholangiopathy. Lab Invest. 2009, 89: 1374-1386. 10.1038/labinvest.2009.105.View ArticlePubMed
- Gressner OA, Lahme B, Gressner AM: Gc-globulin (vitamin D binding protein) is synthesized and secreted by hepatocytes and internalized by hepatic stellate cells through Ca2+−dependent interaction with the megalin/gp330 receptor. Clin Chim Acta. 2008, 390: 28-37. 10.1016/j.cca.2007.12.011.View ArticlePubMed
- Yang W, Maqsodi B, Ma Y, Bui S, Crawford KL, McMaster GK, Witney F, Luo Y: Direct quantification of gene expression in homogenates of formalin-fixed, paraffin-embedded tissues. Biotechniques. 2006, 40: 481-486. 10.2144/000112133.View ArticlePubMed
- Knodell RG, Ishak KG, Black WC, Chen TS, Craig R, Kaplowitz N, Kiernan TW, Wollman J: Formulation and application of a numerical scoring system for assessing histological activity in asymptomatic chronic active hepatitis. Hepatology. 1981, 1: 431-435. 10.1002/hep.1840010511.View ArticlePubMed
- Poynard T, Bedossa P, Opolon P: Natural history of liver fibrosis progression in patients withchronic hepatitis C. The OBSVIRC, METAVIR, CLINIVIR, and DOSVIRC groups. Lancet. 1997, 349: 825-832. 10.1016/S0140-6736(96)07642-8.View ArticlePubMed
- Knudsen BS, Allen AN, McLerran DF, Vessella RL, Karademos J, Davies JE, Maqsodi B, McMaster GK, Kristal AR: Evaluation of the branched-chain DNA assay for measurement of RNA in formalin-fixed tissues. J Mol Diagn. 2008, 10: 169-176. 10.2353/jmoldx.2008.070127.PubMed CentralView ArticlePubMed
- Bushnell S, Budde J, Catino T, Cole J, Derti A, Kelso R, Collins ML, Molino G, Sheridan P, Monahan J, Urdea M: ProbeDesigner: for the design of probesets for branched DNA (bDNA) signal amplification assays. Bioinformatics. 1999, 15: 348-355. 10.1093/bioinformatics/15.5.348.View 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.PubMed CentralView ArticlePubMed
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995, 57: 289-300.
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.