Novel method to estimate the phenotypic variation explained by genome-wide association studies reveals large fraction of the missing heritability

MDVarP: modifier ~ disease-causing variant pairs predictor

BACKGROUND

Modifiers significantly impact disease phenotypes by modulating the effects of disease-causing variants, resulting in varying disease manifestations among individuals. However, identifying genetic interactions between modifier and disease-causing variants is challenging.

RESULTS

We developed MDVarP, an ensemble model comprising 1000 random forest predictors, to identify modifier ~ disease-causing variant combinations. MDVarP achieves high accuracy and precision, as verified using an independent dataset with published evidence of genetic interactions. We identified 25 novel modifier ~ disease-causing variant combinations and obtained supporting evidence for these associations. MDVarP outputs a class label ("Associated-pair" or "Nonrelevant-pair") and two prediction scores indicating the probability of a true association.

CONCLUSIONS

MDVarP prioritizes variant pairs associated with phenotypic modulations, enabling more effective mapping of functional contributions from disease-causing and modifier variants. This framework interprets genetic interactions underlying phenotypic variations in human diseases, with potential applications in personalized medicine and disease prevention.

Neurobiology of numerical learning

Numerical abilities are complex cognitive skills essential for dealing with requirements of the modern world. Although the brain structures and functions underlying numerical cognition in different species have long been appreciated, genetic and molecular techniques have more recently expanded the knowledge about the mechanisms underlying numerical learning. In this review, we discuss the status of the research related to the neurobiological bases of numerical abilities. We consider how genetic factors have been associated with mathematical capacities and how these link to the current knowledge of brain regions underlying these capacities in human and non-human animals. We further discuss the extent to which significant variations in the levels of specific neurotransmitters may be used as potential markers of individual performance and learning difficulties and take into consideration the therapeutic potential of brain stimulation methods to modulate learning and improve interventional outcomes. The implications of this research for formulating a more comprehensive view of the neural basis of mathematical learning are discussed.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Genome-wide association study revealed that the TaGW8 gene was associated with kernel size in Chinese bread wheat

Using Wheat 90 K SNP assay, kernel-related traits of Chinese bread wheat were used to perform association mapping in 14 environments by GWAS. Results indicated that 996 and 953 of 4417 and 3172 significant SNPs for kernel length and thousand-kernel weight were located on the chromosome 7B. Haplotype analysis of these SNPs on 7B generated the block containing the predicted TaGW8-B1 gene. TaGW8-B1 gene was further cloned by sequencing in bread wheat and a 276-bp InDel was found in the first intron. TaGW8-B1 without and with the 276-bp InDel were designated as TaGW8-B1a and TaGW8-B1b, respectively. Analysis of agronomic traits indicated that cultivars with TaGW8-B1a possessed significantly wider kernel width, significantly more kernel number per spike, longer kernel length, higher thousand-kernel weight and more spikelet number per spike than cultivars with TaGW8-B1b. Furthermore, cultivars with TaGW8-B1a possessed significantly higher yield than cultivars with TaGW8-B1b. Therefore, TaGW8-B1a was considered as a potentially superior allele. Meanwhile, TaGW8-B1a possessed a significantly higher expression level than TaGW8-B1b in mature seeds by qRT-PCR. It possibly suggested that the high expression of TaGW8-B1 was positively associated with kernel size in bread wheat. Distribution of TaGW8-B1 allele indicated that TaGW8-B1a has been positively selected in Chinese wheat.

The genetics of adiposity

Genome-wide discovery efforts have identified more than 500 genetic loci associated with adiposity traits. The vast majority of these loci were found through large-scale meta-analyses for body mass index (BMI) and waist-to-hip ratio (WHR), and in European ancestry populations. However, alternative approaches, focusing on non-European ancestry populations, more refined adiposity measures, and low-frequency (minor allele frequency (MAF)<5%) coding variants, identified additional novel loci that had not been identified before. Loci associated with overall obesity implicate pathways that act in the brain, whereas loci associated with fat distribution point to pathways involved in adipocyte biology. Pinpointing the causal gene within each locus remains challenging, but is a critical step towards translation of genome-wide association study (GWAS) loci into new biology. Ultimately, new genes may provide pharmacological targets for the development of weight loss drugs.

