Pre-B-cell leukemia homeobox 1 (PBX1) shows functional and possible genetic association with bone mineral density variation

Comprehensive summary: the role of PBX1 in development and cancers

PBX1 is a transcription factor that can promote the occurrence of various tumors and play a reg-ulatory role in tumor growth, metastasis, invasion, and drug resistance. Furthermore, a variant generated by fusion of E2A and PBX1, E2A-PBX1, has been found in 25% of patients with childhood acute lymphoblastic leukemia. Thus, PBX1 is a potential therapeutic target for many cancers. Here, we describe the structure of PBX1 and E2A-PBX1 as well as the molecular mecha-nisms whereby these proteins promote tumorigenesis to provide future research directions for developing new treatments. We show that PBX1 and E2A-PBX1 induce the development of highly malignant and difficult-to-treat solid and blood tumors. The development of specific drugs against their targets may be a good therapeutic strategy for PBX1-related cancers. Furthermore, we strongly recommend E2A-PBX1 as one of the genes for prenatal screening to reduce the incidence of childhood hematological malignancies.

Epigenetic modifications in spinal ligament aging

Spinal stenosis is a common degenerative spine disorder in the aged population and the spinal ligament aging is a main contributor to this chronic disease. However, the underlying mechanisms of spinal ligament aging remain unclear. Epigenetics is the study of heritable and reversible changes in the function of a gene or genome that occur without any alteration in the primary DNA sequence. Epigenetic alterations have been demonstrated to play crucial roles in age-related diseases and conditions, and they are recently studied as biomarkers and therapeutic targets in the field of cancer research. The main epigenetic modifications, including DNA methylation alteration, histone modifications as well as dysregulated noncoding RNA modulation, have all been implicated in spinal ligament aging diseases. DNA methylation modulates the expression of critical genes including WNT5A, GDNF, ACSM5, miR-497 and miR-195 during spinal ligament degeneration. Histone modifications widely affect gene expression and obvious histone modification abnormalities have been found in spinal ligament aging. MicroRNAs (miRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) exert crucial regulating effects on spinal ligament aging conditions via targeting various osteogenic or fibrogenic differentiation related genes. To our knowledge, there is no systematic review yet to summarize the involvement of epigenetic mechanisms of spinal ligament aging in degenerative spinal diseases. In this study, we systematically discussed the different epigenetic modifications and their potential functions in spinal ligament aging process.

Extensive Interspecific Gene Flow Shaped Complex Evolutionary History and Underestimated Species Diversity in Rapidly Radiated Dolphins

Recently diverged taxa are often characterized by high rates of hybridization, which can complicate phylogenetic reconstruction. For this reason, the phylogenetic relationships and evolutionary history of dolphins are still not very well resolved; the question of whether the genera Tursiops and Stenella are monophyletic is especially controversial. Here, we performed re-sequencing of six dolphin genomes and combined them with eight previously published dolphin SRA datasets and six whole-genome datasets to investigate the phylogenetic relationships of dolphins and test the monophyly hypothesis of Tursiops and Stenella. Phylogenetic reconstruction with the maximum likelihood and Bayesian methods of concatenated loci, as well as with coalescence analyses of sliding window trees, produced a concordant and well-supported tree. Our studies support the non-monophyletic status of Tursiops and Stenella because the species referred these genera do not form exclusive monophyletic clades. This suggests that the current taxonomy of both genera might not reflect their evolutionary history and may underestimate their diversity. A four-taxon D-statistic (ABBA-BABA) test, five-taxon DFOIL test, and tree-based PhyloNet analyses all showed extensive gene flow across dolphin species, which could explain the instability in resolving phylogenetic relationship of oceanic dolphins with different and limited markers. This study could be a good case to demonstrate how genomic data can reveal complex speciation and phylogeny in rapidly radiating animal groups.

Evaluate the effects of serum urate level on bone mineral density: a genome-wide gene-environment interaction analysis in UK Biobank cohort

INTRODUCTION

Serum urate is associated with BMD and may be a protective factor. However, the exact association and mechanism are still unclear. We performed a genome-wide gene-environmental interaction study (GWGEIS) to explore the interaction effects between gene and urate on BMD, using data from the UK Biobank cohort.

