Functional characterization of genetic variation in the Frizzled 1 (FZD1) promoter and association with bone phenotypes: more to the LRP5 story?

Transcriptomics reveals the molecular regulation of Chinese medicine formula on improving bone quality in broiler

Skeletal disorder is of concern to the poultry industry as it affects animal welfare and production performance. Traditional Chinese medicine could improve bone quality and reduce the incidence of bone disease, but the molecular regulation of Chinese medicine formula (CMF) on improving bone quality in broilers is still unclear. This study was performed to research the effects of CMF on skeletal performance of Cobb broilers and reveal the molecular regulation. A total of 120 one-day-old Cobb broilers were randomly allocated into 4 equal groups of 30 chickens, with 5 replicates and 6 chickens in each replicate. The control (CON) group was fed a diet without CMF, while the CMF1, CMF2, and CMF3 groups were supplemented with different CMF at 6,000 mg/kg diet, respectively. The broilers were raised to 60 d of age, then bone tissues were collected for biomechanical properties, micro-CT detection and transcriptomic sequencing analysis. The results showed that CMF3 improved the biomechanical properties of broiler tibia, via increasing the elastic modulus (P < 0.05), yield strength (P > 0.05), maximum stress (P < 0.05) and fracture stress (P < 0.05) of the tibia. Micro-CT analysis indicated that CMF3 increased the bone mineral density (BMD), bone volume/total volume (BV/TV), bone surface density (BS/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), and decreased the trabecular separation (Tb.Sp) of femur cancellous bone (P < 0.05). RNA-seq analysis revealed 2,177 differentially expressed genes (DEGs) (|log2FoldChange| ≥ 1, FDR < 0.05) between the CMF3 group and CON group. Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) analysis showed 13 pathways mostly associated with bone growth and development and bone metabolism, and we identified 39 bone-related DEGs. This study suggests that CMF3 could improve bone strength and bone microstructure of broilers, and showed a positive effect on bone performance. Our research could provide a theoretical reference for the development of pollution-free feed additives to improve the skeletal performance of broilers, which could help promote healthy farming of chickens.

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.

WNT5A and Its Receptors in the Bone-Cancer Dialogue

Wnt signaling is critical for tumorigenesis and skeletal remodeling. However, its contribution to the formation of metastatic bone lesions remains poorly defined. One major challenge of unraveling its role in cancer progression is the high complexity of Wnt signaling, which includes numerous ligands, receptors, and inhibitors, with intricate biological effects and specific signaling pathways depending on the cellular context. In this perspective, we summarize the role of the noncanonical Wnt ligand WNT5A in the development and metastatic process of osteotropic cancer entities. We focus on its tumor-suppressive function in breast cancer, tumor promoting effects in melanoma, and ambiguous role in prostate cancer, and discuss potential challenges and opportunities that may be associated with targeting Wnt signaling for cancer therapy and treatment of bone metastases. © 2016 American Society for Bone and Mineral Research.

Planar cell polarity aligns osteoblast division in response to substrate strain

Exposure of bone to dynamic strain increases the rate of division of osteoblasts and also influences the directional organization of the cellular and molecular structure of the bone tissue that they produce. Here, we report that brief exposure to dynamic substrate strain (sufficient to rapidly stimulate cell division) influences the orientation of osteoblastic cell division. The initial proliferative response to strain involves canonical Wnt signaling and can be blocked by sclerostin. However, the strain-related orientation of cell division is independently influenced through the noncanonical Wnt/planar cell polarity (PCP) pathway. Blockade of Rho-associated coiled kinase (ROCK), a component of the PCP pathway, prevents strain-related orientation of division in osteoblast-like Saos-2 cells. Heterozygous loop-tail mutation of the core PCP component van Gogh-like 2 (Vangl2) in mouse osteoblasts impairs the orientation of division in response to strain. Examination of bones from Vangl2 loop-tail heterozygous mice by µCT and scanning electron microscopy reveals altered bone architecture and disorganized bone-forming surfaces. Hence, in addition to the well-accepted role of PCP involvement in response to developmental cues during skeletal morphogenesis, our data reveal that this pathway also acts postnatally, in parallel with canonical Wnt signaling, to transduce biomechanical cues into skeletal adaptive responses. The simultaneous and independent actions of these two pathways appear to influence both the rate and orientation of osteoblast division, thus fine-tuning bone architecture to meet the structural demands of functional loading.

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.

