Progress and promise of genome-wide association studies for human complex trait genetics

Enormous progress in mapping complex traits in humans has been made in the last 5 yr. There has been early success for prevalent diseases with complex phenotypes. These studies have demonstrated clearly that, while complex traits differ in their underlying genetic architectures, for many common disorders the predominant pattern is that of many loci, individually with small effects on phenotype. For some traits, loci of large effect have been identified. For almost all complex traits studied in humans, the sum of the identified genetic effects comprises only a portion, generally less than half, of the estimated trait heritability. A variety of hypotheses have been proposed to explain why this might be the case, including untested rare variants, and gene-gene and gene-environment interaction. Effort is currently being directed toward implementation of novel analytic approaches and testing rare variants for association with complex traits using imputed variants from the publicly available 1000 Genomes Project resequencing data and from direct resequencing of clinical samples. Through integration with annotations and functional genomic data as well as by in vitro and in vivo experimentation, mapping studies continue to characterize functional variants associated with complex traits and address fundamental issues such as epistasis and pleiotropy. This review focuses primarily on the ways in which genome-wide association studies (GWASs) have revolutionized the field of human quantitative genetics.     

Genephony: a knowledge management tool for genome-wide research

Background  

One of the consequences of the rapid and widespread adoption of high-throughput experimental technologies is an exponential increase of the amount of data produced by genome-wide experiments. Researchers increasingly need to handle very large volumes of heterogeneous data, including both the data generated by their own experiments and the data retrieved from publicly available repositories of genomic knowledge. Integration, exploration, manipulation and interpretation of data and information therefore need to become as automated as possible, since their scale and breadth are, in general, beyond the limits of what individual researchers and the basic data management tools in normal use can handle. This paper describes Genephony, a tool we are developing to address these challenges.     

Bayesian model selection for genome-wide epistatic quantitative trait loci analysis

The problem of identifying complex epistatic quantitative trait loci (QTL) across the entire genome continues to be a formidable challenge for geneticists. The complexity of genome-wide epistatic analysis results mainly from the number of QTL being unknown and the number of possible epistatic effects being huge. In this article, we use a composite model space approach to develop a Bayesian model selection framework for identifying epistatic QTL for complex traits in experimental crosses from two inbred lines. By placing a liberal constraint on the upper bound of the number of detectable QTL we restrict attention to models of fixed dimension, greatly simplifying calculations. Indicators specify which main and epistatic effects of putative QTL are included. We detail how to use prior knowledge to bound the number of detectable QTL and to specify prior distributions for indicators of genetic effects. We develop a computationally efficient Markov chain Monte Carlo (MCMC) algorithm using the Gibbs sampler and Metropolis-Hastings algorithm to explore the posterior distribution. We illustrate the proposed method by detecting new epistatic QTL for obesity in a backcross of CAST/Ei mice onto M16i.     

A rapid generalized least squares model for a genome-wide quantitative trait association analysis in families

Genome-wide association studies (GWAS) using family data involve association analyses between hundreds of thousands of markers and a trait for a large number of related individuals. The correlations among relatives bring statistical and computational challenges when performing these large-scale association analyses. Recently, several rapid methods accounting for both within- and between-family variation have been proposed. However, these techniques mostly model the phenotypic similarities in terms of genetic relatedness. The familial resemblances in many family-based studies such as twin studies are not only due to the genetic relatedness, but also derive from shared environmental effects and assortative mating. In this paper, we propose 2 generalized least squares (GLS) models for rapid association analysis of family-based GWAS, which accommodate both genetic and environmental contributions to familial resemblance. In our first model, we estimated the joint genetic and environmental variations. In our second model, we estimated the genetic and environmental components separately. Through simulation studies, we demonstrated that our proposed approaches are more powerful and computationally efficient than a number of existing methods are. We show that estimating the residual variance-covariance matrix in the GLS models without SNP effects does not lead to an appreciable bias in the p values as long as the SNP effect is small (i.e. accounting for no more than 1% of trait variance).     

Genomic selection and complex trait prediction using a fast EM algorithm applied to genome-wide markers

Background  