METHODS

A total of 4575 participants for femur total BMD, 4561 participants for L1-L4 BMD, and 237799 participants for heel BMD were included in the present study. Linear regression models were used to test for associations between urate and BMD (femur total BMD, L1-L4 BMD, heel BMD) by R software. GWGEIS was conducted by PLINK 2.0 using a generalize linear model, adjusted for age, sex, weight, smoking behavior, drinking behavior, physical activity and 10 principle components for population structure.

RESULTS

Results showed that urate was positively associated with femur total BMD, L1-L4 BMD and heel BMD and similar findings were observed in both the male and female subgroups. GWGEIS identified 261 genome-wide significant (P < 5.00 × 10^-8) SNP × urate interaction effects for femur total BMD (rs8192585 in NOTCH4, rs116080577 in PBX1, rs9409991 in COL5A1), 17 genome-wide significant SNP × urate interaction effects for heel BMD (rs145344540 in PDE11A and rs78485379 in DKK2), 17 suggestive genome-wide SNP × urate interaction effects (P < 1.00 × 10^-5) for L1-L4 BMD (rs10977015 in PTPRD). We also detected genome-wide significant and suggestive SNP × urate interaction effects for BMD in both the male and female subgroups.

CONCLUSIONS

This study reported several novel candidate genes, and strengthen the evidence of the interactive effects between gene and urate on the variations of BMD.

MicroRNA-181 regulates the development of Ossification of Posterior longitudinal ligament via Epigenetic Modulation by targeting PBX1

Objectives: Ossification of the posterior longitudinal ligament (OPLL) presents as the development of heterotopic ossification in the posterior longitudinal ligament of the spine. The etiology of OPLL is genetically linked, as shown by its high prevalence in Asian populations. However, the molecular mechanism of the disease remains obscure. In this study, we explored the function and mechanism of OPLL-specific microRNAs. Methods: The expression levels of the ossification-related OPLL-specific miR-181 family were measured in normal or OPLL ligament tissues. The effect of miR-181a on the ossification of normal or pathogenic ligament cells was tested using real-time polymerase chain reaction (PCR), Western blot, alizarin red staining and alkaline phosphatase (ALP) staining. The candidate targets of miR-181 were screened using a dual luciferase reporter assay and functional analysis. The link between miR-181a and its target PBX1 was investigated using chromatin immunoprecipitation, followed by real-time PCR detection. Histological and immunohistochemical analysis as well as micro-CT scanning were used to evaluate the effects of miR-181 and its antagonist using both tip-toe-walking OPLL mice and in vivo bone formation assays. Results: Using bioinformatic analysis, we found that miR-181a-5p is predicted to play important roles in the development of OPLL. Overexpression of miR-181a-5p significantly increased the expression of ossification-related genes, staining level of alizarin red and ALP activity, while the inhibition of miR-181a-5p by treatment with an antagomir had the opposite effects. Functional analysis identified PBX1 as a direct target of miR-181a-5p, and we determined that PBX1 was responsible for miR-181a-5p's osteogenic phenotype. By chromatin immunoprecipitation assay, we found that miR-181a-5p controls ligament cell ossification by regulating PBX1-mediated modulation of histone methylation and acetylation levels in the promoter region of osteogenesis-related genes. Additionally, using an in vivo model, we confirmed that miR-181a-5p can substantially increase the bone formation ability of posterior ligament cells and cause increased osteophyte formation in the cervical spine of tip-toe-walking mice. Conclusions: Our data unveiled the mechanism by which the miR-181a-5p/PBX1 axis functions in the development of OPLL, and it revealed the therapeutic effects of the miR-181a-5p antagomir in preventing OPLL development both in vivo and in vitro. Our work is the first to demonstrate that microRNA perturbation could modulate the development of OPLL through epigenetic regulation.

NANOG Attenuates Hair Follicle-Derived Mesenchymal Stem Cell Senescence by Upregulating PBX1 and Activating AKT Signaling

