Linkage mapping methods applied to the COGA data set: presentation Group 4 of Genetic Analysis Workshop 14

Presentation Group 4 participants analyzed the Collaborative Study on the Genetics of Alcoholism data provided for Genetic Analysis Workshop 14. This group examined various aspects of linkage analysis and related issues. Seven papers included linkage analyses, while the eighth calculated identity-by-descent (IBD) probabilities. Six papers analyzed linkage to an alcoholism phenotype: ALDX1 (four papers), ALDX2 (one paper), or a combination both (one paper). Methods used included Bayesian variable selection coupled with Haseman-Elston regression, recursive partitioning to identify phenotype and covariate groupings that interact with evidence for linkage, nonparametric linkage regression modeling, affected sib-pair linkage analysis with discordant sib-pair controls, simulation-based homozygosity mapping in a single pedigree, and application of a propensity score to collapse covariates in a general conditional logistic model. Alcoholism linkage was found with > or =2 of these approaches on chromosomes 2, 4, 6, 7, 9, 14, and 21. The remaining linkage paper compared the utility of several single-nucleotide polymorphism (SNP) and microsatellite marker maps for Monte Carlo Markov chain combined oligogenic segregation and linkage analysis, and analyzed one of the electrophysiological endophenotypes, ttth1, on chromosome 7. Linkage was found with all marker sets. The last paper compared the multipoint IBD information content of several SNP sets and the microsatellite set, and found that while all SNP sets examined contained more information than the microsatellite set, most of the information contained in the SNP sets was captured by a subset of the SNP markers with approximately 1-cM marker spacing. From these papers, we highlight three points: a 1-cM SNP map seems to capture most of the linkage information, so denser maps do not appear necessary; careful and appropriate use of covariates can aid linkage analysis; and sources of increased gene-sharing between relatives should be accounted for in analyses.     

Parent-of-origin, imprinting, mitochondrial, and X-linked effects in traits related to alcohol dependence: presentation Group 18 of Genetic Analysis Workshop 14

The participants of Presentation Group 18 of Genetic Analysis Workshop 14 analyzed the Collaborative Study on the Genetics of Alcoholism data set to investigate sex-specific effects for phenotypes related to alcohol dependence. In particular, the participants looked at imprinting (which is also known as parent-of-origin effect), differences between recombination fractions for the two sexes, and mitochondrial and X-chromosomal effects. Five of the seven groups employed newly developed or existing methods that take imprinting into account when testing for linkage, or test for imprinting itself. Single-marker and multipoint analyses were performed for microsatellite as well as single-nucleotide polymorphism markers, and several groups used a sex-specific genetic map in addition to a sex-averaged map. Evidence for paternal imprinting (i.e., maternal expression) was consistently obtained by at least two groups at genetic regions on chromosomes 10, 12, and 21 that possibly harbor genes responsible for alcoholism. Evidence for maternal imprinting (which is equivalent to paternal expression) was consistently found at a locus on chromosome 11. Two groups applied extensions of variance components analysis that model a mitochondrial or X-chromosomal effect to latent class variables and electrophysiological traits employed in the diagnosis of alcoholism. The analysis, without using genetic markers, revealed mitochondrial or X-chromosomal effects for several of these traits. Accounting for sex-specific environmental variances appeared to be crucial for the identification of an X-chromosomal factor. In linkage analysis using marker data, modeling a mitochondrial variance component increased the linkage signals obtained for autosomal loci.     

Multiple testing in high-throughput sequence data: experiences from Group 8 of Genetic Analysis Workshop 17

