WHAP: haplotype-based association analysis

We describe a software tool to perform haplotype-based association analysis, for quantitative and qualitative traits, in population and family samples, using single nucleotide polymorphism or multiallelic marker data. A range of tests is offered: omnibus and haplotype-specific tests; prospective and retrospective likelihoods; covariates and moderators; sliding window analyses; permutation P-values. We focus on the ability to flexibly impose constraints on haplotype effects, which allows for a range of conditional haplotype-based likelihood ratio tests: for example, whether an allele has an effect independent of its haplotypic background, or whether a single variant can explain the overall association at a locus. We illustrate using these tests to dissect a multi-locus association. AVAILABILITY: WHAP is a C/C++ program, freely available from the author's website: http://pngu.mgh.harvard.edu/purcell/whap/     

EFBAT: exact family-based association tests

Background

Previous studies have reported frequent stretches of homozygosity in human subjects but have failed to clarify whether these are due to cytogenetic abnormalities or to autozygosity.

Methods

Trios which had been typed for closely spaced SNPs spanning the genome were studied. Stretches of extended homozygosity were identified in the child members, as were occasions on which the child had been genotyped as not inheriting one parental allele. The number of times such transmission errors occurred within regions of extended homozygosity was compared with the chance expectation.

Results

Transmission errors occurred more rarely in regions of extended homozygosity than would be expected by chance.

Discussion

Regions of extended homozygosity are not generally due to cytogenetic abnormalities such as uniparental isodisomy. They reflect the Mendelian inheritance of haplotypes from a common ancestor. This may have implications for mapping disease genes.     

Population genomics: ageing by association

Mutations in the gene Klotho lead to premature ageing in mice. Recent work on human genetic variation has identified an association between a particular allele of Klotho and human lifespan. A harbinger of things to come, this work illustrates the power and growth potential for association studies of human ageing.     

A Bayesian multilocus association method: allowing for higher-order interaction in association studies

For most common diseases with heritable components, not a single or a few single-nucleotide polymorphisms (SNPs) explain most of the variance for these disorders. Instead, much of the variance may be caused by interactions (epistasis) among multiple SNPs or interactions with environmental conditions. We present a new powerful statistical model for analyzing and interpreting genomic data that influence multifactorial phenotypic traits with a complex and likely polygenic inheritance. The new method is based on Markov chain Monte Carlo (MCMC) and allows for identification of sets of SNPs and environmental factors that when combined increase disease risk or change the distribution of a quantitative trait. Using simulations, we show that the MCMC method can detect disease association when multiple, interacting SNPs are present in the data. When applying the method on real large-scale data from a Danish population-based cohort, multiple interactions are identified that severely affect serum triglyceride levels in the study individuals. The method is designed for quantitative traits but can also be applied on qualitative traits. It is computationally feasible even for a large number of possible interactions and differs fundamentally from most previous approaches by entertaining nonlinear interactions and by directly addressing the multiple-testing problem.     

Protein-protein association: investigation of factors influencing association rates by brownian dynamics simulations

The rate of protein-protein association limits the response time due to protein-protein interactions. The bimolecular association rate may be diffusion-controlled or influenced, and in such cases, Brownian dynamics simulations of protein-protein diffusional association may be used to compute association rates. Here, we report Brownian dynamics simulations of the diffusional association of five different protein-protein pairs: barnase and barstar, acetylcholinesterase and fasciculin-2, cytochrome c peroxidase and cytochrome c, the HyHEL-5 antibody and hen egg lysozyme (HEL), and the HyHEL-10 antibody and HEL. The same protocol was used to compute the diffusional association rates for all the protein pairs in order to assess, by comparison to experimentally measured rates, whether the association of these proteins can be explained solely on the basis of diffusional encounter. The simulation protocol is similar to those previously derived for simulation of the association of barnase and barstar, and of acetylcholinesterase and fasciculin-2; these produced results in excellent agreement with experimental data for these protein pairs, with changes in association rate due to mutations reproduced within the limits of expected computational and modeling errors. Here, we find that for all protein pairs, the effects of mutations can be well reproduced by the simulations, even though the degree of the electrostatic translational and orientational steering varies widely between the cases. However, the absolute values of association rates for the acetylcholinesterase: fasciculin-2 and HyHEL-10 antibody: HEL pairs are overestimated. Comparison of bound and unbound protein structures shows that this may be due to gating resulting from protein flexibility in some of the proteins. This may lower the association rates compared to their bimolecular diffusional encounter rates.     

The multispecies formicary association: the transition from an incomplete to a complete association

