Functional Genomics Complements Quantitative Genetics in Identifying Disease-Gene Associations

An ultimate goal of genetic research is to understand the connection between genotype and phenotype in order to improve the diagnosis and treatment of diseases. The quantitative genetics field has developed a suite of statistical methods to associate genetic loci with diseases and phenotypes, including quantitative trait loci (QTL) linkage mapping and genome-wide association studies (GWAS). However, each of these approaches have technical and biological shortcomings. For example, the amount of heritable variation explained by GWAS is often surprisingly small and the resolution of many QTL linkage mapping studies is poor. The predictive power and interpretation of QTL and GWAS results are consequently limited. In this study, we propose a complementary approach to quantitative genetics by interrogating the vast amount of high-throughput genomic data in model organisms to functionally associate genes with phenotypes and diseases. Our algorithm combines the genome-wide functional relationship network for the laboratory mouse and a state-of-the-art machine learning method. We demonstrate the superior accuracy of this algorithm through predicting genes associated with each of 1157 diverse phenotype ontology terms. Comparison between our prediction results and a meta-analysis of quantitative genetic studies reveals both overlapping candidates and distinct, accurate predictions uniquely identified by our approach. Focusing on bone mineral density (BMD), a phenotype related to osteoporotic fracture, we experimentally validated two of our novel predictions (not observed in any previous GWAS/QTL studies) and found significant bone density defects for both Timp2 and Abcg8 deficient mice. Our results suggest that the integration of functional genomics data into networks, which itself is informative of protein function and interactions, can successfully be utilized as a complementary approach to quantitative genetics to predict disease risks. All supplementary material is available at http://cbfg.jax.org/phenotype

A Polymorphism in a Gene Encoding Perilipin 4 Is Associated with Height but not with Bone Measures in Individuals from the Framingham Osteoporosis Study

There is increasing interest in identifying new pathways and candidate genes that confer susceptibility to osteoporosis. There is evidence that adipogenesis and osteogenesis may be related, including a common bone marrow progenitor cell for both adipocytes and osteoblasts. Perilipin 1 (PLIN1) and Perilipin 4 (PLIN4) are members of the PATS family of genes and are involved in lipolysis of intracellular lipid deposits. A previous study reported gender-specific associations between one polymorphism of PLIN1 and bone mineral density (BMD) in a Japanese population. We hypothesized that polymorphisms in PLIN1 and PLIN4 would be associated with bone measures in adult Caucasian participants of the Framingham Osteoporosis Study (FOS). We genotyped 1,206 male and 1,445 female participants of the FOS for four single-nucleotide polymorphism (SNPs) in PLIN1 and seven SNPs in PLIN4 and tested for associations with measures of BMD, bone ultrasound, hip geometry, and height. We found several gender-specific significant associations with the measured traits. The association of PLIN4 SNP rs8887, G>A with height in females trended toward significance after simulation testing (adjusted P = 0.07) and remained significant after simulation testing in the combined-sex model (adjusted P = 0.033). In a large study sample of men and women, we found a significant association between one SNP in PLIN4 and height but not with bone traits, suggesting that PATS family genes are not important in the regulation of bone. Identification of genes that influence human height may lead to a better understanding of the processes involved in growth and development

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study in 142,487 individuals from the UK Biobank to identify loci associated with BMD estimated by quantitative ultrasound of the heel (“eBMD”). We identified 307 conditionally independent SNPs attaining genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 novel loci, and several rare variants with large effect sizes. To investigate underlying mechanisms we undertook: 1) bioinformatic, functional genomic annotation and human osteoblast expression studies; 2) gene function prediction; 3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and 4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. These studies strongly implicate GPC6 as a novel determinant of BMD and also identify abnormal skeletal phenotypes in knockout mice for a further 100 prioritized genes.This part of the work was supported by Genome Quebec, Genome Canada and the Canadian Institutes of Health Research (CIHR). This work was supported by the Medical Research Council (Programme Grant MC_UU_12013/4 to D.M.E.), the Wellcome Trust (Strategic Award grant number 101123; project grant 094134; to G.R.W., J.H.D.B. and P.I.C.), the Netherlands Organization for Health Research and Development ZonMw VIDI 016.136.367 (funding to F.R., C.M.-G. and K.T.), the mobility stimuli plan of the European Union Erasmus Mundus Action 2: ERAWEB (programme funding to K.T.), NIAMS, NIH (AR060981 and AR060234 to C.L.A.-B.), the National Health and Medical Research Council (Early Career Fellowship APP1104818 to N.M.W.), the Swedish Research Council (funding to E.G.), the Réseau de Médecine Génétique Appliquée (RMGA; J.A.M.), the Fonds de Recherche du Québec–Santé (FRQS; J.A.M. and J.B.R.), the Natural Sciences and Engineering Research Council of Canada (C.M.T.G.), the J. Gibson and the Ernest Heine Family Foundation (P.I.C.), Arthritis Research UK (ref. 20000; to C.L.G.), the Canadian Institutes of Health Research (J.B.R.), the Jewish General Hospital (J.B.R.), and the Australian Research Council (Future Fellowship FT130101709 to D.M.E.). This research was conducted using the UK Biobank Resource (application number 12703). Access to the UK Biobank study data was funded by the University of Queensland (Early Career Researcher Grant 2014002959 to N.M.W.)

Mouse Genome-Wide Association and Systems Genetics Identify Asxl2 As a Regulator of Bone Mineral Density and Osteoclastogenesis

Significant advances have been made in the discovery of genes affecting bone mineral density (BMD); however, our understanding of its genetic basis remains incomplete. In the current study, genome-wide association (GWA) and co-expression network analysis were used in the recently described Hybrid Mouse Diversity Panel (HMDP) to identify and functionally characterize novel BMD genes. In the HMDP, a GWA of total body, spinal, and femoral BMD revealed four significant associations (−log10P>5.39) affecting at least one BMD trait on chromosomes (Chrs.) 7, 11, 12, and 17. The associations implicated a total of 163 genes with each association harboring between 14 and 112 genes. This list was reduced to 26 functional candidates by identifying those genes that were regulated by local eQTL in bone or harbored potentially functional non-synonymous (NS) SNPs. This analysis revealed that the most significant BMD SNP on Chr. 12 was a NS SNP in the additional sex combs like-2 (Asxl2) gene that was predicted to be functional. The involvement of Asxl2 in the regulation of bone mass was confirmed by the observation that Asxl2 knockout mice had reduced BMD. To begin to unravel the mechanism through which Asxl2 influenced BMD, a gene co-expression network was created using cortical bone gene expression microarray data from the HMDP strains. Asxl2 was identified as a member of a co-expression module enriched for genes involved in the differentiation of myeloid cells. In bone, osteoclasts are bone-resorbing cells of myeloid origin, suggesting that Asxl2 may play a role in osteoclast differentiation. In agreement, the knockdown of Asxl2 in bone marrow macrophages impaired their ability to form osteoclasts. This study identifies a new regulator of BMD and osteoclastogenesis and highlights the power of GWA and systems genetics in the mouse for dissecting complex genetic traits