A whole-genome linkage scan suggests several genomic regions potentially containing quantitative trait Loci for osteoporosis

Controversial causal association between IGF family members and osteoporosis: a Mendelian randomization study between UK and FinnGen biobanks

Objectives

Osteoporosis, a prevalent skeletal disorder characterized by reduced bone strength, is closely linked to the IGF system, crucial for skeletal metabolism. However, the precise nature of this relationship remains elusive. In this study, we employed Mendelian randomization (MR) to unravel the associations between genetically predicted serum IGF system member levels and osteoporosis.

Methods

A two-sample MR approach was employed to investigate these causal associations based on two individual datasets. Predictions of 14 serum levels of IGF system members were made using 11,036,163 relevant Single Nucleotide Polymorphisms (SNPs) within a cohort of 4,301 individuals of European descent. Genetic association estimates for osteoporosis were derived from two publicly available GWAS consortia: the Finnish consortium from the FinnGen biobank, comprising 212,778 individuals of Finnish descent (3,203 cases and 209,575 controls), and the UK consortium from the UK Biobank, including 337,159 individuals of European descent (5,266 cases and 331,893 controls).

Results

According to the UK dataset, IGF-1 levels were associated with a reduced risk of osteoporosis, as indicated by the weighted median method (Odds Ratio [OR] = 0.998, 95% CI = 0.997-1.000, P = 0.032). Additionally, higher levels of IGFBP-3 were linked to a decreased risk of osteoporosis using the Inverse-Variance Weighted (IVW) method (OR = 0.999, 95% CI = 0.998-1.000, P = 0.019), and CTGF levels exhibited a negative association with osteoporosis, as determined by the weighted median method (OR = 0.998, 95% CI = 0.996-0.999, P = 0.004). In the FinnGen dataset, IGF-1 and IGFBP-3 were not identified to be associated with osteoporosis. While, IGF-LR1 levels displayed a negative association with osteoporosis, according to the MR-Egger method (OR = 0.886, 95% CI = 0.795-0.987, P = 0.036), while CYR61 was linked to an increased risk of osteoporosis based on both the weighted median and IVW methods (OR = 1.154, 95% CI = 1.009-1.319, P = 0.037, and OR = 1.115, 95% CI = 1.022-1.215, P = 0.014, respectively).

Conclusion

This study provides compelling evidence that certain IGF family members play a role in the pathogenesis of osteoporosis between different datasets, indicating population specific causal effects between IGF family and osteoporosis. Although the results from both datasets demonstrated that IGF family involved in the pathogenesis of osteoporosis, but the responding key molecules might be various among different population. Subsequent research is warranted to evaluate the potential of these biomarkers as targets for osteoporosis prevention and treatment in specific population.

Oxidative-Stress-Related Genes in Osteoporosis: A Systematic Review

Osteoporosis is characterized by a decline in bone mineral density (BMD) and increased fracture risk. Free radicals and antioxidant systems play a central role in bone remodeling. This study was conducted to illustrate the role of oxidative-stress-related genes in BMD and osteoporosis. A systematic review was performed following the PRISMA guidelines. The search was computed in PubMed, Web of Sciences, Scopus, EBSCO, and BVS from inception to November 1st, 2022. The risk of bias was evaluated using the Joanna Briggs Institute Critical Appraisal Checklist tool. A total of 427 potentially eligible articles exploring this search question were detected. After removing duplicates (n = 112) and excluding irrelevant manuscripts based on screenings of their titles and abstracts (n = 317), 19 articles were selected for full-text review. Finally, 14 original articles were included in this systematic review after we applied the exclusion and inclusion criteria. Data analyzed in this systematic review indicated that oxidative-stress-related genetic polymorphisms are associated with BMD at different skeletal sites in diverse populations, influencing the risk of osteoporosis or osteoporotic fracture. However, it is necessary to look deep into their association with bone metabolism to determine if the findings can be translated into the clinical management of osteoporosis and its progression.

Association between ALOX15 gene polymorphism and brick-tea type skeletal fluorosis in Tibetans, Kazaks and Han, China

To evaluate the association between ALOX15 gene polymorphism and skeletal fluorosis (SF), a case-control study was conducted. A total of 1023 individuals, including 308 Tibetans, 290 Kazaks and 425 Han, were enrolled in this study, in which cases and controls were 278 and 745, respectively. SF was diagnosed by X-ray absorptiometry. SNPs were genotyped using the Sequenom Mass ARRAY system. The genotypes of ALOX15 rs7220870, rs2664593 and rs1107852 were not associated with the risk of SF. After reconstructing the haplotype of rs7220870 and rs11078528, the risk effect of haplotype CA was found in Han participants aged ≤45 years or with moderate fluoride intake. Diplotype of CC/CC had a protective effect on SF risk in Han participants; whereas, CA/CC diplotype showed a risk effect on SF risk in participants aged ≥65; Our results provide the first evidence of an association between ALOX15 gene polymorphism and SF risk in Han participants.Abbreviation: SF: Skeletal fluorosis; SNP: Single Nucleotide polymorphism.

Cytosolic proteome profiling of monocytes for male osteoporosis

In male Caucasians with discordant hip bone mineral density (BMD), we applied the subcellular separation and proteome profiling to investigate the monocytic cytosol. Three BMD-associated proteins (ALDOA, MYH14, and Rap1B) were identified based on multiple omics evidence, and they may influence the pathogenic mechanisms of osteoporosis by regulating the activities of monocytes.

INTRODUCTION

Osteoporosis is a serious public health problem, leading to significant mortality not only in aging females but also in males. Peripheral blood monocytes (PBMs) play important roles in bone metabolism by acting as precursors of osteoclasts and producing cytokines important for osteoclast development. The first cytosolic sub-proteome profiling analysis was performed in male PBMs to identify differentially expressed proteins (DEPs) that are associated with BMDs and risk of osteoporosis.

METHODS

Here, we conducted a comparative proteomics analysis in PBMs from Caucasian male subjects with discordant hip BMD (29 low BMD vs. 30 high BMD). To decrease the proteome complexity and expand the coverage range of the cellular proteome, we separated the PBM proteome into several subcellular compartments and focused on the cytosolic fractions, which are involved in a wide range of fundamental biochemical processes.

RESULTS

Of the total of 3796 detected cytosolic proteins, we identified 16 significant (P < 0.05) and an additional 22 suggestive (P < 0.1) DEPs between samples with low vs. high hip BMDs. Some of the genes for DEPs, including ALDOA, MYH14, and Rap1B, showed an association with BMD in multiple omics studies (proteomic, transcriptomic, and genomic). Further bioinformatics analysis revealed the enrichment of DEPs in functional terms for monocyte proliferation, differentiation, and migration.

CONCLUSIONS

The combination strategy of subcellular separation and proteome profiling allows an in-depth and refined investigation into the composition and functions of cytosolic proteome, which may shed light on the monocyte-mediated pathogenic mechanisms of osteoporosis.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

Polymorphisms in inflammation associated genes ALOX15 and IL-6 are associated with bone properties in young women and fracture in elderly

PURPOSE

ALOX12 and ALOX15 encode arachidonate lipoxygenases which produce lipid metabolites involved in inflammatory processes. Metabolites generated by ALOX12 and ALOX15 can activate the expression of the potent pro-inflammatory cytokine IL-6, and produce endogenous ligands for PPARG. In this study, polymorphisms in ALOX12, ALOX15, IL6 and PPARG were investigated for association with bone properties in young and elderly Swedish women.

METHODS

Three SNPs in ALOX12, five in ALOX15, one each in IL6 and PPARG were genotyped in the cohorts PEAK-25 (n=1061 women; all 25y) and OPRA (n=1044 women; all 75y). Bone mineral density (BMD) and quantitative ultrasound (QUS) were analyzed in both cohorts; trabecular bone score (TBS) in PEAK-25; bone loss, fracture incidence and serum C-reactive protein (CRP) were assessed in OPRA.

