Gene-based Genome-wide Association Study Research Articles

Gene-based GWAS analysis for consecutive studies of GEFOS.

Single-nucleotide polymorphism (SNP)-based genome-wide association studies (GWASs) have already identified about 100 loci associated with bone mineral density (BMD), but these loci only explain a small proportion of heritability to osteoporosis risk. In the present study, we performed a gene-based analysis of the largest GWASs in the bone field to identify additional BMD-associated genes. BMD-associated genes were identified by combining the summary statistic P values of SNPs across individual genes in the two consecutive meta-analyses of GWASs from the Genetic Factors for Osteoporosis (GEFOS) studies. The potential functionality of these genes to bone was partially assessed by differential gene expression analysis. Additionally, the consistency of the identification of potential bone mineral density (BMD)-associated variants were evaluated by estimating the correlation of the P values of the same single-nucleotide polymorphisms (SNPs)/genes between the two consecutive Genetic Factors for Osteoporosis Studies (GEFOS) with largely overlapping samples. Compared to the SNP-based analysis, the gene-based strategy identified additional BMD-associated genes with genome-wide significance and increased their mutual replication between the two GEFOS datasets. Among these BMD-associated genes, three novel genes (UBTF, AAAS, and C11orf58) were partially validated at the gene expression level. The correlation analysis presented a moderately high between-study consistency of potential BMD-associated variants. Gene-based analysis as a supplementary strategy to SNP-based genome-wide association studies, when applied here, is shown that it helped identify some novel BMD-associated genes. In addition to its empirically increased statistical power, gene-based analysis also provides a higher testing stability for identification of BMD genes.

A combined analysis of genetically correlated traits identifies 187 loci and a role for neurogenesis and myelination in intelligence

Intelligence, or general cognitive function, is phenotypically and genetically correlated with many traits, including a wide range of physical, and mental health variables. Education is strongly genetically correlated with intelligence (rg = 0.70). We used these findings as foundations for our use of a novel approach—multi-trait analysis of genome-wide association studies (MTAG; Turley et al. 2017)—to combine two large genome-wide association studies (GWASs) of education and intelligence, increasing statistical power and resulting in the largest GWAS of intelligence yet reported. Our study had four goals: first, to facilitate the discovery of new genetic loci associated with intelligence; second, to add to our understanding of the biology of intelligence differences; third, to examine whether combining genetically correlated traits in this way produces results consistent with the primary phenotype of intelligence; and, finally, to test how well this new meta-analytic data sample on intelligence predicts phenotypic intelligence in an independent sample. By combining datasets using MTAG, our functional sample size increased from 199,242 participants to 248,482. We found 187 independent loci associated with intelligence, implicating 538 genes, using both SNP-based and gene-based GWAS. We found evidence that neurogenesis and myelination—as well as genes expressed in the synapse, and those involved in the regulation of the nervous system—may explain some of the biological differences in intelligence. The results of our combined analysis demonstrated the same pattern of genetic correlations as those from previous GWASs of intelligence, providing support for the meta-analysis of these genetically-related phenotypes.

A Bayesian Gene-Based Genome-Wide Association Study Analysis of Osteosarcoma Trio Data Using a Hierarchically Structured Prior.

Osteosarcoma is considered to be the most common primary malignant bone cancer among children and young adults. Previous studies suggest growth spurts and height to be risk factors for osteosarcoma. However, studies on the genetic cause are still limited given the rare occurrence of the disease. In this study, we investigated in a family trio data set that is composed of 209 patients and their unaffected parents and conducted a genome-wide association study (GWAS) to identify genetic risk factors for osteosarcoma. We performed a Bayesian gene-based GWAS based on the single-nucleotide polymorphism (SNP)-level summary statistics obtained from a likelihood ratio test of the trio data, which uses a hierarchically structured prior that incorporates the SNP-gene hierarchical structure. The Bayesian approach has higher power than SNP-level GWAS analysis due to the reduced number of tests and is robust by accounting for the correlations between SNPs so that it borrows information across SNPs within a gene. We identified 217 genes that achieved genome-wide significance. Ingenuity pathway analysis of the gene set indicated that osteosarcoma is potentially related to TP53, estrogen receptor signaling, xenobiotic metabolism signaling, and RANK signaling in osteoclasts.

