# Efficiency clustering for low-density microarrays and its application to QPCR

- Eric F Lock
^{1}, - Ryan Ziemiecki
^{1}, - JS Marron
^{1}and - Dirk P Dittmer
^{2}Email author

**11**:386

**DOI: **10.1186/1471-2105-11-386

© Lock et al; licensee BioMed Central Ltd. 2010

**Received: **31 March 2010

**Accepted: **20 July 2010

**Published: **20 July 2010

## Abstract

### Background

Pathway-targeted or low-density arrays are used more and more frequently in biomedical research, particularly those arrays that are based on quantitative real-time PCR. Typical QPCR arrays contain 96-1024 primer pairs or probes, and they bring with it the promise of being able to reliably measure differences in target levels without the need to establish absolute standard curves for each and every target. To achieve reliable quantification all primer pairs or array probes must perform with the same efficiency.

### Results

Our results indicate that QPCR primer-pairs differ significantly both in reliability and efficiency. They can only be used in an array format if the raw data (so called CT values for real-time QPCR) are transformed to take these differences into account. We developed a novel method to obtain efficiency-adjusted CT values. We introduce transformed confidence intervals as a novel measure to identify unreliable primers. We introduce a robust clustering algorithm to combine efficiencies of groups of probes, and our results indicate that using *n* < 10 cluster-based mean efficiencies is comparable to using individually determined efficiency adjustments for each primer pair (*N* = 96-1024).

### Conclusions

Careful estimation of primer efficiency is necessary to avoid significant measurement inaccuracies. Transformed confidence intervals are a novel method to assess and interprete the reliability of an efficiency estimate in a high throughput format. Efficiency clustering as developed here serves as a compromise between the imprecision in assuming uniform efficiency, and the computational complexity and danger of over-fitting when using individually determined efficiencies.

## Background

Array and microarray based methods are the mainstay of molecular biology. Recently, lower-density or targeted arrays have been introduced. These comprise on the order of 10-1000 targets and represent an intermediate between 1 target assays, e.g. the viral load assay for HIV, and genomewide microarrays, comprising 10,000 -1,000,000 targets, e.g. Affymetrix™ arrays. Analytically they can be treated as either a collection of individual assays or as microarrays.

We are interested in discovering novel statistical approaches to the analysis of these intermediate density arrays. On the one hand high density microarray-based analysis methods do not capture all the information that is available. This results in lower than possible linear range and lower discriminating power. On the other hand, applying analysis methods developed for a single assay to arrays is overly complex and likely to introduce inacurracies due to overfitting problems.

Practically lower density arrays are based on real-time quantitative polymerase chain reaction (QPCR). Real-time QPCR measures the amount of product at each cycle of the reaction either by binding of a uorescent, double strand-specific dye (SYBR™ green) or by hybridization to a third sequence-specific, dual-labeled uorogenic oligonucleotide probe (molecular Beacon, TaqMan™ ). These have been used very successfully to profile messenger RNAs and microRNA [1–5]. Other assay formats use capture oligonucleotides or other means (e.g. Nanostring™ or Luminex™).

The simplest way to compare relative expression ratios is the so-called method, which assumes (i) uniform and (ii) perfect efficiency: the amount of product doubles exactly after each PCR cycle for each primer in the array. However, most QPCR reactions do not attain optimal efficiency [6, 7]. Even small deviations from an assumed efficiency level can lead to multiple-fold inaccuracies in expression measurements [8]. Therefore, differing efficiencies among primer-pairs, as well as the reliability of the efficiency estimates must be considered.

The most accurate and also most complex method is to calculate individual primer pair efficiencies or to use an absolute standard curve. This method does not assume uniform efficiency among different primer pairs and it does not assume perfect efficiency. The well known REST method [9] is an example of this strategy. However, this approach becomes impractical and extremely costly to use even for low density arrays.

More important, using individually adjusted efficiencies brings about the problem of overfitting and it requires that the individual efficiencies were determined with equal confidence. For instance, a 96 primer targeted QPCR array, would require the computation of 96^{2} individual efficiency ratios. These need to be tracked throughout the analysis process and each ratio is associated with its own error, which also needs to be propagated throughout the analysis process.

The methods presented here use the serial dilution method [10] to obtain initial primer pair efficiencies. We expand on this method, introducing transformed normality-based confidence intervals as a novel tool to interpret the reliability of an efficiency estimate. We propose a rule to identify unreliable primers, as well as a robust algorithm to cluster primers with similar efficiencies. Finally, we show how differences in efficiency can be applied for more accurate comparisons of relative gene expression.

