Genome-wide scan identified QTLs underlying femoral neck cross-sectional geometry that are novel studied risk factors of osteoporosis

Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children

CONTEXT

Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied.

OBJECTIVE

We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children.

METHODS

We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry.

RESULTS

Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7).

CONCLUSIONS

We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

Liuwei Dihuang pill cures postmenopausal osteoporosis with kidney-Yin deficiency: Potential therapeutic targets identified based on gene expression profiling

This study aimed to investigate the potential therapeutic targets of Liuwei Dihuang pill (LDP) in the treatment of postmenopausal osteoporosis with kidney-Yin deficiency (PMO-KY).Gene expression data were downloaded from the GEO database, including 4 PMO-KY samples and 3 healthy postmenopausal controls from GSE56116, as well as 3 PMO-KY samples before LDP treatment and 3 PMO-KY samples after three months of LDP treatment from GSE57273. Limma package was used to identify differentially expressed genes (DEGs). Afterwards, the potential target genes of LDP (namely key DEGs) were identified according to the comparison of DEGs in PMO-KY group and the DEGs in LDP treatment groups. Subsequently, iRegulon plugin in Cytoscape software was used to predict potential transcription factors (TFs) that regulated the key DEGs, and Comparative Toxicogenomics Database was utilized to identify known PMO-related genes among the key DEGs.Totally, 202 and 2066 DEGs were identified between PMO-KY and controls, as well as after-treatment and before-treatment groups, respectively. Among them, 52 DEGs were up-regulated in PMO-KY but down-regulated after LDP treatment, and 8 TFs were predicted to these DEGs. Furthermore, 34 DEGs were down-regulated in PMO-KY but up-regulated after treatment, and 7 TFs were predicted to regulate these DEGs. Additionally, 43 of the 86 key DEGs were known PMO-related genes.NCOA3, TCF4, DUSP6, PELI2, and STX7 were predicted to be regulated by HOXA13. In the PMO-KY treatment, NCOA3, TCF4, DUSP6, PELI2, and STX7 might be the potential therapeutic targets of LDP. However, further investigation is required to confirm these genes.

Gender differences in a reference database of age-related femoral neck geometric parameters for Chinese population and their association with femoral neck fractures

Femoral neck geometric parameters (FNGPs) are closely related to the strength of the femoral neck and the risk of fragility fractures. No reference database is available for FNGPs for Chinese population, and gender-related differences in FNGPs as well as their association with the risk of femoral neck fractures are unknown. This investigation aimed to set up reference databases for FNGPs, understand gender-related differences in FNGPs, and examine the association between FNGPs and the risk of osteoporotic fractures of the femoral neck. This study included 5268 females and 2156 males (aged 15-91years) from Chinese population. A total of 384 patients (282 females and 102 males) had sustained femoral neck fractures; 384 age- and sex-matched individuals without any fractures served as controls. Femoral neck DXA images were used to measure bone mineral density (BMD) and eight FNGPs. Our results showed that the age-related trends of FNGPs were fitted with the best goodness-of-fit by applying the cubic regression model. The trends shown by FNGPs were significantly different between male and female subjects, and the fitting curves were significantly higher in male subjects. After adjustments were made for age, height, weight, and body mass index, Cox regression analysis showed that changes in all FNGPs were related to increased hazard ratios (HRs) of femoral neck fractures. After further adjustment was made for BMD of the femoral neck, the HRs related to a cortical thickness (CT) decrease and buckling ratio (BR) increase in females went up by 3.35-folds (95% CI: 2.75-4.07) and 1.86-folds (95% CI: 1.33-2.60), respectively. In males, the HRs related to the decrease in CT and cross-sectional area (CSA) increased by 3.21-folds (95% CI: 2.32-4.45) and 1.88-folds (95% CI: 1.03-3.44), respectively. In conclusions, the reference databases of FNGPs established in this study will assist in the evaluation and prediction of femoral neck fracture risk in the clinic. The decrease in CT and increase in BR of the femoral neck were independent risk factors for osteoporotic fractures of the femoral neck in females from mainland China, while a decrease in CT and CSA were risk factors in male.

Genetic Contribution of Femoral Neck Bone Geometry to the Risk of Developing Osteoporosis: A Family-Based Study

Femoral neck geometry parameters are believed to be as good as bone mineral density as independent factors in predicting hip fracture risk. This study was conducted to analyze the roles of genetic and environmental factors in femoral properties measured in a sample of Spanish families with osteoporotic fractures and extended genealogy. The "Genetic Analysis of Osteoporosis (GAO) Project" involved 11 extended families with a total number of 376 individuals. We studied three categorical phenotypes of particular clinical interest and we used a Hip structural analysis based on DXA to analyze 17 strength and geometrical phenotypes of the hip. All the femoral properties had highly significant heritability, ranging from 0.252 to 0.586. The most significant correlations were observed at the genetic level (ρG). Osteoporotic fracture status (Affected 2) and, particularly, low bone mass and osteoporotic condition (Affected 3) had the highest number of significant genetic correlations with diverse femoral properties. In conclusion, our findings suggest that a relatively simple and easy to use method based on DXA studies can provide useful data on properties of the Hip in clinical practice. Furthermore, our results provide a strong motivation for further studies in order to improve the understanding of the pathophysiological mechanism underlying bone architecture and the genetics of osteoporosis.

Osteoporosis: A Silent Disease with Complex Genetic Contribution

Osteoporosis is the most common multifactorial metabolic bone disorder worldwide with a strong genetic component. In this review, the evidence for a genetic contribution to osteoporosis and related phenotypes is summarized alongside with methods used to identify osteoporosis susceptibility genes. The key biological pathways involved in the skeleton and bone development are discussed with a particular focus on master genes clustered in these pathways and their mode of action. Furthermore, the most studied single nucleotide polymorphisms (SNPs) analyzed for their importance as genetic markers of the disease are presented. New data generated by next-generation sequencing in conjunction with extensive meta-analyses should contribute to a better understanding of the genetic basis of osteoporosis and related phenotype variability. These data could be ultimately used for identifying at-risk patients for disease prevention by both controlling environmental factors and providing possible therapeutic targets.