RESULTS

In the elderly women ALOX15 (rs2619112) was associated with CRP levels (p=0.004) and incident fracture of any type (p=0.014), although not with BMD or ultrasound. In young women, carrying the common T allele (ALOX 15 rs748694) was associated with lower QUS values (p=0.002-0.006). The IL6 SNP was associated with lower BMD in PEAK-25 (femoral neck p=0.034; hip p=0.012). TBS was not associated with variation in any gene. Variants in the ALOX12 and PPARγ were not associated with BMD in either cohort.

CONCLUSIONS

This study suggests that variation in inflammation related genes ALOX15 and IL6 was associated with bone microarchitecture and density in young adult women, but appears to be less important in the elderly, despite an observed association with CRP as a marker of inflammation and incident fracture.

Linkage and association analyses using families identified a locus affecting an osteoporosis-related trait

Osteoporosis is a common disorder characterized by low bone mass and microarchitectural deterioration of bone tissue, resulting in an increase in bone fragility and in susceptibility to fractures. The genetic basis of osteoporosis is complex and involves multiple genes and environmental factors. Here we introduce a family-based study of the genetics of osteoporosis - the Genetic Analysis of Osteoporosis (GAO) Project - to discover genetic variants affecting osteoporosis-related phenotypes. The GAO Project involved 11 extended families from Barcelona, Spain selected through a proband with osteoporosis (N=367). We performed spine, femur and whole body densitometry for all participants and also analyzed strength and geometrical properties of the hip. Our study focused on 23 densitometric phenotypes that we considered of high clinical relevance and four definitions of low bone mass and fracture status. Pedigree validation was carried out through microsatellite genotyping. The same microsatellites were used to interrogate our data (i) for the replication of previous linkage signals and (ii) for the potential discovery of new linkage signals. The linkage analysis identified one region marked by microsatellite D17S787 showing a strong and significant signal of linkage with femoral shaft cross-sectional moment of inertia (CSMI; LOD=3.18; p=6.5×10(-5)). The chromosomal location marked by microsatellite D17S787 includes several genes, among which two are of particular interest: COL1A1 and SOST, coding for collagen alpha-1 (I) chain and sclerostin, respectively. Follow-up association analysis resulted in only one significant result for rs4792909 from the SOST genomic region (p=0.00248). As a result, we provide strong and significant evidence from both linkage and association analyses that the SOST gene may affect the strength of the femoral shaft. Future investigations should study the relationship between bone mass formation and strength properties of the bones.

The transcription factor T-box 3 regulates colony-stimulating factor 1-dependent Jun dimerization protein 2 expression and plays an important role in osteoclastogenesis

Colony-stimulating factor 1 (CSF1) is known to promote osteoclast progenitor survival, but its roles in osteoclast differentiation and mature osteoclast function are less well understood. In a microarray screen, Jun dimerization protein 2 (JDP2) was identified as significantly induced by CSF1. Recent reports indicate that JDP2 is required for normal osteoclastogenesis and skeletal metabolism. Because there are no reports on the transcriptional regulation of this gene, the DNA sequence from -2612 to +682 bp (relative to the transcription start site) of the JDP2 gene was cloned, and promoter activity was analyzed. The T box-binding element (TBE) between -191 and -141 bp was identified as the cis-element responsible for CSF1-dependent JDP2 expression. Using degenerate PCR, Tbx3 was identified as the major isoform binding the TBE. Overexpression of Tbx3 induced JDP2 promoter activity, whereas suppressing Tbx3 expression substantially attenuated CSF1-induced transcription. Suppressing Tbx3 in osteoclast precursors reduced JDP2 expression and significantly impaired RANKL/CSF1-induced osteoclastogenesis. A MEK1/2-specific inhibitor was found to block CSF1-induced JDP2 expression. Consistent with these data, JDP2(-/-) mice were found to have increased bone mass. In summary, CSF1 up-regulates JDP2 expression by inducing Tbx3 binding to the JDP2 promoter. The downstream signaling cascade from activated c-Fms involves the MEK1/2-ERK1/2 pathway. Tbx3 plays an important role in osteoclastogenesis at least in part by regulating CSF1-dependent expression of JDP2.

Replication of Caucasian loci associated with bone mineral density in Koreans

SUMMARY

Most bone mineral density (BMD) loci were reported in Caucasian genome-wide association studies (GWAS). This study investigated the association between 59 known BMD loci (+200 suggestive SNPs) and DXA-derived BMD in East Asian population with respect to sex and site specificity. We also identified four novel BMD candidate loci from the suggestive SNPs.

INTRODUCTION

Most GWAS have reported BMD-related variations in Caucasian populations. This study investigates whether the BMD loci discovered in Caucasian GWAS are also associated with BMD in East Asian ethnic samples.

METHODS

A total of 2,729 unrelated Korean individuals from a population-based cohort were analyzed. We selected 747 single-nucleotide polymorphisms (SNPs). These markers included 547 SNPs from 59 loci with genome-wide significance (GWS, p value less than 5 × 10(-8)) levels and 200 suggestive SNPs that showed weaker BMD association with p value less than 5 × 10(-5). After quality control, 535 GWS SNPs and 182 suggestive SNPs were included in the replication analysis.

RESULTS

Of the 535 GWS SNPs, 276 from 25 loci were replicated (p < 0.05) in the Korean population with 51.6 % replication rate. Of the 182 suggestive variants, 16 were replicated (p < 0.05, 8.8 % of replication rate), and five reached a significant combined p value (less than 7.0 × 10(-5), 0.05/717 SNPs, corrected for multiple testing). Two markers (rs11711157, rs3732477) are for the same signal near the gene CPN2 (carboxypeptidase N, polypeptide 2). The other variants, rs6436440 and rs2291296, were located in the genes AP1S3 (adaptor-related protein complex 1, sigma 3 subunit) and RARB (retinoic acid receptor, beta).

CONCLUSION

Our results illustrate ethnic differences in BMD susceptibility genes and underscore the need for further genetic studies in each ethnic group. We were also able to replicate some SNPs with suggestive associations. These SNPs may be BMD-related genetic markers and should be further investigated.

The insulin-like growth factor system in bone: basic and clinical implications

The insulin-like growth factor (IGF) regulatory system is critical for skeletal growth and maintenance. Initially there was great hope that the recombinant IGFs might be used clinically for disorders ranging from short stature to fracture repair and osteoporosis. Although this potential was not realized, basic and translational studies have continued, providing significant insights into the role of this family of growth factors in skeletal homeostasis and the pathophysiology of several bone disorders. This article reviews the importance of the IGF regulatory system in skeletal growth and maintenance.

Inositol polyphosphate 4-phosphatase B as a regulator of bone mass in mice and humans

Osteoporosis is a multifactorial genetic disease characterized by reduction of bone mass due to dysregulation of osteoclast differentiation or maturation. Herein, we identified a regulator of osteoclastogenesis, the murine homolog of inositol polyphosphate 4-phosphatase type IIα (Inpp4bα). Expression of Inpp4bα is detected from early osteoclast differentiation to activation stage. Targeted expression of native Inpp4bα ex vivo repressed whereas phosphatase-inactive Inpp4bα stimulated osteoclast differentiation. Inpp4bα acts on intracellular calcium level that modulates NFATc1 nuclear translocation and activation. In vivo mice deficient in Inpp4b displayed increased osteoclast differentiation rate and potential resulting in decreased bone mass and osteoporosis. Importantly, INPP4B in human was identified as a susceptibility locus for osteoporosis. This study defined Inpp4b as a major modulator of the osteoclast differentiation and as a gene linked to variability of bone mineral density in mice and humans.

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.

BMD regulation on mouse distal chromosome 1, candidate genes, and response to ovariectomy or dietary fat

