Quantitative trait loci on chromosomes 2p, 4p, and 13q influence bone mineral density of the forearm and hip in Mexican Americans

20(S)-Protopanaxadiol Inhibits Titanium Particle-Induced Inflammatory Osteolysis and RANKL-Mediated Osteoclastogenesis via MAPK and NF-κB Signaling Pathways

Osteolysis is a principal reason for arthroplasty failure like aseptic loosening induced by Titanium (Ti) particle. It is a challenge for orthopedic surgeons. Recent researches show that 20(S)-protopanaxadiol can inhibit inflammatory cytokine release in vitro. This study aims to assess the effect of 20(S)-protopanaxadiol on Ti particle-induced osteolysis and RANKL-mediated osteoclastogenesis. Micro-CT and histological analysis in vivo indicated the inhibitory effects of 20(S)-protopanaxadiol on osteoclastogenesis and the excretion of inflammatory cytokines. Next, we demonstrated that 20(S)-protopanaxadiol inhibited osteoclast differentiation, bone resorption area, and F-actin ring formation in a dose-dependent manner. Moreover, mechanistic studies suggested that the suppression of MAPK and NF-κB signaling pathways were found to mediate the inhibitory effects of 20(S)-protopanaxadiol. In conclusion, 20(S)-protopanaxadiol may suppress osteoclastogenesis in a dose- dependent manner and it could be a potential treatment of Ti particle-induced osteolysis.

Lithium-containing biomaterials inhibit osteoclastogenesis of macrophagesin vitroand osteolysisin vivo

Li-containing bioceramics were promising biomaterials for inhibiting osteoclastogenesis of macrophages and osteolysisin vivo, potentially using for treating osteoporosis.

Heritability of bone mineral density in a multivariate family-based study

There is evidence for a genetic contribution to bone mineral density (BMD×). Different loci affecting BMD have been identified by diverse linkage and genome-wide association studies. We studied the heritability of and the correlations among six densitometric phenotypes and four bone mass/fracture phenotypes. For this purpose, we used a family-based study of the genetics of osteoporosis, the Genetic Analysis of Osteoporosis Project. The primary aim of our study was to examine the roles of genetic and environmental factors in determining osteoporosis-related phenotypes. The project consisted of 11 extended families from Spain. All of them were selected through a proband with osteoporosis. BMD was measured using dual-energy X-ray absorptiometry. The proportion of variance of BMD attributable to significant covariates ranged from 25% (for femoral neck BMD) to 48% (for whole-body total BMD). The vast majority of the densitometric phenotypes had highly significant heritability, ranging from 0.252 (whole-body total BMD) to 0.537 (trochanteric BMD) after correcting for covariate effects. All of the densitometric phenotypes showed high and significant genetic correlations (from -0.772 to -1.000) with a low bone mass/osteopenia condition (Affected 3). Our findings provide additional evidence on the heritability of BMD and a strong genetic correlation between BMD and bone mass/fracture phenotypes in a Spanish population. Our results emphasize the importance of detecting genetic risk factors and the benefit of early diagnosis and especially therapeutic and preventive strategies.

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.

Arguments for the sake of endophenotypes: examining common misconceptions about the use of endophenotypes in psychiatric genetics

Endophenotypes are measurable biomarkers that are correlated with an illness, at least in part, because of shared underlying genetic influences. Endophenotypes may improve our power to detect genes influencing risk of illness by being genetically simpler, closer to the level of gene action, and with larger genetic effect sizes or by providing added statistical power through their ability to quantitatively rank people within diagnostic categories. Furthermore, they also provide insight into the mechanisms underlying illness and will be valuable in developing biologically-based nosologies, through efforts such as RDoC, that seek to explain both the heterogeneity within current diagnostic categories and the overlapping clinical features between them. While neuroimaging, electrophysiological, and cognitive measures are currently most used in psychiatric genetic studies, researchers currently are attempting to identify candidate endophenotypes that are less genetically complex and potentially closer to the level of gene action, such as transcriptomic and proteomic phenotypes. Sifting through tens of thousands of such measures requires automated, high-throughput ways of assessing, and ranking potential endophenotypes, such as the Endophenotype Ranking Value. However, despite the potential utility of endophenotypes for gene characterization and discovery, there is considerable resistance to endophenotypic approaches in psychiatry. In this review, we address and clarify some of the common issues associated with the usage of endophenotypes in the psychiatric genetics community.

Genome-wide approaches for identifying genetic risk factors for osteoporosis

Osteoporosis, the most common type of bone disease worldwide, is clinically characterized by low bone mineral density (BMD) and increased susceptibility to fracture. Multiple genetic and environmental factors and gene-environment interactions have been implicated in its pathogenesis. Osteoporosis has strong genetic determination, with the heritability of BMD estimated to be as high as 60%. More than 80 genes or genetic variants have been implicated in risk of osteoporosis by hypothesis-free genome-wide studies. However, these genes or genetic variants can only explain a small portion of BMD variation, suggesting that many other genes or genetic variants underlying osteoporosis risk await discovery. Here, we review recent progress in genome-wide studies of osteoporosis and discuss their implications for medicine and the major challenges in the field.

Single nucleotide polymorphisms in sFRP4 are associated with bone and body composition related parameters in Danish but not in Belgian men