Severe hypoxia exerts parallel and cell-specific regulation of gene expression and alternative splicing in human mesenchymal stem cells

BACKGROUND

The endosteum of the bone marrow provides a specialized hypoxic niche that may serve to preserve the integrity, pluripotency, longevity and stemness of resident mesenchymal stem cells (MSCs). To explore the molecular genetic consequences of such a niche we subjected human (h) MSCs to a pO2 of 4 mmHg and analyzed global gene expression and alternative splicing (AS) by genome-exon microarray and RT-qPCR, and phenotype by western blot and immunostaining.

RESULTS

Out of 446 genes differentially regulated by >2.5-fold, down-regulated genes outnumbered up-regulated genes by 243:203. Exon analyses revealed 60 hypoxia-regulated AS events with splice indices (SI) >1.0 from 53 genes and a correlation between high SI and degree of transcript regulation. Parallel analyses of a publicly available AS study on human umbilical vein endothelial cells (HUVECs) showed that there was a strong cell-specific component with only 11 genes commonly regulated in hMSCs and HUVECs and 17 common differentially spliced genes. Only 3 genes were differentially responsive to hypoxia at the gene (>2.0) and AS levels in both cell types. Functional assignments revealed unique profiles of gene expression with complex regulation of differentiation, extracellular matrix, intermediate filament and metabolic marker genes. Antioxidant genes, striated muscle genes and insulin/IGF-1 signaling intermediates were down-regulated. There was a coordinate induction of 9 out of 12 acidic keratins that along with other epithelial and cell adhesion markers implies a partial mesenchymal to epithelial transition.

CONCLUSIONS

We conclude that severe hypoxia confers a quiescent phenotype in hMSCs that is reflected by both the transcriptome profile and gene-specific changes of splicosome actions. The results reveal that severe hypoxia imposes markedly different patterns of gene regulation of MSCs compared with more moderate hypoxia. This is the first study to report hypoxia-regulation of AS in stem/progenitor cells and the first molecular genetic characterization of MSC in a hypoxia-induced quiescent immobile state.

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.

Mechanisms of RAS/β-catenin interactions

Signaling through the WNT/β-catenin and the RAS (rat sarcoma)/MAPK (mitogen-activated protein kinase) pathways plays a key role in the regulation of various physiological cellular processes including proliferation, differentiation, and cell death. Aberrant mutational activation of these signaling pathways is closely linked to the development of cancer in many organs, in humans as well as in laboratory animals. Over the past years, more and more evidence for a close linkage of the two oncogenic signaling cascades has accumulated. Using different experimental approaches, model systems, and experimental conditions, a variety of molecular mechanisms have been identified by which signal transduction through WNT/β-catenin and RAS interact, either in a synergistic or an antagonistic manner. Mechanisms of interaction comprise an upstream crosstalk at the level of pathway-activating ligands and their receptors, interrelations of cytosolic kinases involved in either pathways, as well as interaction in the nucleus related to the joint regulation of target gene transcription. Here, we present a comprehensive review of the current knowledge on the interaction of RAS/MAPK- and WNT/β-catenin-driven signal transduction in mammalian cells.

Osteogenic and chondrogenic potential of biomembrane cells from the PMMA-segmental defect rat model

A layer of cells (the "biomembrane") has been identified in large segmental defects between bone and surgically placed methacrylate spacers or antibiotic-impregnated cement beads. We hypothesize that this contains a pluripotent stem cell population with potential valuable applications in orthopedic tissue engineering. Objectives using biomembranes harvested from rat segmental defects were to: (1) Culture biomembrane cells in specialized media to direct progenitor cells along bone or cartilage cell differentiation lineages; (2) evaluate harvested biomembranes for mesenchymal stem cell markers, and (3) define relevant gene expression patterns in harvested biomembranes using microarray analysis. Culture in osteogenic media produced mineralized nodules; culture in chondrogenic media produced masses containing chondroitin sulfate/sulfated proteoglycans. Molecular analysis of biomembrane cells versus control periosteum showed significant upregulation of key genes functioning in mesenchymal stem cell differentiation, development, maintenance, and proliferation. Results identified significant upregulation of WNT receptor signaling pathway genes and significant upregulation of BMP signaling pathway genes. Findings confirm that the biomembrane has a pluripotent stem cell population. The ability to heal large bone defects is clinically challenging, and novel tissue engineering uses of the biomembrane hold great promise in treating non-unions, open fractures with large bone loss and/or infections, and defects associated with tumor resection.