Stem cells derived from elderly donors or harvested by repeated subculture exhibit a marked decrease in proliferative capacity and multipotency, which not only compromises their therapeutic potential but also raises safety concerns for regenerative medicine. NANOG-a well-known core transcription factor-plays an important role in maintaining the self-renewal and pluripotency of stem cells. Unfortunately, the mechanism that NANOG delays mesenchymal stem cell (MSC) senescence is not well-known until now. In our study, we showed that both ectopic NANOG expression and PBX1 overexpression (i) significantly upregulated phosphorylated AKT (p-AKT) and PARP1; (ii) promoted cell proliferation, cell cycle progression, and osteogenesis; (iii) reduced the number of senescence-associated-β-galactosidase- (SA-β-gal-) positive cells; and (iv) downregulated the expression of p16, p53, and p21. Western blotting and dual-luciferase activity assays showed that ectopic NANOG expression significantly upregulated PBX1 expression and increased PBX1 promoter activity. In contrast, PBX1 knockdown by RNA interference in hair follicle- (HF-) derived MSCs that were ectopically expressing NANOG resulted in the significant downregulation of p-AKT and the upregulation of p16 and p21. Moreover, blocking AKT with the PI3K/AKT inhibitor LY294002 or knocking down AKT via RNA interference significantly decreased PBX1 expression, while increasing p16 and p21 expression and the number of SA-β-gal-positive cells. In conclusion, our findings show that NANOG delays HF-MSC senescence by upregulating PBX1 and activating AKT signaling and that a feedback loop likely exists between PBX1 and AKT signaling.

Hematopoietic PBX-interacting protein mediates cartilage degeneration during the pathogenesis of osteoarthritis

Osteoarthritis (OA) has been recognized as the most common chronic age-related disease. Cartilage degeneration influences OA therapy. Here we report that hematopoietic pre-B cell leukemia transcription factor-interacting protein (HPIP) is essential for OA development. Elevated HPIP levels are found in OA patients. Col2a1-CreER^T2/HPIP^f/f mice exhibit obvious skeletal abnormalities compared with their HPIP^f/f littermates. HPIP deficiency in mice protects against developing OA. Moreover, intra-articular injection of adeno-associated virus carrying HPIP-specific short hairpin RNA in vivo attenuates OA histological signs. Notably, in vitro RNA-sequencing and chromatin immunoprecipitation sequencing profiles identify that HPIP modulates OA cartilage degeneration through transcriptional activation of Wnt target genes. Mechanistically, HPIP promotes the transcription of Wnt targets by interacting with lymphoid enhancer binding factor 1 (LEF1). Furthermore, HPIP potentiates the transcriptional activity of LEF1 and acetylates histone H3 lysine 56 in the promoters of Wnt targets, suggesting that HPIP is an attractive target in OA regulatory network.

Muscle-restricted nuclear receptor interaction protein knockout causes motor neuron degeneration through down-regulation of myogenin at the neuromuscular junction

BACKGROUND

Nuclear receptor interaction protein (NRIP) is a calcium/calmodulin (CaM) binding protein. Nuclear receptor interaction protein interacts with CaM to activate calcineurin and CaMKII signalling. The conventional NRIP knockout mice (global knockout) showed muscular abnormality with reduction of muscle oxidative functions and motor function defects.

METHODS

To investigate the role of NRIP on neuromuscular system, we generated muscle-restricted NRIP knockout mice [conditional knockout (cKO)]. The muscle functions (including oxidative muscle markers and muscle strength) and lumbar motor neuron functions [motor neuron number, axon denervation, neuromuscular junction (NMJ)] were tested. The laser-captured microdissection at NMJ of skeletal muscles and adenovirus gene therapy for rescued effects were performed.

RESULTS

The cKO mice showed muscular abnormality with reduction of muscle oxidative functions and impaired motor performances as global knockout mice. To our surprise, cKO mice also displayed motor neuron degeneration with abnormal architecture of NMJ. Specifically, the cKO mice revealed reduced motor neuron number with small neuronal size in lumbar spinal cord as well as denervating change, small motor endplates, and decreased myonuclei number at NMJ in skeletal muscles. To explore the mechanisms, we screened various muscle-derived factors and found that myogenin is a potential candidate that myogenin expression was lower in skeletal muscles of cKO mice than wild-type mice. Because NRIP and myogenin were colocalized around acetylcholine receptors at NMJ, we extracted RNA from synaptic and extrasynaptic regions of muscles using laser capture microdissection and showed that myogenin expression was especially lower at synaptic region in cKO than wild-type mice. Notably, overexpression of myogenin using intramuscular adenovirus encoding myogenin treatment rescued abnormal NMJ architecture and preserved motor neuron death in cKO mice.