The senescence accelerated mouse P6 (SAMP6) has a low bone mass and has previously shown to be a good model for senile osteoporosis in humans. In addition to a reduced bone mass, SAMP6 mice are obese and have hyperlipidemia. Using positional cloning and expression studies, an increased expression of sfrp4 was found in these mice. SFRP4 is a modulator of the Wnt signalling pathway. This pathway has been previously shown to be involved in regulating bone mass. Additional evidence that sFRP4 has an influence on BMD was delivered by linkage and association studies mostly performed in Asian populations. Based on these data we decided to perform an association study between common variants in sFRP4, BMD, hip geometry parameters and body composition parameters in a population consisting of 1383 Danish men (783 aged 20-29 years; 600 aged 60-74 years). Afterwards we tried to replicate the significant results in a population of 994 Belgian men. In the Danish population we found 6 SNPs associated with BMD at the hip and/or femoral neck. Furthermore, all 6 SNPs were associated with several hip geometry parameters. The homozygous presence of the minor allele resulted for all SNPs (except rs4720265) in a decrease in bone density and bone strength. Finally, we observed in the Danish population age specific associations with height and fat mass. In the Belgian population we tried to replicate the results of three SNPs with BMD and body composition parameters. Unfortunately, we were not able to replicate the results found in the Danish cohort but we found one SNP (rs2598116) associated with height. In conclusion, genetic variation in sFRP4 has an influence on hip fracture risk, percentage body fat and height in a Danish male population. However, we were unable to replicate these results in an independent Belgian population.

Investigation of the genetic architecture of a bone carcass weight QTL on BTA6

A previous analysis of an F(2) /Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone-related trait, the sum of the bone weight for all the left fore- and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA-119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone-related QTL.

A functional haplotype in EIF2AK3, an ER stress sensor, is associated with lower bone mineral density

EIF2AK3 is a type I transmembrane protein that functions as an endoplasmic reticulum (ER) stress sensor to regulate global protein synthesis. Rare mutations in EIF2AK3 cause Wolcott-Rallison syndrome (OMIM 226980), an autosomal recessive disorder characterized by diabetes, epiphyseal dysplasia, osteoporosis, and growth retardation. To investigate the role of common genetic variation in EIF2AK3 as a determinant of bone mineral density (BMD) and osteoporosis, we sequenced all exons and flanking regions, then genotyped six potentially functional single nucleotide polymorphisms (SNPs) in this gene in 997 Amish subjects for association analysis, and attempted replication in 887 Mexican Americans. We found that the minor allele of a nonsynonymous SNP rs13045 had borderline associations with decreased forearm BMD in both discovery and replication cohorts (unadjusted p = 0.036 and β = -0.007 for the Amish; unadjusted p = 0.031 and β = -0.008 for Mexican Americans). A meta-analysis indicated this association achieved statistical significance in the combined sample (unadjusted p = 0.003; Bonferroni corrected p = 0.009). Rs13045 and three other potentially functional SNPs, a promoter SNP (rs6547787) and two nonsynonymous SNPs (rs867529 and rs1805165), formed two haplotypes: a low-BMD associated haplotype, denoted haplotype B [minor allele frequency (MAF) = 0.311] and a common haplotype A (MAF = 0.676). There were no differences in mRNA expression in lymphoblastoid cell lines between the two haplotypes. However, after treating lymphoblastoid cell lines with thapsigargin to induce ER stress, cell lines with haplotype B showed increased sensitivity to ER stress (p = 0.014) compared with cell lines with haplotype A. Taken together, our results suggest that common nonsynonymous sequence variants in EIF2AK3 have a modest effect on ER stress response and may contribute to the risk for low BMD through this mechanism.

High dimensional endophenotype ranking in the search for major depression risk genes

BACKGROUND

Despite overwhelming evidence that major depression is highly heritable, recent studies have localized only a single depression-related locus reaching genome-wide significance and have yet to identify a causal gene. Focusing on family-based studies of quantitative intermediate phenotypes or endophenotypes, in tandem with studies of unrelated individuals using categorical diagnoses, should improve the likelihood of identifying major depression genes. However, there is currently no empirically derived statistically rigorous method for selecting optimal endophentypes for mental illnesses. Here, we describe the endophenotype ranking value, a new objective index of the genetic utility of endophenotypes for any heritable illness.

METHODS

Applying endophenotype ranking value analysis to a high-dimensional set of over 11,000 traits drawn from behavioral/neurocognitive, neuroanatomic, and transcriptomic phenotypic domains, we identified a set of objective endophenotypes for recurrent major depression in a sample of Mexican American individuals (n = 1122) from large randomly selected extended pedigrees.

RESULTS

Top-ranked endophenotypes included the Beck Depression Inventory, bilateral ventral diencephalon volume, and expression levels of the RNF123 transcript. To illustrate the utility of endophentypes in this context, each of these traits were utlized along with disease status in bivariate linkage analysis. A genome-wide significant quantitative trait locus was localized on chromsome 4p15 (logarithm of odds = 3.5) exhibiting pleiotropic effects on both the endophenotype (lymphocyte-derived expression levels of the RNF123 gene) and disease risk.

CONCLUSIONS

The wider use of quantitative endophenotypes, combined with unbiased methods for selecting among these measures, should spur new insights into the biological mechanisms that influence mental illnesses like major depression.

Cortical bone health shows significant linkage to chromosomes 2p, 3p, and 17q in 10-year-old children