Molecular genetic studies of gene identification for osteoporosis: the 2009 update

Osteoporosis is a complex human disease that results in increased susceptibility to fragility fractures. It can be phenotypically characterized using several traits, including bone mineral density, bone size, bone strength, and bone turnover markers. The identification of gene variants that contribute to osteoporosis phenotypes, or responses to therapy, can eventually help individualize the prognosis, treatment, and prevention of fractures and their adverse outcomes. Our previously published reviews have comprehensively summarized the progress of molecular genetic studies of gene identification for osteoporosis and have covered the data available to the end of September 2007. This review represents our continuing efforts to summarize the important and representative findings published between October 2007 and November 2009. The topics covered include genetic association and linkage studies in humans, transgenic and knockout mouse models, as well as gene-expression microarray and proteomics studies. Major results are tabulated for comparison and ease of reference. Comments are made on the notable findings and representative studies for their potential influence and implications on our present understanding of the genetics of osteoporosis.

Variation in the PTH gene, hip fracture, and femoral neck geometry in elderly women

Parathyroid hormone (PTH) is a principal regulator of calcium homeostasis. Previously, we studied single-nucleotide polymorphisms present in the major genes in the PTH pathway (PTH, PTHrP, PTHR1, PTHR2) in relation to bone mineral density (BMD) and fracture incidence. We found that haplotypes of the PTH gene were associated with fracture risk independent of BMD. In the present study, we evaluated the relationship between PTH haplotypes and femoral neck bone size. Hip structure analysis and BMD of the femoral neck was assessed by DXA in elderly women from the Malmö Osteoporosis Prospective Risk Assessment study. Data on hip fracture, sustained as a result of low trauma, after the age of 45 years were also analyzed. Haplotypes derived from six polymorphisms in the PTH locus were analyzed in 750 women. Carriers of haplotype 9 had lower values for hip geometry parameters cross-sectional moment of inertia (P = 0.029), femoral neck width (P = 0.049), and section modulous (P = 0.06), suggestive of increased fracture risk at the hip. However, this did not translate into an increased incidence of hip fracture in the studied population. Women who suffered a hip fracture compared to those who had not had longer hip axis length (HAL) (P < 0.001). HAL was not significantly different among haplotypes. Polymorphisms in the PTH gene are associated with differences in aspects of femoral neck geometry in elderly women; however, the major predictor of hip fracture in our population was HAL, to which PTH gene variation does not contribute significantly.

Common genetic variation in the DKK1 gene is associated with hip axis length but not with bone mineral density and bone turnover markers in young adult men: results from the Odense Androgen Study

LRP5 was recently confirmed as an important susceptibility gene for osteoporosis. Our objective was to evaluate the effect of DKK1 polymorphisms on bone mineral density (BMD), hip geometry, and bone turnover. DKK1 is a secreted protein that binds to LRP5/6 receptors and inhibits canonical Wnt signaling. Using HapMap, we selected three SNPs covering the genetic variation in a 13.53-kb region comprising DKK1. The Odense Androgen Study is a population-based study comprising 783 Caucasian men aged 20-29 years. BMD and hip structural parameters were available for study. Bone turnover markers were used as a secondary end point. All analyses were repeated after adjusting for covariables and in subgroups according to physical activity. We found no significant association between DKK1 and BMD or markers of bone turnover; however, a significant association (P = 0.012) was found for rs1569198 with hip axis length (HAL), independent of BMD and height. Moreover, the association seemed to be driven by the non-sedentary subgroup (P = 0.004). Haplotype analysis further confirmed the association of rs1569198 with HAL. Furthermore, we obtained indications for interaction between DKK1 and LRP5 genotypes for different hip geometry parameters. As almost all variance within the DKK1 gene was covered, we conclude that common variation in this gene does not markedly influence BMD or bone turnover markers in young men. In this population, however, a common SNP in DKK1 does have a significant effect on HAL, implying a possible effect on hip fracture risk in the general population. This finding could be of interest but needs replication in independent populations.

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