## Methods

### Experimental Methods

The primer arrays and experimental methods were previously described [11, 12]. Primers were from commercial sources (MWG Inc.) and resuspended at 100 pmol per microliter in 0.1x Tris-EDTA ph 8.0. The QPCR reaction contained 2.5 microliter of primer mix at a starting concentration of 300 nM combined with 7.5 microliter SYBR Green 2× PCR mix (Applied Biosystems, Inc.) and 5 microliter target DNA. It was subjected to real-time QPCR on an Opticon2 cycler (MJR Inc.) using standard cycling conditions [13]. This particular QPCR array was directed against every open reading frame of human herpesvirus 8. This particular virus is made of linear double-stranded DNA. Thus we were able to use the same purified viral DNA as a common target for all primers.

### Quantitative Analysis

The R programming environment for statistical computing and graphics [14] was used for all computation and statistical analysis. Stand-alone R functions were developed to automate the analysis for a given serial dilution table. Functions were created to provide individual efficiency estimates with confidence, identify and remove unreliable primers, cluster amplification efficiencies, and adjust *C*_{
T
}values in a table based on user-defined efficiencies. All functions were included in the R script Primer Efficiency Analysis (Additional file 1 - PEA.r); a reference manual is available with setup instructions, detailed function descriptions and illustrative examples (Additional file 2 - PEA User's Guide).

## Results and Discussion

### Estimating Efficiency with Confidence

The number of cycles required to reach a certain level of product (*C*_{
T
}) is measured for 96 primers, each with initial dilution levels 1, 0.1, 0.01, and 0.001. *C*_{
T
} values above 40 indicate undetectable product and are considered to be missing values (NA) for subsequent calculations. If the dilution curve for a primer contains 2 or more missing values (50% of total), statistical analysis for a linear fit is impossible and the primer is removed from the analysis. Of the 96 original primers in our dataset, 3 were removed due to missing values.

*E*) estimates for each primer-pair are calculated using their serial dilution curve. Theoretically, one expects a linear relationship between

*C*

_{ T }and the logarithm of the initial dilution level. Therefore, we fit a standard linear model of the form

_{2}) rather than the common log

_{10}[15], as the former has the nice interpretation that a unit slope (

*β*

_{1}= -1) corresponds to perfect efficiency (

*E*= 2). Figure 2 shows the fitted linear model for a single primer-pair.

Before continuing, we remove those primers for which a linear fit is not statistically significant, based on the t-statistic for the estimated slope coefficient
. If the corresponding p-value is above 0.05, the primer is considered unreliable and removed from further analysis. Using the t-statistic rather then the regression coefficient *R*^{2} initially serves the following purpose. The t-statistic accounts for the number of dilution steps, whereas a threshhold based on the coefficient of determination *R*^{2} does not (*R*^{2} = 1 for 2 data points and likely decreases with the addition of more measurements). More dilution steps generally result in more accurate efficiency estimates. In our dataset, 8 primers were removed due to an insignificant linear slope.

*β*

_{1}. For each primer a 95% confidence interval for the slope, (

*β*

_{ lower };

*β*

_{ upper }), is given by

where
is the estimated slope, t_{0.025;n-2}is the 2.5% quantile of Student's t-distribution with *n* - 2 degrees of freedom, and SE(
) is the standard error of
. This analysis used 3 or 4 dilution levels for each primer (*n* = 3 or 4), so we use a t-distribution with 1 or 2 degrees of freedom.

*β*

_{1}< 0, confidence intervals for

*β*

_{1}are easily extended to confidence intervals for

*E*. Under the condition

*β*

_{ upper }< 0, a 95% confidence interval for

*E*is given by (

*E*

_{ lower };

*E*

_{ upper }), where

Only keeping those primers with a negative and significant slope ensures that the condition *β*_{
upper
} < 0 is satisfied.

### Identifying Unreliable Primers

Bartlett's test of heteroscadasticity (unequal variability within groups), applied to the residuals of each linear model fit by (1), is highly significant (*p* < 0.0001). This indicates that certain aberrations in efficiency estimation may be due to inherent primer reliability issues, rather than standard experimental and residual error. Purging those primers with missing values and an insignificant slope can be considered a first step to identify primer-pairs yielding unreliable data, but further analysis is warranted. Here we introduce transformed confidence intervals as a novel measure of primer reliability