Genes play an important role in lifelong skeletal health. Genes that influence bone building during childhood have the potential to affect bone health not only throughout childhood but also into adulthood. Given that peak bone mass is a significant predictor of adult fracture risk, it is imperative that the genetic underpinnings of the normal pediatric skeleton are uncovered. In a sample of 600 10-year-old children from 144 families in the Fels Longitudinal Study, we examined radiographic cortical bone measures of the second metacarpal. Morphometic measurements included bone width, medial and lateral cortical thicknesses, and the calculated cortical index representing the amount of cortex relative to bone width. We then conducted genome-wide linkage analysis on these traits in 440 genotyped individuals using the SOLAR analytic platform. Significant quantitative trait loci (QTL) were identified for bone traits on three separate chromosomes. A QTL for medial cortical thickness was localized to chromosome 2p25.2. A QTL for lateral cortical thickness was localized to chromosomal region 3p26.1-3p25.3. Finally, a QTL detected for cortical index was localized to the 17q21.2 chromosomal region. Each region contains plausible candidate genes for pediatric skeletal health, some of which confirm findings from studies of adulthood bone, and for others represent novel candidate genes for skeletal health.

Genetic analysis of cortical thickness and fractional anisotropy of water diffusion in the brain

OBJECTIVES

The thickness of the brain's cortical gray matter (GM) and the fractional anisotropy (FA) of the cerebral white matter (WM) each follow an inverted U-shape trajectory with age. The two measures are positively correlated and may be modulated by common biological mechanisms. We employed four types of genetic analyses to localize individual genes acting pleiotropically upon these phenotypes.

METHODS

Whole-brain and regional GM thickness and FA values were measured from high-resolution anatomical and diffusion tensor MR images collected from 712, Mexican American participants (438 females, age = 47.9 ± 13.2 years) recruited from 73 (9.7 ± 9.3 individuals/family) large families. The significance of the correlation between two traits was estimated using a bivariate genetic correlation analysis. Localization of chromosomal regions that jointly influenced both traits was performed using whole-genome quantitative trait loci (QTL) analysis. Gene localization was performed using SNP genotyping on Illumina 1M chip and correlation with leukocyte-based gene-expression analyses. The gene-expressions were measured using the Illumina BeadChip. These data were available for 371 subjects.

RESULTS

Significant genetic correlation was observed among GM thickness and FA values. Significant logarithm of odds (LOD ≥ 3.0) QTLs were localized within chromosome 15q22-23. More detailed localization reported no significant association (p < 5·10(-5)) for 1565 SNPs located within the QTLs. Post hoc analysis indicated that 40% of the potentially significant (p ≤ 10(-3)) SNPs were localized to the related orphan receptor alpha (RORA) and NARG2 genes. A potentially significant association was observed for the rs2456930 polymorphism reported as a significant GWAS finding in Alzheimer's disease neuroimaging initiative subjects. The expression levels for RORA and ADAM10 genes were significantly (p < 0.05) correlated with both FA and GM thickness. NARG2 expressions were significantly correlated with GM thickness (p < 0.05) but failed to show a significant correlation (p = 0.09) with FA.

DISCUSSION

This study identified a novel, significant QTL at 15q22-23. SNP correlation with gene-expression analyses indicated that RORA, NARG2, and ADAM10 jointly influence GM thickness and WM-FA values.

Whole brain and regional hyperintense white matter volume and blood pressure: overlap of genetic loci produced by bivariate, whole-genome linkage analyses

BACKGROUND AND PURPOSE

The volume of T2-hyperintense white matter (HWM) is an important neuroimaging marker of cerebral integrity with a demonstrated high heritability. Pathophysiology studies have shown that the regional, ependymal, and subcortical HWM lesions are associated with elevated arterial pulse pressure and arterial blood pressure (BP), respectively. We performed bivariate, whole-genome linkage analyses for HWM volumes and BP measurements to identify chromosomal regions that contribute jointly to both traits in a population of healthy Mexican Americans. Our aims were to localize novel quantitative trait loci acting pleiotropically on these phenotypes and to replicate previous genetic findings on whole brain HWM volume and BP measurements.

METHODS

BP measurements and volumes of whole-brain (WB), subcortical, and ependymal HWM lesions, measured from high-resolution (1 mm(3)) 3-dimensional fluid-attenuated inversion recovery images, served as focal quantitative phenotypes. Data were collected from 357 (218 females; mean age=47.9±13.2 years) members of large extended families who participated in the San Antonio Family Heart Study.

RESULTS

Bivariate genomewide linkage analyses localized a significant quantitative trait locus influencing WB and regional (ependymal) HWM volumes and pulse pressure and systolic BP to chromosomal location 1q24 between markers D1S196 and D1S1619. Several other chromosomal regions (1q42, 10q24-q26, and 15q26) exhibited suggestive linkages. The results of the post hoc analyses that excluded 55 subjects taking antihypertensive medication showed no substantive differences from the results obtained in the full cohort.

CONCLUSIONS

This study confirms several previously observed quantitative trait loci influencing BP and cerebral integrity and identifies a novel significant quantitative trait locus at chromosome 1q24. The genetic results strongly support a role for pleiotropically acting genes jointly influencing BP and cerebral white matter integrity.

A common variant in fibroblast growth factor binding protein 1 (FGFBP1) is associated with bone mineral density and influences gene expression in vitro

