Validity, efficiency, and robustness of a family-based test of association

We propose a new test of linkage in the presence of allelic association that uses all available information in a sample of nuclear families, including parental phenotypes, genotypes from both affected and unaffected siblings, and families with homozygous parents. The test is based on the conditional framework developed by Rabinowitz and Laird [2000: Hum Hered 50:211-223] and is thus immune to population stratification and can be applied to families with any pattern of missing information. The test statistic is a conditional likelihood ratio based on a standard two-point linkage model with allelic association, where parameters are estimated from the sample. Through a simulation study, we determined that the proposed test has near optimal power for a wide range of scenarios, outperforming FBAT both when data were complete and when parental genotypes were missing, although differences between the two tests diminish as the genetic effect is reduced. To assess robustness, we also evaluated the performance of the tests under scenarios with population stratification and found that although there is a loss of efficiency, our proposed test remains a strong competitor to FBAT.     

X-LRT: a likelihood approach to estimate genetic risks and test association with X-linked markers using a case-parents design

Recently, there has been interest in family-based tests of association to identify X-chromosome genes. However, none of the approaches allow for estimation of genetic risks. We propose a likelihood approach to estimate disease-related marker relative risks and test genotype association using a case-parents design. The test uses nuclear families with a single affected proband and allows additional siblings and missing parental genotypes. Extension to a haplotype test is based on assumptions of random mating and multiplicative penetrance. We investigate power and type I error rate of the likelihood-based test, using simulated data and apply our method to marker data from the monoamine oxidase A&B genes in families with Parkinson disease. We show how efficiency with missing parental information can be improved with additional sibling genotype information. Our likelihood approach offers great flexibility for testing different penetrance relationships within and between sexes. In addition, estimation of disease-related marker relative risks provides a measure of the magnitude of X-linked genetic effects on complex disorders.     

Tests for genetic association using family data.

We use likelihood-based score statistics to test for association between a disease and a diallelic polymorphism, based on data from arbitrary types of nuclear families. The Nonfounder statistic extends the transmission disequilibrium test (TDT) to accommodate affected and unaffected offspring, missing parental genotypes, phenotypes more general than qualitative traits, such as censored survival data and quantitative traits, and residual correlation of phenotypes within families. The Founder statistic compares observed or inferred parental genotypes to those expected in the general population. Here the genotypes of affected parents and those with many affected offspring are weighted more heavily than unaffected parents and those with few affected offspring. We illustrate the tests by applying them to data on a polymorphism of the SRD5A2 gene in nuclear families with multiple cases of prostate cancer. We also use simulations to compare the power of these family-based statistics to that of the score statistic based on Cox's partial likelihood for censored survival data, and find that the family-based statistics have considerably more power when there are many untyped parents. The software program FGAP for computing test statistics is available at http://www.stanford.edu/dept/HRP/epidemiology/FGAP.     

Ordered-subset analysis (OSA) for family-based association mapping of complex traits

Association analysis provides a powerful tool for complex disease gene mapping. However, in the presence of genetic heterogeneity, the power for association analysis can be low since only a fraction of the collected families may carry a specific disease susceptibility allele. Ordered-subset analysis (OSA) is a linkage test that can be powerful in the presence of genetic heterogeneity. OSA uses trait-related covariates to identify a subset of families that provide the most evidence for linkage. A similar strategy applied to genetic association analysis would likely result in increased power to detect association. Association in the presence of linkage (APL) is a family-based association test (FBAT) for nuclear families with multiple affected siblings that properly infers missing parental genotypes when linkage is present. We propose here APL-OSA, which applies the OSA method to the APL statistic to identify a subset of families that provide the most evidence for association. A permutation procedure is used to approximate the distribution of the APL-OSA statistic under the null hypothesis that there is no relationship between the family-specific covariate and the family-specific evidence for allelic association. We performed a comprehensive simulation study to verify that APL-OSA has the correct type I error rate under the null hypothesis. This simulation study also showed that APL-OSA can increase power relative to other commonly used association tests (APL, FBAT and FBAT with covariate adjustment) in the presence of genetic heterogeneity. Finally, we applied APL-OSA to a family study of age-related macular degeneration, where cigarette smoking was used as a covariate.     

CAPL: a novel association test using case-control and family data and accounting for population stratification