Bivariate genome-wide association analyses identified genes with pleiotropic effects for femoral neck bone geometry and age at menarche

Femoral neck geometric parameters (FNGPs), which include cortical thickness (CT), periosteal diameter (W), buckling ratio (BR), cross-sectional area (CSA), and section modulus (Z), contribute to bone strength and may predict hip fracture risk. Age at menarche (AAM) is an important risk factor for osteoporosis and bone fractures in women. Some FNGPs are genetically correlated with AAM. In this study, we performed a bivariate genome-wide association study (GWAS) to identify new candidate genes responsible for both FNGPs and AAM. In the discovery stage, we tested 760,794 SNPs in 1,728 unrelated Caucasian subject, followed by replication analyses in independent samples of US Caucasians (with 501 subjects) and Chinese (with 826 subjects). We found six SNPs that were associated with FNGPs and AAM. These SNPs are located in three genes (i.e. NRCAM, IDS and LOC148145), suggesting these three genes may co-regulate FNGPs and AAM. Our findings may help improve the understanding of genetic architecture and pathophysiological mechanisms underlying both osteoporosis and AAM.

Insights into the genetics of osteoporosis from recent genome-wide association studies

Osteoporosis, which is characterised by reduced bone mineral density (BMD) and an increased risk of fragility fractures, is the result of a complex interaction between environmental factors and genetic variants that confer susceptibility. Heritability studies have shown that BMD and other osteoporosis-related traits such as ultrasound properties of bone, skeletal geometry and bone turnover have significant inheritable components. Although previous linkage and candidate gene studies have provided few replicated loci for osteoporosis, genome-wide association approaches have produced clear and reproducible findings. To date, 20 genome-wide association studies (GWASs) for osteoporosis and related traits have been conducted, identifying dozens of genes. Further meta-analyses of GWAS data and deep resequencing of rare variants will uncover more novel susceptibility loci and ultimately provide possible therapeutic targets for fracture prevention.

Linkage mapping of principal components for femoral biomechanical performance in a reciprocal HCB-8 × HCB-23 intercross

Studies of bone genetics have addressed an array of related phenotypes, including various measures of biomechanical performance, bone size, bone, shape, and bone mineral density. These phenotypes are not independent, resulting in redundancy of the information they provide. Principal component (PC) analysis transforms multiple phenotype data to a new set of orthogonal "synthetic" phenotypes. We performed PC analysis on 17 femoral biomechanical, anatomic, and body size phenotypes in a reciprocal intercross of HcB-8 and HcB-23, accounting for 80% of the variance in 4 PCs. Three of the 4 PCs were mapped in the cross. The linkage analysis revealed a quantitative trait locus (QTL) with LOD = 4.7 for PC2 at 16 cM on chromosome 19 that was not detected using the directly measured phenotypes. The chromosome 19 QTL falls within a ~10 megabase interval, with Osf1 as a positional candidate gene. PC QTLs were also found on chromosomes 1, 2, 4, 6, and 10 that coincided with those identified for directly measured or calculated material property phenotypes. The novel chromosome 19 QTL illustrates the power advantage that attends use of PC phenotypes for linkage mapping. Constraint of the chromosome 19 candidate interval illustrates an important advantage of experimental crosses between recombinant congenic mouse strains.

Genetic aspects of osteoporosis

The multiple factors contributing to the pathogenesis of osteoporosis include genetic and environmental factors. Because decrease in bone mineral density (BMD) is the major clinical indicator and a useful quantitative trait, many association and linkage studies of BMD have been conducted. Although the series of studies showed apparently significant associations, the genes have not been found that can be utilized in clinical practice. Several genes identified in robust genome-wide association studies will be the new cutting edge in genetic studies of osteoporosis. Our recent reports of functional single nucleotide polymorphism in the tissue-nonspecific alkaline phosphatase gene and gamma-carboxylase gene are presented in this review to discuss the future prospects in the genetic research of osteoporosis from the point of view of genome-nutrition interaction.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterized by reduced bone mass, defects in the microarchitecture of bone tissue, and an increased risk of fragility fractures. Twin and family studies have shown high heritability of bone mineral density (BMD) and other determinants of fracture risk such as ultrasound properties of bone, skeletal geometry, and bone turnover. Osteoporotic fractures also have a heritable component, but this reduces with age as environmental factors such as risk of falling come into play. Susceptibility to osteoporosis is governed by many different genetic variants and their interaction with environmental factors such as diet and exercise. Notable successes in identification of genes that regulate BMD have come from the study of rare Mendelian bone diseases characterized by major abnormalities of bone mass where variants of large effect size are operative. Genome-wide association studies have also identified common genetic variants of small effect size that contribute to regulation of BMD and fracture risk in the general population. In many cases, the loci and genes identified by these studies had not previously been suspected to play a role in bone metabolism. Although there has been extensive progress in identifying the genes and loci that contribute to the regulation of BMD and fracture over the past 15 yr, most of the genetic variants that regulate these phenotypes remain to be discovered.

Refined QTLs of osteoporosis-related traits by linkage analysis with genome-wide SNPs: Framingham SHARe

Genome-wide association studies (GWAS) using high-density array of single-nucleotide polymorphisms (SNPs) offer an unbiased strategy to identify new candidate genes for osteoporosis. We used a subset of autosomal SNPs from the Affymetrix 500K+50K SNP GeneChip marker set to examine genetic linkage with multiple highly heritable osteoporosis-related traits, including BMD of the hip and spine, heel ultrasound (attenuation and speed of sound), and geometric indices of the hip, in two generations from the Framingham Osteoporosis Study. Variance component linkage analysis was performed using normalized residuals (adjusted for age, height, BMI, and estrogen status in women). Multipoint linkage analyses produced LOD scores > or =3.0 for BMD on chromosomes (chr.) 9 and 11 and for ultrasound speed of sound on chr. 5. Hip geometric traits were linked with higher LOD scores, such as with shaft width on chr. 4 (LOD=3.9) and chr. 16 (LOD=3.8) and with shaft section modulus on chr. 22 (LOD=4.0). LOD score > or =5.0 was obtained for femoral neck width on chr. 7. In conclusion, with an SNP-based linkage approach, we identified several novel potential QTLs and confirmed previously identified chromosomal regions linked to bone mass and geometry. Subsequent focus on the spectrum of genetic polymorphisms in these refined regions may contribute to finding variants predisposing to osteoporosis.

