Functional and association analysis of frizzled 1 (FZD1) promoter haplotypes with femoral neck geometry

A systematic review and cross-database analysis of single nucleotide polymorphisms underlying hip morphology and osteoarthritis reveals shared mechanisms

OBJECTIVE

Understanding the mechanisms of hip disease, such as osteoarthritis (OA), is crucial to advance their treatment. Such hip diseases often involve specific morphological changes. Genetic variations, called single nucleotide polymorphisms (SNPs), influence various hip morphological parameters. This study investigated the biological relevance of SNPs correlated to hip morphology in genome-wide association studies (GWAS). The SNP-associated genes were compared to genes associated with OA in other joints, aiming to see if the same genes play a role in both hip development and the risk of OA in other lower limb joints.

METHODOLOGY

A systematic literature review was conducted to identify SNPs correlated with hip morphology, based on the Population, Intervention, Comparison, Outcome, and Study (PICOS) framework. Afterwards, Gene Ontology (GO) analysis was performed, using EnrichR, on the SNP-associated genes and compared with non-hip OA-associated genes, across different databases.

RESULTS

Reviewing 49 GWAS identified 436 SNPs associated with hip joint morphology, encompassing variance in bone size, structure and shape. Among the SNP-associated genes, SOX9 plays a pivotal role in size, GDF5 impacts bone structure, and BMP7 affects shape. Overall, skeletal system development, regulation of cell differentiation, and chondrocyte differentiation emerged as crucial processes influencing hip morphology. Eighteen percent of GWAS-identified genes related to hip morphology were also associated with non-hip OA.

CONCLUSION

Our findings indicate the existence of multiple shared genetic mechanisms across hip morphology and OA, highlighting the necessity for more extensive research in this area, as in contrast to the hip, the genetic background on knee or foot morphology remains largely understudied.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

The role of GPCRs in bone diseases and dysfunctions

The superfamily of G protein-coupled receptors (GPCRs) contains immense structural and functional diversity and mediates a myriad of biological processes upon activation by various extracellular signals. Critical roles of GPCRs have been established in bone development, remodeling, and disease. Multiple human GPCR mutations impair bone development or metabolism, resulting in osteopathologies. Here we summarize the disease phenotypes and dysfunctions caused by GPCR gene mutations in humans as well as by deletion in animals. To date, 92 receptors (5 glutamate family, 67 rhodopsin family, 5 adhesion, 4 frizzled/taste2 family, 5 secretin family, and 6 other 7TM receptors) have been associated with bone diseases and dysfunctions (36 in humans and 72 in animals). By analyzing data from these 92 GPCRs, we found that mutation or deletion of different individual GPCRs could induce similar bone diseases or dysfunctions, and the same individual GPCR mutation or deletion could induce different bone diseases or dysfunctions in different populations or animal models. Data from human diseases or dysfunctions identified 19 genes whose mutation was associated with human BMD: 9 genes each for human height and osteoporosis; 4 genes each for human osteoarthritis (OA) and fracture risk; and 2 genes each for adolescent idiopathic scoliosis (AIS), periodontitis, osteosarcoma growth, and tooth development. Reports from gene knockout animals found 40 GPCRs whose deficiency reduced bone mass, while deficiency of 22 GPCRs increased bone mass and BMD; deficiency of 8 GPCRs reduced body length, while 5 mice had reduced femur size upon GPCR deletion. Furthermore, deficiency in 6 GPCRs induced osteoporosis; 4 induced osteoarthritis; 3 delayed fracture healing; 3 reduced arthritis severity; and reduced bone strength, increased bone strength, and increased cortical thickness were each observed in 2 GPCR-deficiency models. The ever-expanding number of GPCR mutation-associated diseases warrants accelerated molecular analysis, population studies, and investigation of phenotype correlation with SNPs to elucidate GPCR function in human diseases.

Hypermethylation of Frizzled1 is associated with Wnt/β-catenin signaling inactivation in mesenchymal stem cells of patients with steroid-associated osteonecrosis

