# A Pleiotropy-Informed Bayesian False Discovery Rate Adapted to a Shared Control Design Finds New Disease Associations From GWAS Summary Statistics

Many diseases have a significant hereditary component, only part of which has been explained by analysis of genome-wide association studies (GWAS). Shared aetiology, treatment protocols, and overlapping results from existing GWAS suggest similarities in genetic susceptibility between related diseases, which may be exploited to detect more disease-associated SNPs without the need for further data. We extend an existing method for detecting SNPs associated with a given disease by conditioning on association with another disease. Our extension allows GWAS for the two conditions to share control samples, enabling larger overall control groups and application to the common case when GWAS for related diseases pool control samples. We demonstrate that our technique limits the expected overall false discovery rate at a threshold dependent on the two diseases. We apply our technique to genotype data from ten immune mediated diseases. Overall pleiotropy between phenotypes is demonstrated graphically. We are able to declare several SNPs significant at a genome-wide level whilst controlling at a lower false-discovery rate than would be possible using a conventional approach, identifying eight previously unknown disease associations. This technique can improve SNP detection in GWAS by re-analysing existing data, and gives insight into the shared genetic bases of autoimmune diseases.

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

Category:
Research Article

doi: 10.1371/journal.pgen.1004926

## Summary

Many diseases have a significant hereditary component, only part of which has been explained by analysis of genome-wide association studies (GWAS). Shared aetiology, treatment protocols, and overlapping results from existing GWAS suggest similarities in genetic susceptibility between related diseases, which may be exploited to detect more disease-associated SNPs without the need for further data. We extend an existing method for detecting SNPs associated with a given disease by conditioning on association with another disease. Our extension allows GWAS for the two conditions to share control samples, enabling larger overall control groups and application to the common case when GWAS for related diseases pool control samples. We demonstrate that our technique limits the expected overall false discovery rate at a threshold dependent on the two diseases. We apply our technique to genotype data from ten immune mediated diseases. Overall pleiotropy between phenotypes is demonstrated graphically. We are able to declare several SNPs significant at a genome-wide level whilst controlling at a lower false-discovery rate than would be possible using a conventional approach, identifying eight previously unknown disease associations. This technique can improve SNP detection in GWAS by re-analysing existing data, and gives insight into the shared genetic bases of autoimmune diseases.

## Introduction

Genome-wide association studies (GWAS) have enabled identification of genetic variants associated with a wide range of complex phenotypes, but in many cases these variants explain only a proportion of the known heritability [*4*]. There is increasing evidence that this is due to the combined contribution of small effects arising from multiple distinct variants [*5*]. The testing of a large number of potential variants in parallel, with a comparatively low number of samples, mandates a stringent threshold for significance in order to limit false positives (type 1 errors), meaning that discovery of variants responsible for small effects requires very large sample sizes. Detection of such variants by increasing numbers of samples in studies is time-consuming and expensive, particularly for rare phenotypes, but it may be possible to improve detection by re-analysis of existing data [*6*]. One promising strategy is to co-analyse GWAS results from similar phenotypes to exploit potential similarities in genetic aetiology. This has been attempted using several different methods [*2*, *7*, *8*].

The assumption that GWAS for similar diseases may yield overlapping sets of disease-associated variants is based on the phenomenon of pleiotropy, in which a genetic variant is associated with more than one trait or disease [*9*]. Pleiotropy is common in human genes: even when exclusively considering single nucleotide polymorphisms (SNPs) with strong evidence of association, around 15% of those associated with at least one trait are associated with multiple traits [*10*]. Elements of shared genetic aetiology may be suspected in diseases with similar symptomatology, such as bipolar disorder and schizophrenia [*11*] or in diseases with common risk factors, such as type 2 diabetes and obesity [*12*]. If two diseases are known or suspected to share associated genetic variants, a degree of association of a locus with one disease may increase the likelihood of association with the other. Use of external covariates in this way can alleviate some of the effect of multiple testing [*8*], meaning that phenotypic similarity may lead to improved detection of disease-associated variants. Correspondingly, discovery and specification of shared genetic aetiology between two diseases may suggest some shared pathophysiology [*12*].

A technique for improved discovery of disease variants using pleiotropy between pairs of diseases has been successfully developed and applied by Andreasson et al [*3*, *13*, *14*]. The technique extends the empirical Bayesian false discovery rate [*15*] to a two-phenotype scenario, in which association with one phenotype is tested conditional on varying degrees of association with another. We denote the phenotype for which association is being tested as the ‘principal phenotype’ and the other as the ‘conditional phenotype’.

By successively restricting attention to SNPs with a given strength of association in the conditional phenotype, the number of parallel tests to perform for association with the principal phenotype is reduced. If the two phenotypes share common associated variants, this restriction will retain disease-associated SNPs at a higher rate than null SNPs, resulting in a higher proportion of disease-associated SNPs in the restricted group than in the whole. The ‘conditional false discovery rate’ (cFDR), defined as the probability that a SNP is not associated the principal phenotype given its p values for the principal and conditional phenotypes are below some thresholds, exploits this effect. By computing cFDRs for schizophrenia conditioned on bipolar disorder and vice versa, Andreasson et al [*3*] identified multiple previously undiscovered loci for both. In a separate study computing cFDRs for hypertension conditioned on 12 related traits [*13*], 42 new loci associated with hypertension were reported. These constituted considerable improvement on existing results using single GWAS, albeit using a rather relaxed threshold of estimated *cFDR* ≤ 0.01.

A major disadvantage of the algorithm developed and used by Andreasson et al is the requirement that control groups for the two GWAS be distinct, in order to ensure that observed effect sizes are uncorrelated at null SNPs. This requires splitting a pool of potential controls between studies, with the summary statistics for each GWAS computed from only the controls allocated to that study. This may be impractical as it requires access to raw genotype data. More importantly, accuracy of effect size estimates improves with larger control groups, and consequently splitting controls in this way weakens the effect size estimates for individual studies. For this reason, many researchers employ a study design in which controls are pooled into a large group; for example, the Wellcome Trust Case Control and ImmunoChip consortia [*16*, *17*].

Here we extend the cFDR approach to studies with overlapping control groups, exploiting an approach developed by Zaykin et al, following Lin et al [*18*, *19*] to adjust for the effect of shared controls. This allows the strongest available estimates of effect sizes to be used for calculation, and consequently strengthens the power of the technique. Our technique additionally allows cFDR rates to be computed from summary statistics alone, without the need to recalculate effect sizes after re-allocating controls. We demonstrate the improvement arising from sharing controls in a type 1 diabetes data set.

