The Genetics of PPARG and the Skeleton

Osteoporosis is a complex metabolic bone disorder. Recently it has been appreciated that the “obesity in bone” phenomenon occurs at the expense of bone formation, and that is a key component of the pathology of this disease. Mouse models with altered bone expression levels of peroxisome proliferator-activated receptor gamma (PPARG) impact bone formation, but genetic studies connecting PPARG polymorphisms to skeletal phenotypes in humans have proven to be less than satisfactory. One missense polymorphism in exon one has been linked to low bone mineral density (BMD), but the most studied polymorphism, Pro12Ala, has not yet been examined in the context of skeletal phenotype. The studies to date are a promising start in leading to our understanding of the genetic contribution of PPARG to the phenotypes of BMD and fracture risk

Mouse BMD Quantitative Trait Loci Show Improved Concordance With Human Genome-wide Association Loci When Recalculated on a New, Common Mouse Genetic Map

Bone mineral density (BMD) is a heritable trait, and in mice, over 100 quantitative trait loci (QTLs) have been reported, but candidate genes have been identified for only a small percentage. Persistent errors in the mouse genetic map have negatively affected QTL localization, spurring the development of a new, corrected map. In this study, QTLs for BMD were remapped in 11 archival mouse data sets using this new genetic map. Since these QTLs all were mapped in a comparable way, direct comparisons of QTLs for concordance would be valid. We then compared human genome-wide association study (GWAS) BMD loci with the mouse QTLs. We found that 26 of the 28 human GWAS loci examined were located within the confidence interval of a mouse QTL. Furthermore, 14 of the GWAS loci mapped to within 3 cM of a mouse QTL peak. Lastly, we demonstrated that these newly remapped mouse QTLs can substantiate a candidate gene for a human GWAS locus, for which the peak single-nucleotide polymorphism (SNP) fell in an intergenic region. Specifically, we suggest that MEF2C (human chromosome 5, mouse chromosome 13) should be considered a candidate gene for the genetic regulation of BMD. In conclusion, use of the new mouse genetic map has improved the localization of mouse BMD QTLs, and these remapped QTLs show high concordance with human GWAS loci. We believe that this is an opportune time for a renewed effort by the genetics community to identify the causal variants regulating BMD using a synergistic mouse-human approach. © 2010 American Society for Bone and Mineral Research

SEESAW: detecting isoform-level allelic imbalance accounting for inferential uncertainty.

Detecting allelic imbalance at the isoform level requires accounting for inferential uncertainty, caused by multi-mapping of RNA-seq reads. Our proposed method, SEESAW, uses Salmon and Swish to offer analysis at various levels of resolution, including gene, isoform, and aggregating isoforms to groups by transcription start site. The aggregation strategies strengthen the signal for transcripts with high uncertainty. The SEESAW suite of methods is shown to have higher power than other allelic imbalance methods when there is isoform-level allelic imbalance. We also introduce a new test for detecting imbalance that varies across a covariate, such as time

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome [email protected] Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass

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 circadian-regulated gene, Nocturnin, promotes adipogenesis by stimulating PPAR- nuclear translocation

Nocturnin (NOC) is a circadian-regulated protein related to the yeast family of transcription factors involved in the cellular response to nutrient status. In mammals, NOC functions as a deadenylase but lacks a transcriptional activation domain. It is highly expressed in bone-marrow stromal cells (BMSCs), hepatocytes, and adipocytes. In BMSCs exposed to the PPAR-γ (peroxisome proliferator-activated receptor-γ) agonist rosiglitazone, Noc expression was enhanced 30-fold. Previously, we reported that Noc−/− mice had low body temperature, were protected from diet-induced obesity, and most importantly exhibited absence of Pparg circadian rhythmicity on a high-fat diet. Consistent with its role in influencing BMSCs allocation, Noc−/− mice have reduced bone marrow adiposity and high bone mass. In that same vein, NOC overexpression enhances adipogenesis in 3T3-L1 cells but negatively regulates osteogenesis in MC3T3-E1 cells. NOC and a mutated form, which lacks deadenylase activity, bind to PPAR-γ and markedly enhance PPAR-γ transcriptional activity. Both WT and mutant NOC facilitate nuclear translocation of PPAR-γ. Importantly, NOC-mediated nuclear translocation of PPAR-γ is blocked by a short peptide fragment of NOC that inhibits its physical interaction with PPAR-γ. The inhibitory effect of this NOC-peptide was partially reversed by rosiglitazone, suggesting that effect of NOC on PPAR-γ nuclear translocation may be independent of ligand-mediated PPAR-γ activation. In sum, Noc plays a unique role in the regulation of mesenchymal stem-cell lineage allocation by modulating PPAR-γ activity through nuclear translocation. These data illustrate a unique mechanism whereby a nutrient-responsive gene influences BMSCs differentiation, adipogenesis, and ultimately body composition

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women

© 2021, The Author(s). Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256, 523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1p = 4 × 10−17), arthritis (GDF5p = 4 × 10−13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for bone mineral density GWAS associations

Genome-wide association studies (GWASs) for bone mineral density (BMD) have identified over 1,100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant (Bonferroni <= 0.05) using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.First author draf

Canonical A-to-I and C-to-U RNA Editing Is Enriched at 3′UTRs and microRNA Target Sites in Multiple Mouse Tissues

RNA editing is a process that modifies RNA nucleotides and changes the efficiency and fidelity of the central dogma. Enzymes that catalyze RNA editing are required for life, and defects in RNA editing are associated with many diseases. Recent advances in sequencing have enabled the genome-wide identification of RNA editing sites in mammalian transcriptomes. Here, we demonstrate that canonical RNA editing (A-to-I and C-to-U) occurs in liver, white adipose, and bone tissues of the laboratory mouse, and we show that apparent non-canonical editing (all other possible base substitutions) is an artifact of current high-throughput sequencing technology. Further, we report that high-confidence canonical RNA editing sites can cause non-synonymous amino acid changes and are significantly enriched in 3′ UTRs, specifically at microRNA target sites, suggesting both regulatory and functional consequences for RNA editing