The Wnt/β-catenin signaling pathway is associated with the pathogenesis of steroid-induced osteonecrosis. Our investigation studied whether aberrant CpG island hypermethylation of the FZD1 gene was present in patients with osteonecrosis of the femoral head (ONFH), which results in Wnt/β-catenin signaling inactivation and subsequent cell dysfunction. Bone marrow was collected from the proximal femurs of patients with steroid-associated ONFH (n = 21) and patients with new femoral neck fractures (n = 22), and then mesenchymal stem cells (MSCs) were isolated. We investigated cell viability, the transcription and translation levels of Wnt/β-catenin signaling-related genes, the extent of methylation at CpG islands of the FZD1 promoter, and the osteogenic and adipogenic differentiation abilities of MSCs from the control group and from the ONFH group treated with or without 5′-Aza-dC. According to the results, MSCs from the ONFH group showed a reduced proliferation ability, low transcription and translation levels of FZD1, inhibition of the Wnt/β-catenin signaling pathway, weakened osteogenesis and enhanced adipogenesis ability. Aberrant CpG island hypermethylation of FZD1 was observed in the ONFH group. Treatment with 5’-Aza-dC resulted in de novo FZD1 expression, reactivation of the Wnt/β-catenin signaling pathway and promotion of osteogenesis. Taken together, our study not only provides novel insights into the regulation of the Wnt/β-catenin signaling pathway in this disease but also reveals potential for the use of demethylating agents for the treatment of GC-associated ONFH.

Studies of genetic and molecular signaling processes in the bone disease osteonecrosis, when it is associated with steroid use, reveal insights into disease development and suggest new approaches for treatment. Steroid drugs increase the risk of osteonecrosis, in which bone tissue dies due to insufficient blood supply, but the mechanism of this effect is unclear. Researchers in China, led by Zhibo Sun at Wuhan University, investigated a link between the aberrant addition of methyl groups (CH₃) to the DNA of a specific gene and the onset of the disease. They identified an important molecular signaling pathway in cultured bone marrow cells from patients that is inhibited by the gene methylation. Treating these cells with a drug that inhibits methylation led to reactivation of the gene and the associated signalling pathway that promotes healthy bone formation.

Transcriptional Regulation of Frizzled-1 in Human Osteoblasts by Sp1

The wingless pathway has a powerful influence on bone metabolism and is a therapeutic target in skeletal disorders. Wingless signaling is mediated in part through the Frizzled (FZD) receptor family. FZD transcriptional regulation is poorly understood. Herein we tested the hypothesis that Sp1 plays an important role in the transcriptional regulation of FZD1 expression in osteoblasts and osteoblast mineralization. To test this hypothesis, we conducted FZD1 promoter assays in Saos2 cells with and without Sp1 overexpression. We found that Sp1 significantly up-regulates FZD1 promoter activity in Saos2 cells. Chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift (EMSA) assays identified a novel and functional Sp1 binding site at -44 to -40 from the translation start site in the FZD1 promoter. The Sp1-dependent activation of the FZD1 promoter was abolished by mithramycin A (MMA), an antibiotic affecting both Sp1 binding and Sp1 protein levels in Saos2 cells. Similarly, down-regulation of Sp1 in hFOB cells resulted in less FZD1 expression and lower alkaline phosphatase activity. Moreover, over-expression of Sp1 increased FZD1 expression and Saos2 cell mineralization while MMA decreased Sp1 and FZD1 expression and Saos2 cell mineralization. Knockdown of FZD1 prior to Sp1 overexpression partially abolished Sp1 stimulation of osteoblast differentiation markers. Taken together, our results suggest that Sp1 plays a role in human osteoblast differentiation and mineralization, which is at least partially mediated by Sp1-dependent transactivation of FZD1.

AP2 suppresses osteoblast differentiation and mineralization through down-regulation of Frizzled-1

AP2 is an important transcription factor in bone development. AP2 binds to the human Frizzled 1 (FZD1) promoter and down-regulates FZD1 expression in osteoblasts. In addition, AP2 negatively regulates osteoblast differentiation and mineralization in part through down-regulation of FZD1 expression.

SaOS2 Osteosarcoma cells as an in vitro model for studying the transition of human osteoblasts to osteocytes

