# The Genetic Interpretation of Area under the ROC Curve in Genomic Profiling

Genome-wide association studies in human populations have facilitated the creation of genomic profiles which combine the effects of many associated genetic variants to predict risk of disease. The area under the receiver operator characteristic (ROC) curve is a well established measure for determining the efficacy of tests in correctly classifying diseased and non-diseased individuals. We use quantitative genetics theory to provide insight into the genetic interpretation of the area under the ROC curve (AUC) when the test classifier is a predictor of genetic risk. Even when the proportion of genetic variance explained by the test is 100%, there is a maximum value for AUC that depends on the genetic epidemiology of the disease, i.e. either the sibling recurrence risk or heritability and disease prevalence. We derive an equation relating maximum AUC to heritability and disease prevalence. The expression can be reversed to calculate the proportion of genetic variance explained given AUC, disease prevalence, and heritability. We use published estimates of disease prevalence and sibling recurrence risk for 17 complex genetic diseases to calculate the proportion of genetic variance that a test must explain to achieve AUC = 0.75; this varied from 0.10 to 0.74. We provide a genetic interpretation of AUC for use with predictors of genetic risk based on genomic profiles. We provide a strategy to estimate proportion of genetic variance explained on the liability scale from estimates of AUC, disease prevalence, and heritability (or sibling recurrence risk) available as an online calculator.

Published in the journal:
. PLoS Genet 6(2): e32767. doi:10.1371/journal.pgen.1000864

Category:
Research Article

doi: 10.1371/journal.pgen.1000864

## Summary

Genome-wide association studies in human populations have facilitated the creation of genomic profiles which combine the effects of many associated genetic variants to predict risk of disease. The area under the receiver operator characteristic (ROC) curve is a well established measure for determining the efficacy of tests in correctly classifying diseased and non-diseased individuals. We use quantitative genetics theory to provide insight into the genetic interpretation of the area under the ROC curve (AUC) when the test classifier is a predictor of genetic risk. Even when the proportion of genetic variance explained by the test is 100%, there is a maximum value for AUC that depends on the genetic epidemiology of the disease, i.e. either the sibling recurrence risk or heritability and disease prevalence. We derive an equation relating maximum AUC to heritability and disease prevalence. The expression can be reversed to calculate the proportion of genetic variance explained given AUC, disease prevalence, and heritability. We use published estimates of disease prevalence and sibling recurrence risk for 17 complex genetic diseases to calculate the proportion of genetic variance that a test must explain to achieve AUC = 0.75; this varied from 0.10 to 0.74. We provide a genetic interpretation of AUC for use with predictors of genetic risk based on genomic profiles. We provide a strategy to estimate proportion of genetic variance explained on the liability scale from estimates of AUC, disease prevalence, and heritability (or sibling recurrence risk) available as an online calculator.

## Introduction

Genome-wide association studies in human populations have facilitated the creation of genomic profiles which combine the effects of many associated genetic variants to predict risk of disease. Genetic testing has long been available for Mendelian genetic diseases for which variants within one gene are directly responsible for the disease. In contrast, the etiology of complex genetic diseases, such those listed in *Table 1*, comprises both genetic and environmental risk factors. Results from genome-wide association studies have provided empirical evidence that very few associated genetic variants with effect size greater than odds ratio of 1.5 exist *[1]*,*[2]*. Reconciliation of these effect sizes with the, often sizeable, estimates of heritability for many complex diseases (*Table 1*) means that we must expect there to be many (perhaps thousands) of genetic variants underlying complex disease if the effect size of any one variant is very small. It follows that each individual will carry a different, probably unique, portfolio of risk alleles. Whereas common risk variants have size too small to be used individually as risk predictors, profiles based on many associated genetic variants could provide useful predictions of genetic risk *[3]*,*[4]*. We define genetic risk as the risk of disease given an individual's unique multi-locus genotype; genetic risk remains unchanged throughout an individual's lifetime and so could be predicted at birth prior to exposure to many environmental risk factors. Indeed, such risk predictions could be age specific, for example, risk of type 2 diabetes at 10 years, 20 years or 50 years if genomic profile sets based on empirical data were available for these scenarios which have age-specific genetic epidemiologies. As more variants are identified in the coming years, there will be increasing interest in the prospects of genomic profiling. It has been argued that genomic profiles should be assessed in terms of their clinical validity as diagnostic classifiers *[5]*,*[6]*. The receiver operator characteristic (ROC) curve *[7]* is a well established tool for determining the efficacy of clinical diagnostic and prognostic tests in correctly classifying diseased and non-diseased individuals and has been used in the context of genomic profiling e.g., *[6]*,*[8]*,*[9]*. While the area under the ROC curve (AUC) is an important measure for clinical validity it does not tell the whole story as it does not differentiate between the accuracy with which the genomic profile predicts the true genetic risk of individuals and the accuracy with which true genetic risk predicts disease status, which is not under our control. We believe that the ability to differentiate between these components (i.e. the distinction between prediction of genotype and phenotype) is important for interpretation of the value of the genomic profile, particularly as the use of genomic profiles is very much in its infancy at present. Our knowledge of the genetic epidemiology of a disease means that we can know *a priori* that genomic profiles might not, on their own, be accurate diagnostic classifiers. For this reason, genomic profiles should judged in the first instance on the basis of their analytic validity *[10]* as predictors of *genetic* rather than *absolute* risk. Of course, in the long term genomic profiles can be combined with environmental risk factors to predict absolute risk in the context of clinical utility. Genomic profiles should improve upon family history which has long been used as a crude estimate of genetic risk (see *Text S1*).