Towards precision medicine

There is great potential for genome sequencing to enhance patient care through improved diagnostic sensitivity and more precise therapeutic targeting. To maximize this potential, genomics strategies that have been developed for genetic discovery - including DNA-sequencing technologies and analysis algorithms - need to be adapted to fit clinical needs. This will require the optimization of alignment algorithms, attention to quality-coverage metrics, tailored solutions for paralogous or low-complexity areas of the genome, and the adoption of consensus standards for variant calling and interpretation. Global sharing of this more accurate genotypic and phenotypic data will accelerate the determination of causality for novel genes or variants. Thus, a deeper understanding of disease will be realized that will allow its targeting with much greater therapeutic precision.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

A new method for estimating effect size distribution and heritability from genome-wide association summary results

Accurately estimating the distribution and heritability of SNP effects across the genome could help explain the mystery of missing heritability. In this study, we propose a novel statistical method for estimating the distribution and heritability of SNP effects from genome-wide association studies (GWASs), and compare its performance to several existing methods using both simulations and real data. Specifically, we study the full range of GWAS summary results and link observed p values and unobserved effect sizes by (non-central) Chi-square distribution. By modeling the observed full set of association signals using a multinomial distribution, we build a likelihood function of SNP effect sizes using parametric and non-parametric maximum likelihood frameworks. Simulation studies show that the proposed method can accurately estimate effect sizes and the number of associated SNPs. As real applications, we analyze publicly available GWAS summary results for height, body mass index (BMI), and bone mineral density (BMD). Our analyses show that there are over 10,000 SNPs that might be associated with height, and the total heritability attributable to these SNPs exceeds 70 %. The heritabilities for BMI and BMD are ~10 and ~15 %, respectively. The results indicate that the proposed method has the potential to improve the accuracy of estimates of heritability and effect size for common SNPs in large-scale GWAS meta-analyses. These improved estimates may contribute to an enhanced understanding of the genetic basis of complex traits.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Heritability of liver enzyme levels estimated from genome-wide SNP data

Variation in the liver enzyme levels in humans is moderately heritable, as indicated by twin-family studies. At present, genome-wide association studies have traced <2% of the variance back to genome-wide significant single-nucleotide polymorphisms (SNPs). We estimated the SNP-based heritability of levels of three liver enzymes (gamma-glutamyl transferase (GGT); alanine aminotransferase (ALT); and aspartate aminotransferase (AST)) using genome-wide SNP data in a sample of 5421 unrelated Dutch individuals. Two estimation methods for SNP-based heritability were compared, one based on the distant genetic relatedness among all subjects as summarized in a Genetic Relatedness Matrix (GRM), and the other one based on density estimation (DE). The DE method was also applied to meta-analysis results on GGT and ALT. GRM-derived SNP-based heritability estimates were significant for GGT (16%) and AST (11%), but not for ALT (6%). DE estimates in the same sample varied as a function of pruning and were around 23% for all liver enzymes. Application of the DE approach to meta-analysis results for GGT and ALT gave SNP-based heritability estimates of 6 and 3%. The significant results in the Dutch sample indicate that genome-wide SNP platforms contain substantial information regarding the underlying genetic variation in the liver enzyme levels. A major part of this genetic variation remains however undetected. SNP-based heritability estimates, based on meta-analysis results, may point at substantial heterogeneity among cohorts contributing to the meta-analysis. This type of analysis may provide useful information to guide future gene searches.

An atlas of genetic influences on human blood metabolites

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

Genome-wide association study of lung function phenotypes in a founder population

BACKGROUND

Lung function is a long-term predictor of mortality and morbidity.

OBJECTIVE

We sought to identify single nucleotide polymorphisms (SNPs) associated with lung function.

METHODS