Genome-wide association study for femoral neck bone geometry

Poor femoral neck bone geometry at the femur is an important risk factor for hip fracture. We conducted a genome-wide association study (GWAS) of femoral neck bone geometry, examining approximately 379,000 eligible single-nucleotide polymorphisms (SNPs) in 1000 Caucasians. A common genetic variant, rs7430431 in the receptor transporting protein 3 (RTP3) gene, was identified in strong association with the buckling ratio (BR, P = 1.6 x 10(-7)), an index of bone structural instability, and with femoral cortical thickness (CT, P = 1.9 x 10(-6)). The RTP3 gene is located in 3p21.31, a region that we found to be linked with CT (LOD = 2.19, P = 6.0 x 10(-4)) in 3998 individuals from 434 pedigrees. The replication analyses in 1488 independent Caucasians and 2118 Chinese confirmed the association of rs7430431 to BR and CT (combined P = 7.0 x 10(-3) for BR and P = 1.4 x 10(-2) for CT). In addition, 350 hip fracture patients and 350 healthy control individuals were genotyped to assess the association of the RTP3 gene with the risk of hip fracture. Significant association between a nearby common SNP, rs10514713 of the RTP3 gene, and hip fracture (P = 1.0 x 10(-3)) was found. Our observations suggest that RTP3 may be a novel candidate gene for femoral neck bone geometry.

Pathway-based genome-wide association analysis identified the importance of EphrinA-EphR pathway for femoral neck bone geometry

Femoral neck (FN) bone geometry is an important predictor of bone strength with high heritability. Previous studies have revealed certain candidate genes for FN bone geometry. However, the majority of the underlying genetic factors remain to be discovered. In this study, pathway-based genome-wide association analysis was performed to explore the joint effects of genes within biological pathways on FN bone geometry variations in a cohort of 1000 unrelated US whites. Nominal significant associations (nominal p value<0.05) were observed between 76 pathways and a key FN bone geometry variable-section modulus (Z), biomechanically indicative of bone strength subject to bending. Among them, EphrinA-EphR pathway was most significantly associated with FN Z even after multiple testing adjustments (p(FWER) value=0.035). The association of EphrinA-EphR pathway with FN Z was also observed in an independent sample from Framingham Osteoporosis Study. Overall, these results suggest the significant genetic contribution of EphrinA-EphR pathway to femoral neck bone geometry.

Association analyses of vitamin D-binding protein gene with compression strength index variation in Caucasian nuclear families

This study was conducted to test whether there exists an association between vitamin D-binding protein (DBP) gene and compression strength index (CSI) phenotype. Candidate gene association analyses were conducted in total sample, male subgroup, and female subgroup, respectively. Two single-nucleotide polymorphisms (SNPs) with significant association results were found in males, suggesting the importance of DBP gene polymorphisms on the variation in CSI especially in Caucasian males.

INTRODUCTION

CSI of the femoral neck (FN) is a newly developed phenotype integrating information about bone size, body size, and bone mineral density. It is considered to have the potential to improve the performance of risk assessment for hip fractures because it is based on a combination of phenotypic traits influencing hip fractures rather than a single trait. CSI is under moderate genetic determination (with a heritability of approximately 44% found in this study), but the relevant genetic study is still rather scarce.

METHODS

Based on the known physiological role of DBP in bone biology and the relatively high heritability of CSI, we tested 12 SNPs of the DBP gene for association with CSI variation in 405 Caucasian nuclear families comprising 1,873 subjects from the Midwestern US. Association analyses were performed in the total sample, male and female subgroups, respectively.

RESULTS

Significant associations with CSI were found with two SNPs (rs222029, P = 0.0019; rs222020, P = 0.0042) for the male subgroup. Haplotype-based association tests corroborated the single-SNP results.

CONCLUSIONS

Our findings suggest that the DBP gene might be one of the genetic factors influencing CSI phenotype in Caucasians, especially in males.

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.

The contributions of lean tissue mass and fat mass to bone geometric adaptation at the femoral neck in Chinese overweight adults

BACKGROUND

Overweight or obese populations may have lower risk of osteoporotic fractures and higher bone mineral density (BMD), while bone strength is determined not only by bone material but also by bone structural parameters. Thus, the influence of body weight on bone geometry was examined in Chinese overweight adults.

AIM

The purpose of this study was to explore how total body lean mass (TBLM) and total body fat mass (TBFM) contribute to the variation of bone geometry at the femoral neck in Chinese overweight adults.

SUBJECTS AND METHODS

Bone geometric parameters including section modulus (Z), cross-sectional area (CSA), subperiosteal width (W), cortical thickness (CT) and buckling ratio (BR) were compared in 100 overweight (body mass index, BMI >/= 23) vs. 100 underweight subjects (BMI </= 18.5) in Chinese female and male adults aged 20-44 years by multiple regression analyses.

RESULTS

Multiple regression analysis revealed that both TBLM and TBFM were significantly higher in overweight subjects than in underweight subjects. Meanwhile, significant differences in bone geometric parameters (except W) were also detected between the overweight and underweight groups after adjustment for age and height in both sexes (p </= 0.001). Bone bending strength Z and axial strength CSA were 14% and 13% higher in females, as well as 18% and 20% higher in males in the overweight group than in the underweight group, respectively. The significant differences mentioned above were not observed when adjusted for TBLM, age, and height. TBLM seemed to be the strongest significant positive predictor of bone geometric parameters (p < 0.001), with the exception of W in both sexes and BR in females, while TBFM did not contribute significantly to the bone geometric parameters (p > 0.055 for both sexes).

CONCLUSION

Bone geometry may adapt primarily to mechanical load as represented by TBLM, but TBFM seemed to have no independent effect on bone geometry in Chinese overweight subjects.

Bivariate whole-genome linkage scan for bone geometry and total body fat mass

