Volume 8 Supplement 1

## Genetic Analysis Workshop 18: Human sequence data in extended pedigrees

- Proceedings
- Open Access
- Published:

# Family-based association test using normal approximation to gene dropping null distribution

*BMC Proceedings*
**volume 8**, Article number: S18 (2014)

## Abstract

We derive the analytical mean and variance of the score test statistic in gene-dropping simulations and approximate the null distribution of the test statistic by a normal distribution. We provide insights into the gene-dropping test by decomposing the test statistic into two components: the first component provides information about linkage, and the second component provides information about fine mapping under the linkage peak. We demonstrate our theoretical findings by applying the gene-dropping test to the simulated data set from Genetic Analysis Workshop 18 and comparing its performance with existing population and family-based association tests.

## Background

When testing genotype-phenotype association using individuals from extended families, one has to account for correlations in genotypes and/or phenotypes between related individuals. One simple and effective method to account for genotype correlations is to simulate the null genotype distribution by gene dropping [1], which is simulating founder alleles according to estimated allele frequencies and dropping these alleles down the pedigrees according to random segregation of gametes (i.e., Mendel's first law). The gene-dropping method is straightforward to implement (e.g., implemented in by Allen-Brady *et al* [2]) and applies to all pedigree structures, but it is computationally intensive and thus is impractical to use when dealing with millions of single-nucleotide polymorphisms (SNPs).

In this article, we derive the analytical mean and variance of the score test statistic under the gene-dropping setting and approximate the gene-dropping null distribution of the test statistic by a normal distribution with the analytically derived mean and variance. Using this normal approximation, the gene-dropping test becomes computationally efficient and can be easily applied to millions of SNPs.

Furthermore, we provide insights into the gene-dropping test by decomposing the test statistic into two components: the first component resembles a quantity frequently used in variance-component based linkage tests and provides information for linkage, and the second component provides information for fine mapping under the linkage peak. Rabinowitz and Laird [3], among others, have pointed out the subtle distinction between two types of null hypotheses in family-based association analysis: the null hypothesis of no linkage and no association versus the null hypothesis of no association in the presence of linkage. To test the latter, one needs to condition on the inheritance ${S}_{\tau}$ vector at the test locus [3]. Our decomposition provides an explicit separation of linkage and association information in a family-based study.

We compare the performance of the gene-dropping test (using normal approximation) to association tests using only unrelated individuals and to the family-based association test in the software program FBAT [3] by analyzing Genetic Analysis Workshop 18 (GAW18) simulated data set.

## Methods

### Preprocessing of genotype data

We analyzed SNPs from chromosome 3 only. At each of the SNPs, we performed Pearson's chi-squared test for the Hardy-Weinberg equilibrium using 142 unrelated individuals. We excluded SNPs that yielded a *p*-value smaller than 10^{−4} from our analysis. In the gene-dropping test, we excluded SNPs with estimated minor allele frequency (MAF) smaller than 0.001.

### Preprocessing of phenotype data

We focused on the analysis of the quantitative trait systolic blood pressure (SBP) in the simulated data set 1. The true simulation model was known to us [4]. When testing association between genotype doses and trait values (see later discussion), we include factors AGE, SEX, and AGE by SEX interaction as covariates (${Z}_{k}$'s in equation [1]). Including BPMED as a covariate will overcompensate because BPMED is a consequence of SBP level. Instead, we estimated the effect of BPMED from a regression model with only individuals with hypertension. Because BPMED was randomly assigned to individuals with hypertension, the BPMED effect estimated this way will not be biased by its correlation with SBP. We then adjusted the trait values $Y$ by subtracting the estimated BPMED effect.

### Score tests of genotype-phenotype association using unrelated individuals

At locus $\tau $, we consider a quantitative trait model

and test the null hypothesis ${\beta}_{\tau}=0$. In equation (1), $Y$ is the vector of trait values (SBP adjusted for the BPMED effect), $\mu $ is a constant vector of baseline mean trait values, coefficients ${\alpha}_{k}$ represent the effects of the covariates ${Z}_{k},k=1,\dots ,K,$ (e.g., AGE, SEX and AGE by SEX interaction) on trait values, ${X}_{\tau}$ is the vector of genotype doses (the number of minor alleles possessed by each individual) at locus $\tau $, and the coefficient ${\beta}_{\tau}$ represents the effect size of a single allele. The fitted value of ${\beta}_{\tau}$ will reflect the collective effect of all causal SNPs that are in linkage disequilibrium (LD) with the test SNP $\tau $[5].

Let $\u0176$ and ${\stackrel{\u2322}{X}}_{\tau}$ be the vectors of fitted values after regressing the $Y$ and ${X}_{\tau}$ on measured covariates ${Z}_{k}$'s. The score statistic [6, 7] for testing genotype-trait association at a single SNP $\tau $ is $u={X}_{\tau}^{\prime}R$, where $R=Y-\stackrel{\u2322}{Y}$ is the vector of residuals. Under the null hypothesis of no association, the variance of $u$ is estimated by

where $Z=\left(1,{Z}_{1},\dots ,{Z}_{K}\right)$ and ${s}_{YY}$ is the sample variance of the residual trait values ($1$ is a vector of ones) [6]. To test association, ${u}^{2}/v$ is compared with a ${\chi}_{1}^{2}$ distribution.

### Family-based association test by gene dropping

When related individuals are used to compute the score test statistic $u={X}_{\tau}^{\prime}R$, components of ${X}_{\tau}$ can be dependent, and the variance estimator (2) is no longer valid. One can account for correlations between components in ${X}_{\tau}$ by simulating the null distribution of ${X}_{\tau}$ using gene dropping. We now derive the analytical mean and variance of $u$ under the gene-dropping setting. In the score test using unrelated individuals, we treat $R$ as random, and ${X}_{\tau}$ can be viewed as either random or fixed. In a gene-dropping simulation, $R$ is held fixed, and ${X}_{\tau}$ is random.

Let $i,j$ index individuals $\left(i,j=1,\dots ,n\right)$ and let ${X}_{\tau}={\left({X}_{1},\dots ,{X}_{n}\right)}^{\prime}$ and $R={\left({R}_{1},\dots ,{R}_{n}\right)}^{\prime}$. The expected value of $u$ is $\sum _{i=1}^{n}E\left({X}_{i}\right){R}_{i}$ and ${X}_{i}={P}_{i}+{M}_{i}$, where ${P}_{i}$$\left({M}_{i}\right)$ is 1 if the paternal (maternal) allele is the minor allele and 0 otherwise. So $E\left({X}_{i}\right)$ is twice the MAF ${f}_{\tau}$ at SNP $\tau $ and is the same for all individuals and thus $E\left(u\right)=2{f}_{\tau}\sum {R}_{i}=0$ because ${R}_{i}$'s are residuals from a linear regression model with intercept. The variance of $u$ is $E\left({u}^{2}\right)={R}^{\prime}E\mathsf{\text{(}}{X}_{\tau}{X}_{\tau}^{\prime})R.$ The $\left(i,j\right)$ th element in $E\mathsf{\text{(}}{X}_{\tau}{X}_{\tau}^{\prime})$ is $E\left({X}_{i}{X}_{j}\right)=E\left({P}_{i}{P}_{j}+{P}_{i}{M}_{j}+{M}_{i}{P}_{j}+{M}_{i}{M}_{j}\right).$${P}_{i},{M}_{i},{P}_{j},{M}_{j}$ are all Bernoulli random variables with probability ${f}_{\tau}$, and any two of them are identical if the corresponding alleles are identity-by-descent (IBD) and are independent otherwise [8]. Let ${\varphi}_{ij}$ be the number of IBD pairs among the four pairs of alleles ${P}_{i}{P}_{j},{P}_{i}{M}_{j},{M}_{i}{P}_{j},{M}_{i}{M}_{j}$. The value of ${\varphi}_{ij}$ at locus $\tau $ is determined by the inheritance vector ${S}_{\tau}$, which summarizes whether the paternal or the maternal allele is passed from the parent to the child in each meiosis [9]. Given the inheritance vector ${S}_{\tau}$,

(e.g., $E\left({P}_{i}{P}_{j}\right)=E\left({P}_{i}^{2}\right)={f}_{\tau}$ if ${P}_{i}$ and ${P}_{j}$ correspond to IBD alleles and $E\left({P}_{i}{P}_{j}\right)=E\left({P}_{i}\right)E\left({P}_{j}\right)={f}_{\tau}^{2}$ if ${P}_{i}$ and ${P}_{j}$ correspond to non-IBD alleles). In a gene-dropping simulation, the inheritance vector ${S}_{\tau}$ is randomly sampled among all possible inheritance vectors. The expected number of IBD alleles shared between $i$ and $j$, $\mathsf{\text{E}}\left({\varphi}_{ij}\left({S}_{\tau}\right)\right),$ over all possible inheritance vectors is four times the kinship coefficient ${\psi}_{ij}\mathsf{\text{:E}}\left({\varphi}_{ij}\left({S}_{\tau}\right)\right)=4{\psi}_{ij}$. The kinship coefficients are determined by pedigree structures. The expected value of ${X}_{i}{X}_{j}$ in a gene-dropping simulation is thus $E\left(E\left({X}_{i}{X}_{j}|{S}_{\tau}\right)\right)=4{\psi}_{ij}\left({f}_{\tau}-{f}_{\tau}^{2}\right)+4{f}_{\tau}^{2}.$ Letting $\Phi \left({S}_{\tau}\right)=\left({\varphi}_{ij}\right)$ be the matrix of IBD counts and $\Psi $ be the matrix of kinship coefficients, we can rewrite the above as:

where $J$ is a matrix of all ones. Because ${R}^{\prime}JR=0$ for residuals from a linear regression model with an intercept, the variance of $u$ under gene dropping is ${v}_{gd}={R}^{\prime}E\left({X}_{\tau}{X}_{\tau}^{\prime}\right)R=4{R}^{\prime}\Psi R\left({f}_{\tau}-{{f}_{\tau}}^{2}\right)$ if unconditional on the inheritance vector ${S}_{\tau}$, and is ${v}_{\tau}={R}^{\prime}E\left({X}_{\tau}{X}_{\tau}^{\prime}|{S}_{\tau}\right)R={R}^{\prime}\Phi \left({S}_{\tau}\right)R\left({f}_{\tau}-{{f}_{\tau}}^{2}\right)$ if conditional on the inheritance vector ${S}_{\tau}$ (holding ${S}_{\tau}$ fixed). We can approximate the gene-dropping null distribution of $u$ by a normal distribution with mean $0$ and variance ${v}_{gd}$, and compute the gene-dropping *p*-value by comparing $t={u}^{2}/{v}_{gd}$ with a ${\chi}_{1}^{2}$ distribution. To test association in the presence of linkage, one needs to condition on the inheritance ${S}_{\tau}$ vector at $\tau $[3] and use ${v}_{\tau}$. In practice, ${S}_{\tau}$ is not observable, but we estimate ${v}_{\tau}$ by drawing Markov chain Monte Carlo (MCMC) samples of ${S}_{\tau}$ based on observed genotypes in the pedigrees using MORGAN (http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml) [10].

## Results

### Theoretical findings

In a gene-dropping simulation, the analytical mean of the score statistic $u={{X}_{\tau}}^{\prime}R$ is $0$. The variance of the score statistic is ${R}^{\prime}\Phi \left({\mathsf{\text{S}}}_{\tau}\right)R\left({f}_{\tau}-{f}_{\tau}^{2}\right)$ if conditional on the inheritance vector (i.e., holding the inheritance vector fixed during gene-dropping simulation) and is $4{R}^{\prime}\Psi R\left({f}_{\tau}-{f}_{\tau}^{2}\right)$ if unconditional on the inheritance vector. The normal approximation is justified by the central limit theorem because the test statistic is additive over pedigrees. Its performance depends on the number, sizes, and structure of pedigrees and on MAF at the test locus. The approximation may not be accurate for extremely small *p*-values. However, the rankings of the *p*-values will not change.

We can decompose the unconditional gene-dropping test statistic into two components:

The first component can be used as a test statistic for detecting association in the presence of linkage (i.e., fine mapping under a linkage peak) because the denominator is the variance of $u$ conditional upon the observed IBD sharing. The second component provides information about linkage. The kinship coefficients in $\Psi $ are determined by pedigree structure, so ${R}^{\prime}\Psi R$ is a constant in a gene-dropping simulation. ${R}^{\prime}\Phi \left({S}_{\tau}\right)R=\sum _{ij}{r}_{i}{r}_{j}{\varphi}_{ij}\left({S}_{\tau}\right)$ measures the correlation between trait value similarity $\left({r}_{i}{r}_{j}\right)$ and IBD sharing $\left({\varphi}_{ij}\right)$ at locus $\tau $ across all pairs of individuals in a pedigree. This correlation is expected to be stronger if there is stronger linkage between $\tau $ and a true causal locus. Therefore, ${R}^{\prime}\Phi \left({S}_{\tau}\right)R$ can be used as a test statistic to detect linkage, with null distribution obtained by gene-dropping simulations. In a gene-dropping simulation, the inheritance vectors are simulated as if they were from a marker unlinked to any potential causal loci. ${R}^{\prime}\Phi \left({S}_{\tau}\right)R$ resembles similar quantities that are frequently used in linkage analysis methods such as the well-known Haseman-Elston regression [11] as well as many variance components or generalized estimating equation-based methods [12].

### Simulation results

We performed a genome-wide association studies (GWAS) score test using 142 unrelated individuals, the family-based association test using FBAT [3], and the gene-dropping test on SNPs on chromosome 3 (FBAT and the gene-dropping test used 847 individuals from 20 pedigrees). Table 1 summarizes the *p*-value ranks that each test assigns the true causal SNPs. The gene-dropping test for fine mapping (conditional on the inheritance vector) performs very similarly to the unconditional gene-dropping test, so its results are omitted. It is seen that the gene-dropping tests can quickly identify a few true causal SNPs within a short list of top findings. However, if we allow more false positives by considering a greater number of the most significant SNPs, other methods start to pick up true causal SNPs and eventually have a result similar to gene dropping.

Figure 1 shows the physical positions and negative log *p*-values of the top 500 SNPs identified by each of the three tests, as well as the negative log *p*-values of the linkage test based on the linkage component of the gene-dropping test statistic.

We also examined adjusting for population stratification by fitting the first two principal components of genetic variation [13] as covariates in the regression model (1). The *p*-values resulting from this expanded model differed negligibly from the original model. The ranks in Table 1 were essentially unchanged by this adjustment.

## Discussion

### Comparison between genome-wide association studies, FBAT and gene-dropping test

FBAT splits each pedigree into nuclear families. In each nuclear family, FBAT uses information from the offspring while conditioning on the parental marker genotypes. In contrast, GWAS uses information in unrelated individuals. The two methods use almost "orthogonal" sources of information. There is almost no correlation between the log *p*-values from these two methods (Table 2). In contrast, the gene-dropping test applies to multigeneration pedigrees and uses information from all individuals: the gene-dropping test extracts information from founders by resimulating founder genotypes and from offspring by resimulating inheritance vectors.

It is also possible to derive the analytical mean and variance of the test statistic in the gene-dropping test where we permute the founder alleles rather than resimulate the founder alleles. FBAT is more robust to population stratification by conditioning on founder genotypes. The gene-dropping test can gain similar robustness by restricting permutations to founder alleles within each family.

It is somewhat surprising that the gene-dropping test did not outperform GWAS given that it uses more individuals. One possible interpretation is that the effect of LD is stronger when more individuals are used. As we can see in Figure 1, the signals detected by the gene-dropping test come in bigger clusters. In other words, many SNPs ranked high by the gene-dropping test might be in LD with one or more of the causal SNPs.

### Separating linkage and association signals

The gene-dropping test captures both linkage and association signals. One can decompose the test statistic into a linkage component and an association component.

The association component corresponds to testing association in the presence of linkage, which requires one to condition on the true inheritance vector at the test locus. Our results through MCMC approximation show that whether or not to condition on the inheritance vector actually does not make a big difference for this data set because the variance of the test statistic with conditioning only differs slightly from the variance of the test statistic without conditioning. This conclusion might be dependent on the structure of the pedigree.

The linkage component, however, clearly provides valuable information. The linkage signal is stronger in most regions containing causal SNPs. It is obvious that the linkage curve can help eliminate many of the false association signals in this study. It would be interesting to investigate how to use the linkage information more effectively in the future.

## References

- 1.
MacCluer JW, Vandeburg JL, Read B, Ryder OA: Pedigree analysis by computer simulation. Zoo Biol. 1986, 5: 149-160.

- 2.
Allen-Brady K, Wong J, Camp NJ: PedGenie: an analysis approach for genetic association testing in extended pedigrees and genealogies of arbitrary size. BMC Bioinformatics. 2006, 7: 209-10.1186/1471-2105-7-209.

- 3.
Rabinowitz D, Laird N: A unified approach to adjusting association tests for population admixture with arbitrary pedigree structure and arbitrary missing marker information. Hum Hered. 2000, 50: 211-223. 10.1159/000022918.

- 4.
Almasy L, Dyer T, Peralta J, Jun G, Fuchsberger C, Almeida M, Kent JW, Fowler S, Duggirala R, Blangero J: Data for Genetic Analysis Workshop 18: human whole genome sequence, blood pressure, and simulated phenotypes in extended pedigrees. BMC Proc. 2014, 8 (suppl 2): S2-

- 5.
Di Y, Mi G, Sun L, Dong R, Zhu H, Peng L: Power of association tests in the presence of multiple causal variants. BMC Proc. 2011, 5 (suppl 9): S63-10.1186/1753-6561-5-S9-S63.

- 6.
Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA: Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002, 70: 425-434. 10.1086/338688.

- 7.
Clayton D, Chapman J, Cooper J: Use of unphased multilocus genotype data in indirect association studies. Genet Epidemiol. 2004, 27: 415-428. 10.1002/gepi.20032.

- 8.
Lange K, Westlake J, Spence MA: Extensions to pedigree analysis. III. Variance components by the scoring method. Ann Hum Genet. 1976, 39: 485-491. 10.1111/j.1469-1809.1976.tb00156.x.

- 9.
Lander ES, Green P: Construction of multilocus genetic linkage maps in humans. Proc Natl Acad Sci USA. 1987, 84: 2363-2367. 10.1073/pnas.84.8.2363.

- 10.
Thompson EA: MCMC in the analysis of genetic data on pedigrees. Markov Chain Monte Carlo: Innovations and Applications. Edited by: Liang F, Wang J-S, Kendall W. Lecture Note Series of the IMS, National University of Singapore. World Scientific Co Pte Ltd, Singapore, 183-216.

- 11.
Elston RC, Buxbaum S, Jacobs KB, Olson JM: Haseman and Elston revisited. Genet Epidemiol. 2000, 19: 1-17. 10.1002/1098-2272(200007)19:1<1::AID-GEPI1>3.0.CO;2-E.

- 12.
Chen WM, Broman KW, Liang KY: Quantitative trait linkage analysis by generalized estimating equations: unification of variance components and Haseman-Elston regression. Genet Epidemiol. 2004, 26: 265-272. 10.1002/gepi.10315.

- 13.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.

## Acknowledgements

We thank the GAW18 organizers. The GAW18 whole genome sequence data were provided by the Type 2 Diabetes Genetic Exploration by Next-generation sequencing in Ethnic Samples (T2D-GENES) Consortium, which is supported by National Institutes of Health (NIH) grants U01 DK085524, U01 DK085584, U01 DK085501, U01 DK085526, and U01 DK085545. The other genetic and phenotypic data for GAW18 were provided by the San Antonio Family Heart Study and San Antonio Family Diabetes/Gallbladder Study, which are supported by NIH grants P01 HL045222, R01 DK047482, and R01 DK053889. The GAW is supported by NIH grant R01 GM031575.

This article has been published as part of *BMC Proceedings* Volume 8 Supplement 1, 2014: Genetic Analysis Workshop 18. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcproc/supplements/8/S1. Publication charges for this supplement were funded by the Texas Biomedical Research Institute.

## Author information

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

All authors participated in analysis of the data. YJ, SE, and YD conceived of the study and drafted the manuscript. All authors participated in the critical revision of the manuscript and gave final approval of the article.

## Rights and permissions

## About this article

#### Published

#### DOI

### Keywords

- Kinship Coefficient
- Genetic Analysis Workshop
- Linkage Peak
- Causal SNPs
- Pedigree Structure