The distal end of mouse chromosome 1 (Chr 1) harbors quantitative trait loci (QTLs) that regulate bone mineral density (BMD) and share conserved synteny with human chromosome 1q. The objective of this article was to map this mouse distal Chr 1 region and identify gene(s) responsible for BMD regulation in females. We used X-ray densitometry [ie, dual-energy X-ray Absorptiometry (DXA), micro-computed tomography (µCT), and peripheral quantitative computed tomography (pQCT)] to phenotype a set of nested congenic strains constructed from C57BL/6BmJ (B6/Bm) and C3H/HeJ (C3H) mice to map the region associated with the BMD QTL. The critical region has been reduced to an interval of 0.152 Mb that contributes to increased BMD when C3H alleles are present. Histomorphometry and osteoblast cultures indicated that increased osteoblast activity was associated with increased BMD in mouse strains with C3H alleles in this critical region. This region contains two genes, Aim2, which binds with cytoplasmic dsDNA and results in apoptosis, and AC084073.22, a predicted gene of unknown function. Ovariectomy induced bone loss in the B6/Bm progenitor and the three congenic strains regardless of the alleles present in the critical BMD region. High dietary fat treatment (thought to suppress distal Chr 1 QTL for BMD in mice) did not induce bone loss in the congenics carrying C3H alleles in the critical 0.152 Mb carrying the AIM2 and AC084073.22 genes. Gene expression studies in whole bone of key congenics showed differential expression of AC084073.22 for strains carrying B6/Bm versus C3H alleles but not for Aim2. In conclusion, our data suggest that osteoblasts are the cellular target of gene action and that AC084073.22 is the best candidate for female BMD regulation in the distal region of mouse Chr 1.

A genome-wide association analysis implicates SOX6 as a candidate gene for wrist bone mass

Osteoporosis is a highly heritable common bone disease leading to fractures that severely impair the life quality of patients. Wrist fractures caused by osteoporosis are largely due to the scarcity of wrist bone mass. Here we report the results of a genome-wide association study (GWAS) of wrist bone mineral density (BMD). We examined ∼500000 SNP markers in 1000 unrelated homogeneous Caucasian subjects and found a novel allelic association with wrist BMD at rs11023787 in the SOX6 (SRY (sex determining region Y)-box 6) gene (P=9.00×10(-5)). Subjects carrying the C allele of rs11023787 in SOX6 had significantly higher mean wrist BMD values than those with the T allele (0.485:0.462 g cm(-2) for C allele vs. T allele carriers). For validation, we performed SOX6 association for BMD in an independent Chinese sample and found that SNP rs11023787 was significantly associated with wrist BMD in the Chinese sample (P=6.41×10(-3)). Meta-analyses of the GWAS scan and the replication studies yielded P-values of 5.20×10(-6) for rs11023787. Results of this study, together with the functional relevance of SOX6 in cartilage formation, support the SOX6 gene as an important gene for BMD variation.

Clinical review 2: Genetic determinants of bone density and fracture risk--state of the art and future directions

CONTEXT

Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved.

EVIDENCE ACQUISITION AND SYNTHESIS

In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward.

CONCLUSIONS

Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Genes influencing spinal bone mineral density in inbred F344, LEW, COP, and DA rats

Previously, we identified the regions of chromosomes 10q12-q31 and 15p16-q21 harbor quantitative trait loci (QTLs) for lumbar volumetric bone mineral density (vBMD) in female F2 rats derived from Fischer 344 (F344) x Lewis (LEW) and Copenhagen 2331 (COP) x Dark Agouti (DA) crosses. The purpose of this study is to identify the candidate genes within these QTL regions contributing to the variation in lumbar vBMD. RNA was extracted from bone tissue of F344, LEW, COP, and DA rats. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 Arrays. Genes differentially expressed among the rat strains were then ranked based on the strength of the correlation with lumbar vBMD in F2 animals derived from these rats. Quantitative PCR (qPCR) analysis was performed to confirm the prioritized candidate genes. A total of 285 genes were differentially expressed among all strains of rats with a false discovery rate less than 10%. Among these genes, 18 candidate genes were prioritized based on their strong correlation (r (2) > 0.90) with lumbar vBMD. Of these, 14 genes (Akap1, Asgr2, Esd, Fam101b, Irf1, Lcp1, Ltc4s, Mdp-1, Pdhb, Plxdc1, Rabep1, Rhot1, Slc2a4, Xpo4) were confirmed by qPCR. We identified several novel candidate genes influencing spinal vBMD in rats.

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.

Differentially expressed genes strongly correlated with femur strength in rats

The region of chromosome 1q33-q54 harbors quantitative trait loci (QTL) for femur strength in COPxDA and F344xLEW F2 rats. The purpose of this study is to identify the genes within this QTL region that contribute to the variation in femur strength. Microarray analysis was performed using RNA extracted from femurs of COP, DA, F344 and LEW rats. Genes differentially expressed in the 1q33-q54 region among these rat strains were then ranked based on the strength of correlation with femur strength in F2 animals derived from these rats. A total of 214 genes in this QTL region were differentially expressed among all rat strains, and 81 genes were found to be strongly correlated (r(2)>0.50) with femur strength. Of these, 12 candidate genes were prioritized for further validation, and 8 of these genes (Ifit3, Ppp2r5b, Irf7, Mpeg1, Bloc1s2, Pycard, Sec23ip, and Hps6) were confirmed by quantitative PCR (qPCR). Ingenuity Pathway Analysis suggested that these genes were involved in interferon alpha, nuclear factor-kappa B (NFkB), extracellular signal-related kinase (ERK), hepatocyte nuclear factor 4 alpha (HNF4A) and tumor necrosis factor (TNF) pathways.

Bone mineral density variation in men is influenced by sex-specific and non sex-specific quantitative trait loci

INTRODUCTION

A major predictor of age-related osteoporotic fracture is peak areal bone mineral density (aBMD) which is a highly heritable trait. However, few linkage and association studies have been performed in men to identify the genes contributing to normal variation in aBMD. The aim of this study was to perform a genome wide scan in healthy men to identify quantitative trait loci (QTL) that were significantly linked to aBMD and to test whether any of these might be sex-specific.

METHODS

aBMD at the spine and hip were measured in 515 pairs of brothers, aged 18-61 (405 white pairs, 110 black pairs). Linkage analysis in the brother sample was compared with results in a previously published sample of 774 sister pairs to identify sex-specific quantitative trait loci (QTL).

RESULTS

A genome wide scan identified significant QTL (LOD>3.6) for aBMD on chromosomes 4q21 (hip), 7q34 (spine), 14q32 (hip), 19p13 (hip), 21q21 (hip), and 22q13 (hip). Analysis suggested that the QTL on chromosomes 7q34, 14q32, and 21q21 were male-specific whereas the others were not sex-specific.

CONCLUSIONS

This study demonstrates that six QTL were significantly linked with aBMD in men. One was linked to the spine and five were linked to the hip. When compared to published data in women from the same geographical region, the QTL on chromosomes 7, 14 and 21 were male-specific. The occurrence of sex-specific genes in humans for aBMD has important implications for the pathogenesis and treatment of osteoporosis.

Further genetic evidence suggesting a role for the RhoGTPase-RhoGEF pathway in osteoporosis