The length of the transformed confidence interval *E*_{
upper
} - *E*_{
lower
} is used to quantify the precision of the estimated amplification efficiency
. To identify those primers with unreliable data, we first calculate the interquartile range (IQR) of the transformed CI lengths as the difference betweeen the 75th and 25th percentiles (quantiles). Any primer with a transformed CI length higher than 1.5 × *IQR* above the 75th percentile is identified as unreliable. Such a procedure (the "1.5 × *IQR*" rule) is commonly used to identify outliers. In our dataset, 7 primers are identified using this procedure.

*R*

^{2}value of the original linear model for each primer. The coefficient of determination

*R*

^{2}measures the strength of a linear fit, and can be interpreted as the percentage of variation in the response (

*C*

_{ T }) that is explained by

*log*

_{2}(

*Dilution*). Primers identified as outliers by the 1.5 ×

*IQR*rule for transformed CI length are colored red. After removing outliers based on CI length, the remaining primers all show strong linearity with

*R*

^{2}values greater than 0.98 and correlation coefficients greater than 0.99.

Most of those primers identified for large transformed CI's tend to have lower *R*^{2} values, whereas there is no correlation with primer pair efficiency. Hence, a low efficiency primer pair may indeed be very reliable and an apparently good efficiency primer may be the result of a bad fit. This demonstrates that considering a measure of variation for each efficiency estimate adds value over prior methods that consider efficiency as the sole criterion for primer performance.

*C*

_{ T }and log dilution level for every primer (4a), and after removing missing values, insignificant linear slopes and CI length outliers (4b). The purged data is much "cleaner", with non-linear and sporadic relationships removed. This purged set of primers is appropriate for use in subsequent, rank-based, clustering and classification analyses, whereas the original set would produce unacceptable results.

### Comparing Primer Efficiencies

Even after purging the data of unreliable primers, estimates of primer efficiency vary, from *E* = 1.7 to *E* > 2. It has been shown that rather than improving accuracy, individual corrections to efficiency can in fact exaggerate artificial differences in expression [16]. Therefore, we consider the hypothesis that there is no underlying difference in primer efficiencies, and observed differences are due to residual error.

The F-statistic for the interaction between primer and log dilution *log*_{2}(*dil*) is highly significant (*p* = 0.0017). This indicates that differences among slope coefficients, hence the resulting efficiency estimates, are not merely due to residual error.

This significant difference in primer efficiencies shows that simply using the average efficiency (*E* = 1.90) among all primers is inappropriate for comparing expression levels. An alternative approach is to use the individually determined efficiency value for each primer. Both methods generate errors on efficiency, which may be propogated through the calculation of expression quantities [17].

A trade-off exists between the lack of precision in not recognizing differences in primer efficiency, and the computational complexity and danger of over-fitting when using individually determined efficiencies for each primer. Here we introduce a novel, third alternative, which finds a useful balance between the two extreme approaches. We first cluster primer-pairs based on their estimated efficiency values, then use only a measure of the average efficiency within each cluster for subsequent adjustments. Rather than a single efficiency, i.e. no correction, we use 5-10 different efficiencies, based on clustering to obtain adjusted *C*_{
T
} values. Rather than 96-1024 individual, but computational unreliable, efficiencies, we use 5-10 different cluster-derived efficiencies, knowing that within each cluster, there exists no significant difference between individual primer pairs.

*R*

^{2}). If the model with two groups significantly improves over the model with a single slope (the p-value of the corresponding F-statistic is less than 0.05), we divide the two groups into three optimal groups of primers. If three groups represents a significant improvement over two, we divide the primers into four groups, and so on. For this particular dataset, we find that the optimal model with eight slope coefficients does not significantly improve on the model with seven, hence we identify seven clusters. The primer clusters and corresponding efficiencies are shown in Figure 5.

Note that the majority of efficiencies, even across clusters, still are between E = 1.8 and E = 2.0. This is expected since we started with a "good" primerset based on sequence criteria [11]. Still we are able to subdivide the "efficiency space" 1.8 to 2.0 into four clusters and thus improve overall performance. Furthermore, we are able to use even the relatively divergent primers pairs with *E* < 1.8 and *E* > 2 (due to repeats in the target sequence) for further analysis thus maximizing usability. The clustering method described is robust in the sense that significant differences in primer efficiencies should be reflected in the number of clusters. Both the magnitude of observed efficiency differences and precision of those efficiency estimates are accounted for. Arrays where the observed differences in primer efficiency are not significant are expected to result in few clusters, whereas if efficiencies differ greatly and can be measured precisely the method will tend to identify many clusters.