This work describes the transformation of the nest complex of a facultatively dominant species, Formica cunicularia, which was naturally colonized by the families of obligate dominants, Formica pratensis and Lasius fuliginosus. The complex was under observation from 2003 to 2013. During this period, F. pratensis families almost completely replaced F. cunicularia in the tenth part of the territory of this multispecies association. The family L. fuliginosus did not exert such an effect on the settlement of the facultative dominant. The study addresses the interaction of foragers of the species that found themselves in the zone of F. pratensis invasion, the dynamics of expansion of the dominant’s foraging areas, and the changes in the multispecies formicary association on the whole. It is shown that, even with different peaceful coexistence mechanisms at work, the alteration of the association occurs through attacks of the obligate dominant on the neighboring families. The progress of this process and possible key situations are considered.     

Ectomychorrizal DB: a symbiotic association database

Ectomycorrhizal (ECM) fungal species, a "Symbiotic" relationship between tress and fungi in forest has a great ecological and economic importance. Here is an attempt to describe database named "EctomycorrhizalDB", addressing ECM diversity of Central Himalaya (Kumaun region), with special emphasis on their characterization, physical properties and morphological features along with specifications. This database would help the scientific community to draw a better understanding of the environmental factors that affects species diversity. AVAILABILITY: The database is available for free at http://www.kubic.nic.in/ectomychorhiza.     

AssociationDB: web-based exploration of genomic association

Genome-wide association studies use hundreds of thousands of markers making it challenging to present and finally interpret the results. We developed a graphical, web-based solution for an interactive exploration of the results of case-control studies, with a tight integration of related gene information and tissue-specific expression data. Association results are presented as physical position-based vertical bars with known genes included as horizontal bars at their respective physical positions. The interface allows the specification of filtering criteria for the association data and highlights potentially interesting genes with user-specified terms occurring in their reports or with relevant expression patterns. Pop-up windows and hyperlinks provide drill-down capabilities and quick access to relevant data AssociationDB can either be used as a stand-alone solution or as a front-end joining association results obtained by other software with genomic information.     

GWAMA: software for genome-wide association meta-analysis

Background  

Despite the recent success of genome-wide association studies in identifying novel loci contributing effects to complex human traits, such as type 2 diabetes and obesity, much of the genetic component of variation in these phenotypes remains unexplained. One way to improving power to detect further novel loci is through meta-analysis of studies from the same population, increasing the sample size over any individual study. Although statistical software analysis packages incorporate routines for meta-analysis, they are ill equipped to meet the challenges of the scale and complexity of data generated in genome-wide association studies.     

Type VIII collagen: heterotrimeric chain association

Two chains, alpha1(VIII) and alpha2(VIII), have been described for type VIII collagen. Early work suggested that these chains were present in a 2:1 ratio, although recent work has shown that homotrimers can form and predominate in some tissues. In order to address the question of whether the alpha1(VIII) and alpha2(VIII) chains could co-polymerise we made a shortened alpha1(VIII) chain and expressed this with full length alpha2(VIII) chain in an in vitro translation system supplemented with semi-permeabilised cells. Heterotrimers containing either two or one alpha2(VIII) were evident. Interestingly, a point mutation in the NC1 domain of the alpha1(VIII) chain abrogated trimer formation. In addition we were able to demonstrate chain association of the alpha1(X) chain of type X collagen with the shortened alpha1(VIII) chain. Variations in chain association were seen when altered ratios of message were used. These results demonstrate the importance of the NC1 domain in chain association and suggest that gene expression regulates the composition and function of type VIII collagen by varying chain composition.     

Human chromosomology: Random association of acrocentrics

Lymphocytes from one normal male and one normal female were labeled with tritiated thymidine to study the autoradiographic patterns of association of acrocentrics, i.e., D and G groups. Association into groups of two and three acrocentrics was found to be random. These results are supported by a random association found in other cell lines studied without autoradiography. These include: normal persons, leukemic cell lines and several other instances where rearrangements of D and G chromosomes have been found.This paper is dedicated with gratitude to my teacher Dr. MiguelJiménez S. This study was supported by funds from the Universidad Nacional Autónoma de México and the Secretaria de Salubridad y Asistencia.     

Adiponectin concentrations: a genome-wide association study

Adiponectin is associated with obesity and insulin resistance. To date, there has been no genome-wide association study (GWAS) of adiponectin levels in Asians. Here we present a GWAS of a cohort of Korean volunteers. A total of 4,001 subjects were genotyped by using a genome-wide marker panel in a two-stage design (979 subjects initially and 3,022 in a second stage). Another 2,304 subjects were used for follow-up replication studies with selected markers. In the discovery phase, the top SNP associated with mean log adiponectin was rs3865188 in CDH13 on chromosome 16 (p = 1.69 × 10(-15) in the initial sample, p = 6.58 × 10(-39) in the second genome-wide sample, and p = 2.12 × 10(-32) in the replication sample). The meta-analysis p value for rs3865188 in all 6,305 individuals was 2.82 × 10(-83). The association of rs3865188 with high-molecular-weight adiponectin (p = 7.36 × 10(-58)) was even stronger in the third sample. A reporter assay that evaluated the effects of a CDH13 promoter SNP in complete linkage disequilibrium with rs3865188 revealed that the major allele increased expression 2.2-fold. This study clearly shows that genetic variants in CDH13 influence adiponectin levels in Korean adults.     