To quantify the genetic correlations between total body fat mass (TBFM) and femoral neck geometric parameters (FNGPs) and, if possible, to detect the specific genomic regions shared by them, bivariate genetic analysis and bivariate whole-genome linkage scan were carried out in a large Caucasian population. All the phenotypes studied were significantly controlled by genetic factors (P < 0.001) with the heritabilities ranging from 0.45 to 0.68. Significantly genetic correlations were found between TBFM and CSA (cross-section area), W (sub-periosteal diameter), Z (section modulus) and CT (cortical thickness) except between TBFM and BR (buckling ratio). The peak bivariate LOD scores were 3.23 (20q12), 2.47 (20p11), 3.19 (6q27), 1.68 (20p12), and 2.47 (7q11) for the five pairs of TBFM and BR, CSA, CT, W, and Z in the entire sample, respectively. Gender-specific bivariate linkage evidences were also found for the five pairs. 6p25 had complete pleiotropic effects on the variations of TBFM & Z in the female sub-population, and 6q27 and 17q11 had coincident linkages for TBFM & CSA and TBFM & Z in the entire population. We identified moderate genetic correlations and several shared genomic regions between TBFM and FNGPs in a large Caucasian population.

A genome wide linkage scan of metacarpal size and geometry in the Framingham Study

Bone geometry is a significant component of bone strength, and has a clinical utility in predicting fractures and quantifying bone loss. Bone geometry is known to have a substantial genetic component. We performed linkage analysis to identify chromosomal regions governing metacarpal bone geometry. A genome-wide scan (with a set of 615 markers with spacing of approximately 5.7 cM) was performed on 1,702 individuals from 330 extended families of the Framingham Study. Midshaft width was measured and several indices calculated, namely Metacarpal Cortical Thickness (MCT), Cortical Index (MCI), and Section Modulus (MZ), using digitized X-rays of 1,380 participants (men, n = 666, mean age 55.2 yr, women, n = 714, 55.5 yr). Metacarpals 2, 3, and 4 were averaged. Heritability was significant for all indices, ranging from 0.51 to 0.72. Linkage analysis of indices adjusted for age, age(2), and estrogen status in women, identified chromosomal regions 6p21, 9p21, 11q21-q22, and Xq26-Xq27, with LOD scores >2.0. Additional adjustment for smoking, height, and BMI, generally reduced the LOD scores. Finally, bivariate linkage analysis confirmed that a QTL on chr. 6 (51 cM) was shared by midshaft width and MZ (LOD = 2.40, adjusted for all covariates). Neither MCT nor MCI shared linkage loci with width or MZ. In conclusion, we have identified chromosomal regions potentially linked to bone geometry. Genes in these regions may regulate bone geometry via effects on body size. Identification and subsequent characterization of loci for bone geometry can further elucidate the genetic contributions to bone's resistance to stress.

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.

Bivariate genome-wide linkage analysis for traits BMD and AAM: effect of menopause on linkage signals

Osteoporosis is an age-related systemic skeletal disease, characterized by low bone mineral density (BMD). Low BMD is closely associated with late age at menarche (AAM). Our previous bivariate genome-wide linkage analyses (GWLAs) between BMD and AAM identified two shared genomic regions in 2584 Caucasian females including both pre- and post-menopausal females. However, menopause often causes dramatic bone loss in post-menopausal females; this may introduce some confounding effects on the bivariate GWLA for BMD and AAM. To address the effect of menopause on the identification of genetic factors shared by BMD and AAM, we segregated the previously studied population of 2584 females into two separate subgroups consisting of 1462 pre-menopause subjects and 1122 post-menopausal subjects, and performed further bivariate GWLAs. The BMD was measured by Hologic Dual-energy X-ray (DXA) scanners (Hologic Inc., Bedford, MA, USA). Based on the genome-wide thresholds corrected for multiple testing, we found more significant genomic regions in the pre-menopausal group than in total group (including pre- and post-menopausal women), e.g., we found 4, 1, and 2 shared by spine BMD and AAM, femoral neck (FNK) BMD and AAM and ultra distal (UD) BMD and AAM, respectively. We did not found any significant linkage signals in the post-menopausal group. Importantly, the linkage signals at all significant regions were much stronger in pre-menopausal group than in the other groups: post-menopausal females and total females. For example, the linkage LOD score for FNK BMD and AAM is as high as 4.88 in pre-menopausal females, but only 0.24 and 0.31 in post-menopausal and total females, respectively. These results suggest that menopause introduces some noise signals into GWLAs when estimating the shared genetic factors by BMD and AAM. Therefore, it is very important to classify female subjects properly according to their menopause stage when performing such studies.

A bivariate whole genome linkage study identified genomic regions influencing both BMD and bone structure

Areal BMD (aBMD) and areal bone size (ABS) are biologically correlated traits and are each important determinants of bone strength and risk of fractures. Studies showed that aBMD and ABS are genetically correlated, indicating that they may share some common genetic factors, which, however, are largely unknown. To study the genetic factors influencing both aBMD and ABS, bivariate whole genome linkage analyses were conducted for aBMD-ABS at the femoral neck (FN), lumbar spine (LS), and ultradistal (UD)-forearm in a large sample of 451 white pedigrees made up of 4498 individuals. We detected significant linkage on chromosome Xq27 (LOD = 4.89) for LS aBMD-ABS. In addition, we detected suggestive linkages at 20q11 (LOD = 3.65) and Xp11 (LOD = 2.96) for FN aBMD-ABS; at 12p11 (LOD = 3.39) and 17q21 (LOD = 2.94) for LS aBMD-ABS; and at 5q23 (LOD = 3.54), 7p15 (LOD = 3.45), Xq27 (LOD = 2.93), and 12p11 (LOD = 2.92) for UD-forearm aBMD-ABS. Subsequent discrimination analyses indicated that quantitative trait loci (QTLs) at 12p11 and 17q21 may have pleiotropic effects on aBMD and ABS. This study identified several genomic regions that may contain QTLs important for both aBMD and ABS. Further endeavors are necessary to follow these regions to eventually pinpoint the genetic variants affecting bone strength and risk of fractures.

Role of Nonbehavioral Factors in Adjusting Long Bone Diaphyseal Structure in Free-ranging Pan troglodytes