CONCLUSIONS

In summary, we demonstrated that deprivation of NRIP decreases myogenin expression at NMJ, possibly leading to abnormal NMJ formation, denervation of acetylcholine receptor, and subsequent loss of spinal motor neuron. Overexpression of myogenin in cKO mice can partially rescue abnormal NMJ architecture and motor neuron death. Therefore, muscular NRIP is a novel trophic factor supporting spinal motor neuron via stabilization of NMJ by myogenin expression.

Identification and Analysis of Human Sex-biased MicroRNAs

Sex differences are widely observed under various circumstances ranging from physiological processes to therapeutic responses, and a myriad of sex-biased genes have been identified. In recent years, transcriptomic datasets of microRNAs (miRNAs), an important class of non-coding RNAs, become increasingly accessible. However, comprehensive analysis of sex difference in miRNA expression has not been performed. Here, we identified the differentially-expressed miRNAs between males and females by examining the transcriptomic datasets available in public databases and conducted a systemic analysis of their biological characteristics. Consequently, we identified 73 female-biased miRNAs (FmiRs) and 163 male-biased miRNAs (MmiRs) across four tissues including brain, colorectal mucosa, peripheral blood, and cord blood. Our results suggest that compared to FmiRs, MmiRs tend to be clustered in the human genome and exhibit higher evolutionary rate, higher expression tissue specificity, and lower disease spectrum width. In addition, functional enrichment analysis of miRNAs show that FmiR genes are significantly associated with metabolism process and cell cycle process, whereas MmiR genes tend to be enriched for functions like histone modification and circadian rhythm. In all, the identification and analysis of sex-biased miRNAs together could provide new insights into the biological differences between females and males and facilitate the exploration of sex-biased disease susceptibility and therapy.

Polymorphisms in sex steroid receptors: From gene sequence to behavior

Sex steroid receptors have received much interest as potential mediators of human behaviors and mental disorders. Candidate gene association studies have identified about 50 genetic variants of androgen and estrogen receptors that correlate with human behavioral phenotypes. Because most of these polymorphisms lie outside coding regions, discerning their effect on receptor function is not straightforward. Thus, although discoveries of associations improve our ability to predict risk, they have not greatly advanced our understanding of underlying mechanisms. This article is intended to serve as a starting point for psychologists and other behavioral biologists to consider potential mechanisms. Here, I review associations between polymorphisms in sex steroid receptors and human behavioral phenotypes. I then consider ways in which genetic variation can affect processes such as mRNA transcription, splicing, and stability. Finally, I suggest ways that hypotheses about mechanism can be tested, for example using in vitro assays and/or animal models.

Genome-Wide DNA Methylation Study Identifies Significant Epigenomic Changes in Osteoarthritic Subchondral Bone and Similarity to Overlying Cartilage

OBJECTIVE

To perform a genome-wide DNA methylation study to identify differential DNA methylation patterns in subchondral bone underlying eroded and intact cartilage from patients with hip osteoarthritis (OA) and to compare these with DNA methylation patterns in overlying cartilage.

METHODS

Genome-wide DNA methylation profiling using Illumina HumanMethylation 450 arrays was performed on eroded and intact cartilage and subchondral bone from within the same joint of 12 patients undergoing hip arthroplasty. Genes with differentially methylated CpG sites were analyzed to identify shared pathways, upstream regulators, and overrepresented gene ontologies, and these patterns were compared with those of the overlying cartilage. Histopathology was graded by modified Mankin score and assessed for correlation with DNA methylation.

RESULTS

We identified 7,316 differentially methylated CpG sites in subchondral bone underlying eroded cartilage, most of which (∼75%) were hypomethylated, and 1,397 sites in overlying eroded cartilage, 126 of which were shared. Samples clustered into 3 groups with distinct histopathologic scores. We observed differential DNA methylation of genes including the RNA interference-processing gene AGO2, the growth factor TGFB3, the OA suppressor NFATC1, and the epigenetic effector HDAC4. Among known susceptibility genes in OA, 32 were differentially methylated in subchondral bone, 8 were differentially methylated in cartilage, and 5 were shared. Upstream regulator analysis using differentially methylated genes in OA subchondral bone showed a strong transforming growth factor β1 signature (P = 1 × 10(-40) ) and a tumor necrosis factor family signature (P = 3.2 × 10(-28) ), among others.