We performed a genome-wide association study (GWAS) of FEV1, forced vital capacity (FVC), and FEV1/FVC in 1144 Hutterites aged 6 to 89 years, who are members of a founder population of European descent. We performed least absolute shrinkage and selection operation regression to select the minimum set of SNPs that best predict FEV1/FVC in the Hutterites and used the GRAIL algorithm to mine the Gene Ontology database for evidence of functional connections between genes near the predictive SNPs.

RESULTS

Our GWAS identified significant associations between FEV1/FVC and SNPs at the THSD4-UACA-TLE3 locus on chromosome 15q23 (P = 5.7 × 10(-8) to 3.4 × 10(-9)). Nine SNPs at or near 4 additional loci had P < 10(-5) with FEV1/FVC. Only 2 SNPs were found with P < 10(-5) for FEV1 or FVC. We found nominal levels of significance with SNPs at 9 of the 27 previously reported loci associated with lung function measures. Among a predictive set of 80 SNPs, 6 loci were identified that had a significant degree of functional connectivity (GRAIL P < .05), including 3 clusters of β-defensin genes, 2 chemokine genes (CCL18 and CXCL12), and TNFRSF13B.

CONCLUSION

This study identifies genome-wide significant associations and replicates results of previous GWASs. Multimarker modeling implicated for the first time common variation in genes involved in antimicrobial immunity in airway mucosa that influences lung function.

Meta-analysis methods for genome-wide association studies and beyond

Meta-analysis of genome-wide association studies (GWASs) has become a popular method for discovering genetic risk variants. Here, we overview both widely applied and newer statistical methods for GWAS meta-analysis, including issues of interpretation and assessment of sources of heterogeneity. We also discuss extensions of these meta-analysis methods to complex data. Where possible, we provide guidelines for researchers who are planning to use these methods. Furthermore, we address special issues that may arise for meta-analysis of sequencing data and rare variants. Finally, we discuss challenges and solutions surrounding the goals of making meta-analysis data publicly available and building powerful consortia.

A genome-wide association study of brain lesion distribution in multiple sclerosis

Brain magnetic resonance imaging is widely used as a diagnostic and monitoring tool in multiple sclerosis and provides a non-invasive, sensitive and reproducible way to track the disease. Topological characteristics relating to the distribution and shape of lesions are recognized as important neuroradiological markers in the diagnosis of multiple sclerosis, although these have been much less well characterized quantitatively than have traditional measures such as T2 hyperintense or T1 hypointense lesion volumes. Here, we used voxel-level 3 T magnetic resonance imaging T1-weighted scans to reconstruct the 3D topology of lesions in 284 subjects with multiple sclerosis and tested whether this is a heritable phenotype. To this end, we extracted the genotypes from a published genome-wide association study on these same individuals and searched for genetic associations with lesion load, shape and topological distribution. Lesion probability maps were created to identify frequently affected areas and to assess the overall distribution of T1 lesions in the subject population as a whole. We then developed an original algorithm to cluster adjacent lesional voxels (cluxels) in each subject and tested whether cluxel topology was significantly associated with any single-nucleotide polymorphism in our data set. To focus on patterns of lesion distribution, we computed the first 10 principal components. Although principal component 1 correlated with lesion load, none of the remaining orthogonal components correlated with any other known variable. We then conducted genome-wide association studies on each of these and found 31 significant associations (false discovery rate <0.01) with principal component 8, which represents a mode of variation of lesion topology in the population. The majority of the loci can be linked to genes related to immune cell function and to myelin and neural growth; some (SYK, MYT1L, TRAPPC9, SLITKR6 and RIC3) have been previously associated with the distribution of white matter lesions in multiple sclerosis. Finally, we used a bioinformatics approach to identify a network of 48 interacting proteins showing genetic associations (P < 0.01) with cluxel topology in multiple sclerosis. This network also contains proteins expressed in immune cells and is enriched in molecules expressed in the central nervous system that contribute to neural development and regeneration. Our results show how quantitative traits derived from brain magnetic resonance images of patients with multiple sclerosis can be used as dependent variables in a genome-wide association study. With the widespread availability of powerful computing and the availability of genotyped populations, integration of imaging and genetic data sets is likely to become a mainstream tool for understanding the complex biological processes of multiple sclerosis and other brain disorders.