In this paper, we provide insight into the genetic interpretation of AUC. We begin by considering quantitative traits for which the concepts of accuracy of risk prediction are well developed. For disease traits we differentiate between measures on the observed scale of disease versus the underlying scale of disease risk as we believe recognition of scale of measurement is often overlooked. We define *AUC _{max}* as the maximum AUC that could be achieved for a disease when the test classifier is a perfect predictor of genetic risk. We quantify the relationship between

*AUC*and heritability of liability and disease prevalence (lifetime morbidity risk). We show how to interpret AUC (which is a measure on the observed disease scale) of a genomic profile as the proportion of variance explained (or accuracy of prediction squared) on the underlying liability scale. Finally, we benchmark the value of genomic profiles by comparing them to the AUC expected when family history resulting from shared genetic risk factors is used as a predictor of genetic risk.

_{max}## Methods

## Background: quantitative traits

For quantitative traits, in which phenotypic scores are (or can be transformed to be) normally distributed, the efficacy of a genomic profile is naturally expressed as the proportion of the genetic variance explained by the profile. The variance in phenotypes, *V _{P}*, can be partitioned into variance of genetic values,

*V*, so that the proportion of the variance that is genetic is the heritability

_{G}*V*. Genomic profiling provides a direct estimate, , of true genetic values,

_{G}/V_{P}*G*, for individuals in a population and the efficacy of a genomic profile can be expressed as the proportion of the genetic variance explained by the profile /

*V*. We define = /

_{G}*V*, since in selection theory

_{G}*[11]*, used in livestock and plant breeding, the correlation between predicted and true genetic risk () is used as the measure of accuracy of prediction, , and if the predictor is unbiased (the regression of

*G*on is 1), . The ratio

*/V*is estimated as the

_{P}*R*

^{2}from the regression of

*P*on and is interpreted recognising its upper limit to be

*V*or heritability. These measures show that for quantitative traits, the accuracy with which the genomic profile predicts genetic risk is clearly separable from accuracy with which the true genetic risk predicts the phenotype. In contrast, AUC is a measure of the efficacy with which predicts phenotype which, as shown below, has an upper limit constrained by the heritability, and also prevalence, of the disease.

_{G}/V_{P}## Background: disease traits

For disease traits, the phenotype has two possible values, either affected or not affected. On this observed scale, the directly measurable genetic parameters are those of recurrence risks to relatives, *λ _{R}* for relatives of type

*R*, which is the ratio of the prevalence of disease in the relatives of affected individuals (

*K*) compared to the prevalence in the population (

_{R}*K*),where cov(

*X*,

*R*) the covariance in disease status between diseased individuals

*X*and their relatives on the observed disease risk scale

*[12]*. For example, when the relatives are monozygous twins (

*R*=

*MZ*), Cov(

*X*,

*MZ*) = the genetic variance, with the subscript “01” denoting the all-or-none disease risk scale. On this scale, the majority of the genetic variance is non-additive, especially when disease prevalence is low

*[13]*,

*[14]*. The broad sense heritability on this scale is = (

*λ*-1)

_{MZ}*K*/(1-

*K*) where

*λ*is the monozygotic twin recurrence risk, assuming there is no common environmental component to the recurrence risk. is not a normally reported statistic because of its dependence on disease prevalence

_{MZ}*[15]*. If the relatives are siblings (

*R*=

*S*) then

*λ*is the sibling risk ratio and Cov(

_{S}*X*,

*S*) =

*[11]*, where the variance subscripts

*A*and

*D*denote additive and dominance terms, and in combination denote epistatic variance terms. Thus, although

*λ*is an estimable quantity, it is not simply related to the genetic variances on the observed binary scale.