We previously detected strong evidence for linkage of forearm bone mineral density (BMD) to chromosome 4p (lod=4.3) in a set of 29 large Mexican American families. Fibroblast growth factor binding protein 1 (FGFBP1) is a strong candidate gene for bone homeostasis in this region. We sequenced the coding region of FGFBP1 in a subset of our Mexican American study population and performed association studies with BMD on SNPs genotyped in the entire cohort. We then attempted to replicate these findings in an independent study cohort and performed in vitro functional studies on replicated, potentially functional polymorphisms using a luciferase reporter construct to evaluate influence on gene expression. Several SNPs spanning the gene, all in one large block of linkage disequilibrium, were significantly associated with BMD at various skeletal sites (n=872, p=0.001-0.04). The associations were then replicated in an independent population of European ancestry (n=972; p=0.02-0.04). Sex-stratified association analyses in both study populations suggest this association is much stronger in men. Subsequent luciferase reporter gene assays revealed marked differences in FGFBP1 expression among the three common haplotypes. Further experiments revealed that a promoter polymorphism, rs12503796, results in decreased expression of FGFBP1 and inhibits upregulation of the gene by testosterone in vitro. Collectively, these findings suggest that sequence variation in FGFBP1 may contribute to variation in BMD, possibly influencing osteoporosis risk.

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.

Genetics of microstructure of cerebral white matter using diffusion tensor imaging

We analyzed the degree of genetic control over intersubject variability in the microstructure of cerebral white matter (WM) using diffusion tensor imaging (DTI). We performed heritability, genetic correlation and quantitative trait loci (QTL) analyses for the whole-brain and 10 major cerebral WM tracts. Average measurements for fractional anisotropy (FA), radial (L( perpendicular)) and axial (L( vertical line)) diffusivities served as quantitative traits. These analyses were done in 467 healthy individuals (182 males/285 females; average age 47.9+/-13.5 years; age range: 19-85 years), recruited from randomly-ascertained pedigrees of extended families. Significant heritability was observed for FA (h(2)=0.52+/-0.11; p=10(-7)) and L( perpendicular) (h(2)=0.37+/-0.14; p=0.001), while L( vertical line) measurements were not significantly heritable (h(2)=0.09+/-0.12; p=0.20). Genetic correlation analysis indicated that the FA and L( perpendicular) shared 46% of the genetic variance. Tract-wise analysis revealed a regionally diverse pattern of genetic control, which was unrelated to ontogenic factors, such as tract-wise age-of-peak FA values and rates of age-related change in FA. QTL analysis indicated linkages for whole-brain average FA (LOD=2.36) at the marker D15S816 on chromosome 15q25, and for L( perpendicular) (LOD=2.24) near the marker D3S1754 on the chromosome 3q27. These sites have been reported to have significant co-inheritance with two psychiatric disorders (major depression and obsessive-compulsive disorder) in which patients show characteristic alterations in cerebral WM. Our findings suggest that the microstructure of cerebral white matter is under a strong genetic control and further studies in healthy as well as patients with brain-related illnesses are imperative to identify the genes that may influence cerebral white matter.

Analysis of genetic variability and whole genome linkage of whole-brain, subcortical, and ependymal hyperintense white matter volume

BACKGROUND AND PURPOSE

The cerebral volume of T2-hyperintense white matter (HWM) is an important neuroimaging marker of cerebral integrity. Pathophysiology studies identified that subcortical and ependymal HWM are produced by 2 different mechanisms but shared a common risk factor: high arterial pulse pressure. Recent studies have demonstrated high heritability of the whole-brain HMW volume and reported significant and suggestive evidence of genetic linkage. We performed heritability and whole-genome linkage analysis to replicate previous reported findings and to study shared genetic variance, and possible overlap for specific loci, between subcortical and ependymal HWM volumes in a population of healthy Mexican Americans.

METHODS

The volumes of subcortical and ependymal HWM regions were measured from high-resolution (1 mm(3)), 3-dimensional fluid-attenuated inversion recovery images acquired for 459 (283 females, 176 males) active participants in the San Antonio Family Heart Study. Subjects ranged in age from 19 to 85 years of age (47.9+/-13.5 years) and were part of 49 families (9.4+/-8.5 individuals per family).

RESULTS

The volumes of whole-brain, subcortical, and ependymal HWM were highly heritable (h(2)=0.72, 0.66, and 0.73, respectively). The subcortical and ependymal HWM volumes shared 21% of genetic variability indicating significant pleiotropy. Genomewide linkage analysis showed only a suggestive bivariate linkage for subcortical and ependymal HWM volumes (log of odds=2.12) on chromosome 1 at 288 cM.

CONCLUSIONS

We replicated previous findings of high heritability for the whole-brain HWM volume. We also showed that subcortical and ependymal volume shared a significant portion of genetic variability and the bivarate linkage analysis produced a suggestive linkage near the locus previously identified in a study of whole-brain HWM volume and arterial pulse pressure.

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.

New roles for INSL3 in adults