The information provided by dense genome-wide markers using high throughput technology is of considerable potential in human disease studies and livestock breeding programs. Genome-wide association studies relate individual single nucleotide polymorphisms (SNP) from dense SNP panels to individual measurements of complex traits, with the underlying assumption being that any association is caused by linkage disequilibrium (LD) between SNP and quantitative trait loci (QTL) affecting the trait. Often SNP are in genomic regions of no trait variation. Whole genome Bayesian models are an effective way of incorporating this and other important prior information into modelling. However a full Bayesian analysis is often not feasible due to the large computational time involved.     

The mouse as a model for human biology: a resource guide for complex trait analysis

The mouse has been a powerful force in elucidating the genetic basis of human physiology and pathophysiology. From its beginnings as the model organism for cancer research and transplantation biology to the present, when dissection of the genetic basis of complex disease is at the forefront of genomics research, an enormous and remarkable mouse resource infrastructure has accumulated. This review summarizes those resources and provides practical guidelines for their use, particularly in the analysis of quantitative traits.     

Segregation analysis of a complex quantitative trait: approaches for identifying influential data points

BACKGROUND/AIMS: Complex traits pose a particular challenge to standard methods for segregation analysis (SA), and for such traits it is difficult to assess the ability of complex SA (CSA) to approximate the true mode of inheritance. Here we use an oligogenic Bayesian Markov chain Monte Carlo method for SA (OSA) to verify results from a single-locus likelihood-based CSA for data on a quantitative measure of reading ability. METHODS: We compared the profile likelihood from CSA, maximized over the trait allele frequency, to the posterior distribution of genotype effects from OSA to explore differences in the overall parameter estimates from SA on the original phenotype data and the same data Winsorized to reduce the potential influence of three outlying data points. RESULTS: Bayesian OSA revealed two modes of inheritance, one of which coincided with the QTL model from CSA. Winsorizing abolished the model originally estimated by CSA; both CSA and OSA identified only the second OSA model. CONCLUSION: Differences between the results from the two methods alerted us to the presence of influential data points, and identified the QTL model best supported by the data. Thus, the Bayesian OSA proved a valuable tool for assessing and verifying inheritance models from CSA.     

Use of multivariate linkage analysis for dissection of a complex cognitive trait

Replication of linkage results for complex traits has been exceedingly difficult, owing in part to the inability to measure the precise underlying phenotype, small sample sizes, genetic heterogeneity, and statistical methods employed in analysis. Often, in any particular study, multiple correlated traits have been collected, yet these have been analyzed independently or, at most, in bivariate analyses. Theoretical arguments suggest that full multivariate analysis of all available traits should offer more power to detect linkage; however, this has not yet been evaluated on a genomewide scale. Here, we conduct multivariate genomewide analyses of quantitative-trait loci that influence reading- and language-related measures in families affected with developmental dyslexia. The results of these analyses are substantially clearer than those of previous univariate analyses of the same data set, helping to resolve a number of key issues. These outcomes highlight the relevance of multivariate analysis for complex disorders for dissection of linkage results in correlated traits. The approach employed here may aid positional cloning of susceptibility genes in a wide spectrum of complex traits.     

myMIR: a genome-wide microRNA targets identification and annotation tool

miRNA target genes prediction represents a crucial step in miRNAs functional characterization. In this context, the challenging issue remains predictions accuracy and recognition of false positive results. In this article myMIR, a web based system for increasing reliability of miRNAs predicted targets lists, is presented. myMIR implements an integrated pipeline for computing ranked miRNA::target lists and provides annotations for narrowing them down. The system relies on knowledge base data, suitably integrated in order to extend the functional characterization of targeted genes to miRNAs, by highlighting the search on over-represented annotation terms. Validation results show a dramatic reduction in the quantity of predictions and an increase in the sensitivity, when compared to other methods. This improves the predictions accuracy and allows the formulation of novel hypotheses on miRNAs functional involvement.     

ReplicationDomain: a visualization tool and comparative database for genome-wide replication timing data

Background  

Eukaryotic DNA replication is regulated at the level of large chromosomal domains (0.5–5 megabases in mammals) within which replicons are activated relatively synchronously. These domains replicate in a specific temporal order during S-phase and our genome-wide analyses of replication timing have demonstrated that this temporal order of domain replication is a stable property of specific cell types.     

Sex-specific regulation of mitochondrial DNA levels: genome-wide linkage analysis to identify quantitative trait Loci