Osteoporosis is a highly heritable trait that appears to be influenced by multiple genes. Genome-wide linkage studies have highlighted the chromosomal region 3p14-p21 as a quantitative trait locus for BMD. We have previously published evidence suggesting that the ARHGEF3 gene from this region is associated with BMD in women. The product of this gene activates the RHOA GTPase, the gene for which is also located within this region. The aim of this study was to evaluate the influence of genetic polymorphism in RHOA on bone density in women. Sequence variation within the RHOA gene region was determined using 9 single nucleotide polymorphisms (SNPs) in a discovery cohort of 769 female sibs. Of the 9 SNPs, one was found to be monomorphic with the others representing 3 distinct linkage disequilibrium (LD) blocks. Using FBAT software, significant associations were found between two of these LD blocks and BMD Z-score of the spine and hip (P=0.001-0.036). The LD block tagged by the SNP rs17595772 showed maximal association, with the more common G allele at rs17595772 associated with decreased BMD Z-score. Genotyping for rs17595772 in a replication cohort of 780 postmenopausal women confirmed an association with BMD Z-score (P=0.002-0.036). Again, the G allele was found to be associated with a reduced hip and spine BMD Z-score. These results support the implication of the RhoGTPase-RhoGEF pathway in osteoporosis, and suggest that one or more genes in this pathway may be responsible for the linkage observed between 3p14-p21 and BMD.

Insulin-like growth factor-I and bone: lessons from mice and men

Studies of humans and animals have illustrated a strong association between insulin-like growth factor (IGF)-I and skeletal acquisition. However, the precise molecular and cellular mechanisms underlying this effect still largely remain unknown. Recent advances in molecular and genetic techniques for in vivo studies provide excellent tools for us to explore how circulating and skeletal insulin-like growth factor-I (IGF-I) may affect not only peak bone mass but also bone loss. This review highlights recent findings that shed new light on the interaction of the IGF-I signaling pathway with other skeletal networks, and the role of IGF-I in the bone marrow milieu.

The human calcium-sensing receptor and interleukin-6 genes are associated with bone mineral density in Chinese

Calcium sensing receptor (CASR) is a central factor involved in calcium metabolism. Interleukin-6 (IL-6) is a pleiotropic cytokine that plays an important role in osteoclast differentiation. Thus, both CASR and IL-6 are important in bone and mineral metabolism and are prominent candidate genes for osteoporosis. The study aimed to test association and/or linkage between the CASR and IL-6 genes with bone mineral density (BMD) variation in a Chinese population. A cytosine-adenine (CA)n repeat polymorphism in the CASR gene and the IL-6 gene was genotyped, respectively, in 1,263 subjects from 402 Chinese nuclear families. Employing tests implemented in the program QTDT (quantitative transmission disequilibrium tests), a significant total association of the CASR (CA)12 allele (P = 0.006) and (CA)18 allele (P = 0.02) with BMD at the femoral neck was found. For the IL-6 gene, significant within-family associations were found between the (CA)14 allele and BMD at the total hip (P = 0.021), the femoral neck (P = 0.041), and the intertrochanteric region (P = 0.029). A significant linkage was also observed between IL-6 CA repeat polymorphism and BMD at the spine (P = 0.001). The results suggest that the CASR gene and the IL-6 gene may have effects on BMD variation in Chinese.

Comparison of whole genome linkage scans in premenopausal and postmenopausal women: no bone-loss-specific QTLs were implicated

This study was conducted to investigate if there exist bone-loss-specific quantitative trait loci (QTLs) for females. Genome-wide linkage scans were conducted in total, premenopausal, and postmenopausal women, respectively. No QTLs exclusively were found in postmenopausal women, suggesting that no bone-loss-specific QTL was implicated independent of BMD in our sample.

INTRODUCTION

Bone mineral density (BMD) in elderly women is determined jointly by peak bone mass achieved before menopause and by subsequent bone loss upon and after menopause. Peak bone mass is under strong genetic control, but whether bone loss has genetic determination independent of peak BMD is unknown.

MATERIALS AND METHODS

To investigate if there exist bone-loss-specific quantitative trait loci (QTLs) for females, we conducted genome-wide linkage scans in 2,582 Caucasian females from 451 pedigrees including 1,486 premenopausal and 1,096 postmenopausal women. Linkage analyses were performed in the total sample and premenopausal and postmenopausal women subgroups, respectively, and the results were compared.

RESULTS

No linkage evidence was found exclusively in postmenopausal women. Linkage signals identified are largely consistent in the total, premenopausal, and postmenopausal samples. For example, for spine BMD, for the total sample, a significant linkage was obtained on 15q13 (LOD = 3.67), and LOD scores of 1.52 and 2.49 were achieved on 15q13 in premenopausal and postmenopausal women, respectively.

CONCLUSIONS

We did not find any QTLs exclusively in postmenopausal women; hence, no specific QTL for bone loss was implicated independent of BMD in our female sample.

Quantitative trait locus on chromosome 1q influences bone loss in young Mexican American adults

Bone loss occurs as early as the third decade and its cumulative effect throughout adulthood may impact risk for osteoporosis in later life, however, the genes and environmental factors influencing early bone loss are largely unknown. We investigated the role of genes in the change in bone mineral density (BMD) in participants in the San Antonio Family Osteoporosis Study. BMD change in 327 Mexican Americans (ages 25-45 years) from 32 extended pedigrees was calculated from DXA measurements at baseline and follow-up (3.5 to 8.9 years later). Family-based likelihood methods were used to estimate heritability (h(2)) and perform autosome-wide linkage analysis for BMD change of the proximal femur and forearm and to estimate heritability for BMD change of lumbar spine. BMD change was significantly heritable for total hip, ultradistal radius, and 33% radius (h(2) = 0.34, 0.34, and 0.27, respectively; p < 0.03 for all), modestly heritable for femoral neck (h(2) = 0.22; p = 0.06) and not heritable for spine BMD. Covariates associated with BMD change included age, sex, baseline BMD, menopause, body mass index, and interim BMI change, and accounted for 6% to 24% of phenotype variation. A significant quantitative trait locus (LOD = 3.6) for femoral neck BMD change was observed on chromosome 1q23. In conclusion, we observed that change in BMD in young adults is heritable and performed one of the first linkage studies for BMD change. Linkage to chromosome 1q23 suggests that this region may harbor one or more genes involved in regulating early BMD change of the femoral neck.

Genome-wide copy-number-variation study identified a susceptibility gene, UGT2B17, for osteoporosis

Osteoporosis, a highly heritable disease, is characterized mainly by low bone-mineral density (BMD), poor bone geometry, and/or osteoporotic fractures (OF). Copy-number variation (CNV) has been shown to be associated with complex human diseases. The contribution of CNV to osteoporosis has not been determined yet. We conducted case-control genome-wide CNV analyses, using the Affymetrix 500K Array Set, in 700 elderly Chinese individuals comprising 350 cases with homogeneous hip OF and 350 matched controls. We constructed a genomic map containing 727 CNV regions in Chinese individuals. We found that CNV 4q13.2 was strongly associated with OF (p = 2.0 x 10(-4), Bonferroni-corrected p = 0.02, odds ratio = 1.73). Validation experiments using PCR and electrophoresis, as well as real-time PCR, further identified a deletion variant of UGT2B17 in CNV 4q13.2. Importantly, the association between CNV of UGT2B17 and OF was successfully replicated in an independent Chinese sample containing 399 cases with hip OF and 400 controls. We further examined this CNV's relevance to major risk factors for OF (i.e., hip BMD and femoral-neck bone geometry) in both Chinese (689 subjects) and white (1000 subjects) samples and found consistently significant results (p = 5.0 x 10(-4) -0.021). Because UGT2B17 encodes an enzyme catabolizing steroid hormones, we measured the concentrations of serum testosterone and estradiol for 236 young Chinese males and assessed their UGT2B17 copy number. Subjects without UGT2B17 had significantly higher concentrations of testosterone and estradiol. Our findings suggest the important contribution of CNV of UGT2B17 to the pathogenesis of osteoporosis.

A genome-wide scan for quantitative trait loci affecting limb bone lengths and areal bone mineral density of the distal femur in a White Duroc x Erhualian F2 population

Background

Limb bone lengths and bone mineral density (BMD) have been used to assess the bone growth and the risk of bone fractures in pigs, respectively. It has been suggested that limb bone lengths and BMD are under genetic control. However, the knowledge about the genetic basis of the limb bone lengths and mineralisatinon is limited in pigs. The aim of this study was to identify quantitative trait loci (QTL) affecting limb bone lengths and BMD of the distal femur in a White Duroc × Erhualian resource population.