Locus and gene-based GWAS meta-analysis identifies new diabetic nephropathy genes.

Objective Assimilation of SNPs Interacting in Synchrony (OASIS) is a locus-based clustering algorithm recently described that can potentially address false positives and negatives in genome-wide association studies (GWAS) of complex disorders. Diabetic nephropathy (DN) is incompletely understood due to a paucity of genes identified despite several GWAS. OASIS was applied to three DN dbGAP GWAS datasets (4725 subjects; 1.06 million SNPs). OASIS identified 19 DN genes which were verified using single variant replication in a standard association study and gene-based analysis using GATES. CARS and FRMD3 were confirmed as DN genes, and five known diabetes-associated genes, viz. NLRP3, INPPL1, PIK3C2G, NRXN3, and TBC1D4, not previously identified using these datasets were discovered. Furthermore, three additional novel DN genes were found which replicated in two sets of analysis, viz. NTN1, EBF2, and DNAH11. Hence, composite analysis with OASIS, gene-based, and single variant association testing can be universally applied to existing GWAS datasets for the identification of new genes.

Variants in TTC25 affect autistic trait in patients with autism spectrum disorder and general population.

Autism spectrum disorder (ASD) is a highly heritable neurodevelopmental disorder with a complex genetic architecture. To identify genetic variants underlying ASD, we performed single-variant and gene-based genome-wide association studies using a dense genotyping array containing over 2.3 million single-nucleotide variants in a discovery sample of 160 families with at least one child affected with non-syndromic ASD using a binary (ASD yes/no) phenotype and a quantitative autistic trait. Replication of the top findings was performed in Psychiatric Genomics Consortium and Erasmus Rucphen Family (ERF) cohort study. Significant association of quantitative autistic trait was observed with the TTC25 gene at 17q21.2 (effect size=10.2, P-value=3.4 × 10-7) in the gene-based analysis. The gene also showed nominally significant association in the cohort-based ERF study (effect=1.75, P-value=0.05). Meta-analysis of discovery and replication improved the association signal (P-valuemeta=1.5 × 10-8). No genome-wide significant signal was observed in the single-variant analysis of either the binary ASD phenotype or the quantitative autistic trait. Our study has identified a novel gene TTC25 to be associated with quantitative autistic trait in patients with ASD. The replication of association in a cohort-based study and the effect estimate suggest that variants in TTC25 may also be relevant for broader ASD phenotype in the general population. TTC25 is overexpressed in frontal cortex and testis and is known to be involved in cilium movement and thus an interesting candidate gene for autistic trait.

Gene-based genome-wide association study identified 19p13.3 for lean body mass

Lean body mass (LBM) is a complex trait for human health. To identify genomic loci underlying LBM, we performed a gene-based genome-wide association study of lean mass index (LMI) in 1000 unrelated Caucasian subjects, and replicated in 2283 unrelated Caucasians subjects. Gene-based association analyses highlighted the significant associations of three genes UQCR, TCF3 and MBD3 in one single locus 19p13.3 (discovery p = 6.10 × 10−5, 1.65 × 10−4 and 1.10 × 10−4; replication p = 2.21 × 10−3, 1.84 × 10−3 and 6.95 × 10−3; combined p = 2.26 × 10−6, 4.86 × 10−6 and 1.15 × 10−5, respectively). These results, together with the known functional relevance of the three genes to LMI, suggested that the 19p13.3 region containing UQCR, TCF3 and MBD3 genes was a novel locus underlying lean mass variation.

Searching for new loci and candidate genes for economically important traits through gene-based association analysis of Simmental cattle