Frizzleds are receptors for Wnt signaling and are involved in skeletal morphogenesis. Little is known about the transcriptional regulation of frizzleds in bone cells. In the current study, we determined if two common and potentially functional genetic variants (rs2232157, rs2232158) in the frizzled-1 (FZD1) promoter region and their haplotypes influence FZD1 promoter activity in human osteoblast-like cells. We also determined if these variants are associated with femoral neck bone mineral density (BMD) and geometry in 1319 African ancestry men aged > or =40 years. Real-time quantitative PCR and western blot analysis demonstrated FZD1 mRNA and protein expression in the human osteoblast-like cell lines, MG63 and SaOS-2. Promoter activity was next assessed by transient expression of haplotype specific FZD1 promoter reporter plasmids in these cells. In comparison to the common GG haplotype, promoter activity was 3-fold higher for the TC haplotype in both cell lines (p<0.05). We previously demonstrated that rs2232158 is associated with differential FZD1 promoter activity and Egr1 binding and thus focused further functional analyses on the rs2232157 G-to-T polymorphism. Electrophoretic mobility shift assay demonstrated that distinct nuclear protein complexes were associated with rs2232157 in an allele specific manner. Bioinformatics analysis predicted that the G to T transversion creates an E2F1 binding site, further supporting the functional significance of rs2232157 in FZD1 promoter regulation. Individual SNPs and haplotypes were not associated with femoral neck BMD. The TC haplotype was associated with larger subperiosteal width and greater CSMI (p<0.05). These results suggest that FZD1 expression is regulated in a haplotype-dependent manner in osteoblasts and that these same haplotypes may be associated with biomechanical indices of bone strength.

No association of the polymorphisms of the frizzled-related protein gene with peak bone mineral density in Chinese nuclear families

Background

The Wnt/beta-catenin signaling pathway plays an important role in skeletal development. Polymorphisms of frizzled-related protein (FRZB), an antagonist of this pathway, may generate variations in bone mineral density (BMD). In this study, we analyzed the association between FRZB genotypes and peak BMD variation in the spines and hips of two relatively large samples of Chinese female-offspring and male-offspring nuclear families.

Methods

We recruited 1,260 subjects from 401 female-offspring nuclear families and 1,296 subjects from 427 male-offspring nuclear families and genotyped four tagging single nucleotide polymorphisms (tagSNPs) (rs6433993, rs409238, rs288324, and rs4666865) spanning the entire FRZB gene. The SNPs rs288326 and rs7775, which are associated with hip osteoarthritis, were not selected in this study because of their low minor allele frequencies (MAFs) in Chinese people. The quantitative transmission disequilibrium test (QTDT) was used to analyze the association between each SNP and haplotype with peak BMD in female- and male-offspring nuclear families.

Results

In the female-offspring nuclear families, we found no evidence of an association between either single SNPs or haplotypes and peak BMD in the spine or hip. In the male-offspring nuclear families, no within-family association was observed for either SNPs or haplotypes, although a significant total association was found between rs4666865 and spine BMD (P = 0.0299).

Conclusion

Our results suggest that natural variation in FRZB is not a major contributor to the observed variability in peak BMD in either Chinese females or males. Because ethnic differences in the FRZB genotypes may exist, other studies in different population are required to confirm such results.

Genetics of osteoporosis: accelerating pace in gene identification and validation

Osteoporosis is characterized by low bone mineral density and structural deterioration of bone tissue, leading to an increased risk of fractures. It is the most common metabolic bone disorder worldwide, affecting one in three women and one in eight men over the age of 50. In the past 15 years, a large number of genes have been reported as being associated with osteoporosis. However, only in the past 4 years we have witnessed an accelerated pace in identifying and validating osteoporosis susceptibility loci. This increase in pace is mostly due to large-scale association studies, meta-analyses, and genome-wide association studies of both single nucleotide polymorphisms and copy number variations. A comprehensive review of these developments revealed that, to date, at least 15 genes (VDR, ESR1, ESR2, LRP5, LRP4, SOST, GRP177, OPG, RANK, RANKL, COLIA1, SPP1, ITGA1, SP7, and SOX6) can be reasonably assigned as confirmed osteoporosis susceptibility genes, whereas, another >30 genes are promising candidate genes. Notably, confirmed and promising genes are clustered in three biological pathways, the estrogen endocrine pathway, the Wnt/beta-catenin signaling pathway, and the RANKL/RANK/OPG pathway. New biological pathways will certainly emerge when more osteoporosis genes are identified and validated. These genetic findings may provide new routes toward improved therapeutic and preventive interventions of this complex disease.