CONCLUSION

Our data suggest the presence of an epigenetic phenotype associated with eroded OA subchondral bone that is similar to that of overlying eroded OA cartilage.

Gene Expression Profiling in Peripheral Blood Cells and Synovial Membranes of Patients with Psoriatic Arthritis

Background

Psoriatic arthritis (PsA) is an inflammatory arthritis whose pathogenesis is poorly understood; it is characterized by bone erosions and new bone formation. The diagnosis of PsA is mainly clinical and diagnostic biomarkers are not yet available. The aim of this work was to clarify some aspects of the disease pathogenesis and to identify specific gene signatures in paired peripheral blood cells (PBC) and synovial biopsies of patients with PsA. Moreover, we tried to identify biomarkers that can be used in clinical practice.

Methods

PBC and synovial biopsies of 10 patients with PsA were used to study gene expression using Affymetrix arrays. The expression values were validated by Q-PCR, FACS analysis and by the detection of soluble mediators.

Results

Synovial biopsies of patients showed a modulation of approximately 200 genes when compared to the biopsies of healthy donors. Among the differentially expressed genes we observed the upregulation of Th17 related genes and of type I interferon (IFN) inducible genes. FACS analysis confirmed the Th17 polarization. Moreover, the synovial trascriptome shows gene clusters (bone remodeling, angiogenesis and inflammation) involved in the pathogenesis of PsA. Interestingly 90 genes are modulated in both compartments (PBC and synovium) suggesting that signature pathways in PBC mirror those of the inflamed synovium. Finally the osteoactivin gene was upregulared in both PBC and synovial biopsies and this finding was confirmed by the detection of high levels of osteoactivin in PsA sera but not in other inflammatory arthritides.

Conclusions

We describe the first analysis of the trancriptome in paired synovial tissue and PBC of patients with PsA. This study strengthens the hypothesis that PsA is of autoimmune origin since the coactivity of IFN and Th17 pathways is typical of autoimmunity. Finally these findings have allowed the identification of a possible disease biomarker, osteoactivin, easily detectable in PsA serum.

Estrogen Receptor Alpha as a Mediator of Life-History Trade-offs

Trade-offs between competitive and parental strategies often are mediated by sex steroids. The mechanisms underlying steroid signaling and metabolism may therefore serve as targets of disruptive selection that leads to alternative behavioral phenotypes. White-throated sparrows exhibit two color morphs that differ in both competitive and parental behavior; white-striped (WS) birds engage in more territorial singing, whereas tan-striped (TS) birds provision nestlings more often. Although WS birds have higher levels of plasma testosterone (T) and estradiol than do TS birds, experimental equalization of these hormones does not abolish morph differences in singing. Neural sensitivity to sex steroids may differ between the morphs because the gene for estrogen receptor alpha (ERα) has been captured by a chromosomal rearrangement found only in the WS birds. We recently showed that expression of this gene differs between the morphs and may drive the behavioral polymorphism. First, the ERα promoter region contains fixed polymorphisms that affect transcription efficiency in vitro. Second, in a free-living population, local expression of ERα depends strongly on morph and predicts both territorial singing and parental provisioning. Differential ERα expression is particularly striking in the medial amygdala; WS birds have three times more ERα mRNA than do TS birds. This difference persists during the non-breeding season and is unaffected by exogenous T treatment. Finally, preliminary data generated by RNA-seq confirm that ERα expression in MeA is both differentially expressed and correlated with territorial singing. Together, these results suggest that ERα may be a target of disruptive selection that leads to alternative behavioral strategies. Our future directions include a more detailed analysis of the ERα promoter regions to determine the molecular basis of differential expression as well as gene network analyses to identify genes connected to ERα.

NRIP is a novel Z-disc protein to activate calmodulin signaling for skeletal muscle contraction and regeneration