Insulin-like factor 3 (INSL3) is produced primarily by testicular Leydig cells and acts by binding to its specific G-protein-coupled receptor, RXFP2 (relaxin family peptide 2). INSL3 is involved in testicular descent during fetal development, and mutations in the INSL3 and RXFP2 genes cause cryptorchidism. The physiological role of INSL3 in adults is not known, although substantial INSL3 circulating levels are present. After extensive clinical, biochemical, and hormonal investigations, including bone densitometry with dual energy X-ray absorptiometry, on 25 young men (age, 27-41 years) who have the well-characterized T222P mutation in the RXFP2 gene, we found that 16 of them (64%) had significantly reduced bone density. No other cause of osteoporosis was evident in these subjects, whose testosterone and gonadal function were normal. Expression analysis of INSL3 and RXFP2 on human bone biopsy and human and mouse osteoblast cell cultures performed by reverse transcription-PCR and immunohistochemistry showed the presence of RXFP2 in these cells. Real-time cAMP imaging analysis and proliferation assays under the stimulus of INSL3 showed a dose- and time-dependent increase in cAMP and cell proliferation, and specific osteoblast gene activation was observed by real-time PCR after INSL3 stimulation. Lumbar spine and femoral bone of Rxfp2-deficient mice were studied by static and dynamic histomorphometry and micro-computed tomography, respectively, and showed decreased bone mass, mineralizing surface, bone formation, and osteoclast surface compared to wild-type littermates, compatible with a functional osteoblast impairment. This study identified for the first time a role for INSL3 in adults, demonstrating a modulating effect on bone metabolism and linking RXFP2 gene mutations with human osteoporosis.

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.

Quantitative trait loci, genes, and polymorphisms that regulate bone mineral density in mouse

This is an in silico analysis of data available from genome-wide scans. Through analysis of QTL, genes and polymorphisms that regulate BMD, we identified 82 BMD QTL, 191 BMD-associated (BMDA) genes, and 83 genes containing known BMD-associated polymorphisms (BMDAP). The catalogue of all BMDA/BMDAP genes and relevant literatures are provided. In total, there are substantially more BMDA/BMDAP genes in regions of the genome where QTL have been identified than in non-QTL regions. Among 191 BMDA genes and 83 BMDAP genes, 133 and 58 are localized in QTL regions, respectively. The difference was still noticeable for the chromosome distribution of these genes between QTL and non-QTL regions. These results have allowed us to generate an integrative profile of QTL, genes, polymorphisms that determine BMD. These data could facilitate more rapid and comprehensive identification of causal genes underlying the determination of BMD in mouse and provide new insights into how BMD is regulated in humans.

Choice of study phenotype in osteoporosis genetic research

Osteoporosis is a common disease characterized by compromised bone strength predisposing a person to an increased risk of osteoporotic fracture (OF). Recently, extensive efforts have been made to identify candidate genes underlying osteoporosis by the use of surrogate phenotypes, such as bone mineral density (BMD) and bone geometry. Among them, BMD is a suitable choice if we aim to classify the role of biological pathways for bone strength and to understand the bone conditions in the development of osteoporosis. However, evidences show that the genetic correlation between BMD and OF is very limited. In this review, we are mainly concerned with an important issue, i.e., phenotype choice in osteoporosis genetic research. For clarity, we address this issue with several arguments, and comments are made on most representative literature.

Genetic influences on bone loss in the San Antonio Family Osteoporosis study

The genetic contribution to age-related bone loss is not well understood. We estimated that genes accounted for 25-45% of variation in 5-year change in bone mineral density in men and women. An autosome-wide linkage scan yielded no significant evidence for chromosomal regions implicated in bone loss.

INTRODUCTION

The contribution of genetics to acquisition of peak bone mass is well documented, but little is known about the influence of genes on subsequent bone loss with age. We therefore measured 5-year change in bone mineral density (BMD) in 300 Mexican Americans (>45 years of age) from the San Antonio Family Osteoporosis Study to identify genetic factors influencing bone loss.

METHODS

Annualized change in BMD was calculated from measurements taken 5.5 years apart. Heritability (h(2)) of BMD change was estimated using variance components methods and autosome-wide linkage analysis was carried out using 460 microsatellite markers at a mean 7.6 cM interval density.

RESULTS

Rate of BMD change was heritable at the forearm (h(2) = 0.31, p = 0.021), hip (h(2) = 0.44, p = 0.017), spine (h(2) = 0.42, p = 0.005), but not whole body (h(2) = 0.18, p = 0.123). Covariates associated with rapid bone loss (advanced age, baseline BMD, female sex, low baseline weight, postmenopausal status, and interim weight loss) accounted for 10% to 28% of trait variation. No significant evidence of linkage was observed at any skeletal site.

CONCLUSIONS

This is one of the first studies to report significant heritability of BMD change for weight-bearing and non-weight-bearing bones in an unselected population and the first linkage scan for change in BMD.

Human genetics of osteoporosis

A family history of hip fracture carries a twofold increased risk of fracture among descendants. Genetic factors indeed play a major role in the determination of bone mineral density (BMD) and osteoporosis risk. Multiple chromosomal loci have been mapped by linkage approaches which potentially carry hundreds of genes involved in the determination of bone mass and quality. Association studies based on candidate gene polymorphisms and subsequent meta-analyses, and the more recent genome-wide association studies (GWAS), have clearly identified a handful of genes associated with BMD and/or fragility fractures. Among them are genes coding for the LDL-receptor related protein 5 (LRP5), estrogen receptor alpha (ESR1) and osteoprotegerin, OPG (TNFRSf11b). However, the percentage of osteoporosis risk explained by any of these polymorphisms is small, indicating that most genetic risk factors remain to be discovered and/or that interaction with environmental factors needs further consideration.

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.

Vitamin D-binding protein gene microsatellite polymorphism influences BMD and risk of fractures in men