Results

Limb bone lengths and femoral bone mineral density (fBMD) were measured in a total of 1021 and 116 F₂ animals, respectively. There were strong positive correlations among the lengths of limb bones and medium positive correlations between the lengths of limb bones and fBMD. A whole-genome scan involving 183 microsatellite markers across the pig genome revealed 35 QTL for the limb bone lengths and 2 for femoral BMD. The most significant QTL for the lengths of five limb bones were mapped on two chromosomes affecting all 5 limb bones traits. One was detected around 57 cM on pig chromosome (SSC) 7 with the largest F-value of more than 26 and 95% confidence intervals of less than 5 cM, providing a crucial start point to identify the causal genes for these traits. The Erhualian alleles were associated with longer limb bones. The other was located on SSCX with a peak at 50–53 cM, whereas alleles from the White Duroc breed increased the bone length. Many QTL identified are homologous to the human genomic regions containing QTL for bone-related traits and a list of interesting candidate genes.

Conclusion

This study detected the QTL for the lengths of scapula, ulna, humerus and tibia and fBMD in the pig for the first time. Moreover, several new QTL for the pig femoral length were found. As correlated traits, QTL for the lengths of five limb bones were mainly located in the same genomic regions. The most promising QTL for the lengths of five limb bones on SSC7 merits further investigation.

Genetic studies in osteoporosis--the end of the beginning

Osteoporosis and disorders of bone fragility are highly heritable, but despite much effort the identities of few of the genes involved has been established. Recent developments in genetics such as genome-wide association studies are revolutionizing research in this field, and it is likely that further contributions will be made through application of next-generation sequencing technologies, analysis of copy number variation polymorphisms, and high-throughput mouse mutagenesis programs. This article outlines what we know about osteoporosis genetics to date and the probable future directions of research in this field.

Linkage screen for BMD phenotypes in male and female COP and DA rat strains

Because particular inbred strains of experimental animals are informative for only a subset of the genes underlying variability in BMD, we undertook a genome screen to identify quantitative trait loci (QTLs) in 828 F(2) progeny (405 males and 423 females) derived from the Copenhagen 2331 (COP) and dark agouti (DA) strains of rats. This screen was performed to complement our study in female Fischer 344 (F344) and Lewis (LEW) rats and to further delineate the factors underlying the complex genetic architecture of BMD in the rat model. Microsatellite genotyping was performed using markers at an average density of 20 cM. BMD was measured by pQCT and DXA. These data were analyzed in the R/qtl software to detect QTLs acting in both sexes as well as those having sex-specific effects. A QTL was detected in both sexes on chromosome 18 for midfemur volumetric BMD (vBMD; genome-wide, p < 0.01). On distal chromosome 1, a QTL was found for femur and vertebral aBMD as well as distal femur vBMD, and this QTL appears distinct from the proximal chromosome 1 QTL impacting BMD in our F344/LEW cross. Additional aBMD and vBMD QTLs and several sex-specific QTLs were also detected. These included a male-specific QTL (p < 0.01) on chromosome 8 and a female-specific QTL on chromosomes 7 and 14 (p < 0.01). Few of the QTLs identified showed overlap with the significant QTLs from the F344/LEW cross. These results confirm that the genetic influence on BMD in the rat model is quite complex and would seem to be influenced by a number of different genes, some of which have sex-specific effects.

PPARG by dietary fat interaction influences bone mass in mice and humans

Adult BMD, an important risk factor for fracture, is the result of genetic and environmental interactions. A quantitative trait locus (QTL) for the phenotype of volumetric BMD (vBMD), named Bmd8, was found on mid-distal chromosome (Chr) 6 in mice. This region is homologous to human Chr 3p25. The B6.C3H-6T (6T) congenic mouse was previously created to study this QTL. Using block haplotyping of the 6T congenic region, expression analysis in the mouse, and examination of nonsynonymous SNPs, peroxisome proliferator activated receptor gamma (Pparg) was determined to be the most likely candidate gene for the Bmd8 QTL of the 630 genes located in the congenic region. Furthermore, in the C3H/HeJ (C3H) strain, which is the donor strain for the 6T congenic, several polymorphisms were found in the Pparg gene. On challenge with a high-fat diet, we found that the 6T mouse has a lower areal BMD (aBMD) and volume fraction of trabecular bone (BV/TV%) of the distal femur compared with B6 mice. Interactions between SNPs in the PPARG gene and dietary fat for the phenotype of BMD were examined in the Framingham Offspring Cohort. This analysis showed that there was a similar interaction of the PPARG gene and diet (fat intake) on aBMD in both men and women. These findings suggest that dietary fat has a significant influence on BMD that is dependent on the alleles present for the PPARG gene.

Genetic variation in candidate osteoporosis genes, bone mineral density, and fracture risk: the study of osteoporotic fractures

Candidate osteoporosis gene variants were examined for associations with fracture risk and bone mineral density (BMD). A total of 9,704 white women were recruited at four U.S. clinical centers and enrolled into the Study of Osteoporotic Fractures, a longitudinal cohort study. Genotyping of 31 polymorphisms from 18 candidate osteoporosis genes was performed in 6,752 women. Incident radiographic fractures were identified at the third and eighth examinations compared with the baseline examination. BMD was measured at the total hip by dual-energy X-ray absorptiometry. Analyses were adjusted for age, clinic site, and self-reported ethnicity. During a mean follow-up of 14.5 years, a total of 849 hip, 658 vertebral, and 2,496 nonhip/nonvertebral fractures occurred in 6,752 women. Women carrying the ALOX15_G48924T T/T genotype had a higher rate of hip fracture (hazard ratio [HR] = 1.33;95% confidence interval [95% CI] = 1.00-1.77) compared with the G/G genotype. Compared with those carrying the PRL_T228C T/T genotype, women with either the C/C (HR = 0.80; 95% CI = 0.67-0.95) or C/T (HR = 0.81; 95% CI = 0.68-0.97) genotype had a lower rate of nonvertebral/nonhip fractures. Women carrying the BMP2_A125611G G/G genotype had a higher rate of vertebral fracture (odds ratio [OR] = 1.51; 95% CI = 1.03-2.23) compared with the A/A genotype. Women with the ESR1_C1335G G/G genotype had a higher rate of vertebral fracture (OR = 1.64; 95% CI = 1.07-2.50) compared with the C/C genotype. Compared with those with the MMP2_C595T C/C genotype, women with the C/T (OR = 0.79; 95% CI = 0.65-0.96) or T/T (OR = 0.44; 95% CI = 0.27-0.72) genotype had a lower rate of vertebral fracture. In conclusion, polymorphisms in several candidate genes were associated with hip, vertebral, and nonhip/nonvertebral fractures but not with total hip BMD in this large population based cohort study.

Identification of a role for the ARHGEF3 gene in postmenopausal osteoporosis

Osteoporosis is a common and debilitating bone disease characterized by low bone mineral density (BMD), a highly heritable and polygenic trait. Genome-wide linkage studies have identified 3p14-p21 as a quantitative trait locus for BMD. The ARHGEF3 gene is situated within this region and was identified as a strong positional candidate. The aim of this study was to evaluate the role of variation in ARHGEF3 on bone density in women. Sequence variation within ARHGEF3 was analyzed with 17 single-nucleotide polymorphisms (SNPs) in a discovery cohort of 769 female sibs. Significant associations were found with family-based association tests between five SNPs and various measures of age-adjusted BMD (p = 0.0007-0.041) with rs7646054 showing maximal association. Analysis of the data with QPDTPHASE suggested that the more common G allele at rs7646054 is associated with decreased age-adjusted BMD. Significant associations were also demonstrated between 3-SNP haplotypes and age-adjusted spine and femoral-neck BMD (p = 0.002 and 0.003, respectively). rs7646054 was then genotyped in a replication cohort, and significant associations with hip and spine BMD were confirmed (p = 0.003-0.038), as well as an association with fracture rate (p = 0.02). Again, the G allele was associated with a decrease in age-adjusted BMD at each site studied. In conclusion, genetic variation in ARHGEF3 plays a role in the determination of bone density in Caucasian women. This data implicates the RhoGTPase-RhoGEF pathway in osteoporosis.