The use of high-throughput sequence data in genetic epidemiology allows the investigation of common and rare variants in the entire genome, thus increasing the amount of information and the potential number of statistical tests performed within one study. As a consequence, the problem of multiple testing may become even more pressing than in previous studies. As an important challenge, the exact number of statistical tests depends on the actual statistical method used. Furthermore, many statistical approaches for the analysis of sequence data require permutation. Thus it may be difficult to also use permutation to estimate correct type I error levels as in genome-wide association studies. In view of this, a separate group at Genetic Analysis Workshop 17 was formed with a focus on multiple testing. Here, we present the approaches used for the workshop. Apart from tackling the multiple testing problem, the new group focused on different issues. Some contributors developed and investigated modifications of existing collapsing methods. Others aimed at improving the identification of functional variants through a reduction and analysis of the underlying data dimensions. Two research groups investigated the overall accumulation of rare variation across the genome and its value in predicting phenotypes. Finally, other investigators left the path of traditional statistical analyses by reversing null and alternative hypotheses and by proposing a novel resampling method. We describe and discuss all these approaches.     

Genetic analysis of phenotypes derived from longitudinal data: Presentation Group 1 of Genetic Analysis Workshop 13

The participants of Presentation Group 1 used the GAW13 data to derive new phenotypes, which were then analyzed for linkage and, in one case, for association to the genetic markers. Since the trait measurements ranged over longer time periods, the participants looked at the time dependence of particular traits in addition to the trait itself. The phenotypes analyzed with the Framingham data can be roughly divided into 1) body weight-related traits, which also include a type 2 diabetes progression trait, and 2) traits related to systolic blood pressure. Both trait classes are associated with metabolic syndrome. For traits related to body weight, linkage was consistently identified by at least two participating groups to genetic regions on chromosomes 4, 8, 11, and 18. For systolic blood pressure, or its derivatives, at least two groups obtained linkage for regions on chromosomes 4, 6, 8, 11, 14, 16, and 19. Five of the 13 participating groups focused on the simulated data. Due to the rather sparse grid of microsatellite markers, an association analysis for several traits was not successful. Linkage analysis of hypertension and body mass index using LODs and heterogeneity LODs (HLODs) had low power. For the glucose phenotype, a combination of random coefficient regression models and variance component linkage analysis turned out to be strikingly powerful in the identification of a trait locus simulated on chromosome 5. Haseman-Elston regression methods, applied to the same phenotype, had low power, but the above-mentioned chromosome 5 locus was not included in this analysis.     

Genetic association with rheumatoid arthritis-Genetic Analysis Workshop 15: summary of contributions from Group 2

The papers in presentation group 2 of Genetic Analysis Workshop 15 (GAW15) conducted association analyses of rheumatoid arthritis data. The analyses were carried out primarily in the data provided by the North American Rheumatoid Arthritis Consortium (NARAC). One group conducted analyses in the data provided by the Canadian Rheumatoid Arthritis Genetics Study (CRAGS). Analysis strategies included genome-wide scans, the examination of candidate genes, and investigations of a region of interest on chromosome 18q21. Most authors employed relatively new methods, proposed extensions of existing methods, or introduced completely novel methods for aspects of association analysis. There were several common observations; a group of papers using a variety of methods found stronger association, on chromosomes 6 and 18 and in candidate gene PTPN22 among women with early onset. Generally, models that considered haplotypes or multiple markers showed stronger evidence for association than did single marker analyses.     

Comparison of single-nucleotide polymorphisms and microsatellite markers for linkage analysis in the COGA and simulated data sets for Genetic Analysis Workshop 14: Presentation Groups 1, 2, and 3

The papers in presentation groups 1-3 of Genetic Analysis Workshop 14 (GAW14) compared microsatellite (MS) markers and single-nucleotide polymorphism (SNP) markers for a variety of factors, using multiple methods in both data sets provided to GAW participants. Group 1 focused on data provided from the Collaborative Study on the Genetics of Alcoholism (COGA). Group 2 focused on data simulated for the workshop. Group 3 contained analyses of both data sets. Issues examined included: information content, signal strength, localization of the signal, use of haplotype blocks, population structure, power, type I error, control of type I error, the effect of linkage disequilibrium, and computational challenges. There were several broad resulting observations. 1) Information content was higher for dense SNP marker panels than for MS panels, and dense SNP markers sets appeared to provide slightly higher linkage scores and slightly higher power to detect linkage than MS markers. 2) Dense SNP panels also gave higher type I errors, suggesting that increased test thresholds may be needed to maintain the correct error rate. 3) Dense SNP panels provided better trait localization, but only in the COGA data, in which the MS markers were relatively loosely spaced. 4) The strength of linkage signals did not vary with the density of SNP panels, once the marker density was approximately 1 SNP/cM. 5) Analyses with SNPs were computationally challenging, and identified areas where improvements in analysis tools will be necessary to make analysis practical for widespread use.     