Single-marker genome-wide association study (GWAS) is a convenient strategy of genetic analysis that has been successful in detecting the association of a number of single-nucleotide polymorphisms (SNPs) with quantitative traits. However, analysis of individual SNPs can only account for a small proportion of genetic variation and offers only limited knowledge of complex traits. This inadequacy may be overcome by employing a gene-based GWAS analytic approach, which can be considered complementary to the single-SNP association analysis. Here we performed an initial single-SNP GWAS for bone weight (BW) and meat pH value with a total of 770,000 SNPs in 1141 Simmental cattle. Additionally, 21836 cattle genes collected from the Ensembl Genes 83 database were analyzed to find supplementary evidence to support the importance of gene-based association study. Results of the single SNP-based association study showed that there were 11 SNPs significantly associated with bone weight (BW) and two SNPs associated with meat pH value. Interestingly, all of these SNPs were located in genes detected by the gene-based association study.

Multi-strategy genome-wide association studies identify the DCAF16-NCAPG region as a susceptibility locus for average daily gain in cattle

Average daily gain (ADG) is the most economically important trait in beef cattle industry. Using genome-wide association study (GWAS) approaches, previous studies have identified several causal variants within the PLAG1, NCAPG and LCORL genes for ADG in cattle. Multi-strategy GWASs were implemented in this study to improve detection and to explore the causal genes and regions. In this study, we conducted GWASs based on the genotypes of 1,173 Simmental cattle. In the SNP-based GWAS, the most significant SNPs (rs109303784 and rs110058857, P = 1.78 × 10−7) were identified in the NCAPG intron on BTA6 and explained 4.01% of the phenotypic variance, and the independent and significant SNP (rs110406669, P = 5.18 × 10−6) explained 3.32% of the phenotypic variance. Similarly, in the haplotype-based GWAS, the most significant haplotype block, Hap-6-N1416 (P = 2.56 × 10−8), spanned 12.7 kb on BTA6 and explained 4.85% of the phenotypic variance. Also, in the gene-based GWAS, seven significant genes were obtained which included DCAF16 and NCAPG. Moreover, analysis of the transcript levels confirmed that transcripts abundance of NCAPG (P = 0.046) and DCAF16 (P = 0.046) were significantly correlated with the ADG trait. Overall, our results from the multi-strategy GWASs revealed the DCAF16-NCAPG region to be a susceptibility locus for ADG in cattle.

Susceptibility Genes for Multiple Sclerosis Identified in a Gene-Based Genome-Wide Association Study.

Background and PurposeMultiple sclerosis (MS) is a demyelinating and inflammatory disease of the central nervous system. The aim of this study was to identify more genes associated with MS.MethodsBased on the publicly available data of the single-nucleotide polymorphism-based genome-wide association study (GWAS) from the database of Genotypes and Phenotypes, we conducted a powerful gene-based GWAS in an initial sample with 931 family trios, and a replication study sample with 978 cases and 883 controls. For interesting genes, gene expression in MS-related cells between MS cases and controls was examined by using publicly available datasets.ResultsA total of 58 genes was identified, including 20 "novel" genes significantly associated with MS (p<1.40×10-4). In the replication study, 44 of the 58 identified genes had been genotyped and 35 replicated the association. In the gene-expression study, 21 of the 58 identified genes exhibited differential expressions in MS-related cells. Thus, 15 novel genes were supported by replicated association and/or differential expression. In particular, four of the novel genes, those encoding myelin oligodendrocyte glycoprotein (MOG), coiled-coil alpha-helical rod protein 1 (CCHCR1), human leukocyte antigen complex group 22 (HCG22), and major histocompatibility complex, class II, DM alpha (HLA-DMA), were supported by the evidence of both.ConclusionsThe results of this study emphasize the high power of gene-based GWAS in detecting the susceptibility genes of MS. The novel genes identified herein may provide new insights into the molecular genetic mechanisms underlying MS.

Identification of novel risk genes associated with type1 diabetes mellitus using a genome-wide gene-based association analysis.