### Comparison of Methods

*E*

^{ CT }) divided by the initial dilution:

For replications of the same experiment with different initial dilutions, we expect the AFC to be the same. However, an inaccurate assumed efficiency *E* will lead to different AFC values. For each primer, we calculate the coefficient of variation for the AFC values corresponding to the four dilution levels. The coefficient of variation (CoV) is defined as the standard deviation of the four values, divided by their mean.

*E*= 2 yields the worst result, i.e. the largest CoV across all primers (Figure 6, black line). Using instead the experimentally determined single mean efficiency (

*E*= 1.90) across the array improves the CoV significantly (Figure 6, green line). Using the 7 average efficiencies by cluster lowers the CoV even further overall (Figure 6, blue line). As expected it improves most dramatically the performance of the most divergent primers, i.e. the tails (primer 1-10 and 70-80). Notably, the curves for the 7 average efficiencies by cluster and for the individually determined efficiencies (Figure 6, red line) are overlapping. This demonstrates that using just a small number of average cluster efficiencies is as good as using individual primer efficiencies, i.e. the differences in efficiency within each cluster are negligible.

### Comparing Relative Expressions

*C*

_{ T }, which we call . These are defined as follows. The relative expression ratio between a reaction using primer pair A and a second reaction using primer pair B is given by

*E*

_{ A }and

*E*

_{ B }are the efficiencies of primer pair A and B, respectively. If the two primers are assumed to have the same efficiency

*E*(e.g. they belong to the same efficiency cluster), the calculation is simplified:

This insight is captured for group-wise comparison in the REST [9] software package, which also incorporates error estimates and tests of significance, while allowing for manual input of PCR efficiencies.

*C*

_{ T }values after adjusting for amplification efficiency. We propose the following adjustment:

Here,
is the estimated number of cycles required under perfect efficiency, so that the true fold difference is given by
. For example, if a reaction has perfect efficiency (*E* = 2) with observed *C*_{
T
} = 24, while another has sub-perfect efficiency (*E*= 1.8) with observed *C*_{
T
} = 28.5, we calculate
= 24 for both reactions, indicating no difference in expression. A practical advantage of using
is that the data format of the original set-up remains the same. The expression data remain log transformed and can be used directly in any microarray profiling software such as Eisen's original clustering program [18]. Calculating
is a preprocessing step that does not change the runtime or memory requirements of the subsequent analysis programs.

*C*

_{ T }values for seven primer pairs with different efficiencies. Multiple-fold differences are observed between estimated expression under perfect efficiency and expression after correcting for imperfect efficiency. This suggests that the assumption of perfect efficiency can often lead to multiple fold inaccuracies in relative expression.

Adjusted Relative Expression.

## Conclusions

As array-based measurements for DNA, mRNA and microRNA levels migrate into the mainstay of molecular biology, failure to carefully consider the efficiency of each individual reaction or assay can lead to significant measurement inaccuracies. Yet, explicitly calculating and considering individual assay characteristics is not feasible even for low-density arrays. This is an important problem particularly for real-time QPCR based arrays, but it applies to any type of microarray. QPCR primers differ significantly both in reliability and amplification efficiency, hence they need to be experimentally validated. (i) We identified transformed confidence intervals as a useful means to assess and interpret the reliability of an efficiency estimate. Transformed confidence intervals provide a novel, independent measure in addition to calculation of primer efficacy E, with which to assess primer/probe quality. (ii) After purging unreliable estimates we propose a robust clustering algorithm to group efficiencies, reducing computational complexity and potential over-fitting. Our results suggest that use of a limited number of clustering-based efficiencies is comparable to use of individually determined efficiencies for each primer or probe.

## Declarations

### Acknowledgements

This work was supported by funding from the University Cancer Research Fund, the AIDS malignancies clinical trials consortium (AMC), and the NIH (DE018304).

## Authors’ Affiliations

## References