The recent successes of GWAS based on large sample sizes motivate combining independent datasets to obtain larger sample sizes and thereby increase statistical power. Analysis methods that can accommodate different study designs, such as family-based and case-control designs, are of general interest. However, population stratification can cause spurious association for population-based association analyses. For family-based association analysis that infers missing parental genotypes based on the allele frequencies estimated in the entire sample, the parental mating-type probabilities may not be correctly estimated in the presence of population stratification. Therefore, any approach to combining family and case-control data should also properly account for population stratification. Although several methods have been proposed to accommodate family-based and case-control data, all have restrictions. Most of them require sampling a homogeneous population, which may not be a reasonable assumption for data from a large consortium. One of the methods, FamCC, can account for population stratification and uses nuclear families with arbitrary number of siblings but requires parental genotype data, which are often unavailable for late-onset diseases. We extended the family-based test, Association in the Presence of Linkage (APL), to combine family and case-control data (CAPL). CAPL can accommodate case-control data and families with multiple affected siblings and missing parents in the presence of population stratification. We used simulations to demonstrate that CAPL is a valid test either in a homogeneous population or in the presence of population stratification. We also showed that CAPL can have more power than other methods that combine family and case-control data.     

Method for using complete and incomplete trios to identify genes related to a quantitative trait

A number of tests for linkage and association with qualitative traits have been developed, with the most well-known being the transmission/disequilibrium test (TDT). For quantitative traits, varying extensions of the TDT have been suggested. The quantitative trait approach we propose is based on extending the log-linear model for case-parent trio data (Weinberg et al. [1998] Am. J. Hum. Genet. 62:969-978). Like the log-linear approach for qualitative traits, our proposed polytomous logistic approach for quantitative traits allows for population admixture by conditioning on parental genotypes. Compared to other methods, simulations demonstrate good power and robustness of the proposed test under various scenarios of the genotype effect, distribution of the quantitative trait, and population stratification. In addition, missing parental genotype data can be accommodated through an expectation-maximization (EM) algorithm approach. The EM approach allows recovery of most of the lost power due to incomplete trios.     

Informative-transmission disequilibrium test (i-TDT): combined linkage and association mapping that includes unaffected offspring as well as affected offspring

