A genomewide scan for quantitative trait loci underlying areal bone size variation in 451 Caucasian families

Genome-wide association study identifies HMGN3 locus for spine bone size variation in Chinese

Bone size (BS) is one of the major risk factors for osteoporotic fractures. BS variation is genetically determined to a substantial degree with heritability over 50%, but specific genes underlying variation of BS are still largely unknown. To identify specific genes for BS in Chinese, initial genome-wide association scan (GWAS) study and follow-up replication study were performed. In initial GWAS study, a group of 12 contiguous single-nucleotide polymorphism (SNP)s, which span a region of ~25 kb and locate at the upstream of HMGN3 gene (high-mobility group nucleosomal binding domain 3), achieved moderate association signals for spine BS, with P values ranging from 6.2E-05 to 1.8E-06. In the follow-up replication study, eight of the 12 SNPs were detected suggestive replicate associations with BS in 1,728 unrelated female Caucasians, which have well-known differences from Chinese in ethnic genetic background. The SNPs in the region of HMGN3 gene formed a tightly combined haplotype block in both Chinese and Caucasians. The results suggest that the genomic region containing HMGN3 gene may be associated with spine BS in Chinese.

PLCL1 rs7595412 variation is not associated with hip bone size variation in postmenopausal Danish women

BACKGROUND

Bone size (BS) variation is under strong genetic control and plays an important role in determining bone strength and fracture risk. Recently, a genome-wide association study identified polymorphisms associated with hip BS variation in the PLCL1 (phospholipase c-like 1) locus. Carriers of the major A allele of the most significant polymorphism, rs7595412, have around 17% larger hip BS than non-carriers. We therefore hypothesized that this polymorphism may also influence postmenopausal complications.

METHODS

The effects of rs7595412 on hip BS, bone mineral density (BMD), vertebral fractures, serum Crosslaps and osteocalcin levels were analyzed in 1,191 postmenopausal Danish women.

RESULTS

This polymorphism had no influence on hip and spine BS as well as on femur and spine BMD. Women carrying at least one copy of the A allele had lower levels of serum osteocalcin as compared with those homozygous for the G allele (p = 0.03) whereas no effect on serum Crosslaps was detected. Furthermore, women homozygous for the A allele were more affected by vertebral fractures than those carrying at least one copy of the G allele (p = 0.04).

CONCLUSIONS

In postmenopausal women, our results suggest that the PLCL1 rs7595412 polymorphism has no obvious effect on hip BS or BMD but may be nominally associated with increased proportion of vertebral fracture and increased levels of osteocalcin.

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.

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 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.

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.

Evidence for major pleiotropic effects on bone size variation from a principal component analysis of 451 Caucasian families

AIM

To identify pleiotropic quantitative trait loci (QTL) influencing bone size (BS) at different skeletal sites in Caucasians.

METHODS

In a sample containing 3899 Caucasians from 451 pedigrees, 410 microsatellite markers spaced approximately 8.9 cM apart across the human genome were genotyped. Phenotypical and genetic correlations of BS at lumbar spine, hip (femoral neck, trochanter, and intertrochanter regions), and wrist (ultradistal, mid-distal, and one-third distal sites) were determined using bivariate quantitative genetic analysis. A principal component analysis (PCA) was performed to obtain principal component (PC) factors that were then subjected to variance components linkage analysis to identify regions linked to the PC.

RESULTS

Genetic correlations of BS at different skeletal sites ranged from 0.40 to 0.79 (P<0.001). The PCA yielded a PC named PCtotal, which explained up to 76% of the total (co)variation of all the BS at the 7 skeletal sites for the whole sample. We identified a QTL influencing the BS of multiple skeletal sites on chromosome 7 at 140 cM [logarithm of odds (LOD)=2.85] in the overall sample. Sex-specific evidence for linkage was observed on chromosome 11 at 53 cM (LOD =2.82) in the male-only data subset.

CONCLUSION

Our study identified several genomic regions that may have pleiotropic effects on different skeletal sites. These regions may contain genes that play a critical role in overall bone development and osteoporosis at multiple skeletal sites, hence are biologically and clinically important.

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.

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.