Data mining and computationally intensive methods: summary of Group 7 contributions to Genetic Analysis Workshop 13

The Framingham Heart Study data, as well as a related simulated data set, were generously provided to the participants of the Genetic Analysis Workshop 13 in order that newly developed and emerging statistical methodologies could be tested on that well-characterized data set. The impetus driving the development of novel methods is to elucidate the contributions of genes, environment, and interactions between and among them, as well as to allow comparison between and validation of methods. The seven papers that comprise this group used data-mining methodologies (tree-based methods, neural networks, discriminant analysis, and Bayesian variable selection) in an attempt to identify the underlying genetics of cardiovascular disease and related traits in the presence of environmental and genetic covariates. Data-mining strategies are gaining popularity because they are extremely flexible and may have greater efficiency and potential in identifying the factors involved in complex disorders. While the methods grouped together here constitute a diverse collection, some papers asked similar questions with very different methods, while others used the same underlying methodology to ask very different questions. This paper briefly describes the data-mining methodologies applied to the Genetic Analysis Workshop 13 data sets and the results of those investigations.     

Summary of Genetic Analysis Workshop 15: Group 9 linkage analysis of the CEPH expression data

Group 9 participants carried out linkage analysis of the Centre d'Etude de Polymorphism Humain (CEPH) expression data, using strategies that ranged from focused investigation of a small number of traits to full genome scans of all available traits. Results from five key areas encompass the most important results within and across the 17 participating groups. First, both extensive genetic heterogeneity and poor predictability of mapping results based on heritability have key implications for study design. Second, choice of the map used for linkage analysis is influential, with the implication that meiotic maps are preferable to physical maps. Third, performance of different analytic methods was in general fairly consistent, with the exception of one variance-component method that uses marker allele sharing as the dependent rather than independent variable. Fourth, multivariate analysis approaches did not generally appear to provide advantages over univariate approaches for linkage detection. Finally, there were computational and analytic challenges in working with a large public data set, along with need for more data documentation.     

Approaches to detecting gene x gene interaction in Genetic Analysis Workshop 14 pedigrees

Whether driven by the general lack of success in finding single-gene contributions to complex disease, by increased knowledge about the potential involvement of specific biological interactions in complex disease, or by recent dramatic increases in computational power, a large number of approaches to detect locus x locus interactions were recently proposed and implemented. The six Genetic Analysis Workshop 14 (GAW14) papers summarized here each applied either existing or refined approaches with the goal of detecting gene x gene, or locus x locus, interactions in the GAW14 data. Five of six papers analyzed the simulated data; the other analyzed the Collaborative Study on the Genetics of Alcoholism data. The analytic strategies implemented for detecting interactions included multifactor dimensionality reduction, conditional linkage analysis, nonparametric linkage correlation, two-locus parametric linkage analysis, and a joint test of linkage and association. Overall, most of the groups found limited success in consistently detecting all of the simulated interactions due, in large part, to the nature of the generating model.     

Data mining of RNA expression and DNA genotype data: presentation group 5 contributions to Genetic Analysis Workshop 15