Nuclear receptor interaction protein (NRIP, also known as DCAF6 and IQWD1) is a calcium-dependent calmodulin binding protein (Ca2+/CaM). In this study, we found that NRIP is a novel Z-disc protein in skeletal muscle. NRIP knockout mice (NRIP KO) were generated and found to have reduced muscle strength, susceptibility to fatigue and impaired adaptive exercise performance. The mechanisms of NRIP-regulated muscle contraction depend on NRIP being downstream of calcium signaling, where it stimulates phosphorylation of both calcineurin-nuclear factor of activated T-cells, cytoplasmic 1 (CaN-NFATc1) and calmodulin-dependent protein kinase II (CaMKII) through interaction with CaM, resulting in the induction of slow myosin gene expression and mitochondrial activity, and balancing of Ca2+ homeostasis of the internally stored Ca2+ of the sarcoplasmic reticulum. Moreover, NRIP KO mice have delayed regenerative capacity. The amount of NRIP can be enhanced after muscle injury and is responsible for muscle regeneration, coupled with the increased expression of myogenin, desmin and embryonic myosin heavy chain for myogenesis, as well as myotube formation. In conclusion, NRIP is a novel Z-disc protein important for skeletal muscle strength and regenerative capacity.

Estrogen receptor α polymorphism in a species with alternative behavioral phenotypes

The evolution of behavior relies on changes at the level of the genome; yet the ability to attribute a behavioral change to a specific, naturally occurring genetic change is rare in vertebrates. In the white-throated sparrow (Zonotrichia albicollis), a chromosomal polymorphism (ZAL2/2(m)) is known to segregate with a behavioral phenotype. Individuals with the ZAL2(m) haplotype engage in more territorial aggression and less parental behavior than individuals without it. These behaviors are thought to be mediated by sensitivity to sex steroids, and the chromosomal rearrangement underlying the polymorphism has captured a prime candidate gene: estrogen receptor 1 (ESR1), which encodes estrogen receptor α (ERα). We therefore hypothesized that the behavioral effects of the ZAL2(m) rearrangement are mediated by polymorphism in ESR1. We report here that (i) the ESR1 promoter region contains fixed polymorphisms distinguishing the ZAL2(m) and ZAL2 alleles; (ii); those polymorphisms regulate transcription efficiency in vitro and therefore potentially do the same in vivo (iii); the local expression of ERα in the brain depends strongly on genotype in a free-living population; and (iv) ERα expression in the medial amygdala and medial preoptic area may fully mediate the effects of genotype on territorial aggression and parenting, respectively. Thus, our study provides a rare glimpse of how a chromosomal polymorphism has affected the brain and social behavior in a vertebrate. Our results suggest that in this species, differentiation of ESR1 has played a causal role in the evolution of phenotypes with alternative life-history strategies.

Association of CDX1 binding site of periostin gene with bone mineral density and vertebral fracture risk

SUMMARY

Periostin (POSTN) as a regulator of osteoblast differentiation and bone formation may affect susceptibility to osteoporosis. This study suggests POSTN as a candidate gene for bone mineral density (BMD) variation and vertebral fracture risk, which could better our understanding about the genetic pathogenesis of osteoporosis and will be useful in clinic in the future.

INTRODUCTION

The genetic determination of osteoporosis is complex and ill-defined. Periostin (POSTN), an extracellular matrix secreted by osteoblasts and a regulator of osteoblast differentiation and bone formation, may affect susceptibility to osteoporosis.

METHODS

We adopted a tag-single nucleotide polymorphism (SNP) based association method followed by imputation-based verification and identification of a causal variant. The association was investigated in 1,572 subjects with extreme-BMD and replicated in an independent population of 2,509 subjects. BMD was measured by dual X-ray absorptiometry. Vertebral fractures were identified by assessing vertebral height from X-rays of the thoracolumbar spine. Association analyses were performed with PLINK toolset and imputation analyses with MACH software. The top imputation finding was subsequently validated by genotyping. Interactions between POSTN and another BMD-related candidate gene sclerostin (SOST) were analyzed using MDR program and validated by logistical regression analyses. The putative transcription factor binding with target sequence was confirmed by electrophoretic mobility shift assay (EMSA).

RESULTS