We also identify a previously undiscussed difficulty with the technique potentially leading to a falsely low estimate of false discovery rate amongst SNPs declared non-null. Multiple overlapping sets of SNPs may be defined each of which has *cFDR* ≤ *α*. However, the union of these sets does not necessarily have an expected false-discovery rate less than *α* and is generally higher. An implication of this is that if we declare non-null all SNPs for which estimated cFDR is less than *α*, the expected overall false-discovery rate amongst SNPs declared non-null is greater than *α*. We describe an upper bound on the false discovery rate amongst such SNPs based on areas of regions of the unit square.

We apply our method to summary SNP association statistics for ten phenotypically distinct autoimmune diseases: type 1 diabetes (T1D) [*20*], autoimmune thyroid disease (ATD) [*21*], coeliac disease (CEL) [*22*], multiple sclerosis (MS) [*23*], narcolepsy (NAR) [*24*], primary biliary cirrhosis (PBC) [*25*], psoriasis (PS) [*26*], rheumatoid arthritis (RA) [*27*], ulcerative colitis [*28*], and Crohn’s disease [*28*]. All were genotyped using a common SNP array: the ImmunoChip, designed to provide dense genotype coverage of regions associated with autoimmune disease. Many autoimmune traits are known to have significant heritability, much of which remains unexplained [*29*]. We hypothesised that our method can improve detection of disease-associated variants in these diseases without the need for distinct control groups.

## Results

## Overview of method

The unconditional false discovery rate for a set of SNPs with p values < *p*_{i} is defined as the probability that a random SNP from this set is null. We denote this as *uFDR*(*p*_{i}), and our estimate as u F D R ̂ ( p i ).

The conditional false discovery rate (*cFDR*) is defined [*3*, *14*] as the probability that a random SNP is null for a phenotype *i* given that the observed p values at that SNP for phenotypes *i* and *j* are less than (*p*_{i},*p*_{j}); that is, P r ( H 0 ( i ) ∣ P i ≤ p i,P j ≤ p j ), where H 0 ( i ) is the null hypothesis that the SNP is not associated with phenotype *i*. We denote this quantity as *cFDR*(*p*_{i}∣*p*_{j}), and call phenotype *i* the ‘principal phenotype’ and phenotype *j* the ‘conditional phenotype’.

We first apply genomic control to allow the assumption that, globally, P values for null SNPs are uniformly distributed on [0, 1]. We compute an estimate of the *cFDR*, which we denote c F D R ̂ ( p i ∣ p j ), in a similar manner to that proposed by Andreasson et al, but incorporating expected non-uniformity in the distribution of *P*_{i} due to the sharing of controls. As u F D R ̂ ( p i ) is monotonically related to *p*_{i}, we set a significance cutoff at the maximum value of u F D R ̂ ( p i ) with *p*_{i} < 5 × 10^{−8}. Correspondingly, we set a significance cutoff for c F D R ̂ ( p i,p j ) at the maximum c F D R ̂ ( p i,p j ) with *p*_{i} < 5 × 10^{−8}. Implementation of these steps in R is available from https://github.com/jamesliley/cFDR-common-controls.

## Sharing of control subjects

If no controls are shared between studies, it is reasonable to assume that observed effect sizes for the two phenotypes are independent under a null hypothesis for the principal phenotype. This implies that the expected quantile of a given SNP’s p value for the principal phenotype is simply the p value itself regardless of its p value for the conditional phenotype. However, when control samples are shared, this assumption is invalid. Shared controls induce a positive correlation on estimated effect sizes for the principal and conditional phenotype [*18*, *19*], meaning that when attention is restricted to SNPs with a given degree of association with the conditional phenotype, the p values for the principal phenotype will be falsely low; that is, the probability P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ) will not in general be equal to *p*_{i} (= P r ( P i ≤ p i ∣ H 0 ( i ) )); in fact it will usually be higher.

When controls are shared, the distribution of p values for the principal phenotype given p values for the conditional phenotype depends on the underlying effect of each SNP on the conditional phenotype. For any given SNP, this underlying effect size, which we denote *η*, is not known. However, across all SNPs, *η* may be considered to be realisations of a random variable *H* whose distribution is mirrored by the distribution of observed effect sizes for the conditional phenotype. By integrating over this unknown true effect size for the conditional phenotype, allowance can be made for shared controls, and the ‘expected quantile’ of a p value for the principal phenotype, defined as P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ), can be calculated, as detailed in the Methods section.

We assume that *H* has a mixture distribution defined by two parameters (*π*_{0},*σ*^{2}), such that *H* = 0 with probability *π*_{0} and *H* ∼ *N*(0,*σ*_{2}) with probability 1−*π*_{0}. The parameters (*π*_{0},*σ*) are estimated from the observed distribution of effect sizes for the conditional phenotype. In order to show the effect of our p value adjustment, we simulated p values for 20,000 SNPs for a principal and conditional phenotype, with controls shared between simulated studies. All SNPs were null for the principal phenotype, and were variably null or non-null for the conditional phenotype with probability 0.9, 0.1 respectively. *Z* scores at non-null SNPs for the conditional phenotype were distributed as *N*(0,*σ*^{2}), as per our assumption. A value of *σ* = 3 was used, which was similar to the values of *σ* in real data estimated by our E-M algorithm.

We considered the set of simulated SNPs with p values for the conditional phenotype less than 0.05 (*Fig. 1*). In the absence of shared controls, we expect the distribution of *p*_{i} amongst this set to be uniform, and hence expect the black dots to lie along the x-y line. However, we see the principal p values are biased downward in this set (black dots, *Fig. 1*). Our computed expected quantile (blue dots) agrees closely with the observed quantile. In a sense, this constitutes ‘adjusting’ the p values for the principal phenotype so that the expected distribution is uniform under the null hypothesis. Software to generate this simulation is available at https://github.com/jamesliley/cFDR-common-controls.

Our formula can easily be adapted to arbitrary distributions of *H* at the cost of increased computational time, but the form of the distribution of *H* is not generally known. We show in *S1 Text* (section A) that, even for distributions of *H* which differ markedly from our assumption of normality, the error in the estimate is not large, and generally translates to a negligible difference in the set of SNPs declared non-null using the c F D R ̂ method. While our assumption has the potential to be anti-conservative if *H* is bimodal, nonparametric estimates for distributions of effect sizes suggest they have a uni-modal distribution centred on zero[*30*]. Reassuringly, our assumption is conservative if *H* has heavier tails than a normal.