To date, there is no test valid for the composite null hypothesis of no linkage or no association that utilizes transmission information from heterozygous parents to their unaffected offspring as well as the affected offspring from ascertained nuclear families. Since the unaffected siblings also provide information about linkage and association, we introduce a new strategy called the informative-transmission disequilibrium test (i-TDT), which uses transmission information from heterozygous parents to all of the affected and unaffected offspring in ascertained nuclear families and provides a valid chi-square test for both linkage and association. The i-TDT can be used in various study designs and can accommodate all types of independent nuclear families with at least one affected offspring. We show that the transmission/disequilibrium test (TDT) (Spielman et al. [1993] Am. J. Hum. Genet. 52:506-516) is a special case of the i-TDT, if the study sample contains only case-parent trios. If the sample contains only affected and unaffected offspring without parental genotypes, the i-TDT is equivalent to the sibship disequilibrium test (SDT) (Horvath and Laird [1998] Am. J. Hum. Genet. 63:1886-1897. In addition, the test statistic of i-TDT is simple, explicit and can be implemented easily without intensive computing. Through computer simulations, we demonstrate that power of the i-TDT can be higher in many circumstances compared to a method that uses affected offspring only. Applying the i-TDT to the Framingham Heart Study data, we found that the apolipoprotein E (APOE) gene is significantly linked and associated with cross-sectional measures and longitudinal changes in total cholesterol.     

A method for identifying genes related to a quantitative trait, incorporating multiple siblings and missing parents

When studying either qualitative or quantitative traits, tests of association in the presence of linkage are necessary for fine-mapping. In a previous report, we suggested a polytomous logistic approach to testing linkage and association between a di-allelic marker and a quantitative trait locus, using genotyped triads, consisting of an individual whose quantitative trait has been measured and his or her two parents. Here we extend that approach to incorporate marker information from entire nuclear families. By computing a weighted score function instead of a maximum likelihood test, we allow for both an unspecified correlation structure between siblings and "informative" family size. Both this approach and our original approach allow for population admixture by conditioning on parental genotypes. The proposed method allows for missing parental genotype data through a multiple imputation procedure. We use simulations based on a population with admixture to compare our method to a popular non-parametric family-based association test (FBAT), testing the null of no association in the presence of linkage.     

Family-based tests for associating haplotypes with general phenotype data: application to asthma genetics

We provide a general purpose family-based testing strategy for associating disease phenotypes with haplotypes when phase may be ambiguous and parental genotype data may be missing. These tests for linkage and association can be used in candidate gene studies with tightly linked markers. Our proposed weighted conditional approach extends the method described in Rabinowitz and Laird to multiple markers. It is attractive because it provides haplotype tests for family-based studies that are efficient and robust to population admixture, phenotype distribution specification, and ascertainment based on phenotypes. It can handle missing parental genotypes and/or missing phase in both offspring and parents. It yields either haplotype-specific (univariate) tests or multi-haplotype (global) tests. This extension has been implemented in the freely available software haplotype FBAT. We used the haplotype FBAT program to test for associations between asthma phenotypes and single nucleotide polymorphisms (SNPs) in the beta-2 adrenergic receptor gene. Whereas no single SNP showed significant association with asthma diagnosis or bronchodilator responsiveness (quantitative trait), a haplotype-based global test found a highly significant association with asthma diagnosis (P value <0.00005) and the measure of bronchodilator responsiveness (P value =0.016).     

Testing association in the presence of linkage using the GRE and multiple markers

It has recently been shown that testing for association in the presence of linkage using a score test based on a gamma random effects (GRE) model is substantially more powerful than using the Family-Based Association Test. A reason for the increased power lies in better specification of the within family correlation structure, induced by linkage. The GRE, as presented in (Jonasdottir et al. 2007 Genet Epidemiol. 31:528-540), only considers one marker at a time and does not readily handle missing parental information. Here we extend the GRE to incorporate information from more than one marker. This extension leads to a haplotype GRE test and also to efficient handling of missing data on parental genotypes. We show that the haplotype GRE, the H-GRE, is substantially more powerful than the haplotype FBAT, the Haplotype-Based-Association Test. We demonstrate the usefulness of the extended GRE, by reanalyzing the collaborative study on the genetics of alcoholism data, allowing for missing parental information.     

Genotype-based association test for general pedigrees: the genotype-PDT

Many family-based tests of linkage disequilibrium (LD) are based on counts of alleles rather than genotypes. However, allele-based tests may not detect interactions among alleles at a single locus that are apparent when examining associations with genotypes. Family-based tests of LD based on genotypes have been developed, but they are typically valid as tests of association only in families with a single affected individual. To take advantage of families with multiple affected individuals, we propose the genotype-pedigree disequilibrium test (geno-PDT) to test for LD between marker locus genotypes and disease. Unlike previous tests for genotypic association, the geno-PDT is valid in general pedigrees. Simulations to compare the power of the allele-based PDT and geno-PDT reveal that under an additive model, the allele-based PDT is more powerful, but that the geno-PDT can have greater power when the genetic model is recessive or dominant. Perhaps the most important property of the geno-PDT is the ability to test for association with particular genotypes, which can reveal underlying patterns of association at the genotypic level. These genotype-specific tests can be used to suggest possible underlying genetic models that are consistent with the pattern of genotypic association. This is illustrated through an application to a candidate gene analysis of the MLLT3 gene in families with Alzheimer disease. The geno-PDT approach for testing genotypes in general family data provides a useful tool for identifying genes in complex disease, and partitioning individual genotype contributions will help to dissect the influence of genotype on risk.     

Gene‐Based Association Analysis for Censored Traits Via Fixed Effect Functional Regressions

Genetic studies of survival outcomes have been proposed and conducted recently, but statistical methods for identifying genetic variants that affect disease progression are rarely developed. Motivated by our ongoing real studies, here we develop Cox proportional hazard models using functional regression (FR) to perform gene‐based association analysis of survival traits while adjusting for covariates. The proposed Cox models are fixed effect models where the genetic effects of multiple genetic variants are assumed to be fixed. We introduce likelihood ratio test (LRT) statistics to test for associations between the survival traits and multiple genetic variants in a genetic region. Extensive simulation studies demonstrate that the proposed Cox RF LRT statistics have well‐controlled type I error rates. To evaluate power, we compare the Cox FR LRT with the previously developed burden test (BT) in a Cox model and sequence kernel association test (SKAT), which is based on mixed effect Cox models. The Cox FR LRT statistics have higher power than or similar power as Cox SKAT LRT except when 50%/50% causal variants had negative/positive effects and all causal variants are rare. In addition, the Cox FR LRT statistics have higher power than Cox BT LRT. The models and related test statistics can be useful in the whole genome and whole exome association studies. An age‐related macular degeneration dataset was analyzed as an example.     

Genetic association tests with age at onset

Many diseases or traits exhibit a varying age at onset. Recent data examples of prostate cancer and childhood diabetes show that compared to simply treating the disease outcome as affected vs. unaffected, incorporation of age-at-onset information into the transmission/disequilibrium type of test (TDT) does not appear to change the results much. In this paper, we evaluate the power of TDT as a function of age at onset, and show that age-at-onset information is most useful when the disease is common, or the relative risk associated with the high-risk genotype varies with age. Moreover, an extremely old unaffected subject can contribute substantially to the power of the TDT, sometimes as much as old-onset subjects. We propose a modified test statistic for testing no association between the marker at the candidate locus and age at onset. The simulation study was conducted to evaluate the finite sample properties of proposed and the TDT test statistics under various sampling schemes for trios of parents and offspring, as well as for sibling clusters where unaffected siblings were used as controls.     

Testing for association with a case-parents design in the presence of genotyping errors

Genotyping errors can create a problem for the analysis of case-parents data because some families will exhibit genotypes that are inconsistent with Mendelian inheritance. The problem with correcting Mendelian inconsistent genotype errors by regenotyping or removing families in which they occur is that the remaining unidentified genotype errors can produce excess type I (false positive) error for some family-based tests for association. We address this problem by developing a likelihood ratio test (LRT) for association in a case-parents design that incorporates nuisance parameters for a general genotype error model. We extend the likelihood approach for a single SNP to include short haplotypes consisting of 2 or 3 SNPs. The extension to haplotypes is based on assumptions of random mating, multiplicative penetrances, and at most a single genotype error per family. For a single SNP, we found, using Monte Carlo simulation, that type I error rate can be controlled for a number of genotype error models at different error rates. Simulation results suggest the same is true for 2 and 3 SNPs. In all cases, power declined with increasing genotyping error rates. In the absence of genotyping errors, power was similar whether nuisance parameters for genotype error were included in the LRT or not. The LRT developed here does not require prior specification of a particular model for genotype errors and it can be readily computed using the EM algorithm. Consequently, this test may be generally useful as a test of association with case-parents data in which Mendelian inconsistent families are observed.     

The v-MFG test: investigating maternal, offspring and maternal-fetal genetic incompatibility effects on disease and viability

The MFG test is a family-based association test that detects genetic effects contributing to disease in offspring, including offspring allelic effects, maternal allelic effects and MFG incompatibility effects. Like many other family-based association tests, it assumes that the offspring survival and the offspring-parent genotypes are conditionally independent provided the offspring is affected. However, when the putative disease-increasing locus can affect another competing phenotype, for example, offspring viability, the conditional independence assumption fails and these tests could lead to incorrect conclusions regarding the role of the gene in disease. We propose the v-MFG test to adjust for the genetic effects on one phenotype, e.g., viability, when testing the effects of that locus on another phenotype, e.g., disease. Using genotype data from nuclear families containing parents and at least one affected offspring, the v-MFG test models the distribution of family genotypes conditional on offspring phenotypes. It simultaneously estimates genetic effects on two phenotypes, viability and disease. Simulations show that the v-MFG test produces accurate genetic effect estimates on disease as well as on viability under several different scenarios. It generates accurate type-I error rates and provides adequate power with moderate sample sizes to detect genetic effects on disease risk when viability is reduced. We demonstrate the v-MFG test with HLA-DRB1 data from study participants with rheumatoid arthritis (RA) and their parents, we show that the v-MFG test successfully detects an MFG incompatibility effect on RA while simultaneously adjusting for a possible viability loss.     

Combining unrelated family studies to improve the power of genetic association test

Family-based studies provide powerful inferences regarding associations between genetic variants and risks, but have limitations. Since very often, the availability of the parental genotypes can pose a problem for using family-based design, especially when the disease of interest has a late age of onset. To improve the efficiency of the studies, a popular approach is to reconstruct the missing genotypes from the genotypes of their offspring and correct the biases resulting from the reconstruction. In this paper, the author shows that two or more unrelated family studies, for the same candidate marker but different diseases, can also be combined to construct a more efficient test for association analysis. The usual case-control study with parental genotypes is a special case of the data discussed here. The author used a simulation study to compare the performance of the new method with other well-known methods. The results showed that the new test has an advantage of having larger power when there is no effect of population stratification between two study samples. However, if there is effect of population stratification between the two samples, the new test still maintains the expected type I error rate and has comparable power performance. Since the unrelated family studies not for the disease of interest are often readily accessible with minimal cost, the proposed method has practical value. The new approach can also be easily modified to allow for missing parental data.     

Testing association in the presence of linkage--a powerful score for binary traits

We present a score for testing association in the presence of linkage for binary traits. The score is robust to varying degrees of linkage, and it is valid under any ascertainment scheme based on trait values as well as under population stratification. The score test is derived from a mixed effects model where population level association is modeled using a fixed effect and where correlation among related individuals is allowed for by using log-gamma random effects. The score, as presented in this paper, does not assume full information about the inheritance pattern in families or parental genotypes. We compare the score to the semi-parametric family-based association test (FBAT), which has won ground because of its flexible and simple form. We show that a random effects formulation of co-inheritance can improve the power substantially. We apply the method to data from the Collaborative Study on the Genetics of Alcoholism. We compare our findings to previously published results.     

Efficient use of siblings in testing for linkage and association

Tests of linkage and association between a disease and either a candidate gene or marker allele can be based on sibships with at least one affected and one unaffected sibling. However, specialized techniques are required to account for within-sibship correlation if some sibships contain more than one affected or more than one unaffected sib. In this paper, we propose Within Sibship Paired Resampling (WSPR), a technique that is designed to test the null hypothesis of no linkage or no association, even when sibships contain variable numbers of sibs. One repeatedly generates data subsets based on randomly sampling one affected and one unaffected sibling from each sibship, and each subset is analyzed individually. Then, evidence is combined by averaging results across these resampled data sets, applying a variance expression that implicitly accounts for the correlation among siblings. While the general WSPR procedure allows for numerous testing strategies, we describe two in detail. Simulation results for scenarios with varying degrees of population stratification demonstrate good power for the WSPR testing methods compared to the sib TDT (S-TDT) and the sibship disequilibrium test (SDT).     

A powerful association test of multiple genetic variants using a random‐effects model

There is an emerging interest in sequencing‐based association studies of multiple rare variants. Most association tests suggested in the literature involve collapsing rare variants with or without weighting. Recently, a variance‐component score test [sequence kernel association test (SKAT)] was proposed to address the limitations of collapsing method. Although SKAT was shown to outperform most of the alternative tests, its applications and power might be restricted and influenced by missing genotypes. In this paper, we suggest a new method based on testing whether the fraction of causal variants in a region is zero. The new association test, T _REM, is derived from a random‐effects model and allows for missing genotypes, and the choice of weighting function is not required when common and rare variants are analyzed simultaneously. We performed simulations to study the type I error rates and power of four competing tests under various conditions on the sample size, genotype missing rate, variant frequency, effect directionality, and the number of non‐causal rare variant and/or causal common variant. The simulation results showed that T _REM was a valid test and less sensitive to the inclusion of non‐causal rare variants and/or low effect common variants or to the presence of missing genotypes. When the effects were more consistent in the same direction, T _REM also had better power performance. Finally, an application to the Shanghai Breast Cancer Study showed that rare causal variants at the FGFR2 gene were detected by T _REM and SKAT, but T _REM produced more consistent results for different sets of rare and common variants. Copyright © 2013 John Wiley & Sons, Ltd.     

Genotype relative-risks and association tests for nuclear families with missing parental data

The development of a new method for testing the association of genetic markers with disease is presented. This approach is applicable when sampling nuclear families with one or more affected siblings and when neither, one, or both parents are missing marker genotype data. All siblings, affected and not affected, are used to probabilistically infer the missing parental marker data. A likelihood ratio statistic, which treats marker allele frequencies as nuisance parameters, is presented to test whether all marker relative risks are equal to one (i.e., no marker association). This approach offers a solution to test for marker associations when parents are difficult to obtain. © 1997 Wiley-Liss, Inc.