Here we report the results of a vitamin D-binding protein gene microsatellite polymorphism study in 170 men, comprising healthy male subjects and men with osteoporosis-related symptomatic vertebral fractures. We confirm the results of an earlier study in a different cohort, showing relationship between certain genotypes of (TAAAn)-Alu repeats and reduced BMD and vertebral fractures.

INTRODUCTION

Vitamin D-binding protein (DBP) plays a critical role in the transport and metabolism of metabolites of vitamin D, including the key calciotropic hormone 1alpha,25-dihydroxyvitamin D3 (1,25(OH)2D3).

METHODS

We have investigated intra-intronic variable tandem (TAAA)n-Alu repeat expansion in the DBP gene in 170 men, comprising healthy male subjects and men with idiopathic osteoporosis and low trauma fractures.

RESULTS AND CONCLUSIONS

The predominant DBP-Alu genotype in the control subjects was 10/10 (frequency 0.421), whereas the frequency of this genotype in men with osteoporosis was 0.089. DBP-Alu alleles *10, *8 and *9, respectively, were the three commonest in both healthy subjects and men with osteoporosis. Allele *10 was associated with a lower risk of osteoporosis (OR 0.39, 95% CI 0.25-0.64; p < 0.0005), as was allele *11 (odds ratio 0.09, 95% CI 0.01-0.67; p < 0.007). Logistic regression gave similar results, showing that individuals with genotype 10/10 and 19-20 repeats (genotypes 9/10, 9/11, 10/10,) are protected from fracture or osteoporosis. Overall, there was a relationship between DBP Alu genotype and BMD, suggesting that DBP-Alu genotype may influence fracture risk. This effect may be mediated by changes in the circulating concentrations of DBP which influences free concentrations of vitamin D.

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.

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.

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.

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.

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.

The association of common polymorphisms in the QPCT gene with bone mineral density in the Chinese population

Evidence of the linkage of chromosome 2p to bone mineral density (BMD) has previously been reported in multiple populations. However, the identification of the BMD quantitative trait loci (QTL) gene at chromosome 2p remains a challenge. We performed a gene-wide and tag single nucleotide polymorphism (SNP)-based association study of four positional and functional candidate genes (CALM2, CYP1B1, QPCT, and POMC) in a sample of 1,243 cases and matched controls. Thirteen HapMap tag SNPs were selected and genotyped by using the high-throughput Sequenom genotyping platform. Binary logistic regression analyses were performed to test for associations between each SNP genotype and BMD. Haplotype association analyses were performed by WHAP. The rs3770748 within the QPCT gene showed a significant association with spine BMD in both single-marker (P = 0.002) and haplotype association analyses (P = 0.0482 for the global test; P = 0.00092 for the haplotype-specific test). Subgroup analysis revealed that the effect was primarily driven by an association in the postmenopausal women, presumably suggesting that the rs3770748 affects postmenopausal bone loss rather than peak bone mass. Our results suggest that QPCT may be the QTL gene at chromosome 2p for spine BMD variation in the Chinese population.

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.

Identification of mouse Duffy antigen receptor for chemokines (Darc) as a BMD QTL gene

It is now well known that bone mineral density (BMD) variance is determined by both genetic and environmental factors. Accordingly, studies in human and animal models have revealed evidence for the presence of several quantitative trait loci (QTL) that contribute to BMD variations. However, the identification of BMD QTL genes remains a big challenge. In the current study, we focused our efforts to identify the BMD candidate gene in chromosome 1 (Chr 1) QTL that was detected from a cross involving high BMD CAST/EiJ (CAST) and low BMD C57BL/6J (B6) mice. To this end, we have combined several approaches including: (1) fine mapping the BMD QTL in Chr 1 of the B6.CAST F2 female mice using a large number of polymorphic markers; (2) the generation of congenic sublines of mice by repeated backcrossing of CAST with B6 mice and phenotype characterization; (3) expression profiling genes in the QTL region; and (4) SNP analyses to identify the mouse Duffy Antigen Receptor for Chemokines (Darc) as a candidate gene for Chr 1 BMD QTL2. We verified the involvement of the Darc protein in BMD variation by evaluating the skeletal phenotype of Darc-knockout mice and congenic sublines of mice carrying small chromosomal segments from CAST BMD QTL. Based on the findings that Darc-antibody blocked formation of multinucleated osteoclasts in vitro and that Darc from CAST binds chemokines, known to regulate osteoclast formation, with reduced affinity compared with Darc from B6 mice, we conclude that Darc regulates BMD negatively by increasing osteoclast formation, and that the genetic association between Darc gene polymorphism and BMD variations in humans merits investigation.

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.

Linkage to chromosome 11p12 in two Maltese families with a highly penetrant form of osteoporosis

Osteoporosis is a metabolic bone disease with a strong genetic component. Family-based linkage studies were performed by a number of investigators to try to identify loci that might contain genes responsible for an increased susceptibility to osteoporosis. A whole-genome linkage scan using 400 microsatellite markers was performed in 27 members from two Maltese families with a highly penetrant form of osteoporosis. The phenotype was defined by lumbar and femoral z-scores calculated after measurement of bone mineral density by DEXA. Both males and females were among the affected individuals. Multipoint parametric and non-parametric linkage analyses were performed by EasyLinkage v4.01 using GENEHUNTER v2.1, assuming dominant and recessive modes of inheritance with variable penetrance. Evidence of linkage was observed to a marker at 11p12 where a non-parametric LOD score of 5.77 (P=0.0006) was obtained. A maximum heterogeneity LOD score of 2.55 for this region was obtained for the dominant mode of inheritance with 90% penetrance and a phenocopy rate of 1%. Following fine mapping, the critical interval was narrowed to a region that is 52.94 cM from 11p-telomere. In this region, the gene for tumour necrosis factor receptor-associated factor 6 (TRAF6) is located approximately 1 cM away from the indicated marker. Sequencing of the promoter region and exons of the TRAF6 gene revealed three sequence variants, one of which was found in three affected members within one family.