## Comparison to split control approach

We compared Andreasson’s approach to SNP discovery which advocated splitting controls into non-overlapping subsets to our extended shared-control approach using a type 1 diabetes dataset with a total 12,175 cases and 15,171 controls. Controls and cases were each split into two sets (control sets had size 7,585 and 7,586, cases 6087 and 6088). ‘Split’ p values were computed using one set of controls and one set of cases and corresponding ‘shared’ p values were computed using the complete set of controls. As expected, more shared p values reached genome-wide significance than did split p values (*Fig. 2*).

We computed *cFDR* values by labelling one set of cases ‘conditional’ and the other ‘principal’ using the split-control p values using Andreasson’s approach and using the shared control p values using our method. For reference, we compared these to a naive application of Andreasson’s method on the shared-control p values (*Fig. 2B*). More SNPs can be declared significant according to *cFDR* using the shared-control than split-control approach at all reasonable thresholds, and naive application of Andreasson’s approach to shared-control p values again increases the number declared significant.

Because the quantity P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ) is systematically underestimated when using this naive method (by assuming it is equal to *p*_{i}) as shown in *S1 Text* (section C), it leads to a falsely low c F D R ̂. The increase in observed number of SNPs declared significant when using the naive method shows that it can indeed lead to false discoveries.

For principal phenotype p values in the range 5 × 10^{−6} - 5 × 10^{−8} - effectively the region from which ‘new’ SNPs may be discovered by cFDR rather than p value alone—the naive *cFDR* is frequently underestimated by 2–3 fold (*S1 Fig.*, left panel). For lower p values, the naive *cFDR* may underestimate by hundreds- or thousand-fold, with the potential fold underestimation increasing with decreasing p value (*S1 Fig.*, right panel). Because of the relatively high ratio of number of controls to number of cases, the correlation between effect sizes is lower in this constructed case (c. 0.22) than between most phenotypes in our study (c. 0.5). The underestimation of *cFDR* using the ‘naive’ method worsens with higher correlation, so we would expect that the fold-underestimate we see here is less severe than that which would be observed if applying this to other studies.

## An upper bound on the false discovery rate of all declared SNPs

An important property of our method is the control of the expected false discovery rate (FDR): the expected proportion of false positives among the SNPs found by our method. The p values at a SNP for the principal and conditional phenotype correspond to a point in the unit square. In this sense, we can define the expected FDR of a region *R* of the unit square as the ratio of the expected number of null SNPs whose p values are in *R* divided by the expected total number of SNPs whose p values are in *R*. From a result of Benjamini and Hochberg [*31*], the expected FDR when *R* is a rectangle with vertices at (0,0),(*p*_{i},0),(*p*_{i},*p*_{j}),(0,*p*_{j}) is at most α = cFDR^
( p i ∣ p j ). If we denote by *L* the closed region defined by the set of p value pairs *p*_{i},*p*_{j} such that cFDR^
( p i ∣ p j ) ≤ α, then *L* has the property that the FDR of any rectangle of this form contained within *L* is less than *α*.

However, the expected FDR over *L* is not necessarily bounded by *α* (*Fig. 3*). This can be seen most easily in the extreme scenario in which all non-null SNPs are concentrated in the lower left corner of the unit square, and all null SNPs are also null for the conditional phenotype. In this case, the expected number of null SNPs in a rectangle is proportional to its area, so L is the union of all rectangles of the form above of a given area containing the lower-left corner of the unit square; that is, a hyperbola. Clearly the area of L is larger than the area of these constituent rectangles, yet it contains the same number of non-null SNPs, so it has a lower FDR.

In the original method [*3*], SNPs were declared significant if they were contained within any rectangular regions with a cFDR^ value of less than 0.01. Our reasoning demonstrates that the expected false-discovery rate amongst all such SNPs was higher than 0.01. We can derive an upper bound for the expected FDR of *L* by considering *M*^{*}, the largest rectangle in *L*. We show in *S1 Text* (section B) that the bound may be expressed simply as v ( L ) v ( M * ) α *, where *v*() denotes the expected number of null SNPs contained within *L* or *M*^{*} (approximately the area of *L* and *M*^{*}) and *α*^{*} is the cFDR at the upper-right vertex of *M*^{*} (*Fig. 3*).

## Application to ten immune mediated diseases

We obtained summary statistics in the form of p values for ten immune mediated diseases from ImmunoBase (www.immunobase.org, accessed 19/3/14). For each pair of diseases, the number of shared controls was estimated according to the description of the control samples in each paper. The numbers of cases, controls and our estimated numbers of shared controls for each study are shown in *Table 1*. Uniform quality control criteria were applied to all SNPs, and the MHC region, which exhibits both strong LD and strong association with immune mediated diseases was excluded. P values were corrected within each trait for genomic inflation using a standard algorithm [*32*] applied to SNPs included on the ImmunoChip to replicate a GWAS study of reading and maths ability (Steve Eyre and Cathryn Lewis, personal communication), unlikely to be related to any immune mediated disease studied here.

P values for each principal phenotype were adjusted to *p*^{′} as described above in order to account for the effect of shared controls. For each ordered pair of phenotypes, a Q-Q plot was generated as per Andreasson et al [*3*]. A Q-Q plot is a graph of the observed distribution of a random variable against the expected distribution. We overlaid Q-Q plots for *log*_{10}(*p*^{′}) values for the principal phenotype for subsets of SNPs exhibiting successively smaller p values for the conditional phenotype. *Fig. 4* shows QQ plots for T1D conditional on RA and PSO; plots for all other pairwise comparisons may be found in *S4*–*S13* Figs. Notably, if lines shift further left with more stringent cutoffs on association with the conditional phenotype, then SNPs which are associated with the conditional phenotype are more likely to be associated with the principal phenotype, indicating pleiotropic effects of SNPs on the two phenotypes. In many cases, the Q-Q plots demonstrate considerable leftward shift with conditioning on association with a second disease, and we see strong evidence for pleiotropy for T1D conditioned on RA and little or no evidence for pleiotropy for T1D conditioned on PSO.

