Several methods have been proposed to impute genotypes at untyped markers using observed genotypes and genetic data from a reference panel. We used the Genetic Analysis Workshop 16 rheumatoid arthritis case-control dataset to compare the performance of four of these imputation methods: IMPUTE, MACH, PLINK, and fastPHASE. We compared the methods' imputation error rates and performance of association tests using the imputed data, in the context of imputing completely untyped markers as well as imputing missing genotypes to combine two datasets genotyped at different sets of markers. As expected, all methods performed better for single-nucleotide polymorphisms (SNPs) in high linkage disequilibrium with genotyped SNPs. However, MACH and IMPUTE generated lower imputation error rates than fastPHASE and PLINK. Association tests based on allele "dosage" from MACH and tests based on the posterior probabilities from IMPUTE provided results closest to those based on complete data. However, in both situations, none of the imputation-based tests provide the same level of evidence of association as the complete data at SNPs strongly associated with disease.

The nuclear factor-kappaB (NF-kappaB) family of transcription factors regulates the expression of a variety of genes involved in apoptosis and immune response. We examined relationships between genotypes at five NF-kappaB subunits (NFKB1, NFKB2, REL, RELA, and RELB) and variable expression levels of 15 NF-kappaB regulated proteins with heritability greater than 0.40: BCL2A1, BIRC2, CD40, CD44, CD80, CFLAR, CR2, FAS, ICAM1, IL15, IRF1, JUNB, MYC, SLC2A5, and VCAM1. SNP genotypes and expression phenotypes from pedigrees of Utah residents with ancestry from northern and western Europe were provided by Genetic Analysis Workshop 15 and supplemented with additional genotype data from the International HapMap Consortium. We conducted association, linkage, and family-based association analyses between each candidate gene and the 15 heritable expression phenotypes. We observed consistent results in association and linkage analyses of the NFKB1 region (encoding p50) and levels of FAS and IRF1 expression. FAS is a cell surface protein that also belongs to the TNF-receptor family; signals through FAS are able to induce apoptosis. IRF1 is a member of the interferon regulatory transcription factor family, which has been shown to regulate apoptosis and tumor-suppression. Analyses in the REL region (encoding c-Rel) revealed linkage and association with CD40 phenotype. CD40 proteins belong to the tumor necrosis factor (TNF)-receptor family, which mediates a broad variety of immune and inflammatory responses. We conclude that variation in the genes encoding p50 and c-Rel may play a role in NF-kappaB-related transcription of FAS, IRF1, and CD40.

Several methods to identify tagging single-nucleotide polymorphisms (SNPs) are in common use for genetic epidemiologic studies; however, there may be loss of information when using only a subset of SNPs. We sought to compare the ability of commonly used pairwise, multimarker, and haplotype-based tagging SNP selection methods to detect known associations with quantitative expression phenotypes. Using data from HapMap release 21 on unrelated Utah residents with ancestors from northern and western Europe (CEPH-Utah, CEU), we selected tagging SNPs in five chromosomal regions using ldSelect, Tagger, and TagSNPs. We found that SNP subsets did not substantially overlap, and that the use of trio data did not greatly impact SNP selection. We then tested associations between HapMap genotypes and expression phenotypes on 28 CEU individuals as part of Genetic Analysis Workshop 15. Relative to the use of all SNPs (n = 210 SNPs across all regions), most subset methods were able to detect single-SNP and haplotype associations. Generally, pairwise selection approaches worked extremely well, relative to use of all SNPs, with marked reductions in the number of SNPs required. Haplotype-based approaches, which had identified smaller SNP subsets, missed associations in some regions. We conclude that the optimal tagging SNP method depends on the true model of the genetic association (i.e., whether a SNP or haplotype is responsible); unfortunately, this is often unknown at the time of SNP selection. Additional evaluations using empirical and simulated data are needed.

No abstract available

Due to the growing need to combine data across multiple studies and to impute untyped markers based on a reference sample, several analytical tools for imputation and analysis of missing genotypes have been developed. Current imputation methods rely on single imputation, which ignores the variation in estimation due to imputation. An alternative to single imputation is multiple imputation. In this paper, we assess the variation in imputation by completing both single and multiple imputations of genotypic data using MACH, a commonly used hidden Markov model imputation method. Using data from the North American Rheumatoid Arthritis Consortium genome-wide study, the use of single and multiple imputation was assessed in four regions of chromosome 1 with varying levels of linkage disequilibrium and association signals. Two scenarios for missing genotypic data were assessed: imputation of untyped markers and combination of genotypic data from two studies. This limited study involving four regions indicates that, contrary to expectations, multiple imputations may not be necessary.

No abstract available

The incorporation of disease-associated covariates into studies aiming to identify susceptibility genes for complex human traits is a challenging problem. Accounting for such covariates in genetic linkage and association analyses may help reduce the genetic heterogeneity inherent in these complex phenotypes. For Genetic Analysis Workshop 15 (GAW15) Problem 3 simulated data, our goal was to compare the power of several two-stage study designs to identify rheumatoid arthritis-related genes on chromosome 9 (disease severity), 11 (IgM), and 18 (anti-cyclic citrinullated protein), with knowledge of the answers. Five study designs incorporating an initial linkage step, followed by a case-selection scheme and case-control association analysis by logistic regression, were considered. The linkage step was either qualitative-trait linkage analysis as implemented in MERLIN-nonparametric linkage (NPL), or quantitative-trait locus analysis as implemented in MERLIN-REGRESS. A set of cases representing either one case from each available family, one case per linked family (NPL >/= 0), or one case from each family identified by ordered-subset analysis was chosen for comparison with the full set of 2000 simulated controls. As expected, the performance of these study designs depended on the disease model used to generate the data, especially the simulated allele frequency difference between cases and controls. The quantitative trait loci analysis performed well in identifying these loci, and the power to identify disease-associated alleles was increased by using ordered-subset analysis as a case selection tool.

After genetic linkage has been identified for a complex disease, the next step is often fine-mapping by association analysis, using single-nucleotide polymorphisms (SNPs) within a linkage region. If a SNP shows evidence of association, it is useful to know whether the linkage result can be explained in part or in full by the candidate SNP. The genotype identity-by-descent sharing test (GIST) and linkage and association modeling in pedigrees (LAMP) are two methods that were specifically proposed to address this question. GIST determines whether there is significant correlation between family-specific weights, defined by the presence of a tentatively associated allele in affected siblings, and family-specific nonparametric linkage scores. LAMP constructs a pedigree likelihood function of the marker data conditional on the trait data, and implements three likelihood ratio tests to characterize the relationship between the candidate SNP and the disease locus. The goal of our study was to compare the two approaches and evaluate their ability to identify disease-associated SNPs in the Genetic Analysis Workshop 15 (GAW15) simulated data. Our results can be summarized as follows: 1) GIST is simple and fast but, as a test of association, did not perform well in the GAW15 data, especially with adjustment for multiple testing; 2) as a test of association, the LAMP-LE test performs best when the linkage evidence is strong, or when there is at least moderate linkage disequilibrium between the candidate SNP and the trait locus. We conclude that LAMP is more flexible and reliable to use in practice.