Meta-analysis of genome-wide scans provides evidence for sex- and site-specific regulation of bone mass

Several genome-wide scans have been performed to detect loci that regulate BMD, but these have yielded inconsistent results, with limited replication of linkage peaks in different studies. In an effort to improve statistical power for detection of these loci, we performed a meta-analysis of genome-wide scans in which spine or hip BMD were studied. Evidence was gained to suggest that several chromosomal loci regulate BMD in a site-specific and sex-specific manner.

INTRODUCTION

BMD is a heritable trait and an important predictor of osteoporotic fracture risk. Several genome-wide scans have been performed in an attempt to detect loci that regulate BMD, but there has been limited replication of linkage peaks between studies. In an attempt to resolve these inconsistencies, we conducted a collaborative meta-analysis of genome-wide linkage scans in which femoral neck BMD (FN-BMD) or lumbar spine BMD (LS-BMD) had been studied.

MATERIALS AND METHODS

Data were accumulated from nine genome-wide scans involving 11,842 subjects. Data were analyzed separately for LS-BMD and FN-BMD and by sex. For each study, genomic bins of 30 cM were defined and ranked according to the maximum LOD score they contained. While various densitometers were used in different studies, the ranking approach that we used means that the results are not confounded by the fact that different measurement devices were used. Significance for high average rank and heterogeneity was obtained through Monte Carlo testing.

RESULTS

For LS-BMD, the quantitative trait locus (QTL) with greatest significance was on chromosome 1p13.3-q23.3 (p = 0.004), but this exhibited high heterogeneity and the effect was specific for women. Other significant LS-BMD QTLs were on chromosomes 12q24.31-qter, 3p25.3-p22.1, 11p12-q13.3, and 1q32-q42.3, including one on 18p11-q12.3 that had not been detected by individual studies. For FN-BMD, the strongest QTL was on chromosome 9q31.1-q33.3 (p = 0.002). Other significant QTLs were identified on chromosomes 17p12-q21.33, 14q13.1-q24.1, 9q21.32-q31.1, and 5q14.3-q23.2. There was no correlation in average ranks of bins between men and women and the loci that regulated BMD in men and women and at different sites were largely distinct.

CONCLUSIONS

This large-scale meta-analysis provided evidence for replication of several QTLs identified in previous studies and also identified a QTL on chromosome 18p11-q12.3, which had not been detected by individual studies. However, despite the large sample size, none of the individual loci identified reached genome-wide significance.

Large-scale genome-wide linkage analysis for loci linked to BMD at different skeletal sites in extreme selected sibships

Few genome-wide linkage studies of osteoporosis have been conducted in the Asian population. We performed a genome-wide scan involving 3093 adult siblings with at least one sib-pair extremely concordant or discordant for hip BMD. Our results indicated four genome-wide significant QTLs for BMD. In comparison with 12 previous reported linkage studies, we reveal novel linkage regions that have reaching global significance.

INTRODUCTION

The genetic basis for osteoporosis has been firmly established, but efforts to identify genes associated with this complex trait have been incomplete, especially in Asian populations. The purpose of this study was to identify quantitative trait loci (QTLs) for BMD in a Chinese population.

MATERIALS AND METHODS

We performed a genome-wide scan involving 3093 siblings 25-64 years of age from 941 families, with at least one sib-pair extreme concordant or discordant for total hip BMD from a large community-based cohort (n = 23,327) in Anhui, China. Linkage analysis was performed on BMD residuals adjusted for age, height, weight, occupation, cigarette smoking, physical activity, and alcohol consumption using the revised Haseman-Elston regression-based linkage model.

RESULTS

Our results revealed significant QTLs on chromosome 7p21.2 for femoral neck BMD (LOD = 3.68) and on chromosome 2q24.3 for total hip BMD (LOD = 3.65). Suggestive linkage regions were found to overlap among different skeletal sites on chromosomes 2q, 7p, and 16q. Sex-specific linkage analysis further revealed a significant QTL for lumbar spine BMD on chromosome 13q21.1 (LOD = 3.62) in women only. When performing multivariate linkage analysis by combining BMDs at four skeletal sites (i.e., whole body, total hip, femoral neck, and lumbar spine BMD), an additional significant QTL was found at chromosome 5q21.2 (LOD = 4.56). None of these significant QTLs found in our study overlapped with major QTLs reported by other studies.

CONCLUSIONS

This study reveals four novel QTLs in a Chinese population and suggests that BMD at different skeletal sites may also share common genetic determinants.

Quantitative trait loci modulate vertebral morphology and mechanical properties in a population of 18-month-old genetically heterogeneous mice