Limb bones deform during locomotion and can resist the deformations by adjusting their shapes. For example, a tubular-shaped diaphysis best resists variably-oriented deformations. As behavioral profiles change during adulthood, patterns of bone deformation may exhibit age trends. Habitat characteristics, e.g., annual rainfall, tree density, and elevation changes, may influence bone deformations by eliciting individual components of behavioral repertoires and suppressing others, or by influencing movements during particular components. Habituated chimpanzee communities provide a unique opportunity to examine these factors because of the availability of morphological data and behavioral observations from known-age individuals inhabiting natural habitats. We evaluated adult femora and humeri of 18 female and 10 male free-ranging chimpanzees (Pan troglodytes) from communities in Gombe (Tanzania), Mahale Mountains (Tanzania), and Taï Forest (Côte d'Ivoire) National Parks. We compare cross sections at several locations (35%, 50%, 65% diaphyseal lengths). Community comparisons highlight different diaphyseal shapes of Taï females relative to Mahale and Gombe females, particularly in humeral diaphyses. Age trends in diaphyseal shapes are consistent with reduced activity levels in general, not only reduced arboreal activity. Age-related bone loss is apparent among community females, but is less striking among males. Community trends in diaphyseal shape are qualitatively consistent with ranked annual rainfall at localities, tree density, and elevation change or ruggedness of terrain. Habitat characteristics may contribute to variation in diaphyseal shape among chimpanzee communities, much like among modern human groups, but verification awaits further rigorous experimental and comparative analyses.

Identification of PLCL1 gene for hip bone size variation in females in a genome-wide association study

Osteoporosis, the most prevalent metabolic bone disease among older people, increases risk for low trauma hip fractures (HF) that are associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of the key measurable risk factors for HF. Although hip BS is highly genetically determined, genetic factors underlying the trait are still poorly defined. Here, we performed the first genome-wide association study (GWAS) of hip BS interrogating ∼380,000 SNPs on the Affymetrix platform in 1,000 homogeneous unrelated Caucasian subjects, including 501 females and 499 males. We identified a gene, PLCL1 (phospholipase c-like 1), that had four SNPs associated with hip BS at, or approaching, a genome-wide significance level in our female subjects; the most significant SNP, rs7595412, achieved a p value of 3.72×10⁻⁷. The gene's importance to hip BS was replicated using the Illumina genotyping platform in an independent UK cohort containing 1,216 Caucasian females. Two SNPs of the PLCL1 gene, rs892515 and rs9789480, surrounded by the four SNPs identified in our GWAS, achieved p values of 8.62×10⁻³ and 2.44×10⁻³, respectively, for association with hip BS. Imputation analyses on our GWAS and the UK samples further confirmed the replication signals; eight SNPs of the gene achieved combined imputed p values<10⁻⁵ in the two samples. The PLCL1 gene's relevance to HF was also observed in a Chinese sample containing 403 females, including 266 with HF and 177 control subjects. A SNP of the PLCL1 gene, rs3771362 that is only ∼0.6 kb apart from the most significant SNP detected in our GWAS (rs7595412), achieved a p value of 7.66×10⁻³ (odds ratio = 0.26) for association with HF. Additional biological support for the role of PLCL1 in BS comes from previous demonstrations that the PLCL1 protein inhibits IP3 (inositol 1,4,5-trisphosphate)-mediated calcium signaling, an important pathway regulating mechanical sensing of bone cells. Our findings suggest that PLCL1 is a novel gene associated with variation in hip BS, and provide new insights into the pathogenesis of HF.

Autosome-wide linkage analysis of hip structural phenotypes in the Old Order Amish

INTRODUCTION

Fracture risk is associated with bone mineral density (BMD) and with other indices of bone strength, including hip geometry. While the heritability and associated fracture risk of BMD are well described, less is known about genetic influences of bone geometry. We derived hip structural phenotypes using the Hip Structural Analysis program (HSA) and performed autosome-wide linkage analysis of hip geometric structural phenotypes.

MATERIALS AND METHODS

The Amish Family Osteoporosis Study was designed to identify genes affecting bone health. BMD was measured at the hip using dual X-ray absorptiometry (DXA) in 879 participants (mean age+/-SD=49.8+/-16.1 years, range 18-91 years) from large multigenerational families. From DXA scans, we computed structural measures of hip geometry at the femoral neck (NN) and shaft (S) by HSA, including cross-sectional area (CSA), endocortical or inner diameter (ID), outer diameter (OD) buckling ratio (BR) and section modulus (Z). Genotyping of 731 highly polymorphic microsatellite markers (average spacing of 5.4 cM) and autosome-wide multipoint linkage analysis was performed.

RESULTS

The heritability of HSA-derived hip phenotypes ranged from 40 to 84%. In the group as a whole, autosome-wide linkage analysis suggested evidence of linkage for QTLs related to NN_Z on chromosome 1p36 (LOD=2.36). In subgroup analysis, ten additional suggestive regions of linkage were found on chromosomes 1, 2, 5, 6, 11, 12, 14, 15 and 17, all with LOD>2.3 except for our linkage at 17q11.2-13 for men and women age 50 and under for NN_CSA, which had a lower LOD of 2.16, but confirmed a previous linkage report.

CONCLUSIONS

We found HSA-derived measures of hip structure to be highly heritable independent of BMD. No strong evidence of linkage was found for any phenotype. Confirmatory evidence of linkage was found on chromosome 17q11.2-12 for NN_CSA. Modest evidence was found for genes affecting hip structural phenotypes at ten other chromosomal locations.

Genetic variation in the patterns of skeletal progenitor cell differentiation and progression during endochondral bone formation affects the rate of fracture healing

These studies examined how genetic differences that regulate architectural and bone material properties would be expressed during fracture healing and determine whether any of these features would affect rates of healing as defined by regain of strength. Controlled fractures were generated in three inbred strains of mice: A/J, C57Bl/6J (B6), and C3H/HeJ (C3H). Both the A/J and B6 strains showed faster healing than the C3H strain based on regains in strength and stiffness. Strain-specific architectural features such as moment of inertia, cross-sectional area, and cortical thickness were all recapitulated during the development of the callus tissues. None of these traits were directly relatable to rates of fracture healing. However, rates of healing were related to variations in the temporal patterns of chondrogenic and osteogenic lineage development. The B6 strain expressed the highest percentage of cartilage gene products and had the longest period of chondrocyte maturation and hypertrophy. The slowest healing strain (C3H) had the shortest period of chondrogenic development and earliest initiation of osteogenic development. Although the A/J strain showed an almost identical pattern of chondrogenic development as the C3H strain, A/J initiated osteogenic development several days later than C3H during fracture healing. Long bone growth plates at 28 days after birth showed similar strain-specific variation in cartilage tissue development as seen in fracture healing. Thus, the B6 strain had the largest growth plate heights, cell numbers per column, and the largest cell size, whereas the C3H columns were the shortest, had the smallest number of cells per column, and showed the smallest cell sizes. These results show that (1) different strains of mice express variations of skeletal stem cell lineage differentiation and (2) that these variations affect the rate of fracture healing.

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.

Genetics of the musculoskeletal system: a pleiotropic approach

The risk of osteoporotic fracture can be viewed as a function of loading conditions and the ability of the bone to withstand the load. Skeletal loads are dominated by muscle action. Recently, it has become clear that bone and muscle share genetic determinants. Involution of the musculoskeletal system manifests as bone loss (osteoporosis) and muscle wasting (sarcopenia). Therefore, the consideration of pleiotropy is an important aspect in the study of the genetics of osteoporosis and sarcopenia. This Perspective will provide the evidence for a shared genetic influence on bone and muscle. We will start with an overview of accumulating evidence that physical exercise produces effects on the adult skeleton, seeking to unravel some of the contradictory findings published thus far. We will provide indications that there are pleiotropic relationships between bone structure/mass and muscle mass/function. Finally, we will offer some insights and practical recommendations as to the value of studying shared genetic factors and will explore possible directions for future research. We consider several related questions that together comprise the general paradigm of bone responses to mechanical loading and the relationship between muscle strength and bone parameters, including the genetic factors that modulate these responses. We believe that further progress in understanding the common genetic etiology of osteoporosis and sarcopenia will provide valuable insight into important biological underpinnings for both conditions and may translate into new approaches to reduce the burdens of both conditions through improved diagnosis, prevention, and early targeted treatment.

Contribution of gender-specific genetic factors to osteoporosis risk

Common diseases result from the complex relationship between genetic and environmental factors. The aim of this review is to provide perspective for a conceptual framework aimed at studying the interplay of gender-specific genetic and environmental factors in the etiology of complex disease, using osteoporosis as an example. In recent years, gender differences in the heritability of the osteoporosis-related phenotypes have been reported and sex-specific quantitative-trait loci were discovered by linkage studies in humans and mice. Results of numerous allelic association studies also differed by gender. In most cases, it was not clear whether or not this phenomenon should be attributed to the effect of sex-chromosomes, sex hormones, or other intrinsic or extrinsic differences between the genders, such as the level of bioavailable estrogen and of physical activity. We conclude that there is need to consider gender-specific genetic and environmental factors in the planning of future association studies on the etiology of osteoporosis and other complex diseases prevalent in the general population.

Molecular genetic studies of gene identification for osteoporosis

This review comprehensively summarizes the most important and representative molecular genetics studies of gene identification for osteoporosis published up to the end of September 2007. It is intended to constitute a sequential update of our previously published reviews covering the available data up to the end of 2004. Evidence from candidate gene-association studies, genome-wide linkage and association studies, as well as functional genomic studies (including gene-expression microarray and proteomics) on osteogenesis and osteoporosis, are reviewed separately. Studies of transgenic and knockout mice models relevant to osteoporosis are summarized. The major results of all studies are tabulated for comparison and ease of reference. Comments are made on the most notable findings and representative studies for their potential influence and implications on our present understanding of genetics of osteoporosis. The format adopted by this review should be ideal for accommodating future new advances and studies.

Genome-wide association with bone mass and geometry in the Framingham Heart Study

BACKGROUND

Osteoporosis is characterized by low bone mass and compromised bone structure, heritable traits that contribute to fracture risk. There have been no genome-wide association and linkage studies for these traits using high-density genotyping platforms.

METHODS

We used the Affymetrix 100K SNP GeneChip marker set in the Framingham Heart Study (FHS) to examine genetic associations with ten primary quantitative traits: bone mineral density (BMD), calcaneal ultrasound, and geometric indices of the hip. To test associations with multivariable-adjusted residual trait values, we used additive generalized estimating equation (GEE) and family-based association tests (FBAT) models within each sex as well as sexes combined. We evaluated 70,987 autosomal SNPs with genotypic call rates > or =80%, HWE p > or = 0.001, and MAF > or =10% in up to 1141 phenotyped individuals (495 men and 646 women, mean age 62.5 yrs). Variance component linkage analysis was performed using 11,200 markers.

RESULTS

Heritability estimates for all bone phenotypes were 30-66%. LOD scores > or =3.0 were found on chromosomes 15 (1.5 LOD confidence interval: 51,336,679-58,934,236 bp) and 22 (35,890,398-48,603,847 bp) for femoral shaft section modulus. The ten primary phenotypes had 12 associations with 100K SNPs in GEE models at p < 0.000001 and 2 associations in FBAT models at p < 0.000001. The 25 most significant p-values for GEE and FBAT were all less than 3.5 x 10(-6) and 2.5 x 10(-5), respectively. Of the 40 top SNPs with the greatest numbers of significantly associated BMD traits (including femoral neck, trochanter, and lumbar spine), one half to two-thirds were in or near genes that have not previously been studied for osteoporosis. Notably, pleiotropic associations between BMD and bone geometric traits were uncommon. Evidence for association (FBAT or GEE p < 0.05) was observed for several SNPs in candidate genes for osteoporosis, such as rs1801133 in MTHFR; rs1884052 and rs3778099 in ESR1; rs4988300 in LRP5; rs2189480 in VDR; rs2075555 in COLIA1; rs10519297 and rs2008691 in CYP19, as well as SNPs in PPARG (rs10510418 and rs2938392) and ANKH (rs2454873 and rs379016). All GEE, FBAT and linkage results are provided as an open-access results resource at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007 webcite.

CONCLUSION

The FHS 100K SNP project offers an unbiased genome-wide strategy to identify new candidate loci and to replicate previously suggested candidate genes for osteoporosis.

Bivariate linkage study of proximal hip geometry and body size indices: the Framingham study

Femoral geometry and body size are both characterized by substantial heritability. The purpose of this study was to discern whether hip geometry and body size (height and body mass index, BMI) share quantitative trait loci (QTL). Dual-energy X-ray absorptiometric scans of the proximal femur from 1,473 members in 323 pedigrees (ages 31-96 years) from the Framingham Osteoporosis Study were studied. We measured femoral neck length, neck-shaft angle, subperiosteal width (outer diameter), cross-sectional bone area, and section modulus, at the narrowest section of the femoral neck (NN), intertrochanteric (IT), and femoral shaft (S) regions. In variance component analyses, genetic correlations (rho ( G )) between hip geometry traits and height ranged 0.30-0.59 and between hip geometry and BMI ranged 0.11-0.47. In a genomewide linkage scan with 636 markers, we obtained nominally suggestive linkages (bivariate LOD scores > or =1.9) for geometric traits and either height or BMI at several chromosomes (4, 6, 9, 15, and 21). Two loci, on chr. 2 (80 cM, BMI/shaft section modulus) and chr. X (height/shaft outer diameter), yielded bivariate LOD scores > or =3.0; although these loci were linked in univariate analyses with a geometric trait, neither was linked with either height or BMI. In conclusion, substantial genetic correlations were found between the femoral geometric traits, height and BMI. Linkage signals from bivariate linkage analyses of bone geometric indices and body size were similar to those obtained in univariate linkage analyses of femoral geometric traits, suggesting that most of the detected QTL primarily influence geometry of the hip.

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.

Bivariate whole genome linkage analysis for femoral neck geometric parameters and total body lean mass

A genome-wide bivariate analysis was conducted for femoral neck GPs and TBLM in a large white sample. We found QTLs shared by GPs and TBLM in the total sample and the sex-specific samples. QTLs with potential pleiotropy were also disclosed.

INTRODUCTION

Previous studies have suggested that femoral neck cross-section geometric parameters (FNCS-GPs), including periosteal diameter (W), cross-sectional area (CSA), cortical thickness (CT), buckling ratio (BR), and section modulus (Z), are genetically correlated with total body lean mass (TBLM). However, the shared genetic factors between them are unknown.

MATERIALS AND METHODS

To identify the specific QTLs shared by FNCS-GPs and TBLM, we performed bivariate whole genome linkage analysis (WGLA) in a large sample of 451 white families made up of 4498 subjects.

RESULTS

Multipoint bivariate linkage analyses for 22 autosomes showed evidence of suggestive or significant linkages (thresholds of LOD = 2.3 and 3.7, respectively) to chromosomes 3q12 and 20q13 in the entire sample, 6p25 and 10q24 in women, and 4p15, 5q34-35 and 7q21 in men. Two-point linkage analyses for chromosome X showed strong linkage to Xp22.13, Xp11.4, Xq22.3, Xq23-24, and Xq25. Complete pleiotropy was identified on 10q24 and 5q35 for TBLM and BR in women and for TBLM and CT in men, respectively. Furthermore, chromosomes 5q34-35, 7q21, 10q24, 20q13, Xp22.13, Xp11.4, and Xq25 are also of importance because of their linkage to multiple trait pairs. For example, linkage to chromosome 10q24 was found for TBLM x W (LOD = 2.31), TBLM x CT (LOD = 2.51), TBLM x CSA (LOD = 2.51), TBLM x BR (LOD = 2.64), and TBLM x Z (LOD = 2.55) in women.

CONCLUSIONS

In this study, we identified several genomic regions (e.g., 3q12 and 20q13) that seem to be linked to both FNCS-GPs and TBLM. These regions are of interesting because they may harbor genes that may contribute to variation in both FNCS-GPs and TBLM.

Proximal hip geometry is linked to several chromosomal regions: genome-wide linkage results from the Framingham Osteoporosis Study

INTRODUCTION

Femoral geometry contributes to bone strength and predicts hip fracture risk. The purpose of this study was to evaluate heritability (h(2)) of geometric indices of the proximal hip and to perform whole-genome linkage analyses of these traits, adjusted for body size.

METHODS

DXA scans of the proximal femur from 1473 members of 323 pedigrees (age range 31-96 years) from the population-based Framingham Osteoporosis Study were obtained. Using the hip structural analysis program, we measured femoral neck length (FNL, cm) and neck-shaft angle (NSA); subperiosteal width (WID, cm), cross-sectional area (CSA, cm(2)); and section modulus (Z, cm(3)) at the narrowest section of the neck (NN), intertrochanteric (IT) and femoral shaft (S) regions. Linkage analyses were performed for the above indices with a set of 636 markers using variance components maximum likelihood method.

RESULTS

Substantial genetic influences were found for all geometric phenotypes, with h(2) values between 0.28 (NSA) and 0.70 (IT_WID). Adjustment for height and BMI did not alter h(2) of NSA and FNL but decreased h(2) of the cross-sectional indices. We obtained substantial linkage (multipoint LOD >3.0) for S_Z at 2p21 and 21q11 and S_WID at Xq25-q26. Inclusion of height and BMI as covariates resulted in much lower LOD scores for S_Z, whereas linkage signals for S_Z at 4q25, S_CSA at 4q32 and S_CSA and S_Z at 15q21 increased after the adjustment. Linkage of FNL at 1q and 13q, NSA at 2q and NN_WID at 16q did not change after the adjustment.

CONCLUSION

Suggestive linkages of bone geometric indices were found at 1q, 2p, 4q, 13q, 15q and Xq. The identification of significant linkage regions after adjustment for BMI and height may point to QTLs influencing femoral bone geometry independent of body size.

Association and linkage analysis of COL1A1 and AHSG gene polymorphisms with femoral neck bone geometric parameters in both Caucasian and Chinese nuclear families

AIM

To simultaneously investigate the contribution of the alpha 1 chain of collagen type 1 (COL1A1) and alpha2-HS-glycoprotein (AHSG) genes to the variation of bone geometric parameters in both Caucasians and Chinese.

METHODS

Six hundred and five Caucasian individuals from 157 nuclear families and 1228 Chinese subjects from 400 nuclear families were genotyped at the AHSG-SacI, COL1A1- PCOL2 and Sp1 polymorphisms using polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP). 5 FN bone geometric parameters were calculated based on bone mineral density and bone area of femoral neck (FN) measured by dual energy X-ray absorptiometry. Population stratification, total family association, within-family association, and linkage tests were performed by the quantitative transmission disequilibrium test program.