We estimated the unconditional and conditional false discovery rates, uFDR^
( p i ) and cFDR^
( p i ∣ p j ), at each SNP for each phenotype and each ordered pair of phenotypes respectively. *Fig. 5* shows cFDR^ for T1D conditioned on RA. The advantage gained by cFDR^ can be seen in the left-shift of the region in which a SNP can be declared significant (blue dots), corresponding to a higher p-value cutoff for significance for T1D among SNPs with low p values for RA. Indeed, if only SNPs with a p value for RA less than some threshold *ζ* are considered, a p value cutoff for significance for T1D is given by the leftmost border of the blue dots on the line *P*_{j} = *ζ*.

The degree of leftward shift in the Q-Q plots clearly contains information about the degree of pleiotropy between diseases. We defined a statistic summarizing some aspects of this evidence for pleiotropy and used it to visualise the set of pairwise relationships between diseases as a network (*Fig. 6*). The network encouragingly reflects several pathophysiological associations: UC is linked to CRO, and T1D to ATD. Strong linkage is also seen both ways for MS and PBC, and between T1D and RA, findings which can also be seen in the Q-Q plots (*S4*–*S13* Figs.). One way relationships suggest the presence of a larger total number of associated SNPs for the disease at the start of the arrow than at the end.

## Discovery of novel associations

The numbers of SNPs deemed significant for each phenotype by analysis using unconditional and conditional approaches are shown in *table 2*, with details in *S2*–*S11* Tables. cFDR^ allows certain SNPs with p values as high as 3 × 10^{−6} to be declared significant while controlling the false discovery rate at a relatively low value. Fifty-one of the 59 SNPs we identify uniquely through cFDR^ have previously been reported to be associated with the relevant disease through use of alternative significance thresholds, other genomic control procedures, other GWAS or additional samples not genotyped by ImmunoChip, a useful verification of our technique. Eight of the SNPs we discover uniquely through cFDR were in regions not previously known to be associated with the corresponding disease (*table 3*). These will require replication in independent samples to be declared truly associated, but they contain some potentially interesting signals, such as an association for RA at SNP rs72928038 near existing MS, ATD and T1D associations in *BACH2*, a transcriptional regulator involved in transcription repression and activation by *MAFK* [*33*]