The central importance of osteocytes in regulating bone homeostasis is becoming increasingly apparent. However, the study of these cells has been restricted by the relative paucity of cell line models, especially those of human origin. Therefore, we investigated the extent to which SaOS2 human osteosarcoma cells can differentiate into osteocyte-like cells. During culture under the appropriate mineralising conditions, SaOS2 cells reproducibly synthesised a bone-like mineralised matrix and temporally expressed the mature osteocyte marker genes SOST, DMP1, PHEX and MEPE and down-regulated expression of RUNX2 and COL1A1. SaOS2 cells cultured in 3D collagen gels acquired a dendritic morphology, characteristic of osteocytes, with multiple interconnecting cell processes. These findings suggest that SaOS2 cells have the capacity to differentiate into mature osteocyte-like cells under mineralising conditions. PTH treatment of SaOS2 cells resulted in strong down-regulation of SOST mRNA expression at all time points tested. Interestingly, PTH treatment resulted in the up-regulation of RANKL mRNA expression only at earlier stages of differentiation. These findings suggest that the response to PTH is dependent on the differentiation stage of the osteoblast/osteocyte. Together, our results demonstrate that SaOS2 cells can be used as a human model to investigate responses to osteotropic stimuli throughout differentiation to a mature osteocyte-like stage.

E2F1 effects on osteoblast differentiation and mineralization are mediated through up-regulation of frizzled-1

Frizzled homolog 1 (FZD1) is a transmembrane receptor that mediates Wnt signaling. The transcriptional regulation of FZD1 and the role of FZD1 in osteoblast biology are not well understood. We examined the role of E2F1 in FZD1 promoter activation and osteoblast differentiation and mineralization. A putative E2F1 binding site in the FZD1 promoter region was initially identified in silico and characterized further in Saos2 cells in vitro by chromatin immunoprecipitation (ChIP), electrophoretic mobility shift (EMSA) and promoter reporter assays. Over-expression of E2F1 transactivated the FZD1 promoter and increased endogenous FZD1 mRNA and protein levels in Saos2 cells. Over-expression of E2F1 in Saos2 cells up-regulated osteoblast differentiation markers alkaline phosphatase (ALP), type I collagen α (COL1A), and osteocalcin (OCN). Furthermore, E2F1 over-expression enhanced mineralization of differentiated Saos2 cells, whereas siRNA knockdown of FZD1 diminished the effects of E2F1 on osteoblast mineralization. The effects of E2F1 on FZD1 expression and osteoblast mineralization were further confirmed in normal human FOB osteoblasts. Taken together, our experiments demonstrate a role of E2F1 in osteoblast differentiation and mineralization and suggest that FZD1 is required, in part, for E2F1 regulation of osteoblast mineralization.

The genetic pleiotropy of musculoskeletal aging

Musculoskeletal aging is detrimental to multiple bodily functions and starts early, probably in the fourth decade of an individual's life. Sarcopenia is a health problem that is expected to only increase as a greater portion of the population lives longer; prevalence of the related musculoskeletal diseases is similarly expected to increase. Unraveling the biological and biomechanical associations and molecular mechanisms underlying these diseases represents a formidable challenge. There are two major problems making disentangling the biological complexity of musculoskeletal aging difficult: (a) it is a systemic, rather than "compartmental," problem, which should be approached accordingly, and (b) the aging per se is neither well defined nor reliably measurable. A unique challenge of studying any age-related condition is a need of distinguishing between the "norm" and "pathology," which are interwoven throughout the aging organism. We argue that detecting genes with pleiotropic functions in musculoskeletal aging is needed to provide insights into the potential biological mechanisms underlying inter-individual differences insusceptibility to the musculoskeletal diseases. However, exploring pleiotropic relationships among the system's components is challenging both methodologically and conceptually. We aimed to focus on genetic aspects of the cross-talk between muscle and its "neighboring" tissues and organs (tendon, bone, and cartilage), and to explore the role of genetics to find the new molecular links between skeletal muscle and other parts of the "musculoskeleton." Identification of significant genetic variants underlying the musculoskeletal system's aging is now possible more than ever due to the currently available advanced genomic technologies. In summary, a "holistic" genetic approach is needed to study the systems's normal functioning and the disease predisposition in order to improve musculoskeletal health.