Volume 3 Supplement 7
Representation of genetic association via attributable familial relative risks in order to identify polymorphisms functionally relevant to rheumatoid arthritis
© Bermejo et al; licensee BioMed Central Ltd. 2009
Published: 15 December 2009
The results from association studies are usually summarized by a measure of evidence of association (frequentist or Bayesian probability values) that does not directly reflect the impact of the detected signals on familial aggregation. This article investigates the possible advantage of a two-dimensional representation of genetic association in order to identify polymorphisms relevant to disease: a measure of evidence of association (the Bayes factor, BF) combined with the estimated contribution to familiality (the attributable sibling relative risk, λs). Simulation and data from the North American Rheumatoid Consortium (NARAC) were used to assess the possible benefit under several scenarios. Simulation indicated that the allele frequencies to reach the maximum BF and the maximum attributable λs diverged as the size of the genetic effect increased. The representation of BF versus attributable λs for selected regions of NARAC data revealed that SNPs involved in replicated associations clearly departed from the bulk of SNPs in these regions. In the 12 investigated regions, and particularly in the low-recombination major histocompatibility region, the ranking of SNPs according to BF differed from the ranking of SNPs according to attributable λs. The present results should be generalized using more extensive simulations and additional real data, but they suggest that a characterization of genetic association by both BF and attributable λs may result in an improved ranking of variants for further biological analyses.
Susceptibility to rheumatoid arthritis (RA) is determined by both genetic and environmental factors, with an estimated sibling relative risk of 5-10. The HLA-DRB1 and PTPN22 genes explain approximately 50% of the familial aggregation. A recent genome-wide association (GWA) study carried out by the Wellcome Trust Case Control Consortium (WTCCC) found nine single-nucleotide polymorphisms (SNPs) with associated probability values of less than 10-5 that did not map to loci previously related to RA . In agreement with the common disease-common variant hypothesis, the novel variants were common (minor allele frequencies from 0.07 to 0.41) and they showed genotype relative risks (GRRs) of less than 2.3. The authors of the WTCCC article recognized a disparity between the large population-attributable fraction (PAF) and the modest sibling relative risk (λs) explained by the newly identified variants . This disparity is related to the technical characteristics of the platforms used in GWA studies, designed to scrutinize polymorphisms (common variants).
Let us assume that a polymorphism is a marker of a rarer causal variant, and the marker-specific PAF equals the PAF of the causal variant. Under this assumption, we have demonstrated that the causal variant has a higher attributable λs than the marker . In other words: for a fixed PAF, causal variants contribute more to familial aggregation (λs) than the polymorphisms they are represented by.
There is much debate concerning the representation of statistical evidence in GWA studies. The frequentist p-value depends on the power of the study and it does not provide a consistent measure of evidence because identical p-values do not imply identical evidence of genetic association . In contrast, the Bayes factor (BF) statistic provides a single measure of the strength of evidence of association.
We hypothesize here that the representation of genetic association by BF together with the attributable λs may help to identify polymorphisms relevant to RA, and use simulation and data from the North American Rheumatoid Consortium (NARAC) to evaluate this assumption.
Derivation of attributable risks and Bayes factors
The sibling relative risk (λs) attributable to a gene reflects the extent to which a particular genetic locus contributes to familial aggregation. The calculation of attributable λs values was first described in an early paper of James . If the frequency of the high-risk allele is denoted by RAF, the relative risk of RA for high-risk allele homozygotes compared to low-risk allele homozygotes by GRRhom, the relative risk of heterozygotes compared to low-risk allele homozygotes by GRRhet and the disease prevalence among low-risk allele homozygotes by κ0, the attributable λs is given by: λs = 1+(1/2Va+1/4Vd)/K2, where Va (additive genetic variance divided by κ02) equals 2RAF(1-RAF)[(1-RAF)(1-GRRhet)-RAF(GRRhom-GRRhet)]2, Vd (dominance genetic variance divided by κ02) equals RAF2(1-RAF)2[1+GRRhom-2GRRhet]2 and K (population prevalence divided by κ0) equals RAF2GRRHom+2RAF(1-RAF)GRRhet+ (1-RAF)2. Note that if the overall disease prevalence is 50%, λs cannot exceed two by definition. However, an important property of the attributable λs is that it only depends on RAF, GRRhom, and GRRhet. That is, the attributable λs is independent of the disease prevalence among low-risk homozygotes κ0.
The BF statistic is the ratio of the probability of the observed data under the assumption that there is a true association to its probability under the null hypothesis (absence of association). A small BF provides evidence in favor of a true association. To investigate the relationship between attributable λs and the BF, we calculated the expected distribution of genotypes in 1000 cases and 1000 controls when GRRhom = 1.5, 2, 3, 4, 5 and 10. Calculations were carried out for the dominant model, which assumes that GRRhet = GRRhom, the recessive model (GRRhet = 1), and for the additive model (GRRhet = (1+GRRhom)/2). Under Hardy-Weinberg equilibrium, the expected distribution of genotypes in controls (D = 0) is: Pr(G = aa|D = 0) = (1-RAF)2; Pr(G = Aa|D = 0) = 2RAF(1-RAF); Pr(G = AA|D = 0) = RAF2 The expected distribution of genotypes in cases (D = 1) satisfies:
GRRhom = (Pr(D = 1 | G = AA)/Pr(D = 0 | G = AA))/(Pr(D = 1 | G = aa)/Pr(D = 0 | G = aa)) and GRRhet = (Pr(D = 1 | G = Aa)/Pr(D = 0 | G = Aa))/(Pr(D = 1 | G = aa)/Pr(D = 0 | G = aa)). For example, if RAF = 0.2 and GRRhom = GRRhet = 2 (dominance), the expected distribution is (aa = 640, Aa = 320, AA = 40) in controls and (aa = 471, Aa = 471, AA = 59) in cases.
The expected distributions of genotypes were investigated by Bayesian logistic regression. We considered a general three-genotype model of association. Calculation of BFs requires assumptions about effect sizes. We assumed N(0,1) priors on the mean and the two genetic effects. The function MCMClogit (R package MCMCpack ) was used to generate a 10,000 sample dataset from the posterior distribution using a random walk Metropolis algorithm. The BayesFactor function in the same package was used to compute log marginal likelihoods for a model 'with' compared to a model 'without' genotypic information using a Laplace approximation.
Description and results from the 12 investigated regions
MAD of log10
Outlying SNPs, gene regions
λs Median (5th-95th)
Calculation of attributable risks and Bayes factors for NARAC data
SNPs in the 12 selected intervals were extracted from NARAC data, and the association between RA and the retrieved SNPs was represented by 1) the base 10 logarithm of the BF (log10BF) and 2) the attributable λs. RA status (affected or unaffected) was modelled by logistic regression, taking into account the covariates gender, number of HLA-DRB1 shared-epitope alleles (NN = 0, SN = 1 or SS = 2) and genotype (three levels). Individuals with missing information were excluded from the calculations. Only SNPs with the three genotypes represented in both cases and controls were considered.
As already mentioned, calculation of BFs requires assumptions on effect sizes. We calculated frequentist estimates of logistic regression coefficients over the corresponding entire chromosomes using NARAC data. Density plots and tests of goodness of fit of frequentist estimates of log10GRRHom and log10GRRHet for the seven investigated chromosomes indicated that GRR variation was better represented by the median absolute deviation (MAD) than by the standard deviation (data not shown). Therefore, we assumed N(0, MAD2) priors on the logarithms of genetic effects. MCMClogit and BayesFactor  were used to calculate BFs for each SNP in the selected regions based on NARAC data. The number of individuals with complete data and the observed allele frequencies in controls were used to draw 10,000 samples from a binomial distribution, which were combined with the posterior estimates of genetic effects to simulate the distribution of attributable λs for each SNP.
A bagplot is a bivariate generalization of the non-parametric univariate boxplot . In a bagplot, the region which contains 50% of the observations with greatest bivariate depth is denominated the bag. Observations outside the bag expanded three times are considered outliers, they are too far away from the data's central bulk. These outliers are indicated by closed dots in univariate boxplots. Bagplots have been shown to be equivariant for linear transformations and not limited to elliptical distributions. Due to these properties and the characteristics of our data (unknown shape of the underlying distributions), we used non-parametric bagplots to identify deviating SNPs. Bagplots were determined using the function compute.bagplot in the R package aplpack .
Theoretical relationship between Bayes factor and attributable risk
Allele frequencies at maximum λs and maximum log10(BF)
Bayes factors and attributable sibling relative risks in the investigated regions
The analysis of complete chromosomes by frequentist logistic regression resulted in 2(GRRhom, GRRhet) × 8(chromosomes 1, 4, 6, 7, 10, 13, 21, 22) sets of genetic effects. Table 1 shows the MADs of log10GRRhom and log10GRRhet used as scale parameters in the normal prior distributions. As expected, the a priori variance of genetic effects was highest for chromosome 6 (MHC region).
This study explored the advantage of combining BF, a measure of statistical evidence of association, and the attributable λs, a measure of contribution to familial aggregation, in order to characterize the effect of a particular polymorphism on disease risk. Both simulation and NARAC-based results showed that BF usually correlates with the attributable λs. Most markers identified in GWA studies confer a slightly increased risk of disease. For this kind of modest genetic effects (GRR ≤ 2), simulation indicated that the maximum log10BF and the maximum attributable λs were expected to be found for variants with similar frequencies: around 20% (dominant), 40% (additive) and 70% (recessive penetrances). The larger the size of the genetic effect, the lower the RAF to reach the maxima, and the larger the difference between RAFs. Causal variants generally show stronger effects than markers from GWA studies. For example, a recent GWA study identified five novel breast cancer susceptibility loci. The highest GRR was 1.6, which can be compared with GRR = 6 to 30 for carriers of BRCA1 mutations. Simulation suggested that, for this kind of strong genetic effects, attributable λs and BF are two different dimensions.
The representation of λs versus log10BF based on NARAC data confirmed that SNPs which have shown replicated associations clearly departed from the majority of polymorphisms in the selected regions. For example, the functionally relevant SNP rs2476601 was unequivocally separated from the rest of SNPs in the 1p13 region. Other genes located in the proximity of the most extreme SNPs were ZMYM2, MX1, MX2, IL2RB, and CSF2RB. This list includes SNPs which show strong evidence of association (log10BF higher than three) and outlying λs values. SNPs with high attributable λs values but relatively low BFs may also be of interest. Simulation showed that one BF value corresponds to two different λs values. In the 12 regions investigated using NARAC data, the ranking of SNPs according to BF differed from the ranking according to attributable λs. The difference was particularly clear in the low-recombination MHC region, where rs7765379 showed the maximum λs and rs2395175 showed the maximum log10BF. This difference in rankings suggests that, in addition to a measure of association evidence, the extensive fine-mapping needed to characterize etiological variants may benefit from the consideration of the attributable λs. Future studies based on more extensive simulations and additional real data should investigate this possibility.
The present study made use of Bayesian and robust statistics. The Bayesian approach has been overlooked in the analysis of GWA studies, where the adjustment for multiple testing, the relationship between power and statistical significance, and the selection of disease models are important issues [1, 3]. We only considered SNPs with the three genotypes represented in both cases and controls. This strategy would initially hamper the identification of rare causal variants but, on the other hand, logistic regression has not been robust to particular types of outliers, and a three-genotype model with well represented categories would result in more robust estimates. Robust methods aim to describe the structure best fitting to the bulk of the data, and to identify deviating observations. Robust statistics are used a lot less than they ought to be in genetic epidemiology. The use of MADs and bagplots in the present article illustrates their advantage over classical methods in practical situations.
The association results from GWA studies are usually summarized by a measure of evidence of association (frequentist or Bayesian probability values), which do not reflect the contribution of the identified signals to familial aggregation. We propose here a two-dimensional characterization of genetic association consisting of the attributable λs and the BF.
List of abbreviations used
Genotype relative risk
- λ s :
sibling relative risk
Median absolute deviation
Major histocompatibility complex
Wellcome trust case control consortium
The Genetic Analysis Workshops are supported by NIH grant R01 GM031575 from the National Institute of General Medical Sciences. JLB was supported by Deutsche Krebshilfe and the Swedish Cancer Society. RH was supported by NGFN 2 (SMP-GEM, PGE-S30T09), funded by the BMBF.
This article has been published as part of BMC Proceedings Volume 3 Supplement 7, 2009: Genetic Analysis Workshop 16. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/3?issue=S7.
- Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447: 661-678. 10.1038/nature05911.View ArticleGoogle Scholar
- Hemminki K, Försti A, Bermejo JL: The 'common disease-common variant' hypothesis and familial risks. PLoS ONE. 2008, 3: e2504-10.1371/journal.pone.0002504.PubMed CentralView ArticlePubMedGoogle Scholar
- Wakefield J: Bayes factors for genome-wide association studies: comparison with P-values. Genet Epidemiol. 2009, 33: 79-86. 10.1002/gepi.20359.View ArticlePubMedGoogle Scholar
- James JW: Frequency in relatives for an all-or-none trait. Ann Hum Genet. 1971, 35 (1): 47-49. 10.1111/j.1469-1809.1956.tb01377.x.View ArticlePubMedGoogle Scholar
- Martin AD, Quinn KM: MCMCpack: (MCMC) Package. [http://mcmcpack.wustl.edu/wiki/index.php/Main_Page]
- Rousseeuw PJ, Ruts I, Tukey JW: The Bagplot: a Bivariate boxplot. Am Stat. 1999, 53: 382-387. 10.2307/2686061.Google Scholar
- Wolf P: aplpack: Another Plot PACKage. [http://www.wiwi.unibielefeld.de/~wolf/software/aplpack/]
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.