The SNP rs1034290 in region 1p13.1, which we found to be associated with PBC, is in intron three of *CD58*, which is a surface receptor involved in binding and activation of T-lymphocytes. The protective effect of the MS-associated allele is postulated to arise from upregulation of the transcription factor *FOXP3* [*34*] and the patterns of association in the region suggest the two diseases may share a causal variant here (http://www.immunobase.org).

## Discussion

We have extended a technique for computing conditional Bayesian False Discovery Rates to GWAS for independent diseases with shared control groups. This technique enables improved detection of disease-associated SNPs compared to conventional methods. By enabling larger control groups for each study, our method uses data more efficiently than in corresponding study designs in which control groups are independent, and is applicable to a wider range of GWAS datasets for which only summary statistics are available.

Combination of GWAS by analysis of pleiotropy in this sense has several attractive advantages over single-phenotype analysis. The most obvious advantage is improved detection of disease-associated SNPs using GWAS without the need for additional samples. A secondary advantage arises from understanding of the pleiotropic structure between phenotypes: if a SNP is known to exhibit pleiotropy between two conditions, it may be causative for a shared risk factor or pre-disease state. Analysis of such SNPs has the potential to yield information on disease aetiology, with implications for preventative medicine and development of treatment.

A further potential use for this technique could be the genomic analysis of diseases with complex phenotypes. In many cases, distinction between two diseases may be difficult; for instance, Crohn’s disease and Ulcerative Colitis [*35*]. Additionally, many diseases, including narcolepsy (http://www.uptodate.com/contents/clinical-features-and-diagnosis-of-narcolepsy, accessed 20/6/14), are definitively diagnosed on clinical grounds. This implies that these diseases may constitute a range of biochemical and genetic states. Inclusion criteria based on objective biochemical grounds, such as that used for narcolepsy in the context of this paper [*24*] are unlikely to characterise all patients with these diseases, and conclusions drawn from studies will not necessarily be medically applicable to the whole patient population. Given this, diseases defined phenotypically with potential genomic diversity may be better analysed by separate consideration of biochemically-defined subtypes, with a collective analysis performed by a method such as cFDR^, avoiding the assumption that the genomic bases of disease subtypes are identical.

We identify a counter intuitive property that the FDR in the union of all regions with cFDR^ less than a given *α* may be greater than *α*, and propose a method to overcome this problem. Our methods for adjusting cutoffs to control FDR and account for multiple testing demonstrate the geometrical elegance of the theory of these techniques, with the possibility for further improvements and understanding. They are complex to apply, but could be much simplified if interest was directed to SNPs with conditional p values less than some threshold *p*_{0}. Our method would ensure that the expected false discovery rate at SNPs with cFDR^
( p i ∣ p 0 ) ≤ α would indeed be controlled at *α*. Our more complicated method to control FDR is necessary if the variable *p*_{j} is used in place of the constant *p*_{0}.

An important consideration in both our method and the original Andreasson method is that a cFDR^
( p i ∣ p j ) value which reaches significance does not constitute genome-wide evidence of association with the conditional phenotype *j*; indeed, the probability of association with the conditional phenotype relates to cFDR^
( p j ∣ p i ) and in general cFDR^
( p i ∣ p j ) ≠ cFDR^
( p j ∣ p i ). In some cases, where the principal p value is very close to genome-wide significance, even conditioning on *p*_{j} ≤ 0.5 can theoretically be enough to reach the relevant cFDR^ threshold. This is not a weakness of the cFDR^ method as such, but a consequence of using a discrete technique (a significance cutoff) on a variable which essentially continuous in two dimensions (cFDR^). Principal p values greater than 5 × 10^{−8} which can be declared significant conditioning on large conditional p value cutoffs correspond to an increase in the area of the region *L* (see *results* section), which is accounted for by our FDR-controlling method.

Our method enables improved detection of SNPs compared to analysis of unconditional FDR (principal p value alone). However, the improvement is smaller than that reported by Andreasson et al [*3*, *13*, *14*], who detected almost twice as many SNPs using *cFDR* as they would have detected with *uFDR*. This is expected for two reasons. Firstly, the gain in power from *cFDR* essentially comes from an increase in the total number of controls and the effective number of cases. If controls are shared, the only information gain can come from increasing the number of effective cases. Consequently, the difference in power between *cFDR* and *uFDR* will not be as large when controls are shared, although both outperform their counterparts when controls are split. Secondly, we were careful to use stringent cutoffs for FDR which were chosen to mirror the established genomewide significance threshold of *p* ≤ 5 × 10^{−8}, generally equivalent to a false discovery rate around 5 × 10^{−6} to 5 × 10^{−5}, compared to Andreasson et al who declared non-null all SNPs with cFDR^
< 0 . 01.

One alternative way to exploit pleiotropic relationships is by meta-analysing two related diseases together, as though the diseases were the same. Our method confers several advantages over this approach. The most important of these is that our method borrows strength from other SNPs according to the level of genome wide pleiotropy between diseases; that is, if the two GWAS suggest extensive pleiotropy (such as *Fig. 4* for T1D — RA), a low p value for a conditional phenotype will ‘sway’ our judgement of association with the principal phenotype more than the same p value for a conditional phenotype with poor pleiotropy (such as *Fig. 4*, for T1D — PSO). A meta-analysis would not distinguish these two scenarios. A secondary advantage of our technique is that SNP detection is not systematically weakened if the two diseases do not exhibit pleiotropy, as would be the case in meta-analysis; this arises because we are testing association with only one of the two phenotypes at a time.

## Methods

## Ethics statement

This paper re-analyses previously published datasets. All patient data were handled in accordance with the policies and procedures of the participating organisations.

## Datasets

We obtained SNP summary statistics from ten studies on autoimmune diseases from ImmunoBase (www.immunobase.org). Inclusion and exclusion criteria for the studies are described in detail in the original publications ([*20*–*28*, *28*]. Generally, some or all controls from different studies were obtained from common data sources, resulting in overlapping control groups. All studies used the ImmunoChip array [*17*].

P values for type 1 diabetes were from a meta-analysis of a case-control study and familial study using the transmission disequilibrium test (TDT). In order to calculate the correlation between *p* values for different diseases, we needed to calculate effective numbers of cases and controls for the combined T1D study. For a case control study, under the assumptions of Hardy-Weinberg and the null hypothesis, the variance of the log odds ratio may be expressed as

where

*n*

_{0}and

*n*

_{1}are the numbers of cases and controls and

*f*is the minor allele frequency in controls.

Given the standard error of a log OR for the TDT study, *σ̂*, and a minor allele frequency, we estimated *M* = *σ̂*^{2} *f*(1−*f*) for all ImmunoChip SNPs which did not show deviation from the null hypothesis (*p* > 0.5). The distribution of log(*M*) is shown in *S3 Fig*. By equating the median of *M* with n 0 + n 1 n 0 n 1, and assuming that each TDT family contributed the equivalent information to one control in a case-control study, ie *n*_{0} = 2943, we estimated an equivalent number of cases as 4126. This seemed reasonable, given that there are a total of 5505 (dependent) cases across those families.

SNPs were excluded on the basis of QC summaries calculated on 12,888 common controls: call rate less than 99%, minor allele frequency less than 0.02, or deviation from Hardy-Weinberg equilibrium (∣*Z*∣ > 5). Given the strong association of immune mediated diseases with the MHC and the extended LD in the region, we were concerned that MHC SNPs might cause inaccurate estimation of pleiotropy. We therefore excluded SNPs in a wide band around the MHC region on chromosome 6 (co-ordinates 24500000: 34800000, build NCBI36). After quality control, genotype data was available for at least one phenotype at a total of 110677 SNPs.

## Genomic control

P values were corrected for genomic inflation using a genomic control algorithm [*32*]. A set of SNPs known to be unassociated with autoimmune disease was obtained from the Wellcome Trust Case Control Consortium (WTCCC) study on reading and mathematics ability. These SNPs were pruned so that none were in LD with *r*^{2} > 0.2, and any SNPs within 500 kb of known autoimmune-associated regions were removed. The average degree of inflation was computed for each disease at the remaining 1761 SNPs, and all effect sizes and p values were adjusted accordingly.

## Procedures for uFDR and cFDR

We assume that the p-values for a phenotype *i* across all SNPs are instances of a random variable *P*_{i}. If *p*_{i} is an instance of this random variable corresponding to a SNP of interest, the unconditional false discovery rate *uFDR*(*p*_{i}) is defined as

where H 0 ( i ) is the null hypothesis that the SNP of interest is not associated with phenotype

*i*. Given a set of observed p values { p i 1,p i 2…p i N } for a phenotype

*i*at

*N*different SNPs, and an observed p value

*p*

_{i}for a SNP of interest, we estimate this quantity as

Because we make the approximation P r ( H 0 i ) = 1, the estimate uFDR^ is a an upwards-biased estimate of

*uFDR*; that is, its expected value is greater than the true

*uFDR*, making it a conservative estimator.

We compute the quantity (*1*) for each SNP at each phenotype, declaring any SNP for which uFDR^
( p i ) ≤ α as non-null for phenotype *i*. Defining *V* as the number of SNPs falsely declared non-null, *R* as the total number of SNPs declared non-null, and *Q* = *V*/*R*, a theorem of Benjamini and Hochberg [*31*] shows the expected false discovery rate *E*(*Q*) among SNPs with uFDR^
≤ α is less than *α*.

The *cFDR* constitutes a natural extension of this idea. We assume that the p-values for two phenotypes *i* and *j* across all SNPs are instances of a pair of random variables *P*_{i},*P*_{j}. If *p*_{i} and *p*_{j} are instances of these variables corresponding to a SNP of interest then the conditional false discovery rate *cFDR* is defined for the set of SNPs with p values for each phenotype less than or equal to those at this SNP (as per Andreasson et al [*3*]) as

The estimation of this quantity proceeds in a similar way to

*uFDR*. Given a set of observed p value pairs { ( p i 1,p j 1 ),( p i 2,p j 2 )… ( p i N,p j N ) } for two phenotypes

*i*and

*j*at

*N*different SNPs, and an observed p value pair (

*p*

_{i},

*p*

_{j}) for a SNP of interest, we define

*N*

_{1}as the number of p value pairs with

*P*

_{j}≤

*p*

_{j}, and estimate the

*cFDR*as

Again, this estimate is conservative, due to the approximation P r ( H 0 ( i ) ∣ P j ≤ p j ) = 1.

We compute the quantity (*2*) for each SNP at each pair of phenotypes, declaring any SNP for which cFDR^
( p i ∣ p j ) ≤ α ′ as non-null for phenotype *i*. However, as noted earlier, this does not guarantee that the expected false discovery rate amongst such SNPs is less than *α*^{′}. We show that the FDR is controlled at a higher level dependent on the region of the unit square defined by rectangles for which cFDR^
( x ∣ y ) ≤ α ′.

Our method here diverges from the original method proposed by Andreasson et al, in the use of the expected quantile P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ) in place of the p-value *p*_{i}. If studies share no controls, it can be reasonably assumed that, for a SNP which is null for phenotype *i*, the p values (*p*_{i},*p*_{j}) are independent, so p i ′ = P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ) = p i. This is the approach taken by Andreasson et al [*3*]. We propose a method for computing p i ′ when controls are shared between studies, and the independence assumption above is not valid.

Our approach is to compute the related quantity P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ),H η ( j ) ), where *η* is the (unobserved) effect size we would observe for a given SNP for phenotype *j* if the observed MAFs agreed exactly with the population MAFs for that SNP, and H η ( j ) is the hypothesis that *Z*_{j} ∼ *N*(*η*,1) for that SNP. This quantity can be thought of as the ‘expected quantile’ of *p*_{i}; that is, the proportion of p values we expect to be less than *p*_{i}.

## Computation of expected quantile

From the first part of (*2*), we have:

As per Andreasson et al [

*3*], the quantity P r ( H 0 ( i ) ∣ P j ≤ p j ) is set conservatively at 1, and the quantity

*Pr*(

*P*

_{i}≤

*p*

_{i}∣

*P*

_{j}≤

*p*

_{j}) is estimated empirically as the proportion of pairs of observed p values (p′i,p′j) with p′j≤pj which also satisfy p′i≤pi .

For a given SNP, let *η* denote the standardised mean allele frequency (MAF) difference; that is, the *Z* value we would compute if the observed MAFs agreed exactly with the population MAFs. We consider *η* for a random SNP as being an instance of a random variable *H*, and that the observed *z* value for that SNP *Z*∣*H* = *η* is distributed as

We further assume that

*H*follows a mixture distribution taking the value 0 with probability π 0 ( j ) and a normal

This implies

Thus, given the observed distribution of

*Z*

_{j}, the parameters π 0 ( j ) and

*σ*

_{j}may be estimated by an expectation - maximisation algorithm (https://gist.github.com/chr1swallace/11421212).

We assume as per Zaykin [*18*] that the distribution of pairs of observed z values (*Z*_{i},*Z*_{j}) for a single given SNP is bivariate normal. Denote by H η ( j ) the event that, for a given SNP, the values *Z*_{j} are distributed as *N*(*η*,1), with *η* depending on the SNP.

Under our assumption of the null hypothesis H 0 ( i ) for the principal phenotype and a population MAF difference corresponding to *η* for the conditional phenotype, we have

The correlation *ρ* arises from the shared controls between groups [*18*, *19*] and is asymptotically equal to

where

*N*

_{i}and

*N*

_{j}are the numbers of cases,

*N*

_{0i}and

*N*

_{0j}are the numbers of non-shared controls, and

*N*

_{0}is the number of shared controls for the original GWAS for the principal and conditional phenotypes respectively. There is good agreement with the asymptotic correlation when group sizes are greater than 100 [

*18*].

Given equations (*5*)–(*8*), the joint distribution of *Z*_{i} and *Z*_{j} can be computed under only the assumption H 0 ( i ). The value of the partial PDF of ( Z i,Z j ∣ H 0 ( i ) ) at (*x*,*y*) can be derived in a similar way to (*6*):

We now compute the final probability in equation (*3*). Define

as the probability of observing events

*X*for a particular SNP with true effect size

*η*(which may be 0, corresponding to the general null). Then,

If the distribution of *H* is estimable by other means, quantity (*11*) can be calculated numerically without the assumption that the non-null component of *H* be normally distributed, at the cost of higher computation time. Under our assumptions, equations (*6*) and (*9*) enable the fast computation of quantity (*11*) by normal CDFs; writing

we have

## Point expected quantile

Because the formula for P r ( P i ≤ p i ∣ P j ≤ p j,H 0 ( i ) ) is differentiable on the unit square, an expression for the expected quantile of *p*_{i} given an exact value for *p*_{j} can be computed by taking the partial derivative with respect to *p*_{j}:

where

where

*N*

_{μ,σ2}(

*x*) denotes the value of the normal

*μ*and variance

*σ*

^{2}at

*x*.

## Significance thresholds

Because uFDR^ values are monotonically related to p values, the widely accepted GWAS p value cutoff of 5 × 10^{−8} corresponds naturally to a cutoff for uFDR^. For each phenotype *i*, we set a significance threshold *β*^{i} for uFDR^
( p i ) as the lowest possible value of *γ* for which uFDR^
( p i ) ≤ γ ⇔ p i ≤ 5 × 10 − 8.

We then applied an analagous approach to cFDR^. For each pair of phenotypes (*i*,*j*), we set a significance threshold αji
as the lowest possible value of *γ*^{′} for which cFDR^
( p i ∣ p j ) ≤ γ ′ ⇔ p i ≤ 5 × 10 − 8. Given the distribution of *P*_{j}, it is possible that this could lead to declaring SNPs with *p*_{i} > 5 × 10^{−8},*p*_{j} ≈ 1 as significant. To avoid this, if αji
was larger (less stringent) than *β*^{i}, we set αji
= β i.

For each ordered pair of phenotypes (*i*,*j*), we declared all SNPs with cFDR^
( p i ∣ p j ) ≤ αji as non-null for phenotype *i*. This included all SNPs with uFDR^
( p i ) ≤ β i. We then used a technique described in *S1 Text* (section B) to compute upper bounds c j ( i ) on the false discovery rate amongst SNPs for which cFDR^
( p i,p j ) ≤ αji. For each phenotype, this gave nine upper bounds, corresponding to each of the nine conditional phenotypes.

## Network and heatmap representation of pleiotropy

We compared the degree of pleiotropy between diseases by considering how much the p-value threshold for significance for the principal phenotype changed when conditioning on a small p-value threshold for the conditional phenotype. We used the cFDR^ algorithm to compute the number p i j * such that P ( H 0 ( i ) ∣ P i ≤ p i j *,p j ≤ 5 × 10 − 6 ) = P ( H 0 ( i ) ∣ P i ≤ 5 × 10 − 8 ); that is, cFDR^ ( p i j * ∣ 5 × 10 − 6 ) = uFDR^ ( 5 × 10 − 8 ). We then considered the ratio p i j * / 5 × 10 − 8; that is, the fold increase in significance cutoff after conditioning.

We note that because of the fixed value of *p*_{j} = 5 × 10^{−6}, the expected false discovery rate amongst the set of SNPs which satisfy P ( H 0 ( i ) ∣ P i ≤ p i j *,p j ≤ 5 × 10 − 6 ) ≤ P ( H 0 ( i ) ∣ P i ≤ 5 × 10 − 8 ) is bounded above by P ( H 0 ( i ) ∣ P i ≤ 5 × 10 − 8 ), by the Benjamini-Hochberg result. Thus the expected false discovery rate amongst SNPs with P i ≤ p i j * and *P*_{j} ≤ 5 × 10^{−6} is bounded by the same value as the expected false discovery rate amongst SNPs with *P*_{i} ≤ 5 × 10^{−8}.

We visualised the ratio p i j * / 5 × 10 − 8 as a heatmap (*S2 Fig.*). We also produced a network (*Fig. 6*), with an edge from vertex *i* to vertex *j* if and only if, by conditioning on *P*_{j} ≤ 5 × 10^{−6}, the cutoff for significance for *P*_{i} could be increased from 5 × 10^{−8} to 4 × 10^{−7}. This cutoff was chosen as the smallest value such that the network was weakly connected; that is, each vertex had an arrow either to it or from it.

## Discovery of novel SNP associations

Cutoffs *β*^{i} were chosen such that uFDR^
( p i ) ≤ β i ⇔ p i ≤ 5 × 10 − 8. SNPs were deemed significant for a principal phenotype *i* if cFDR^
( p i ∣ p j ) ≤ α j i for any conditional phenotype *j* and α j i ≤ β i.

The expected false discovery rate amongst SNPs for which uFDR^
( p i ) ≤ β i is less than *β*^{i} due to a theorem of Benjamini and Hochberg. However, as discussed above and in *S1 Text* (section B), the expected false discovery rate amongst SNPs for which cFDR^
( p i ∣ p j ) ≤ α j i is not necessarily lower than α j i. For each ordered pair of phenotypes (*i*,*j*), an upper bound c j i was computed for the expected false discovery rates amongst SNPs with cFDR^
( p i ∣ p j ) ≤ α j i. The list of SNPs declared non-null for phenotype *i* was pruned to allow for linkage disequilibrium (LD) by listing all SNPs in increasing order of minj(cFDR^(pi|pj) and stepping through the list from left to right, at each stage removing all SNPs in LD with *r*^{2} ≥ 0.1 to the right of the current SNP. This ideally leads to the inclusion of at most one SNP from each LD block.

## Multiple testing

A multiple testing problem arises from considering p values for one disease conditioned separately on nine others. Specifically, if the criterion for declaring a SNP non-null for phenotype *i* is that cFDR^
( p i ∣ p j ) ≤αji for at least one of the nine possible values of *j*, then the FDR for all SNPs declared non-null will be greater than the FDR among the smaller set of SNPs for which cFDR^
( p i ∣ p j ) ≤αji for only one value of *j*, due to multiple testing.

However, this excess FDR is not enough to warrant a Bonferroni (Sidak) correction; the cFDR^
( p i ∣ p j ) values for a phenotype *i* are highly correlated, as all are in turn highly correlated with *p*_{i}. A Bonferroni correction tends to remove any advantage in SNP detection gained from cFDR^, though an advantage may still be seen when only considering one conditional phenotype *j*.

We opted to use a method proposed by Nyholt [*36*] to correct for multiple testing in SNPs with high LD. We estimated a correlation matrix Ω for potentially non-null cFDR values using Spearman’s rank correlation. The variance of the eigenvalues of Ω, *Var*(*λ*_{obs}), was computed and used to estimate the effective number of variables *M*_{eff} according to the equation

Note that *Var*(*λ*_{obs}) is between 9 (completely correlated variables, effectively one test) and 0 (completely uncorrelated variables, essentially a Bonferroni correction).

Denoting by n j i the number of SNPs with cFDR^
( p i ∣ p j ) ≤ α j i, corresponding to an upper bound on the FDR of c j i, an upper bound for the FDR among all SNPs declared significant for phenotype *i* was then computed as

intuitively, multiplying the expected average number of false discoveries across conditional phenotypes (c j i n j i) by the effective number of tests. Values of

*M*

_{eff}and c 0 i are shown in

*S1 Table*.

## Supporting Information

##### Zdroje

1. Speed D, Hemani G, Johnson MR, Balding DJ (2012) Improved heritability estimation from genome-wide SNPs. American Journal of Human Genetics 91: 1011–1021. doi: 10.1016/j.ajhg.2012.10.010 23217325

2. Cotsapas C, Voight BF, Rossin E, Lage K, Neale BM, et al. (2011) Pervasive sharing of genetic effects in autoimmune disease. PLOS Genetics 7: 1–8. doi: 10.1371/journal.pgen.1002254 21852963

3. Andreassen OA, Thompson WK, Schork AJ, Ripke S, Mattingsdal M, et al. (2013) Improved detection of common variants associated with schizophrenia and bipolar disorder using pleiotropy-informed conditional false discovery rate. PLOS Genetics 9(4). doi: 10.1371/journal.pgen.1003455 23637625

4. Hirschhorn JN, Daly MJ (2005) Genome-wide association studies for common diseases and complex traits. Nature Reviews Genetics 6: 95–108. doi: 10.1038/nrg1521 15716906

5. Glazier AM, Nadeau JH, Aitman TJ (2002) Finding genes that underlie complex traits. Science 298: 2345–2349. doi: 10.1126/science.1076641 12493905

6. Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, et al. (2011) Genome partitioning of genetic variation for complex traits using common snps. Nature Genetics 43(6): 519–525. doi: 10.1038/ng.823 21552263

7. Smyth D, Plagnol V, Walker NM, Cooper JD, Downes K, et al. (2008) Shared and distinct genetic variants in type 1 diabetes and celiac disease. New England Journal of Medicine 359: 2767–2777. doi: 10.1056/NEJMoa0807917 19073967

8. Ferkingstad E, Frigessi A, Rue H, Thorleifsson G, Kong A (2008) Unsupervised empirical Bayesian multiple testing with external covariates. Annals of Applied Statistics 2: 714–735. doi: 10.1214/08-AOAS158

9. Paaby AB, Rockman MV (2013) The many faces of pleiotropy. Trends in Genetics 29(2): 66–73. doi: 10.1016/j.tig.2012.10.010 23140989

10. Sivakumaran S, Agakov F, Theodoratou E, Theodoratou E, Prendergast JG, et al. (2011) Abundant pleiotropy in human complex diseases and traits. American Journal of Human Genetics 89: 607–618. doi: 10.1016/j.ajhg.2011.10.004 22077970

11. Lichtenstein P, Yip BH, Bjöork C, Pawitan Y, Cannon TD, et al. (2009) Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet 373(9659): 234–239. doi: 10.1016/S0140-6736(09)60072-6 19150704

12. Hasstedt SJ, Hanis CL, Das SK, Elbein SC, The American Diabetes Association GENNID Study Group (2011) Pleiotropy of type 2 diabetes with obesity. Journal of Human Genetics 56(7): 491–495. doi: 10.1038/jhg.2011.46 21525879

13. Andreassen OA, McEvoy LK, Thompson WK, Wang Y, Reppe S, et al. (2014) Identifying common genetic variants in blood pressure due to poly genic pleiotropy with associated phenotypes. Hypertension 63. doi: 10.1161/HYPERTENSIONAHA.113.02077 24396023

14. Andreasson OA, Harbo HF, Wang Y, Thompson WK, Schork AJ, et al. (2014) Genetic pleiotropy between multiple sclerosis and schizophrenia but not bipolar disorder: differential involvement of immune-related gene loci. Molecular psychiatry: 1–8. doi: 10.1038/mp.2013.195

15. Efron B, Tibshirani R (2002) Empirical Bayes methods and false discovery rates for microarrays. Genetic Epidemiology 23: 70–86. doi: 10.1002/gepi.1124 12112249

16. The Wellcome trust case control consortium (2007) Genome-wide association study of 14000 cases of seven common diseases and 3000 shared controls. Nature 447: 661–678. doi: 10.1038/nature05911 17554300

17. Cortes A, Brown MA (2011) Promise and pitfalls of the ImmunoChip. Arthritis Research and Therapy 13. doi: 10.1186/ar3204 21345260

18. Zaykin DV, Kozbur DO (2010) P-value based analysis for shared controls design in genome-wide association studies. Genetic Epidemiology 34: 725–738. doi: 10.1002/gepi.20536 20976797

19. Lin D, Sullivan PF (2009) Meta-analysis of genome-wide association studies with overlapping subjects. American Journal of Human Genetics 85: 862–872. doi: 10.1016/j.ajhg.2009.11.001 20004761

20. Onengut-Gumuscu S, Chen WM, Burren O, Cooper NJ, Quinlan AR, et al. (2014) Fine Mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nature Genetics In Press.

21. Cooper JD, Simmonds MJ, Walker NM, Burren O, Brand OJ, et al. (2012) Seven newly identified loci for autoimmune thyroid disease. Human Molecular Genetics 21: 5202–5208. doi: 10.1093/hmg/dds357 22922229

22. Trynka G, Hunt KA, Bockett NA, Romanos J, Mistry V, et al. (2012) Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nature Genetics 43: 1193–1201. doi: 10.1038/ng.998 22057235

23. International Multiple Sclerosis Genetics Consortium (IMSGC), Beecham AH, Patsopoulos NA, Xifara DK, Davis MF, et al. (2013) Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nature Genetics 45. doi: 10.1038/ng.2770 24076602

24. Faraco J, Lin L, Kornum BR, Kenny EE, Trynka G, et al. (2013) Immunochip study implicates antigen presentation to T cells in narcolepsy. PLOS Genetics 9. doi: 10.1371/journal.pgen.1003270

25. Liu JZ, Almarri MA, Gaffney DJ, Mells GF, Jostins L, et al. (2012) Dense fine-mapping study identifies new susceptibility loci for primary biliary cirrhosis. Nature Genetics 44: 1137–1141. doi: 10.1038/ng.2395 22961000

26. Tsoi LC, Spain SL, Knight J, Ellinghaus E, Stuart PE, et al. (2012) Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity. Nature Genetics 44: 1341–1350. doi: 10.1038/ng.2467 23143594

27. Eyre S, Bowes J, Diogo D, Lee A, Barton A, et al. (2012) High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nature Genetics 44: 1336–1342. doi: 10.1038/ng.2462 23143596

28. Jostins L, Ripke S, Weersma RK, Duerr RH, McGovern DPB, et al. (2012) Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 491: 119–124. doi: 10.1038/nature11582 23128233

29. Cotsapas C, Hafler DA (2013) Immune-mediated disease genetics: the shared basis of pathogenesis. Trends in Immunology 34(1): 22–26. doi: 10.1016/j.it.2012.09.001 23031829

30. Park JH, Wacholder S, Gail MH, Peters U, Jacobs KB, et al. (2010) Estimation of effect size distribution from genome-wide association studies and implications for future discoveries. Nat Genet 42: 570–575. doi: 10.1038/ng.610 20562874

31. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, series B (methodological) 57: 289–300.

32. Devlin B, Roeder K, Wasserman L (2001) Genomic control, a new approach to genetic-based association studies. Theoretical Population Biology 60: 155–166. doi: 10.1006/tpbi.2001.1542 11855950

33. Dubois PCA, Trynka G, Franke L, Hunt KA, Romanos J, et al. (2010) Multiple common variants for celiac disease influencing immune gene expression. Nature Genetics 42: 295–302. doi: 10.1038/ng.543 20190752

34. De Jager PL, Beacher-Allan C, Maier LM, Arthur AT, Ottoboni L, et al. (2009) The role of the CD58 locus in multiple sclerosis. Proceedings of the National Academy of Sciences USA 106: 5264–5269. doi: 10.1073/pnas.0813310106

35. Farmer M, Petras RE, Hunt LE, Janosky JE, Galandiuk S (2000) The importance of diagnostic accuracy in colonic inflammatory bowel disease. The American Journal of Gastroenterology 95: 3184–3188. doi: 10.1111/j.1572-0241.2000.03199.x 11095339

36. Nyholt D (2004) A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. American Journal of Human Genetics 74: 765–769. doi: 10.1086/383251 14997420

##### Štítky

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

### PLOS Genetics

2015 Číslo 2

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

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

- Genomic Selection and Association Mapping in Rice (): Effect of Trait Genetic Architecture, Training Population Composition, Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite, Tropical Rice Breeding Lines
- Discovery of Transcription Factors and Regulatory Regions Driving Tumor Development by ATAC-seq and FAIRE-seq Open Chromatin Profiling
- Evolutionary Signatures amongst Disease Genes Permit Novel Methods for Gene Prioritization and Construction of Informative Gene-Based Networks
- Membrane Recognition and Dynamics of the RNA Degradosome