- Schmittgen T, Jiang J, Liu Q, Yang L: A high-throughput method to monitor the expression of microRNA precursors.
*Nucleic Acids Res*2004, 32: e43. 10.1093/nar/gnh040View ArticlePubMedPubMed CentralGoogle Scholar - Dittmer D: Transcription profile of Kaposi's sarcoma-associated herpesvirus in primary Kaposi's sarcoma lesions as determined by real-time PCR arrays.
*Cancer Res*2003, 63: 2010–5.PubMedGoogle Scholar - O'Hara A, Vahrson W, Dittmer D: Gene alteration and precursor and mature microRNA transcription changes contribute to the miRNA signature of primary effusion lymphoma.
*Blood*2008, 111: 2347–53. 10.1182/blood-2007-08-104463View ArticlePubMedPubMed CentralGoogle Scholar - O'Hara A, Chugh P, Wang L, Netto E, Luz E, Harrington W, Dezube B, Damania B, Dittmer D: Pre-micro RNA signatures delineate stages of endothelial cell transformation in Kaposi sarcoma.
*PLoS Pathog*2009, 5: e1000389. 10.1371/journal.ppat.1000389View ArticlePubMedPubMed CentralGoogle Scholar - Gaur A, Jewell D, Liang Y, Ridzon D, Moore J, Chen C, Ambros V, Israel M: Characterization of microRNA expression levels and their biological correlates in human cancer cell lines.
*Cancer Res*2007, 67: 2456–68. 10.1158/0008-5472.CAN-06-2698View ArticlePubMedGoogle Scholar - Yuan J, Reed A, Chen F, Stewart C: Statistical analysis of real-time PCR data.
*BMC Bioinformatics*2006, 7: 85. 10.1186/1471-2105-7-85View ArticlePubMedPubMed CentralGoogle Scholar - Karlen Y, McNair A, Persequers S, Mazza C, N M: Statistical significance of quantitative PCR.
*BMC Bioinformatics*2007, 8: 131. 10.1186/1471-2105-8-131View ArticlePubMedPubMed CentralGoogle Scholar - Bar T, Stahlberg A, Muszta A, M K: Kinetic Outlier Detection (KOD) in real-time PCR.
*Nucleic Acids Res*2003, 31: e105. 10.1093/nar/gng106View ArticlePubMedPubMed CentralGoogle Scholar - Pfaffl M, Graham W, Dempe L: Relative expression software tool (REST ©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.
*Nucleic Acids Res*2002, 30: e36. 10.1093/nar/30.9.e36View ArticlePubMedPubMed CentralGoogle Scholar - Pfaffl M: A new mathematical model for relative quantification in real-time RT-PCR.
*Nucleic Acids Res*2001, 29: e45. 10.1093/nar/29.9.e45View ArticlePubMedPubMed CentralGoogle Scholar - Papin J, Vahrson W, Hines-Boykin R, Dittmer D: Real-time quantitative PCR analysis of viral transcription.
*Methods Mol Biol*2005, 292: 449–80.PubMedGoogle Scholar - Fakhari F, Dittmer D: Charting latency transcripts in Kaposi's sarcoma-associated herpesvirus by whole-genome real-time quantitative PCR.
*J Virol*2002, 76: 6213–23. 10.1128/JVI.76.12.6213-6223.2002View ArticlePubMedPubMed CentralGoogle Scholar - Hilscher C, Vahrson W, DP D: Faster quantitative real-time PCR protocols may lose sensitivity and show increased variability.
*Nucleic Acids Res*2005, 33: e182. 10.1093/nar/gni181View ArticlePubMedPubMed CentralGoogle Scholar - R Development Core Team: R Foundation for Statistical Computing, Vienna, Austria.
*R: A Language and Environment for Statistical Computing*2009. [ISBN 3–900051–07–0] [http://www.R-project.org] [ISBN 3-900051-07-0]Google Scholar - Schefe J, Lehmann K: Quantitative real-time RT-PCR data analysis: current concepts and the novel "gene expression's Ct difference" formula.
*J Mol Med*2006, 84: 901–910. 10.1007/s00109-006-0097-6View ArticlePubMedGoogle Scholar - Peirson S, Butler B, Foster R: Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis.
*Nucleic Acids Res*2003, 31: e73. 10.1093/nar/gng073View ArticlePubMedPubMed CentralGoogle Scholar - Hellemans J, Mortier G, Paepe A, Speleman F, J V: qBase relative quantification framework and software for management and automated of real-time quantitative PCR data.
*Genome Biology*2007, 8: R19. 10.1186/gb-2007-8-2-r19View ArticlePubMedPubMed CentralGoogle Scholar - Eisen M, Spellman P, Brown P, Botstein D: Cluster analysis and display of genome-wide expression patterns.
*Proc Natl Acad Sci*1998, 95: 14863–8. 10.1073/pnas.95.25.14863View ArticlePubMedPubMed CentralGoogle Scholar

## Copyright

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.