RESULTS

The t-test showed the significant differences of all bone geometric phenotypes (except ED) between Caucasians and Chinese in the offspring using both unadjusted and adjusted (by age, height, weight, and gender) data. In Caucasians, we found significant within-family association results between the COL1A1-Sp1 polymorphism (rs1800012) and cross sectional area (CSA), cortical thickness (CT), endocortical diameter (ED), buckling ratio (BR) (P=0.018, 0.002, 0.023, and 0.001, respectively); the COL1A1-Sp1 polymorphism also detected significant linkage with BR (P=0.039). In the population of China, the within-family associations between the COL1A1-PCOL2 polymorphism (rs1107946) and CT, BR were significant (P=0.012 and 0.008, respectively). Furthermore, evidence of linkage were observed between the AHSG-SacI polymorphism (rs4918) and CT, BR (P=0.042 and 0.014, respectively) in Caucasians, but not in Chinese.

CONCLUSION

Our results suggest that the COL1A1 gene may have significantly association with bone geometry in both Caucasians and Chinese, and the AHSG gene may be linked to bone geometry in Caucasians, but not in Chinese. This study represents our first efforts on investigating the importance of the COL1A1 and AHSG genes on bone geometry in both Caucasians and Chinese.

Is a gene important for bone resorption a candidate for obesity? An association and linkage study on the RANK (receptor activator of nuclear factor-kappaB) gene in a large Caucasian sample