_{S}The genetic properties of disease are much more easily understood by using the threshold liability model *[11]*, in which risk of disease is transformed to a normally distributed liability scale *P ∼N*(0, 1) and *P* = *A* + *E*, where *A*∼*N*(0, ) are the genetic effects on the liability scale. On this scale the genetic effects combine in an additive way; is the narrow sense heritability on the liability scale (or heritability of liability) and on this scale broad sense and narrow sense heritability are equal. *E* are independent environmental effects, *E∼N*(0,1-). The biological plausibility of an underlying normally distributed liability to disease is based on the assumption that complex traits are influenced by many variables; the central limit theorem states that the distribution of the sum of independent random variables approaches normality as the number of variables increases. Under the threshold liability model individuals are affected when *P* >*T*, where *T* is the threshold on the normal distribution which truncates the proportion of affected individuals or disease prevalence (i.e., *K*), *T* = Φ^{−1}(1-*K*), Φ(*T*) = 1-*K*, where Φ(*T*) is the cumulative density function of the normal distribution up to values of *T*, e.g. if *K* = 0.05, *T* = 1.645. The threshold liability of risk scale has much nicer properties than the observed disease scale and provides a framework for comparison of scenarios independent of disease prevalence. The relationship between heritability of liability and the directly estimable parameters of *K* and *λ _{S} is*(1)

*[16]*with and

*z*the height of the standard normal curve and

*T*

_{1}= Φ

^{−1}(1-

*λ*), i.e. the threshold

_{S}K*T*

_{1}<

*T*when

*λ*>1, reflecting that the prevalence amongst sibs of affected individuals,

_{S}*K*is greater than the prevalence in the population as a whole (e.g. if K = 0.05 and

_{S}*λ*= 2,

_{S}*z*= 0.103,

*T*

_{1}= 1.282, = 0.371).

## Area under the ROC curve

The AUC is a statistic calculated on the observed disease scale and is a measure of the efficacy of prediction of phenotype using a test classifier. The ROC plots the true positive rate (TPR or sensitivity) against the false-positive rate (FPR or 1-specificity). TPR = probability (positive test result|diseased) and FPR = probability (positive test result|not diseased). Since these probabilities are conditional, they are not dependent on the number of cases or controls tested, except through the sampling variance associated with them. In genomic profiling the ROC is obtained by ranking a set of individuals with known disease status by their genomic profile from lowest estimated risk (i.e., profile score) to highest estimated risk and then assessing sensitivity and specificity assuming a cut-off after each rank (starting with the highest ranked individual). If *n _{d}* and

*n*are the numbers of diseased and not diseased individuals, and if the individual with the highest predicted genetic risk has rank