Altered mitochondrial DNA (mtDNA) levels have been associated with common diseases in humans. We investigated the genetic mechanism that controls mtDNA levels using genome-wide linkage analyses in families from the Genetic Analysis of Idiopathic Thrombophilia Project (GAIT). We measure mtDNA levels by quantitative real-time PCR in 386 subjects from 21 extended Spanish families. A variance component linkage method using 485 microsatellites was conducted to evaluate linkage and to detect quantitative trait loci (QTLs) involved in the control of mtDNA levels. The heritalibility of mtDNA levels was 0.33 (p = 1.82e-05). We identified a QTL on Chromosome 2 (LOD = 2.21) using all of the subjects, independently on their sex. When females and males were analysed separately, three QTLs were identified. Females showed the same QTL on Chromosome 2 (LOD = 3.09), indicating that the QTL identified in the analysis using all of the subjects was a strong female QTL, and another one on Chromosome 3 (LOD = 2.67), whereas in males a QTL was identified on Chromosome 1 (LOD = 2.81). These QTLs were fine-mapped to find associations with mtDNA levels. The most significant SNP association was for the rs10888838 on Chromosome 1 in males. This SNP mapped to the gene MRPL37, involved in mitochondrial protein translation. The rs2140855 on Chromosome 2 showed association in the analysis using all of the subjects. It was near the gene CMPK2, which encodes a mitochondrial enzyme of the salvage pathway of deoxyribonucleotide synthesis. Our results provide evidence of a sex-specific genetic mechanism for the control of mtDNA levels and provide a framework to identify new genes that influence mtDNA levels.     

ACE-it: a tool for genome-wide integration of gene dosage and RNA expression data

SUMMARY: We describe a tool, called ACE-it (Array CGH Expression integration tool). ACE-it links the chromosomal position of the gene dosage measured by array CGH to the genes measured by the expression array. ACE-it uses this link to statistically test whether gene dosage affects RNA expression. AVAILABILITY: ACE-it is freely available at http://ibivu.cs.vu.nl/programs/acewww/.     

Quantitative trait locus analysis of longitudinal quantitative trait data in complex pedigrees

There is currently considerable interest in genetic analysis of quantitative traits such as blood pressure and body mass index. Despite the fact that these traits change throughout life they are commonly analyzed only at a single time point. The genetic basis of such traits can be better understood by collecting and effectively analyzing longitudinal data. Analyses of these data are complicated by the need to incorporate information from complex pedigree structures and genetic markers. We propose conducting longitudinal quantitative trait locus (QTL) analyses on such data sets by using a flexible random regression estimation technique. The relationship between genetic effects at different ages is efficiently modeled using covariance functions (CFs). Using simulated data we show that the change in genetic effects over time can be well characterized using CFs and that including parameters to model the change in effect with age can provide substantial increases in power to detect QTL compared with repeated measure or univariate techniques. The asymptotic distributions of the methods used are investigated and methods for overcoming the practical difficulties in fitting CFs are discussed. The CF-based techniques should allow efficient multivariate analyses of many data sets in human and natural population genetics.     

Immunogenetic Management Software: a new tool for visualization and analysis of complex immunogenetic datasets

Here we describe the Immunogenetic Management Software (IMS) system, a novel web-based application that permits multiplexed analysis of complex immunogenetic traits that are necessary for the accurate planning and execution of experiments involving large animal models, including nonhuman primates. IMS is capable of housing complex pedigree relationships, microsatellite-based MHC typing data, as well as MHC pyrosequencing expression analysis of class I alleles. It includes a novel, automated MHC haplotype naming algorithm and has accomplished an innovative visualization protocol that allows users to view multiple familial and MHC haplotype relationships through a single, interactive graphical interface. Detailed DNA and RNA-based data can also be queried and analyzed in a highly accessible fashion, and flexible search capabilities allow experimental choices to be made based on multiple, individualized and expandable immunogenetic factors. This web application is implemented in Java, MySQL, Tomcat, and Apache, with supported browsers including Internet Explorer and Firefox on Windows and Safari on Mac OS. The software is freely available for distribution to noncommercial users by contacting Leslie.kean@emory.edu. A demonstration site for the software is available at http://typing.emory.edu/typing_demo , user name: imsdemo7@gmail.com and password: imsdemo.