Several SNPs of POSTN were associated with BMD or vertebral fractures. The most significant polymorphism was rs9547970, located at the -2,327 bp upstream (P = 6.8 × 10(-4)) of POSTN. Carriers of the minor allele G per copy of rs9547970 had 1.33 higher risk of vertebral fracture (P = 0.007). An interactive effect between POSTN and SOST upon BMD variation was suggested (P < 0.01). A specific binding of CDX1 to the sequence of POSTN with the major allele A of rs9547970 but not the variant G allele was confirmed by EMSA.

CONCLUSIONS

Our results suggest POSTN as a candidate gene for BMD variation and vertebral fracture risk.

Genetic epidemiology of age-related osteoporosis and its clinical applications

Osteoporosis is an important and complex disorder that is highly prevalent worldwide. This disease poses a major challenge to modern medicine and its treatment is associated with high costs. Numerous studies have endeavored to decipher the pathogenesis of this disease. The clinical assessment of patients often incorporates information about a family history of osteoporotic fractures. Indeed, the observation of an increased risk of fracture in an individual with a positive parental history of hip fracture provides strong evidence for the heritability of osteoporosis. The onset and progression of osteoporosis are generally controlled by multiple genetic and environmental factors, as well as interactions between them, with rare cases determined by a single gene. In an attempt to identify the genetic markers of complex diseases such as osteoporosis, there has been a move away from traditional linkage mapping studies and candidate gene association studies to higher-density genome-wide association studies. The advent of high-throughput technology enables genotyping of millions of DNA markers in the human genome, and consequently the identification and characterization of causal variants and loci that underlie osteoporosis. This Review presents an overview of the major findings since 2007 and clinical applications of these genome-wide linkage and association studies.

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.

Hip geometry variation is associated with bone mineralization pathway gene variants: The Framingham Study

Mineralization of bone matrix is an important process in bone formation; thus defects in mineralization have been implicated in bone mineral density (BMD) and bone structure alterations. Three central regulators of phosphate balance, ALPL, ANKH, and ENPP1, are central in the matrix mineralization process; therefore, the genes encoding them are considered important candidates genes for BMD and bone geometry. To test for an association between these three candidate genes and BMD and bone geometry traits, 124 informative singlenucleotide polymorphisms (SNPs) were selected and genotyped in 1513 unrelated subjects from the Framingham offspring cohort. Initial results showed that SNP rs1974201 in the gene ENPP1 was a susceptibility variant associated with several hip geometric indices, with the strongest p value of 3.8 × 10(7) being observed for femoral neck width. A few modest associations were observed between SNPs in or near ALPL and several bone traits, but no association was observed with ANKH. The association signals observed for SNPs around rs1974201 were attenuated after conditional analysis on rs1974201. Transcription factor binding-site prediction revealed that the HOXA7 binding site was present in the reference sequence with the major allele, whereas this potential binding site is lost in the sequence with the minor allele of rs1974201. In conclusion, we found evidence for association of bone geometry variation with an SNP in ENPP1, a gene in the mineralization pathway. The alteration of a binding site of the deregulator of extracellular matrix HOXA7 warrants further investigation.

Pbx1 represses osteoblastogenesis by blocking Hoxa10-mediated recruitment of chromatin remodeling factors

Abdominal-class homeodomain-containing (Hox) factors form multimeric complexes with TALE-class homeodomain proteins (Pbx, Meis) to regulate tissue morphogenesis and skeletal development. Here we have established that Pbx1 negatively regulates Hoxa10-mediated gene transcription in mesenchymal cells and identified components of a Pbx1 complex associated with genes in osteoblasts. Expression of Pbx1 impaired osteogenic commitment of C3H10T1/2 multipotent cells and differentiation of MC3T3-E1 preosteoblasts. Conversely, targeted depletion of Pbx1 by short hairpin RNA (shRNA) increased expression of osteoblast-related genes. Studies using wild-type and mutated osteocalcin and Bsp promoters revealed that Pbx1 acts through a Pbx-binding site that is required to attenuate gene activation by Hoxa10. Chromatin-associated Pbx1 and Hoxa10 were present at osteoblast-related gene promoters preceding gene expression, but only Hoxa10 was associated with these promoters during transcription. Our results show that Pbx1 is associated with histone deacetylases normally linked with chromatin inactivation. Loss of Pbx1 from osteoblast promoters in differentiated osteoblasts was associated with increased histone acetylation and CBP/p300 recruitment, as well as decreased H3K9 methylation. We propose that Pbx1 plays a central role in attenuating the ability of Hoxa10 to activate osteoblast-related genes in order to establish temporal regulation of gene expression during osteogenesis.