'Smoking genes': a genetic association study

Some controversy exists on the specific genetic variants that are associated with nicotine dependence and smoking-related phenotypes. The purpose of this study was to analyse the association of smoking status and smoking-related phenotypes (included nicotine dependence) with 17 candidate genetic variants: CYP2A6*1×2, CYP2A6*2 (1799T>A) [rs1801272], CYP2A6*9 (-48T>G) [rs28399433], CYP2A6*12, CYP2A13*2 (3375C>T) [rs8192789], CYP2A13*3 (7520C>G), CYP2A13*4 (579G>A), CYP2A13*7 (578C>T) [rs72552266], CYP2B6*4 (785A>G), CYP2B6*9 (516G>T), CHRNA3 546C>T [rs578776], CHRNA5 1192G>A [rs16969968], CNR1 3764C>G [rs6928499], DRD2-ANKK1 2137G>A (Taq1A) [rs1800497], 5HTT LPR, HTR2A -1438A>G [rs6311] and OPRM1 118A>G [rs1799971]. We studied the genotypes of the aforementioned polymorphisms in a cohort of Spanish smokers (cases, N = 126) and ethnically matched never smokers (controls, N = 80). The results showed significant between-group differences for CYP2A6*2 and CYP2A6*12 (both P<0.001). Compared with carriers of variant alleles, the odds ratio (OR) for being a non-smoker in individuals with the wild-type genotype of CYP2A6*12 and DRD2-ANKK1 2137G>A (Taq1A) polymorphisms was 3.60 (95%CI: 1.75, 7.44) and 2.63 (95%CI: 1.41, 4.89) respectively. Compared with the wild-type genotype, the OR for being a non-smoker in carriers of the minor CYP2A6*2 allele was 1.80 (95%CI: 1.24, 2.65). We found a significant genotype effect (all P≤0.017) for the following smoking-related phenotypes: (i) cigarettes smoked per day and CYP2A13*3; (ii) pack years smoked and CYP2A6*2, CYP2A6*1×2, CYP2A13*7, CYP2B6*4 and DRD2-ANKK1 2137G>A (Taq1A); (iii) nicotine dependence (assessed with the Fagestrom test) and CYP2A6*9. Overall, our results suggest that genetic variants potentially involved in nicotine metabolization (mainly, CYP2A6 polymorphisms) are those showing the strongest association with smoking-related phenotypes, as opposed to genetic variants influencing the brain effects of nicotine, e.g., through nicotinic acetylcholine (CHRNA5), serotoninergic (HTR2A), opioid (OPRM1) or cannabinoid receptors (CNR1).     

Issues in association analysis: error control in case-control association studies for disease gene discovery

Several sources of errors are discussed. While genotyping errors have little effect on power in case-control association studies, they tend to strongly increase false positive results in TDT type tests unless occurrence of errors is allowed for in the analysis (e.g., TDTae test). Disregarding non-genetic risk factors is shown to lead to a form of hidden heterogeneity, which can strongly reduce power. Stratification of data into more homogeneous subgroups is advocated as a simple solution to allowing for non-genetic risk factors such as socio-economic status and food preferences.     

Family-based association tests for genomewide association scans

With millions of single-nucleotide polymorphisms (SNPs) identified and characterized, genomewide association studies have begun to identify susceptibility genes for complex traits and diseases. These studies involve the characterization and analysis of very-high-resolution SNP genotype data for hundreds or thousands of individuals. We describe a computationally efficient approach to testing association between SNPs and quantitative phenotypes, which can be applied to whole-genome association scans. In addition to observed genotypes, our approach allows estimation of missing genotypes, resulting in substantial increases in power when genotyping resources are limited. We estimate missing genotypes probabilistically using the Lander-Green or Elston-Stewart algorithms and combine high-resolution SNP genotypes for a subset of individuals in each pedigree with sparser marker data for the remaining individuals. We show that power is increased whenever phenotype information for ungenotyped individuals is included in analyses and that high-density genotyping of just three carefully selected individuals in a nuclear family can recover >90% of the information available if every individual were genotyped, for a fraction of the cost and experimental effort. To aid in study design, we evaluate the power of strategies that genotype different subsets of individuals in each pedigree and make recommendations about which individuals should be genotyped at a high density. To illustrate our method, we performed genomewide association analysis for 27 gene-expression phenotypes in 3-generation families (Centre d'Etude du Polymorphisme Humain pedigrees), in which genotypes for ~860,000 SNPs in 90 grandparents and parents are complemented by genotypes for ~6,700 SNPs in a total of 168 individuals. In addition to increasing the evidence of association at 15 previously identified cis-acting associated alleles, our genotype-inference algorithm allowed us to identify associated alleles at 4 cis-acting loci that were missed when analysis was restricted to individuals with the high-density SNP data. Our genotype-inference algorithm and the proposed association tests are implemented in software that is available for free.