Sex-specific genetic loci for femoral neck bone mass and strength identified in inbred COP and DA rats

INTRODUCTION

Hip fracture is the most devastating osteoporotic fracture type with significant morbidity and mortality. Several studies in humans identified chromosomal regions linked to hip size and bone mass. Animal models, particularly the inbred rat, serve as complementary approaches for studying the genetic influence on hip fragility. The purpose of this study is to identify sex-independent and sex-specific quantitative trait loci (QTLs) for femoral neck density, structure, and strength in inbred Copenhagen 2331 (COP) and Dark Agouti (DA) rats.

MATERIALS AND METHODS

A total of 828 (405 males and 423 females) F(2) progeny derived from the inbred COP and DA strains of rats were phenotyped for femoral neck volumetric BMD (vBMD), cross-sectional area, polar moment of inertia (Ip), neck width, ultimate force, and energy to break. A whole genome screen was performed using 93 microsatellite markers with an average intermarker distance of 20 cM. Recombination-based marker maps were generated using MAPMAKER/EXP from the COP x DA F(2) data and compared with published Rat Genome Database (RGD) maps. These maps were used for genome-wide linkage analyses to detect sex-independent and sex-specific QTLs.

RESULTS

Significant evidence of linkage (p < 0.01) for sex-independent QTLs were detected for (1) femoral neck vBMD on chromosomes (Chrs) 1, 6, 10, and 12, (2) femoral neck structure on Chrs 5, 7, 10, and 18, and (3) biomechanical properties on Chrs 1 and 4. Male-specific QTLs were discovered on Chrs 2, 9, and 18 for total vBMD, on Chr 17 for trabecular vBMD, on Chr 9 for total bone area, and on Chr 15 for ultimate force. A female-specific QTL was discovered on Chr 2 for ultimate force. The effect size of the individual QTL varied between 1% and 4%.

CONCLUSIONS

We detected evidence that sex-independent and sex-specific QTLs contribute to hip fragility in the inbred rat. Several QTLs regions identified in this study are homologous to human chromosomal regions previously linked to QTLs contributing to femoral neck and related phenotypes.

How is mRNA expression predictive for protein expression? A correlation study on human circulating monocytes

A key assumption in studying mRNA expression is that it is informative in the prediction of protein expression. However, only limited studies have explored the mRNA-protein expression correlation in yeast or human tissues and the results have been relatively inconsistent. We carried out correlation analyses on mRNA-protein expressions in freshly isolated human circulating monocytes from 30 unrelated women. The expressed proteins for 71 genes were quantified and identified by 2-D electrophoresis coupled with mass spectrometry. The corresponding mRNA expressions were quantified by Affymetrix gene chips. Significant correlation (r=0.235, P<0.0001) was observed for the whole dataset including all studied genes and all samples. The correlations varied in different biological categories of gene ontology. For example, the highest correlation was achieved for genes of the extracellular region in terms of cellular component (r=0.643, P<0.0001) and the lowest correlation was obtained for genes of regulation (r=0.099, P=0.213) in terms of biological process. In the genome, half of the samples showed significant positive correlation for the 71 genes and significant correlation was found between the average mRNA and the average protein expression levels in all samples (r=0.296, P<0.01). However, at the study group level, only five studied genes had significant positive correlation across all the samples. Our results showed an overall positive correlation between mRNA and protein expression levels. However, the moderate and varied correlations suggest that mRNA expression might be sometimes useful, but certainly far from perfect, in predicting protein expression levels.

Mutations in the insulin-like factor 3 receptor are associated with osteoporosis

INTRODUCTION

Insulin-like factor 3 (INSL3) is produced primarily by testicular Leydig cells. It acts by binding to its specific G protein-coupled receptor RXFP2 (relaxin family peptide 2) and is involved in testicular descent during fetal development. The physiological role of INSL3 in adults is not known, although substantial INSL3 circulating levels are present. The aim of this study was to verify whether reduced INSL3 activity could cause or contribute to some signs of hypogonadism, such as reduced BMD, currently attributed to testosterone deficiency.

MATERIALS AND METHODS

Extensive clinical, biochemical, and hormonal study, including bone densitometry by DXA, was performed on 25 young men (age, 27-41 yr) with the well-characterized T222P mutation in the RXFP2 gene. Expression analysis of INSL3 and RXFP2 on human bone biopsy and human and mouse osteoblast cell cultures was performed by RT-PCR, quantitative RT-PCR, and immunohistochemistry. Real-time cAMP imaging analysis and proliferation assay under the stimulus of INSL3 was performed on these cells. Lumbar spine and femoral bone of Rxfp2-deficient mice were studied by static and dynamic histomorphometry and muCT, respectively.

RESULTS

Sixteen of 25 (64%) young men with RXFP2 mutations had significantly reduced BMD. No other apparent cause of osteoporosis was evident in these subjects, whose testosterone levels and gonadal function were normal. Expression analyses showed the presence of RXFP2 in human and mouse osteoblasts. Stimulation of these cells with INSL3 produced a dose- and time-dependent increase in cAMP and cell proliferation, confirming the functionality of the RXFP2/INSL3 receptor-ligand complex. Consistent with the human phenotype, bone histomorphometric and muCT analyses of Rxfp2(-/-) mice showed decreased bone mass, mineralizing surface, bone formation, and osteoclast surface compared with wildtype littermates.

CONCLUSIONS

This study suggests for the first time a role for INSL3/RXFP2 signaling in bone metabolism and links RXFP2 gene mutations with human osteoporosis.

Chromosomal regions 22q13 and 3p25 may harbor quantitative trait loci influencing both age at menarche and bone mineral density

Late age at menarche (AAM), an important type of endocrinopathy in females, is associated with lower bone mineral density (BMD), a major risk factor for osteoporosis. The correlation is mainly mediated through common genetic factors, which are largely unknown. A bivariate genome-wide linkage scan was conducted on 2,522 females from 414 Caucasian pedigrees to identify quantitative trait loci influencing both AAM and BMD. The strongest linkage signal was detected on chromosome 22q13. Other regions such as the 3q13, 3p25, 7p15, and 15q13 were also suggested. The inferred promising candidate genes in the linkage regions may contribute to our understanding of pathogenesis of endocrinopathy and osteoporosis in females.

A whole genome linkage scan for QTLs underlying peak bone mineral density

We conducted a whole genome linkage scan for quantitative trait loci (QTLs) underlying peak bone mineral density (PBMD). Our efforts identified several potential genomic regions for PBMD and highlighted the importance of epistatic interaction and sex-specific analyses in identifying genetic regions underlying PBMD variation.

INTRODUCTION

Peak bone mineral density (PBMD) is an important clinical risk predictor of osteoporosis and explains a large part of bone mineral density (BMD) variation.

METHODS

To detect susceptive quantitative trait loci (QTLs) for PBMD variation including consideration of epistatic and sex-specific effects, we conducted a whole genome linkage scan (WGLS) for PBMD using 2,200 Caucasians from 207 pedigrees, aged 20-50 years. All the individuals were genotyped with 410 microsatellite markers. In addition to WGLS in the total combined sample of males and females, we conducted epistatic interaction analyses, and sex-specific subgroup linkage analyses.

RESULTS