In light of findings that osteoporosis and obesity may share some common genetic determination and previous reports that RANK (receptor activator of nuclear factor-kappaB) is expressed in skeletal muscles which are important for energy metabolism, we hypothesize that RANK, a gene essential for osteoclastogenesis, is also important for obesity. In order to test the hypothesis with solid data we first performed a linkage analysis around the RANK gene in 4,102 Caucasian subjects from 434 pedigrees, then we genotyped 19 SNPs in or around the RANK gene. A family-based association test (FBAT) was performed with both a quantitative measure of obesity [fat mass, lean mass, body mass index (BMI), and percentage fat mass (PFM)] and a dichotomously defined obesity phenotype-OB (OB if BMI > or = 30 kg/m(2)). In the linkage analysis, an empirical P = 0.004 was achieved at the location of the RANK gene for BMI. Family-based association analysis revealed significant associations of eight SNPs with at least one obesity-related phenotype (P < 0.05). Evidence of association was obtained at SNP10 (P = 0.002) and SNP16 (P = 0.001) with OB; SNP1 with fat mass (P = 0.003); SNP1 (P = 0.003) and SNP7 (P = 0.003) with lean mass; SNP1 (P = 0.002) and SNP7 (P = 0.002) with BMI; SNP1 (P = 0.003), SNP4 (P = 0.007), and SNP7 (P = 0.002) with PFM. In order to deal with the complex multiple testing issues, we performed FBAT multi-marker test (FBAT-MM) to evaluate the association between all the 18 SNPs and each obesity phenotype. The P value is 0.126 for OB, 0.033 for fat mass, 0.021 for lean mass, 0.016 for BMI, and 0.006 for PFM. The haplotype data analyses provide further association evidence. In conclusion, for the first time, our results suggest that RANK is a novel candidate for determination of obesity.