Extracting actionable information from genome scans

Genome-wide association studies discovered numerous genetic variants significantly associated with various phenotypes. However, significant signals explain only a small portion of the variation in many traits. One explanation is that missing variation is found in "suggestive signals," i.e., variants with reasonably small P-values. However, it is not clear how to capture this information and use it optimally to design and analyze future studies. We propose to extract the available information from a genome scan by accurately estimating the means of univariate statistics. The means are estimated by: (i) computing the sum of squares (SS) of a genome scan's univariate statistics, (ii) using SS to estimate the expected SS for the means (SSM) of univariate statistics, and (iii) constructing accurate soft threshold (ST) estimators for means of univariate statistics by requiring that the SS of these estimators equals the SSM. When compared to competitors, ST estimators explain a substantially higher fraction of the variability in true means. The accuracy of proposed estimators can be used to design two-tier follow-up studies in which regions close to variants having ST-estimated means above a certain threshold are sequenced at high coverage and the rest of the genome is sequenced at low coverage. This follow-up approach reduces the sequencing burden by at least an order of magnitude when compared to a high coverage sequencing of the whole genome. Finally, we suggest ways in which ST methodology can be used to improve signal detection in future sequencing studies and to perform general statistical model selection.

A genome-wide association study reveals loci influencing height and other conformation traits in horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.

The combination of a genome-wide association study of lymphocyte count and analysis of gene expression data reveals novel asthma candidate genes

Recent genome-wide association studies (GWAS) have identified a number of novel genetic associations with complex human diseases. In spite of these successes, results from GWAS generally explain only a small proportion of disease heritability, an observation termed the ‘missing heritability problem’. Several sources for the missing heritability have been proposed, including the contribution of many common variants with small individual effect sizes, which cannot be reliably found using the standard GWAS approach. The goal of our study was to explore a complimentary approach, which combines GWAS results with functional data in order to identify novel genetic associations with small effect sizes. To do so, we conducted a GWAS for lymphocyte count, a physiologic quantitative trait associated with asthma, in 462 Hutterites. In parallel, we performed a genome-wide gene expression study in lymphoblastoid cell lines from 96 Hutterites. We found significant support for genetic associations using the GWAS data when we considered variants near the 193 genes whose expression levels across individuals were most correlated with lymphocyte counts. Interestingly, these variants are also enriched with signatures of an association with asthma susceptibility, an observation we were able to replicate. The associated loci include genes previously implicated in asthma susceptibility as well as novel candidate genes enriched for functions related to T cell receptor signaling and adenosine triphosphate synthesis. Our results, therefore, establish a new set of asthma susceptibility candidate genes. More generally, our observations support the notion that many loci of small effects influence variation in lymphocyte count and asthma susceptibility.

Sensitivity of genome-wide-association signals to phenotyping strategy: the PROP-TAS2R38 taste association as a benchmark

Natural genetic variation can have a pronounced influence on human taste perception, which in turn may influence food preference and dietary choice. Genome-wide association studies represent a powerful tool to understand this influence. To help optimize the design of future genome-wide-association studies on human taste perception we have used the well-known TAS2R38-PROP association as a tool to determine the relative power and efficiency of different phenotyping and data-analysis strategies. The results show that the choice of both data collection and data processing schemes can have a very substantial impact on the power to detect genotypic variation that affects chemosensory perception. Based on these results we provide practical guidelines for the design of future GWAS studies on chemosensory phenotypes. Moreover, in addition to the TAS2R38 gene past studies have implicated a number of other genetic loci to affect taste sensitivity to PROP and the related bitter compound PTC. None of these other locations showed genome-wide significant associations in our study. To facilitate further, target-gene driven, studies on PROP taste perception we provide the genome-wide list of p-values for all SNPs genotyped in the current study.