We identified several potential genomic regions that met the criteria for suggestive linkage. The most impressing region is 12p12 for hip PBMD (LOD = 2.79) in the total sample. Epistatic interaction analyses found a significant epistatic interaction between 12p12 and 22q13 (p = 0.0021) for hip PBMD. Additionally, we detected suggestive linkage evidence at 15q26 (LOD = 2.93), 2p13 (LOD = 2.64), and Xq27 (LOD = 2.64). Sex-specific analyses suggested the presence of sex-specific QTLs for PBMD variation.

CONCLUSIONS

Our efforts identified several potential regions for PBMD and highlighted the importance of epistatic interaction and sex-specific analyses in identifying genetic regions underlying PBMD variation.

Genetic loci affecting bone structure and strength in inbred COP and DA rats

Previous studies have shown that the Copenhagen 2331 (COP) and Dark Agouti (DA) rats have significant differences in bone structure and strength despite their similar body mass. Thus, these inbred rat strains may provide a unique resource to identify the genetics underlying the phenotypic variation in bone fragility. A sample of 828 (405 males and 423 females) COPxDA F2 progeny had extensive phenotyping for bone structure measures including cortical bone area and polar moment of inertia at the femur midshaft and total, cortical and trabecular bone areas, for the lumbar vertebra 5 (L5). Bone strength phenotypes included ultimate force, stiffness and work to failure of femur and L5. These skeletal phenotypes were measured using peripheral quantitative computed tomography (pQCT) and mechanical testing. A whole-genome screen was conducted in the F2 rats, using microsatellite markers spaced at approximately 20 cM intervals. Genetic marker maps were generated from the F2 data and used for genome-wide linkage analyses to detect linkage to the bone structure and strength phenotypes. Permutation testing was employed to obtain the thresholds for genome-wide significance (p<0.01). Significant QTL for femur structure and strength were identified on chromosome (Chr) 1 with a maximum LOD score of 33.5; evidence of linkage was found in both the male and female rats. In addition, Chrs 6, 7, 10, 13, 15 and 18 were linked to femur midshaft structure. QTL linked to femur strength were identified on Chrs 5 and 10. For L5 vertebrae, Chrs 2, 16, and 18 harbored QTL for cortical structure and trabecular structure for L5 was linked to Chrs 1, 7, 12, and 18. One female-specific QTL for femur ultimate force was identified on Chr 5, and two male-specific QTL for L5 cortical area were found on Chrs 2 and 18. Our study demonstrates strong evidence of linkage for bone structure and strength to multiple rat chromosomes.

Polymorphisms in the low-density lipoprotein receptor-related protein 5 (LRP5) gene are associated with peak bone mass in non-sedentary men: results from the Odense androgen study

PURPOSE

To investigate the impact of the Ala1330Val (rs3736228, exon 18) and Val667Met (rs4988321, exon 9) polymorphisms of the low-density lipoprotein receptor-related protein 5 (LRP5) gene on peak bone mass in young men.

METHODS

The Odense Androgen Study (OAS) is a population-based study comprising 783 Caucasian men aged 20-30 years. Genotyping was performed using real-time polymerase chain reaction (PCR) or fluorescence polarization. Bone mineral density (BMD) measurements were performed using dual-energy X-ray absorptiometry.

RESULTS

The CC, CT, and TT genotypes in Ala1330Val were found in 75.6%, 21.8%, and 2.6% of the participants, respectively. Similarly, the GG, GA, and AA genotypes of Val667Met were found in 89.7%, 9.8%, and 0.5%, respectively. For the Ala1330Val polymorphism, no significant differences between the genotypes were found regarding BMD in the overall study population. However, when analysis was restricted to non-sedentary men (n = 589), a significant association between the number of T-alleles and BMD in the spine and whole body were found. Each copy of the T-allele changed the Z-score of the spine by (median and 95% confidence interval) -0.21 [95% CI: -0.40; -0.03] (p < 0.02). Analysis suggested an association between the AA genotype in the Val667Met polymorphism and increased body height and decreased BMD of the femoral neck; however, no significant gene-dose effect of the A-allele could be demonstrated in the whole population. When the analysis was restricted to non-sedentary subjects, however, each number of A-alleles was associated with a change in Z-score of -0.26 [95% CI: -0.51; -0.01] (p = 0.04). No further significant results emerged with haplotype analysis.

CONCLUSION

The Ala1330Val and Val667Met polymorphisms in the LRP5 gene are significantly associated with peak bone mass in physically active men.

Chromosome 2q32 may harbor a QTL affecting BMD variation at different skeletal sites

BMDs at different skeletal sites share some common genetic determinants. Using PCA and bivariate linkage analysis, we identified a QTL on chromosome 2q32 with significant pleiotropic effects on BMDs at different skeletal sites.

INTRODUCTION

BMDs at the hip, spine, and forearm are genetically correlated, suggesting the existence of quantitative trait loci (QTLs) with concurrent effects on BMDs at these three skeletal sites. Consequently, it is important to identify these QTLs in the human genome and, for those implicated QTLs, it is important to differentiate between pleiotropic effects, caused by a single gene that concurrently effects these traits, and co-incident linkage, caused by multiple, closely linked, genes that independently effect these traits.

MATERIALS AND METHODS

For a sample of 451 American white pedigrees made up of 4,498 individuals, we evaluated the correlations between BMDs at the three skeletal sites. We carried out principal component analysis (PCA) for the three correlated traits and obtained a major component, PC1, which accounts for >75% of the co-variation of BMDs at the three sites. We subsequently conducted a whole genome linkage scan for PC1 and performed bivariate linkage analysis for pairs of the three traits (i.e., forearm/spine BMD, hip/forearm BMD, and hip/spine BMD).

RESULTS

Chromosome region 2q32, near the marker GATA65C03M, showed strong linkage to PC1 (LOD = 3.35). Subsequent bivariate linkage analysis substantiated linkage at 2q32 for each trait pair (LOD scores were 2.65, 2.42, and 2.13 for forearm/spine BMD, hip/forearm BMD, and hip/spine BMD, respectively). Further analyses rejected the hypothesis of co-incident linkage (p(0)[forearm/spine] = 0.0005, p(0)[hip/forearm] = 0.004, p(0)(hip/spine] = 0.001) but failed to reject the hypothesis of pleiotropy (p(1)[forearm/spine] = 0.35, p(1)[hip/forearm] = 0.07, p(1)[hip/spine] = 0.15).

CONCLUSIONS

Our results strongly support the conclusion that chromosome region 2q32 may harbor a QTL with pleiotropic effects on BMDs at different skeletal sites.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Genetic dissection of mouse distal chromosome 1 reveals three linked BMD QTLs with sex-dependent regulation of bone phenotypes

Genetic analyses with mouse congenic strains for distal Chr1 have identified three closely linked QTLs regulating femoral vBMD, mid-diaphyseal cortical thickness, and trabecular microstructure in a sex-dependent fashion. The homologous relationship between distal mouse Chr 1 and human 1q21-24 offers the possibility of finding common regulatory genes for cortical and trabecular bone.

INTRODUCTION

The distal third of mouse chromosome 1 (Chr 1) has been shown to carry a major quantitative trait locus (QTL) for BMD from several inbred mouse strain crosses. Genetic and functional analyses are essential to identify genes and cellular mechanisms for acquisition of peak bone mass.

MATERIALS AND METHODS

Nested congenic sublines of mice were developed with a C57BL/6J (B6) background carrying <1- to 9-Mbp-sized segments donated from C3H/HeJ (C3H). Isolated femurs from 16-wk-old female and male mice were measured by pQCT and microCT40 for volumetric (v)BMD, mid-diaphyseal cortical thickness, and distal trabecular phenotypes. Static and dynamic histomorphologic data were obtained on selected females and males at 16 wk.

RESULTS AND CONCLUSIONS