_{d'}*r*=

_{1}*n*+

_{d}*n*=

_{d'}*n*, AUC can be calculated directly from the mean rank of the diseased individuals (),(2)(see example in

*Figure S1*). Equally, AUC can be calculated as AUC = 0.5(1 +

*D*) where

*D*is the Somers' rank correlation

*[17]*between risk profile and disease status (1 = diseased, 0 = not diseased). Another equivalent definition of AUC is the probability that a randomly selected pair of diseased (

*d*) and non-diseased (

*d'*) individuals are accurately classified

*[18]*. The probability is the same as the probability that difference between the genetic liability of the

*d*and

*d'*individuals is greater than zero. This difference is approximately normally distributed with mean

*μ*-

_{d}*μ*and variance . Using the liability threshold model and results of standard genetic selection theory

_{d'}*[11]*the means (

*μ*) and variances (

*σ*

^{2}) of the genetic liability of

*d*and

*d*' individuals arewhere

*v*= -

*iK/*(1 –

*K*). The genetic liabilities of the

*d*and

*d'*groups are each approximately normally distributed, the approximation being less accurate for high heritabilities.

Therefore,(3)

## Using AUC measured by a genomic profile to estimate the proportion of genetic variance explained

A useful property of AUC (as discussed above) is that for a given disease the estimated AUC is independent of the relative proportions of cases and controls in the sample being classified *[7]*, i.e. the mean rank is approximately the same if the proportion of cases: controls is *K*: (1-*K*) or 1∶1. Or equally, the probability of a randomly selected case and control being correctly ranked is independent (except for sampling) of the number of cases and controls measured. We can use equation 3 to estimate the variance on the liability scale explained by a genomic profile, *x*, by making the subject of the equation, but renaming it as , recognising that it represents the proportion of variance explained by the profile. Then, from two measurable parameters, *K* and *AUC*, we can calculate ,(4)

Where *Q* = Φ^{−1}(*AUC*). From this, we can calculate the proportion of the known genetic variance explained by the genomic profile(5)using the estimates of *K* and *λ _{S}* to calculate (equation 1). We can also calculate the proportion of the sibling risk explained by the profile, (

*λ*

_{S}_{[x]}– 1)/(

*λ*– 1), where

_{S}*λ*

_{S}_{[x]}= (1-Φ(

*T*

_{1[x]}))/

*K*and(6)

*[19]*. and (

*λ*

_{S}_{[x]}– 1)/(

*λ*– 1) measure the same concept but in different ways and on different scales; both are useful criteria for assessing the extent to which the genomic profile accounts for the known genetic component of disease. We consider family history as a predictor of genetic risk in the

_{S}*Text S1*.

## Simulation

We used simulation under the liability threshold model *[11]*,*[14]* to check our derivations. We simulated 100,000 nuclear families sampling risk on the liability scale, *P* = *A* + *E*, *A ∼ N*(0, ) for parents, and *A* = ½*A _{dad}*+½

*A*+

_{mum}*A*for children, where the Mendelian segregation terms were random numbers sampled as

_{mend}*A*∼

_{mend}*N*(0, ½); E ∼

*N*(0,1 - ). Individuals were considered affected,

*P*

_{01}= 1, if

*P >*Φ

^{−1}(1-

*K*) =

*T*, otherwise individuals were not affected and

*P*

_{01}= 0. Genetic values on the observed scale,

*G*, were calculated as the normal probability, . From this we could calculate , , (using the

_{01}*G*and

_{01}*P*of the first child from each family) and sibling recurrence risk.

_{01}*AUC*was calculated from the mean rank of diseased individuals using equation 2 when ranked on

_{max}*A*.

## Results

## The maximum value of AUC when the test classifier is a genetic predictor depends on heritability and disease prevalence

In *Figure 1A* we consider two diseases both with heritability of liability, = 0.2, plotting probability of disease (*i.e. G _{01}*) vs genetic liability (

*i.e. A*). To allow an extreme comparison, one of the diseases has prevalence

*K*= 0.5 and the other,

*K*= 0.01.

*Figure 1B*also considers two diseases with prevalences

*K*= 0.5 and 0.01, but in this case both have = 0.8. In

*Figure 1A and 1B*, the position of the rise in probability of disease along the x-axis reflects the disease prevalence and the steepness of the rise reflects the heritability of the disease. In

*Figure 1A*the distribution of genetic liabilities on the underlying scale is exactly the same for these two diseases, but when

*K*= 0.01 higher genetic liabilities are needed before probability of disease rises above virtual zero (virtual because it is not exactly zero, but very close to zero); similarly for the diseases in

*Figure 1B*.

*Figure 1C and 1D*plot the ROC curves for the diseases considered in

*Figure 1A and 1B*, respectively. These graphs demonstrate firstly (not unexpectedly), that for diseases with the same prevalence, genetic liability is a better predictor of disease status for diseases with higher heritability and secondly, that for diseases with the same heritability, genetic liability is a better predictor of disease status for rarer diseases, because a higher proportion of those with high genetic liability are actually diseased. For example, if we used genetic liability of ≥1 as our predictor of disease, then the TPR is 0.26 and the FPR = 0.00, when

*K*= 0.5, compared to TPR = 0.99 and the FPR = 0.12, when

*K*= 0.01. These graphs demonstrate that maximum value of AUC (i.e.

*AUC*) when the test classifier is a genetic predictor is dependent on both and

_{max}*K*.

## Prediction of *AUC*_{max} from (or *λ*_{S}) and *K*

_{max}

_{S}

*Figure 2* plots *AUC _{max}* vs , for

*K*= 0.001, 0.01, 0.1, 0.3 from simulation (dashed line) and from equation 3 (solid line) and shows that

*AUC*is particularly constrained for more common or low heritability diseases. Jannsens et al

_{max}*[3]*, in their Fig. 4, have shown the relationship between AUC and the proportion of variance on the disease scale explained by the genomic profile; since their genomic profile assumed all genetic variants were known without error their graph represents the relationship between

*AUC*and . Our simulation results provided the same relationship when plotted on this scale (

_{max}*Figure 3*, solid line). In

*Figure 3*we show the relationship of

*AUC*with and (for each simulation combination of

_{max}*K*and , the

*AUC*and are calculated).

_{max}

## Complex genetic diseases

*Table 1* lists *AUC _{max}* for a range of complex genetic diseases calculated using equation 3, with calculated using equation 1 from published estimates of

*K*and

*λ*. Despite being observable, the parameters

_{S}*K*and

*λ*are subject to considerable sampling variance; we have tried, where possible, to take estimates from reviews or large studies, but large study samples simply do not exist for some low prevalence disorders. The values of

_{S}*AUC*show that it should be possible for a genomic profile for complex diseases to exceed 0.75, the threshold regarded

_{max}*[20]*as making a diagnostic classifier clinically useful when applied to a sample considered to be at increased risk. However, based on the results in

*Table 1*only the diseases with high heritability and low prevalence, such as Type I diabetes, Crohn's Disease and Lupus, can achieve an AUC, by genomic profiling alone, above the 0.99 threshold regarded

*[20]*as being required for a diagnostic classifier to be applied in the general population. In

*Table 1*, we also consider the AUC expected under scenarios where a genomic profile accounts for only a half (

*AUC*) or a quarter (

_{half}*AUC*) of the known genetic variance. These results show that for rare diseases genomic profiles can be useful classifiers of disease (

_{quar}*AUC*>0.8 when

*K*<0.01), when the profile explains only a quarter of the genetic variance.

Using equations (4) and (5) we calculate for the diseases listed in *Table 1* when AUC = 0.75. The results (*Table 1*) show that the same AUC can represent quite different successes of the genomic profile in representing the known genetic variance, ranging from 0.10 to 0.74. If we are able to explain half of the known genetic variance with identified risk variants then genomic profiles for most complex genetic disease (*AUC _{half},*

*Table 1*) will achieve some clinical validity as AUC is >0.75 for all but bladder cancer, for the examples provided.

## Example: age related macular degeneration

Consider the first listed example in *Table 1*, age related macular degeneration (AMD).

Based on the review of Scholl et al *[21]* and the large twin study of Seddon et al *[22]* we have used a prevalence after 80 years age of advanced AMD *K* = 11.8% and a sibling recurrence risk representing the genetic contribution of *λ _{S}* = 2.2, which correspond to heritability on the liability scale of = 0.68 (equation 1). If the genetic test explains all the genetic variance ( = 1), the maximum AUC that could be achieved by a genomic profile is

*AUC*= 0.92. If only half or a quarter of the genetic variance can be detected by genomic markers then the maximum AUC that can achieved are

_{max}*AUC*= 0.81 and

_{half}*AUC*= 0.72, respectively, values that exceed the prediction of genetic risk based of the most optimistic scenario from a prediction based on family history (

_{quar}*Text S1*). If complete disease status is known for all siblings, parents, grandparents, aunts, uncles and cousins then the maximum AUC that could be achieved is 0.71, translating to a genomic profile that explains 0.21 of the genetic variance (

*Table S1*). In practice, the AUC for a risk predictor based on rs1061170 a single nucleotide polymorphism in the complement factor H (

*CFH*) gene was 0.69

*[23]*(and was approximately equal for advanced AMD cases vs controls and all AMD cases vs controls). From equations 4–6, = 0.12,

*λ*

_{S}_{[x]}= 1.17, = 0.17 and (

*λ*

_{S}_{[x]}– 1)/(

*λ*– 1) = 0.15.

_{S}## Discussion

## Relationship of *AUC*_{max} to heritability and disease prevalence when the disease classifier is a genetic risk predictor

_{max}

The AUC is a widely used statistic that summarises the clinical validity of a diagnostic or prognostic test. However, the AUC statistic of a genomic profile alone has an upper limit (*i.e. AUC _{max}*) which depends on the genetic epidemiology of the disease, namely the disease prevalence and heritability. It is important that in the first instance, particularly when genomic profiling is in its infancy, that genomic profiles are judged on their ability to predict genetic risk (their analytic validity) rather than on the basis of clinical validity

*[10]*. Since AUC is estimated as a function of a rank correlation its genetic interpretation is not immediately obvious. Here we provide a genetic interpretation of the AUC expressed in terms of it genetic epidemiology parameters (equation 3). A relationship between

*AUC*and heritability was first demonstrated graphically by Janssens et al

_{max}*[3]*(see solid line

*Figure 3*). However, their representation was of broad sense heritability on the observed scale (i.e. ) which is a little used measure of heritability because of its dependence on disease prevalence

*[13]*. Here we show (

*Figure 2*and equation 3) the relationship between

*AUC*and the more commonly used measure of heritability, the heritability of liability (i.e., ) We show that

_{max}*AUC*is dependent on both and disease prevalence (

_{max}*i.e. K*).

Initially, it may seem counter-intuitive that AUC depends on disease prevalence since for an individual disease TPR and FPR are independent of the proportion of cases and controls measured and therefore of the sample prevalence. However, as we have clearly shown (*Figure 1A and 1B*) the dependence on disease prevalence results from our ability to generalise across diseases in the context of a test classifier being a genomic profile.

In contrast to our results and those of Janssens et al *[3]*, Clayton *[24]* provided an expression for ROC under a polygenic model which is independent of population disease prevalence. His derivation assumes that the effect of each locus is additive on the log risk scale *[25]*. Slatkin *[26]* and we *[27]* have found that this model allows probabilities of disease that exceed one, which although they occur with low frequency can have substantial impact on the estimates of recurrence risk and genetic variance. Under this model there is a relationship between recurrence risk to monozygotic twins and to siblings of *λ _{MZ}*/ = 1; this ratio is not achieved when probabilities of disease are constrained to their natural parameter space of a maximum of 1. Furthermore, empirical estimates of the ratio of

*λ*/ from the studies listed in

_{MZ}*Table 1*that provide estimates of

*λ*and

_{MZ}*λ*are mostly less than 1.0

_{S}*[27]*, particularly for low prevalence diseases. Recognising that these estimates are subject to sampling variance, the estimates of

*λ*/ are 1.1 (AMD), 0.4 (coronary artery disease), 0.8 (breast cancer), 0.7 (schizophrenia

_{MZ}*[25]*), 0.9 (rheumatoid arthritis) and 0.4 (Type I diabetes). Therefore, we believe the model used by Clayton to derive the relationship between AUC and heritability (or sibling recurrence risk) independent of disease prevalence is not valid.

## AUC and accuracy of genetic profiles

AUC is a useful measure because of its independence of the numbers of diseased and diseased individuals tested, but we advocate the reporting of an estimate of the proportion of the known genetic variance on the liability scale () or the proportion of sibling risk accounted for by the profile and we provide a method to do this using the estimated AUC, disease prevalence and heritability on the liability scale or sibling recurrence risk (equation 5). An AUC of 0.75 can imply anything from 0.10 to 0.74 of the genetic variance explained by the genomic profile for the complex diseases listed in *Table 1*. The correlation has long been the benchmark in non-human genetics of accuracy of genetic risk predictors. can be calculated from three measurable statistics, disease prevalence, sibling recurrence risk and AUC of the profile (using equations (1) and (4)). In this way, estimates of AUC can provide direct estimates of the proportion of ‘missing heritability’ *[28]* which takes into account the interdependence of identified associated variants.

Currently, the derivation of genomic profiles is very much in its infancy. As the sample size of genome-wide association studies increase, we can expect genomic profiles to include more and more validated associated variants. However, is constrained by the variance that could be detected by the markers that are genotyped recognising that the current generation of genome-wide chips explain at most ∼80% of the known variance in single nucleotide polymorphisms across the Caucasian genome *[29]*. This, in turn, may only be a fraction of the total genomic variance once structural variants such as copy number variants are included *[30]*. The actual variance explained by the profile depends on the sample size (i.e., power) of the studies from which associated genetic variants have been detected. It is likely that there are many variants which have such a small effect size that they will be impossible to detect even with very large samples. Although each such variant makes only a very small contribution to the genetic variance, there may be so many that a sizeable proportion of the variance will go undetected. Even if only quarter of the genetic variance is detectable by our future genotyping technology, the AUC is still greater for the genomic profile than for family history (ignoring shared environmental risks of family members, *Text S1*).

## Limitations

In our derivations we have assumed the liability threshold model *[11]*,*[14]*. Slatkin *[26]* demonstrated that the threshold model was one of several genetic models that provided the necessary steep increase in probability of disease with increasing load of genetic risk alleles *[26]*. The main assumption of the liability threshold model is that the distribution of liability scores is unimodal which should be achieved as long as there is no single unidentified genetic or environmental of very large effect *[11]*. The model accommodates any distribution of risk allele effect sizes and risk allele frequencies as long as there are sufficient (“more than one or a few” *[11]*) risk alleles in the population to create an approximately normal distribution of genetic liability scores. Since our simulation results of *AUC _{max}* vs (

*Figure 3*) based on the liability threshold model agree with those of Janssens et al

*[3]*who used a logit model to combine genetic risks from individual genetic variants, it is clear that the dependence of

*AUC*on heritability and disease prevalence is not a function of the threshold model.

_{max}We have also assumed that a genetic profile is applied in the same “average” environment as the genetic risks were estimated and we have assumed that all familiality is of genetic origin. The *AUC _{max}* will be lower than those derived here if any part of the sibling recurrence risk reflects co-variation of non-genetic origin. Using recurrence risks from different types of relatives, the importance of common environmental factors can be assessed and a

*λ*which reflects the genetic contribution of sibling recurrence can be used in our calculations. We have also assumed that the genomic profile consists of genetic markers associated with disease that are passed on according to the rules of Mendelian inheritance. In the future, a genomic profile might include non-heritable genetic variants, for example recurrent

_{S}*de novo*copy number variants or perhaps methylation status variants (for which the inheritance pattern, if any, is currently unclear

*[31]*). Such variants, although genetic, do not contribute to the similarity between relatives, and so would be included in the environmental component when partitioning variance. Under these circumstances it is possible that a genomic profile could exceed the

*AUC*based on sibling recurrence ratio. Our calculations assume that we know the population parameters

_{max}*K*and

*λ*(and therefore ). Estimates of these parameters are sometimes based on small sample size and are subject to sampling bias or different definitions of the disease. In particular, prevalence rates can depend on the age distribution of the population in which they are measured. In addition, recurrence risk ratios of relatives have a maximum possible value which is dependent on the disease prevalence, so that higher risk ratios are achievable when disease prevalence is lower; and estimates of sibling risk ratio and disease prevalence calculated in different studies sometimes reflect this dependence. In

_{S}*Table 1*, we included two different estimates for both schizophrenia and bipolar disorder, but for these examples the estimates of

*AUC*are robust to the magnitude of differences reported in genetic epidemiology parameters for individual diseases. At present, genomic profiles based on validated associated variants do not come anywhere close to the maximum implied by their

_{max}*AUC*; Jakobsdottir et al

_{max}*[6]*have reported AUC of 0.80 for risk of cardiovascular events, 0.64 for type 2 diabetes, 0.56 for prostate cancer, 0.66 for Crohn's Disease and 0.79 for age related macular degeneration. This is not surprising given the effect size of individual associated variants discovered in genome-wide association studies, which imply that much larger sample sizes will be needed to discover the majority of the variants that explain the genetic variance

*[4]*. However, already these genomic profiles outperform family history (resulting from shared genetic risk only) for four out of five of these diseases. Although the AUC is a useful summary statistic for clinical validity, in practice clinical utility depends on many other factors such as the benefits versus risks of the intervention strategies that follow from the risk prediction

*[5]*,

*[32]*; these important factors are not considered here.

## Conclusion

We have provided a genetic interpretation of and insight into the AUC statistic calculated under a genomic profile. Time will tell if genetic variants amenable to genotyping are able to reconstruct the known genetic variance in its totality. Even if it is possible to explain only a quarter of the known genetic variance, the genomic profile will be a more useful predictor of genetic risk than self-reported family history (in the absence of shared environmental risk factors) which is a commonly used measure for targeted screening programmes for complex genetic diseases. In practice, predictions of risk to disease will incorporate both genetic and environmental risk factors to produce the best predictions of absolute risk to disease. Here we provide a benchmark for the expected contribution from the genetic component of the prediction illustrating that the same AUC estimated for different diseases can imply quite different proportions of genetic variance explained by the genomic profile, which is often overlooked (e.g. *[5]*). Ultimately, genomic profiles may be used without contributions from environmental risk factors, since the contribution from the genomic profile can be estimated perinatally, prior to exposure by many environmental risk factors and when limited family history of disease is available. Indeed, one purpose of a genetic risk predictor is to allow individuals to choose to modify their exposure to environmental risks. We provide a simple online calculator (http://gump.qimr.edu.au/genroc) to calculate i) the maximum AUC for a genomic profile of a disease given estimates of disease prevalence and sibling recurrence risk or heritability of liability, ii) the proportion of variance explained on the liability scale given an estimate of AUC from a risk predictor and disease prevalence and iii) proportion of genetic variance or of sibling risk explained given an estimate AUC, disease prevalence and sibling recurrence risk *[2]*.

## Supporting Information

##### Zdroje

1. McCarthyMI

AbecasisGR

CardonLR

GoldsteinDB

LittleJ

2008 Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nature Reviews Genetics 9 356 369

2. IlesMM

2008 What can genome-wide association studies tell us about the genetics of common disease? PLoS Genet 4 e33 doi:10.1371/journal.pgen.0040033

3. JanssensAC

AulchenkoYS

ElefanteS

BorsboomGJ

SteyerbergEW

2006 Predictive testing for complex diseases using multiple genes: fact or fiction? Genet Med 8 395 400

4. WrayNR

GoddardME

VisscherPM

2007 Prediction of individual genetic risk to disease from genome-wide association studies. Genome Res 17 1520 1528

5. KraftP

WacholderS

CornelisMC

HuFB

HayesRB

2009 OPINION Beyond odds - ratios communicating disease risk based on genetic profiles. Nature Reviews Genetics 10 264 269

6. JakobsdottirJ

GorinMB

ConleyYP

FerrellRE

WeeksDE

2009 Interpretation of genetic association studies: markers with replicated highly significant odds ratios may be poor classifiers. PLoS Genet 5 e1000337 doi:10.1371/journal.pgen.1000337

7. MetzCE

1978 Basic principles of ROC analysis. Seminars in Nuclear Medicine 8 283 298

8. LuQ

ElstonRC

2008 Using the optimal receiver operating characteristic curve to design a predictive genetic test, exemplified with type 2 diabetes. American Journal of Human Genetics 82 641 651

9. van der NetJB

JanssensA

DefescheJC

KasteleinJJP

SijbrandsEJG

2009 Usefulness of Genetic Polymorphisms and Conventional Risk Factors to Predict Coronary Heart Disease in Patients With Familial Hypercholesterolemia. American Journal of Cardiology 103 375 380

10. GrosseSD

KhouryMJ

2006 What is the clinical utility of genetic testing? Genet Med 8 448 450

11. FalconerD

MackayT

1996 Introduction to Quantitative Genetics. England Longman 464

12. JamesJW

1971 Frequency in relatives for an all-or-none trait. Ann Hum Genet 35 47 49

13. DempsterER

LernerIM

1950 Heritability of Threshold Characters. Genetics 35 212 236

14. LynchM

WalshB

1998 Genetics and Analysis of Quantitative Traits. Sunderland, Massachusetts Sinauer Associates, Inc

15. RobertsonA

LernerIM

1949 The heritability of all-or-none traits - viability of poultry. Genetics 34 395 411

16. ReichT

JamesJW

MorrisCA

1972 The use of multiple thresholds in determining the mode of transmission of semi-continuous traits. Ann Hum Genet 36 163 184

17. SomersRH

1962 A new asymmetric measure of association for ordinal variables. American Sociological Review 27 799 811

18. HanleyJ

McNeilB

1982 The Meaning and Use of the Area under a Receiver Operating Characteristic (ROC) Curve. Radiology 143

19. YangJ

VisscherPM

WrayNR

2009 Sporadic cases are the norm for common disease. European Journal of Human Genetics 2009 Oct 14. [Epub ahead of print]

20. JanssensAC

MoonesingheR

YangQ

SteyerbergEW

van DuijnCM

2007 The impact of genotype frequencies on the clinical validity of genomic profiling for predicting common chronic diseases. Genet Med 9 528 535

21. SchollHPN

FleckensteinM

IssaPC

KeilhauerC

HolzFG

2007 An update on the genetics of age-related macular degeneration. Molecular Vision 13 196 205

22. SeddonJM

CoteJ

PageWF

AggenSH

NealeMC

2005 The US twin study of age-related macular degeneration - Relative roles of genetic and einivironmental influences. Archives of Ophthalmology 123 321 327

23. GuJ

PauerGJ

YueX

NarendraU

SturgillGM

2009 Assessing susceptibility to age-related macular degeneration with proteomic and genomic biomarkers. Mol Cell Proteomics 8 1338 1349

24. ClaytonDG

2009 Prediction and interaction in complex disease genetics: experience in type 1 diabetes. PLoS Genet 5 e1000540 doi:10.1371/journal.pgen.1000540

25. RischN

1990 Linkage strategies for genetically complex traits. I. Multilocus models. Am J Hum Genet 46 222 228

26. SlatkinM

2008 Exchangeable models of complex inherited diseases. Genetics 179 2253 2261

27. WrayNR

GoddardME

2010 Multi-locus models of genetic risk of disease. Genome Medicine In press

28. MaherB

2008 Personal genomes: The case of the missing heritability. Nature 456 18 21

29. BhangaleTR

RiederMJ

NickersonDA

2008 Estimating coverage and power for genetic association studies using near-complete variation data. Nature Genetics 40 841 843

30. RedonR

IshikawaS

FitchKR

FeukL

PerryGH

2006 Global variation in copy number in the human genome. Nature 444 444 454

31. YoungsonNA

WhitelawE

2008 Transgenerational epigenetic effects. Annual Review of Genomics and Human Genetics 9 233 257

32. BakerSG

CookNR

VickersA

KramerBS

2009 Using relative utility curves to evaluate risk prediction. Journal of the Royal Statistical Society 172 729 748

33. LevinsonDF

2006 The genetics of depression: A review. Biological Psychiatry 60 84 92

34. SullivanPF

NealeMC

KendlerKS

2000 Genetic epidemiology of major depression: Review and meta-analysis. American Journal of Psychiatry 157 1552 1562

35. MarenbergME

RischN

BerkmanLF

FloderusB

DefaireU

1994 Genetic susceptibility to death from coronary heart disease in a study of twins. New England Journal of Medicine 330 1041 1046

36. RischN

2001 The genetic epidemiology of cancer: interpreting family and twin studies and their implications for molecular genetic approaches. Cancer Epidemiol Biomarkers Prev 10 733 741

37. DasSK

ElbeinSC

2006 The Genetic Basis of Type 2 Diabetes. Cellscience 2 100 131

38. HemminkiK

LiX

SundquistK

SundquistJ

2007 Familial risks for asthma among twins and other siblings based on hospitalizations in Sweden. Clinical and Experimental Allergy 37 1320 1325

39. CraddockN

KhodelV

Van EerdeweghP

ReichT

1995 Mathematical limits of multilocus models: the genetic transmission of bipolar disorder. Am J Hum Genet 57 690 702

40. LichtensteinP

YipBH

BjorkC

PawitanY

CannonTD

2009 Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet 373 234 239

41. McGueM

GottesmanII

RaoDC

1983 The transmission of schizophrenia under a multifactorial threshold model. American Journal of Human Genetics 35 1161 1178

42. HarneyS

WordsworthBP

2002 Genetic epidemiology of rheumatoid arthritis. Tissue Antigens 60 465 473

43. HyttinenV

KaprioJ

KinnunenL

KoskenvuoM

TuomilehtoJ

2003 Genetic liability of type 1 diabetes and the onset age among 22,650 young Finnish twin pairs - A nationwide follow-up study. Diabetes 52 1052 1055

44. WTCCC 2007 Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447 661 678

45. HarleyJB

Alarcon-RiquelmeME

CriswellLA

JacobCO

KimberlyRP

2008 Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci. Nat Genet 40 204 210

46. SingT

SanderO

BeerenwinkelN

LengauerT

2005 ROCR: visualizing classifier performance in R. Bioinformatics 21 3940 3941

##### Štítky

Genetika Reprodukční medicínaČlánek vyšel v časopise

### PLOS Genetics

2010 Číslo 2

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…

- Population Genomics of Parallel Adaptation in Threespine Stickleback using Sequenced RAD Tags
- Genome-Wide Identification of Binding Sites Defines Distinct Functions for PHA-4/FOXA in Development and Environmental Response
- Ku Regulates the Non-Homologous End Joining Pathway Choice of DNA Double-Strand Break Repair in Human Somatic Cells
- Nuclear Pore Proteins Nup153 and Megator Define Transcriptionally Active Regions in the Genome