HDC gene polymorphisms are associated with age at natural menopause in Caucasian women

Histidine decarboxylase gene (HDC) encodes histidine decarboxylase which is the crucial enzyme for the biosynthesis of histidine. Studies have shown that histamine is likely to be involved in the regulation of reproduction system. To find the possible correlation between HDC gene and AANM (age at natural menopause), we selected 265 postmenopausal women from 131 nuclear families and performed a transmission disequilibrium test. Significant within-family associations with AANM for SNP rs854163 and SNP rs854158 of HDC gene were observed (P values=0.0018 and 0.0197, respectively). After 1000 permutations, SNP rs854163 still remained significant within-family association with AANM. Consistently, we also detected a significant within-family association between haplotype block 2 (defined by SNP rs854163 and rs860526) and AANM in the haplotype analyses (P value=0.0397). Our results suggest that the HDC gene polymorphisms are significantly associated with AANM in Caucasian women.

Genetic and environmental correlations between bone geometric parameters and body compositions

The purpose of this study was to determine the genetic and environmental correlations between weight, lean mass and bone geometric parameters (sub-periosteal diameter, W; cross-sectional area, CSA; cortical thickness, CT; section modulus, Z; and buckling ratio, BR) of femoral neck. The sample was composed of 512 Caucasian pedigrees, including 2667 females and 1822 males. Bivariate quantitative genetic analyses were performed to evaluate the genetic (rho(G)), environmental (rho(E)) and phenotypic (rho(P)) correlations between the study traits. Univariate genetic analyses showed that the heritabilities (h(2)) for bone geometric parameters were significant (P < 0.001) ranging from 0.50 to 0.60. The significant common household effects indicated the common environment shared by household members for W, CSA, CT, Z and BR (P < 0.05), but the magnitude was small compared with heritabilities. rho(E), rho(G) and rho(P) between bone geometric parameters and weight, lean mass were generally significant. Interestingly, lean mass showed both stronger genetic and environmental correlations with the bone geometric parameters than weight. In addition, according to the magnitude of correlation coefficients, the rho(G) between body compositions and bone geometric parameters were generally stronger than rho(E) (except for that between BR and body compositions). These data suggested that the geometric parameters of femoral neck are under strong genetic control. Furthermore, some common genetic and environmental factors are shared by bone geometric parameters and weight, lean mass. The results may help understand the intertwined relationships between bone metabolisms, mechanical loading and body compositions.