Aims/IntroductionType 1 diabetes mellitus is a serious disorder characterized by destruction of pancreatic β-cells, culminating in absolute insulin deficiency. Genetic factors contribute to the susceptibility of type 1 diabetes mellitus. The aim of the present study was to identify more susceptibility genes of type 1 diabetes mellitus.Materials and MethodsWe carried out an initial gene-based genome-wide association study in a total of 4,075 type 1 diabetes mellitus cases and 2,604 controls by using the Gene-based Association Test using Extended Simes procedure. Furthermore, we carried out replication studies, differential expression analysis and functional annotation clustering analysis to support the significance of the identified susceptibility genes.ResultsWe identified 452 genes associated with type 1 diabetes mellitus, even after adapting the genome-wide threshold for significance (P < 9.05E-04). Among these genes, 171 were newly identified for type 1 diabetes mellitus, which were ignored in single-nucleotide polymorphism-based association analysis and were not previously reported. We found that 53 genes have supportive evidence from replication studies and/or differential expression studies. In particular, seven genes including four non-human leukocyte antigen (HLA) genes (RASIP1, STRN4, BCAR1 and MYL2) are replicated in at least one independent population and also differentially expressed in peripheral blood mononuclear cells or monocytes. Furthermore, the associated genes tend to enrich in immune-related pathways or Gene Ontology project terms.ConclusionsThe present results suggest the high power of gene-based association analysis in detecting disease-susceptibility genes. Our findings provide more insights into the genetic basis of type 1 diabetes mellitus.

Gene-based GWAS and biological pathway analysis of the resilience of executive functioning.

Resilience in executive functioning (EF) is characterized by high EF measured by neuropsychological test performance despite structural brain damage from neurodegenerative conditions. We previously reported single nucleotide polymorphism (SNP) genome-wide association study (GWAS) results for EF resilience. Here, we report gene- and pathway-based analyses of the same resilience phenotype, using an optimal SNP-set (Sequence) Kernel Association Test (SKAT) for gene-based analyses (conservative threshold for genome-wide significance = 0.05/18,123 = 2.8 × 10(-6)) and the gene-set enrichment package GSA-SNP for biological pathway analyses (False discovery rate (FDR) < 0.05). Gene-based analyses found a genome-wide significant association between RNASE13 and EF resilience (p = 1.33 × 10(-7)). Genetic pathways involved with dendritic/neuron spine, presynaptic membrane, postsynaptic density, etc., were enriched with association to EF resilience. Although replication of these results is necessary, our findings indicate the potential value of gene- and pathway-based analyses in research on determinants of cognitive resilience.

A Rapid Gene-Based Genome-Wide Association Test with Multivariate Traits

Objectives: A gene-based genome-wide association study (GWAS) provides a powerful alternative to the traditional single single nucleotide polymorphism (SNP) association analysis due to its substantial reduction in the multiple testing burden and possible gain in power due to modeling multiple SNPs within a gene. A gene-based association analysis on multivariate traits is often of interest, but it imposes substantial analytical as well as computational challenges to implement it at a genome-wide level. Methods: We propose a rapid implementation of the multivariate multiple linear regression (RMMLR) approach in unrelated individuals as well as in families. Our approach allows for covariates. Moreover, the asymptotic distribution of the test statistic is not heavily influenced by the linkage disequilibrium (LD) among the SNPs and hence can be used efficiently to perform a gene-based GWAS. We have developed a corresponding R package to implement such multivariate gene-based GWAS with this RMMLR approach. Results: Through extensive simulation, we compared several approaches for both single and multivariate traits. Our RMMLR approach maintained a correct type I error level even for sets of SNPs in strong LD. It also demonstrated a substantial gain in power to detect a gene when it is associated with a subset of the traits. We also studied performances of the approaches on the Minnesota Center for Twin Family Research dataset. Conclusions: In our overall comparison, our RMMLR approach provides an efficient and powerful tool to perform a gene-based GWAS with single or multivariate traits and maintains the type I error appropriately.

Meta-analysis of gene-based genome-wide association studies of bone mineral density in Chinese and European subjects