We found that the original BMD QTL, Bmd5, mapped to distal Chr 1 consists of three QTLs with different effects on vBMD and trabecular bone in both sexes. Compared with B6 controls, femoral vBMD, BMD, and cortical thickness (p < 0.0001) were significantly increased in congenic subline females, but not in males, carrying C3H alleles at QTL-1. Both females and males carrying C3H alleles at QTL-1 showed marked increases in BV/TV by microCT compared with B6 mice (p < 0.0001). Females increased BV/TV by increasing trabecular thickness, whereas males increased trabecular number. In addition, the microCT40 data showed two unique QTLs for male trabecular bone, QTL-2 and QTL-3, which may interact to regulate trabecular thickness and number. These QTLs are closely linked with and proximal to QTL-1. The histomorphometric data revealed sex-specific differences in cellular and bone formation parameters. Mice and humans share genetic homology between distal mouse Chr 1 and human Chr 1q20-24 that is associated with adult human skeletal regulation. Sex- and compartment-specific regulatory QTLs in the mouse suggest the need to partition human data by sex to improve accuracy of mapping and genetic loci identification.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterised by reduced bone mass and an increased risk of fragility fractures. Twin and family studies have shown that genetic factors contribute to osteoporosis by influencing bone mineral density (BMD), and other phenotypes that are associated with fracture risk, although the heritability of fracture itself is modest. Linkage studies have identified several quantitative trait loci that regulate BMD but most causal genes remain to be identified. In contrast, linkage studies in monogenic bone diseases have been successful in gene identification, and polymorphisms in many of these genes have been found to contribute to the regulation of bone mass in the normal population. Population-based studies have identified polymorphisms in several candidate genes that have been associated with bone mass or osteoporotic fracture, although individually these polymorphisms only account for a small amount of the genetic contribution to BMD regulation. Environmental factors such as diet and physical activity are also important determinants of BMD, and in some cases specific nutrients have been found to interact with genetic polymorphisms to regulate BMD. From a clinical standpoint, advances in knowledge about the genetic basis of osteoporosis are likely to be important in increasing the understanding of the pathophysiology of the disease; providing new genetic markers with which to assess fracture risk and in identifying genes and pathways that form molecular targets for the design of the next generation of drug treatments.

Polymorphisms in ALOX12, but not ALOX15, are significantly associated with BMD in postmenopausal women

The murine arachidonate 15-lipoxygenase gene (Alox15) has recently been identified as a negative regulator of peak bone mineral density (BMD). The human ALOX15 gene shares significant sequence homology with the murine Alox15 gene; however, the human arachidonate 12-lipoxygenase gene (ALOX12) is functionally more similar to the mouse gene. Multiple single-nucleotide polymorphisms (SNPs) in the human ALOX15 and ALOX12 genes have previously been reported to be significantly associated with BMD in humans. On the basis of these data, we carried out our own investigation of the human ALOX15 and ALOX12 genes and their relationship with hip and spine BMD parameters. The study population consisted of 779 postmenopausal women with a mean (+/- standard deviation) age of 62.5 +/- 5.9 years at BMD measurement and was recruited from a single large general practice in Chingford, northeast London. Three SNPs from ALOX15 and five from ALOX12 were analyzed. None of the SNPs that we analyzed in ALOX15 were significantly associated with any of the BMD parameters or fracture data. However, we found that three SNPs from ALOX12, all previously associated with spine BMD in women, were significantly associated with spine and various hip BMD parameters in our cohort (P = 0.029-0.049). In conclusion, we found no association between polymorphism in ALOX15 and BMD phenotypes but were able to replicate previous findings that genetic variation in ALOX12 seems to play a role in determining bone structure in Caucasian women.

A bivariate whole-genome linkage scan suggests several shared genomic regions for obesity and osteoporosis

CONTEXT

A genome-wide bivariate analysis was conducted for body fat mass (BFM) and bone mineral density (BMD) in a large Caucasian sample. We found some quantitative trait loci shared by BFM and BMD in the total sample and the gender-specific subgroups, and quantitative trait loci with potential pleiotropy were disclosed. BFM and BMD, as the respective measure for obesity and osteoporosis, are phenotypically and genetically correlated. However, specific genomic regions accounting for their genetic correlation are unknown.

OBJECTIVE

To identify systemically the shared genomic regions for BFM and BMD, we performed a bivariate whole-genome linkage scan in 4498 Caucasian individuals from 451 families for BFM and BMD at the hip, spine, and wrist, respectively. Linkage analyses were performed in the total sample and the male and female subgroups, respectively.

RESULTS

In the entire sample, suggestive linkages were detected at 7p22-p21 (LOD 2.69) for BFM and spine BMD, 6q27 (LOD 2.30) for BFM and hip BMD, and 11q13 (LOD 2.64) for BFM and wrist BMD. Male-specific suggestive linkages were found at 13q12 (LOD 3.23) for BFM and spine BMD and at 7q21 (LOD 2.59) for BFM and hip BMD. Female-specific suggestive LOD scores were 3.32 at 15q13 for BFM and spine BMD and 3.15 at 6p25-24 for BFM and wrist BMD.

CONCLUSIONS

Several shared genomic regions for BFM and BMD were identified here. Our data may benefit further positional and functional studies, aimed at eventually uncovering the complex mechanism underlying the shared genetic determination of obesity and osteoporosis.

Genetics and osteoporosis

Over the past 10 years, many advances have been made in understanding the mechanisms by which genetic factors regulate susceptibility to osteoporosis. It has become clear from studies in man and experimental animals that different genes regulate BMD at different skeletal sites and in men and women. Linkage studies have identified several chromosomal regions that regulate BMD, but only a few causative genes have been discovered so far using this approach. In contrast, significant advances have been made in identifying the genes that cause monogenic bone diseases, and polymorphic variation is some of these genes has been found to contribute to the genetic regulation of BMD in the normal population. Other genes that have been investigated as possible candidates for susceptibility to osteoporosis because of their role in bone biology, such as vitamin D, have yielded mixed results. Many candidate gene association studies have been underpowered, and meta-analysis has been used to try to confirm or refute potential associations and gain a better estimate of their true effect size in the population. Most of the genetic variants that confer susceptibility to osteoporosis remain to be discovered. It is likely that new techniques such as whole-genome association will provide new insights into the genetic determinants of osteoporosis and will help to identify genes of modest effect size. From a clinical standpoint, genetic variants that are found to predispose to osteoporosis will advance our understanding of the pathophysiology of the disease. They could be developed as diagnostic genetic tests or form molecular targets for design of new drugs for the prevention and treatment of osteoporosis and other bone diseases.

The genetic, environmental and phenotypic correlations of bone phenotypes at the spine and hip in Chinese

BACKGROUND

Bone mineral density (BMD), bone mineral content (BMC), and bone size have been widely studied individually as important risk factors for osteoporotic fracture, but little is known about the correlation and the degree of sharing genetic and environmental factors between the pairs of the three phenotypes.

AIM

The study investigated genetic correlation (rhoG), environmental correlation (rhoE) and phenotypic correlation (rhoP) between BMD, BMC and bone size.

SUBJECTS AND METHODS

Bivariate variance decomposition analyses were performed in 904 subjects from 287 Chinese nuclear families.

RESULTS

Significant rhoE, rhoG and rhoP were detected between BMD, BMC and bone size, except for rhoE between BMD and bone size at the hip (rhoE = 0.121, p = 0.361). Common shared genetic factors explained 86.1% and 60% of BMD and BMC genetic variations at the spine and hip, respectively. However, the genetic and environmental correlations between BMD and bone size were limited. rhoE and rhoG at the spine were 0.392 and 0.381, and at the hip were 0.121 and -0.205, respectively. Only 14.5% and 4.2% of variations between BMD and bone size at the spine and hip may be due to the shared genetic factors.

CONCLUSION

The obtained results suggested that bone size may be used as another surrogate phenotype independently of BMD for eventual elucidation of the pathogenesis of osteoporosis because of the limited correlations between BMD and bone size.