The aim of this study was to examine effects of polymorphic genes on vertebral bone morphology and mechanical properties. Genotypes from 525 18-month-old female mice were compared to geometric traits obtained from micro-computed tomography and mechanical properties from compression testing. Genetic markers were associated with traits on at least 13 different chromosomes, demonstrating the complexity of genetic control over vertebral form, function and aging.

Identification of quantitative trait loci affecting body composition in a mouse intercross

Gravimetric analysis and dual energy x-ray absorptiometry densitometry were used to determine lean, fat, and bone tissue traits in a F(2) mouse population from a C57BL/6J and CASA/Rk intercross (B6CASAF2). These traits were used in a linkage analysis to identify quantitative trait loci that affect body composition. Linkage mapping showed that body weight (BW) loci on proximal chromosome 2 occurred in the same region as body length, lean tissue mass, and bone mineral content and on chromosome 13 in the same region as lean tissue mass, bone mineral density, and bone mineral content. Fat-related loci occurring on mid-chromosome 2 near 60 cM, proximal chromosome 6, and mid-chromosome 10 were distinct from BW, lean tissue, and bone tissue loci. In B6CASAF2 females, heterozygotes and CASA/Rk homozygotes at the chromosome 6 locus marker had higher body fat percentages, and this locus was responsible for 11% of the variance for body fat percentage. Female heterozygotes and C57BL/6J homozygotes at the chromosome 15 locus marker had higher bone mineral densities, and this locus could explain 8% of that trait's variance. A survey of the literature did not reveal any previous reports of fat-specific loci in the chromosomal 10 region near 42 cM reported in this study. The results of this study indicate that BW and BMI have limited usefulness as phenotypes in linkage or association studies when used as obesity phenotypes.

Secreted frizzled-related protein 4 is a negative regulator of peak BMD in SAMP6 mice

We segregated a QTL for peak BMD on Chr 13 by generating congenic sublines of the senescence-accelerated mouse SAMP6. Sfrp 4 within this locus was responsible for lower BMD of SAMP6.

INTRODUCTION

Our genome-wide linkage study using SAMP6 and SAMP2 showed a significant quantitative trait locus (QTL) for peak BMD on chromosome (Chr) 13. To verify the gene that regulates peak BMD, we generated a congenic strain, P6.P2-Pbd2(b), which carried a 15-cM SAMP2 interval on an osteoporotic SAMP6 background, and showed that this Pbd2 locus increased peak BMD in SAMP6.

MATERIALS AND METHODS

To narrow down this interval, we generated a new congenic subline P6.P2-13. We studied the effect of this locus on morphological and histomorphological features in vivo and on osteoblasts in vitro. The levels of expression of all genes in the segregated interval were examined, and we clarified the effect of the candidate gene, secreted frizzled-related protein (Sfrp4), on osteoblasts in vitro.

RESULTS

The new congenic strain, P6.P2-13, retained the 2.4-Mb SAMP2 interval on the SAMP6 background, and 11 genes existed in this interval. In morphometrical analysis, P6.P2-13 increased the bone area fraction (BA/TA) by 6.6% at the diaphysial cortex (p < 0.001) and increased the trabecular bone volume (BV/TV) by 54.2% at the distal metaphysis (p < 0.05) in the femora compared with those of SAMP6. The bone formation rate of P6.P2-13 was markedly increased at the periosteal surface of femoral cortex and that was caused by a higher proliferation rate of osteoblasts in P6.P2-13 compared with those in SAMP6. Quantitative RT-PCR analysis of calvaria tissue showed approximately 40-fold higher levels of expression of Sfrp4 in SAMP6 than in P6.P2-13. Taken together with the result that recombinant Sfrp4 suppressed the proliferation of osteoblasts, we hypothesized that Sfrp4 inhibited the proliferation of osteoblasts through its antagonistic effect on Wnt signaling. TCF/beta-catenin-dependent reporter activity in osteoblasts derived from SAMP6 showed lower responsiveness for the Wnt ligand, Wnt3A, than that in osteoblasts from P6.P2-13.

CONCLUSIONS

In SAMP6 mice, Sfrp4 negatively regulates bone formation and decreases BMD through the inhibition of Wnt signaling.

Genomic regions identified for BMD in a large sample including epistatic interactions and gender-specific effects

A genome-wide linkage scan was conducted using a large white sample to identify QTLs for BMD. We found QTLs in the total sample and the gender-specific subgroups, as well as significant epistatic interactions underlying BMD variations.

INTRODUCTION

Low BMD is an important risk factor for osteoporosis and under strong genetic control.

MATERIALS AND METHODS

To identify quantitative trait loci (QTLs) for regulation of BMD, we performed a large-scale whole genome linkage scan (WGS) involving 4126 individuals from 451 families. In addition to the conventional linkage analyses in the total combined sample of males and females, we conducted epistatic interaction analyses and gender-specific linkage analyses.

RESULTS

Significant linkage was detected on 5q23 for wrist BMD (LOD = 3.39) and 15q13 for female spine BMD (LOD = 4.49). For spine BMD, we revealed significant epistatic interactions between 3p25 and 2q32 (p = 0.0022) and between 3p25 and 11q23 (p = 0.0007). We replicated several genomic regions that showed linkage with BMD in previous studies by others and ours, such as 3p21, 1p36, and Xq27.

CONCLUSIONS

This study highlights the importance of large sample size, incorporation of epistatic interaction, and consideration of gender-specific effects in identifying QTLs for BMD variation. The results of this study provide a foundation for the future fine mapping and gene identification in our population.