SummaryGene-based association approach could be regarded as a complementary analysis to the single SNP association analysis. We meta-analyzed the findings from the gene-based association approach using the genome-wide association studies (GWAS) data from Chinese and European subjects, confirmed several well established bone mineral density (BMD) genes, and suggested several novel BMD genes.IntroductionThe introduction of GWAS has greatly increased the number of genes that are known to be associated with common diseases. Nonetheless, such a single SNP GWAS has a lower power to detect genes with multiple causal variants. We aimed to assess the association of each gene with BMD variation at the spine and hip using gene-based GWAS approach.MethodsWe studied 778 Hong Kong Southern Chinese (HKSC) women and 5,858 Northern Europeans (dCG); age, sex, and weight were adjusted in the model. The main outcome measure was BMD at the spine and hip.ResultsNine genes showed suggestive p value in HKSC, while 4 and 17 genes showed significant and suggestive p values respectively in dCG. Meta-analysis using weighted Z-transformed test confirmed several known BMD genes and suggested some novel ones at 1q21.3, 9q22, 9q33.2, 20p13, and 20q12. Top BMD genes were significantly associated with connective tissue, skeletal, and muscular system development and function (p < 0.05). Gene network inference revealed that a large number of these genes were significantly connected with each other to form a functional gene network, and several signaling pathways were strongly connected with these gene networks.ConclusionOur gene-based GWAS confirmed several BMD genes and suggested several novel BMD genes. Genetic contribution to BMD variation may operate through multiple genes identified in this study in functional gene networks. This finding may be useful in identifying and prioritizing candidate genes/loci for further study.Electronic supplementary materialThe online version of this article (doi:10.1007/s00198-011-1779-7) contains supplementary material, which is available to authorized users.

Gene-Centric Genomewide Association Study via Entropy

Genes are the functional units in most organisms. Compared to genetic variants located outside genes, genic variants are more likely to affect disease risk. The development of the human HapMap project provides an unprecedented opportunity for genetic association studies at the genomewide level for elucidating disease etiology. Currently, most association studies at the single-nucleotide polymorphism (SNP) or the haplotype level rely on the linkage information between SNP markers and disease variants, with which association findings are difficult to replicate. Moreover, variants in genes might not be sufficiently covered by currently available methods. In this article, we present a gene-centric approach via entropy statistics for a genomewide association study to identify disease genes. The new entropy-based approach considers genic variants within one gene simultaneously and is developed on the basis of a joint genotype distribution among genetic variants for an association test. A grouping algorithm based on a penalized entropy measure is proposed to reduce the dimension of the test statistic. Type I error rates and power of the entropy test are evaluated through extensive simulation studies. The results indicate that the entropy test has stable power under different disease models with a reasonable sample size. Compared to single SNP-based analysis, the gene-centric approach has greater power, especially when there is more than one disease variant in a gene. As the genomewide genic SNPs become available, our entropy-based gene-centric approach would provide a robust and computationally efficient way for gene-based genomewide association study.

Melting point determination

To address the lack of statistical power and interpretability of genome-wide association studies (GWAS), gene-level analyses combine the p-values of individual single nucleotide polymorphisms (SNPs) into gene statistics. However, using all SNPs mapped to a gene, including those with low association scores, can mask the association signal of a gene. We therefore propose a new two-step strategy, consisting in first selecting the SNPs most associated with the phenotype within a given gene, before testing their joint effect on the phenotype. The recently proposed kernelPSI framework for kernel-based post-selection inference makes it possible to model non-linear relationships between features, as well as to obtain valid p-values that account for the selection step. In this paper, we show how we adapted kernelPSI to the setting of quantitative GWAS, using kernels to model epistatic interactions between neighboring SNPs, and post-selection inference to determine the joint effect of selected blocks of SNPs on a phenotype. We illustrate this tool on the study of two continuous phenotypes from the UKBiobank. We show that kernelPSI can be successfully used to study GWAS data and detect genes associated with a phenotype through the signal carried by the most strongly associated regions of these genes. In particular, we show that kernelPSI enjoys more statistical power than other gene-based GWAS tools, such as SKAT or MAGMA. kernelPSI is an effective tool to combine SNP-based and gene-based analyses of GWAS data, and can be used successfully to improve both statistical performance and interpretability of GWAS.