Association of JAG1 with bone mineral density and osteoporotic fractures: a genome-wide association study and follow-up replication studies

Bone mineral density (BMD), a diagnostic parameter for osteoporosis and a clinical predictor of fracture, is a polygenic trait with high heritability. To identify genetic variants that influence BMD in different ethnic groups, we performed a genome-wide association study (GWAS) on 800 unrelated Southern Chinese women with extreme BMD and carried out follow-up replication studies in six independent study populations of European descent and Asian populations including 18,098 subjects. In the meta-analysis, rs2273061 of the Jagged1 (JAG1) gene was associated with high BMD (p = 5.27 x 10(-8) for lumbar spine [LS] and p = 4.15 x 10(-5) for femoral neck [FN], n = 18,898). This SNP was further found to be associated with the low risk of osteoporotic fracture (p = 0.009, OR = 0.7, 95% CI 0.57-0.93, n = 1881). Region-wide and haplotype analysis showed that the strongest association evidence was from the linkage disequilibrium block 5, which included rs2273061 of the JAG1 gene (p = 8.52 x 10(-9) for LS and 3.47 x 10(-5) at FN). To assess the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of c-Myc to the "G" but not "A" allele of rs2273061. A mRNA expression study in both human bone-derived cells and peripheral blood mononuclear cells confirmed association of the high BMD-related allele G of rs2273061 with higher JAG1 expression. Our results identify the JAG1 gene as a candidate for BMD regulation in different ethnic groups, and it is a potential key factor for fracture pathogenesis.

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.

Collaborative meta-analysis: associations of 150 candidate genes with osteoporosis and osteoporotic fracture

BACKGROUND

Osteoporosis is a highly heritable trait. Many candidate genes have been proposed as being involved in regulating bone mineral density (BMD). Few of these findings have been replicated in independent studies.

OBJECTIVE

To assess the relationship between BMD and fracture and all common single-nucleotide polymorphisms (SNPs) in previously proposed osteoporosis candidate genes.

DESIGN

Large-scale meta-analysis of genome-wide association data.

SETTING

5 international, multicenter, population-based studies.

PARTICIPANTS

Data on BMD were obtained from 19 195 participants (14 277 women) from 5 populations of European origin. Data on fracture were obtained from a prospective cohort (n = 5974) from the Netherlands.

MEASUREMENTS

Systematic literature review using the Human Genome Epidemiology Navigator identified autosomal genes previously evaluated for association with osteoporosis. We explored the common SNPs arising from the haplotype map of the human genome (HapMap) across all these genes. BMD at the femoral neck and lumbar spine was measured by dual-energy x-ray absorptiometry. Fractures were defined as clinically apparent, site-specific, validated nonvertebral and vertebral low-energy fractures.

RESULTS

150 candidate genes were identified and 36 016 SNPs in these loci were assessed. SNPs from 9 gene loci (ESR1, LRP4, ITGA1, LRP5, SOST, SPP1, TNFRSF11A, TNFRSF11B, and TNFSF11) were associated with BMD at either site. For most genes, no SNP was statistically significant. For statistically significant SNPs (n = 241), effect sizes ranged from 0.04 to 0.18 SD per allele. SNPs from the LRP5, SOST, SPP1, and TNFRSF11A loci were significantly associated with fracture risk; odds ratios ranged from 1.13 to 1.43 per allele. These effects on fracture were partially independent of BMD at SPP1 and SOST.

LIMITATION

Only common polymorphisms in linkage disequilibrium with SNPs in HapMap could be assessed, and previously reported associations for SNPs in some candidate genes could not be excluded.

CONCLUSION

In this large-scale collaborative genome-wide meta-analysis, 9 of 150 candidate genes were associated with regulation of BMD, 4 of which also significantly affected risk for fracture. However, most candidate genes had no consistent association with BMD.