The complexity of data available in human genetics continues to grow at an explosive rate. With that growth, the challenges to understanding the meaning of the underlying information also grow. A currently popular approach to dissecting such information falls under the broad category of data mining. This can apply to any approach that tries to extract relevant information from large amounts of data, but often refers to methods that deal, in a non-linear fashion, with very large numbers of variables that cannot be simultaneously handled by more conventional statistical methods. To explore the usefulness of some of these approaches, 13 groups applied a variety of strategies to the first dataset provided to GAW 15 participants. With the extensive microarray and SNP data provided for 14 CEPH families, these groups explored multistage analyses, machine learning methods, network construction, and other techniques to try to answer questions about gene-gene interaction, functional similarities, co-regulated gene expression and the mapping of gene expression determinants, among others. In general, the methods offered strategies to provide a better understanding of the complex pathways involved in gene expression and function. These are still "works in progress," often exploratory in nature, but they provide insights into ways in which the data might be interpreted. Despite the still preliminary nature of some of these methods and the diversity of the approaches, some common themes emerged. The collection of papers and methods offer a starting point for further exploration of complex interactions in human genetic data now readily available.     

Genetic Analysis Workshop II: segregation and linkage analysis

Familial segregation and linkage analyses were performed on two sets of the Genetic Analyses Workshop II data. The salient features of the mode of inheritance of the disease trait and its linkage/association with polymorphic markers and also marker-marker linkages were delineated using statistical-genetic techniques.     

Multiple testing in the genomics era: findings from Genetic Analysis Workshop 15, Group 15

Recent advances in molecular technologies have resulted in the ability to screen hundreds of thousands of single nucleotide polymorphisms and tens of thousands of gene expression profiles. While these data have the potential to inform investigations into disease etiologies and advance medicine, the question of how to adequately control both type I and type II error rates remains. Genetic Analysis Workshop 15 datasets provided a unique opportunity for participants to evaluate multiple testing strategies applicable to microarray and single nucleotide polymorphism data. The Genetic Analysis Workshop 15 multiple testing and false discovery rate group (Group 15) investigated three general categories for multiple testing corrections, which are summarized in this review: statistical independence, error rate adjustment, and data reduction. We show that while each approach may have certain advantages, adequate error control is largely dependent upon the question under consideration and often requires the use of multiple analytic strategies.     

Genetic Analysis Workshop II: segregation and three-locus linkage analysis

Data stimulated for Genetic Analysis Workshop II were analyzed using PAP. Segregation analysis showed a simple recessive mode of inheritance for data set 2 while no conclusions could be made about the mode of inheritance for data set 3. Pairwise linkage analysis suggested three linkage groups, but three-locus analysis did not provide strong evidence for the gene order within these groups. For three of the four three-locus comparisons performed, three-locus analysis suggested the simulated order. In only one case did the pairwise analysis suggest the simulated order, indicating the necessity for multi-locus analysis for gene order.     

Genetic Analysis Workshop II: multiple-locus segregation analysis incorporating linkage markers

The workshop data were examined using a newly developed methodology (MILINK, Risch, 1984) for combined segregation, linkage, and association analysis of a complex disease trait in pedigree data. Results from problems two and three suggest that the method is powerful both for determining mode of disease inheritance and for resolution of linkage disequilibrium versus pleiotrophy (with epistasis) of marker alleles.     

Genetic Analysis Workshop II: sib pair screening tests for linkage

For each marker locus and for every pair of sibs with data available in the 1983 workshop data, the proportion of genes identical by descent was estimated. The mean proportions were compared between concordant and discordant sib pairs, and the mean proportion for concordantly affected pairs was compared with one half. Together with standard tests of association, these found to be sensitive screening tests for linkage.     

Genetic Analysis Workshop II: results of incorporating a linkage disequilibrium parameter

The corporation of a linkage disequilibrium parameter, delta, into linkage analysis is illustrated for data from Genetic Analysis Workshop II. Points from a joint likelihood surface are calculated and displayed on a recombination fraction-linkage disequilibrium grid using a simple modification of LIPED. The approach is shown to increase the power of linkage analysis and the power of tests for heterogeneity of linkage for the simulated examples.     

Genetic Analysis Workshop II: combined segregation, linkage, and association analysis

A combined segregation, linkage, and association analysis using the program COMBIN was performed on the simulated pedigree data prepared for the Second Genetic Analysis Workshop. The model used in COMBIN is described and the presented results illustrate its effectiveness in the analysis of such data. Linkage analysis was performed and maps for each linkage group are presented.