|
|
|
REVIEW
1 Rheumatic Diseases Unit, Molecular Medicine Centre, Western General Hospital, Edinburgh EH4 2XU, United Kingdom; 2 Department of Molecular Genetics, M.D. Anderson Cancer Center, University of Texas, Houston, Texas 77030, USA
| Abstract |
|---|
|
|
|---|
[Keywords: Human molecular genetics; linkage; osteoporosis; association studies]
| Genetic influences on osteoporosis |
|---|
|
|
|---|
It is currently believed that in the general population, genetic influences on BMD and the other phenotypes mentioned above are polygenic in nature and are mediated by the influence of several genetic variants, each of modest effect size, and their interaction with environmental factors. It is important to note, however, that severe osteoporosis and fragility fractures or unusually high bone mass (HBM) can also be an important feature of rare diseases that are primarily genetic in nature and are inherited in a classical Mendelian manner. Such diseases include osteopetrosis, sclerosing bone dysplasias, osteogenesis imperfecta, osteoporosis-pseudoglioma syndrome (OPPS), and osteoporotic syndromes associated with inactivating mutations of the estrogen receptor
and aromatase genes (Bilezikian et al. 1998
; Janssens and Van Hul 2002
; Van Wesenbeeck et al. 2003
; Balemans et al. 2005
). Although these diseases are rare, increasing evidence suggests that subtle polymorphic variations in genes that are mutated in these disorders also regulate BMD in the general population.
| Identifying osteoporosis genes |
|---|
|
|
|---|
|
Association studies have been widely used in the field of osteoporosis genetics. These studies typically involve identifying polymorphisms of a particular candidate gene and relating alleles to BMD or osteoporotic fractures in a population-based study or a case-control study. Recently, however, it has also become possible to perform association studies on a genome-wide basis by analyzing a large number of closely spaced single-nucleotide polymorphisms (SNPs) spread randomly across the genome (Cardon and Abecasis 2003
). The rationale for these studies is that these SNPs will be in linkage disequilibrium with causal variants in genes that predispose to the disease under study. Classical association studies typically involve choosing a candidate gene that is known to have effects on bone metabolism on the basis of knowledge about skeletal physiology or bone biology. Genes that are mutated in rare monogenic bone diseases have also proved to be a rich source of polymorphisms that regulate BMD in the normal population (Grant et al. 1996
; Ferrari et al. 2004
; Sobacchi et al. 2004
; Pettersson et al. 2005
). Association studies are relatively easy to perform and can be powered to detect small effects, but spurious results are often reported due to confounding factors, small sample size, and population stratification (Ioannidis 2003
). In view of this, meta-analysis has been increasingly used to estimate the true size of the effect of polymorphisms that have been implicated in the pathogenesis of complex diseases (Lohmueller et al. 2003
).
Linkage analysis in humans
Linkage studies are a tried and tested way of identifying genes responsible for monogenic diseases. Linkage studies have been spectacularly successful in identifying the genes that are responsible for rare, monogenic bone diseases characterized by abnormalities of bone mass such as osteopetrosis and sclerosing bone dysplasias (Janssens and Van Hul 2002
; Balemans et al. 2005
). Linkage analysis has also been applied to the identification of chromosomal regions that harbor genes that regulate quantitative traits such as bone mass, bone geometry, and ultra-sound properties of bone in the normal population (Koller et al. 2003
; Wilson et al. 2003
, 2004
; Ralston 2005
). Linkage studies in osteoporosis classically involved genotyping a series of polymorphic microsatellite markers, typically spread at 10-cM intervals throughout the genome, and relating inheritance of alleles to the inheritance of bone mass within the family members. More recently, however, much denser maps of SNP have been used for linkage studies of complex disease, and initial experience suggests that these may be more powerful than classical microsatellite-based methods (Sawcer et al. 2004
). A disadvantage of linkage studies in complex disease is that they have a relatively low statistical power to detect genes that have modest effects on the trait of interest, even when sample sizes of several thousand are used (Altmuller et al. 2001
). Reflecting this fact, linkage studies of outbred human populations have so far met with limited success in terms of identification of genes that predispose to osteoporosis, whereas studies of isolated populations have been more successful in identifying genetic variants that predispose to osteoporosis (Styrkarsdottir et al. 2003
) and other complex diseases (Gianfrancesco et al. 2003
; Gretarsdottir et al. 2003
).
Model organisms
Genetic linkage studies in experimental animals have long been used in the identification of genes responsible for complex traits (Rogers et al. 1997
) and, over recent years, many such studies have been performed in an attempt to identify loci involved in the regulation of BMD in mice (Benes et al. 2000
; Beamer et al. 2001
; Orwoll et al. 2001
). The approach used has been to cross inbred laboratory strains of mice with low and high bone density. By interbreeding offspring from the first generation (F1), a second generation (F2) of mice is established with varying levels of BMD, because of segregation of the alleles that regulate BMD in the F2 offspring. A genome-wide search is then performed in the F2 generation. This approach has several advantages: Environment can be carefully controlled, thus minimizing the influence of confounding factors, and large numbers of progeny can be generated, giving excellent statistical power. Genetic mapping studies in mice have identified a large number of loci that regulate BMD and have shown that the loci that regulate BMD are gender specific (Orwoll et al. 2001
). Only one gene has been identified by positional cloning in inbred strains of mice to date (Klein et al. 2004
), reflecting the fact that moving from QTL identification to gene identification is challenging. The reason for this is that the typical size of a QTL identified by genetic mapping within inbred strains of mice is 2040 cM, but within this region, the genetic variants are essentially in complete linkage disequilibrium with each other, making identification of the causal variant difficult. The region of interest can be narrowed by generating successive backcrosses against one of the parental strains, but this is time consuming and several generations of backcrossing are required to narrow the interval appreciably (Turner et al. 2003
). Microarray analysis to define differentially regulated genes offers a potential approach to identifying causal genes, and this was sucessfully used by Klein et al. (2004)
to identify a candidate gene (Alox15) within a mouse QTL that regulates BMD. Recent studies have indicated that polymorphic variation in a human homolog of this gene (Alox12) also regulates BMD in man (Ichikawa et al. 2006
).
| Candidate genes for osteoporosis |
|---|
|
|
|---|
|
|
The LRP5 and LRP6 receptors are transmembrane proteins that function as coreceptors for canonical Wnt signaling. Wnt signaling plays an important role in several key developmental processes, including cell fate decisions, limb patterning, and osteoblast and chondrocyte differentiation, and development of the central nervous system and other organs (Johnson et al. 2004
). Transduction of the Wnt signal involves both LRP5 and LRP6, which form a receptor complex with Frizzled (Fz) to activate the transcriptional activity of
-catenin, the downstream effector of canonical Wnt signaling (Cong et al. 2004
). Regulation of signaling through the WntLRP
-catenin pathway occurs at several levels, and key components of the pathway are illustrated in Figure 3. Extra-cellular inhibitors of Wnt signaling include the secreted Wnt antagonist; Dickkopf (Dkk) proteins, which bind to LRP5 and LRP6 and inhibit Wnt signaling (Bafico et al. 2001
; Mao and Niehrs 2003
); and a family of secreted frizzled-related proteins (Sfrps), which bind Wnt polypeptides and function as decoy receptors (Kawano and Kypta 2003
). Dickkopf proteins can also down-regulate Wnt signaling by interacting with the transmembrane protein Kremen. As a result of this interaction, a complex is formed containing Dkk, Kremen, and the LRP5/6 proteins that is internalized, thereby down-regulating Wnt signaling (Mao et al. 2002
).
|
-propeller motif of the molecule (Little et al. 2002
-propeller motif of LRP5 (Van Wesenbeeck et al. 2003
The mutations that cause OPPS produce a truncated or nonfunctional LRP5 protein (Gong et al. 2001
), and in accordance with this, complete inactivation of LRP5 by gene targeting in mice causes a low bone mass phenotype, which phenocopies the human disease OPPS (Kato et al. 2002
). Analysis of bone histomorphometry in LRP5 knockout mice has shown that the low bone mass is a consequence of decreased osteoblast proliferation and reduced bone matrix deposition rather than an effect of bone resorption (Kato et al. 2002
). The G171V mutation that was associated with HBM in the family studied by Johnson and colleagues (Johnson et al. 1997
; Little et al. 2002
) was found to cause increased bone mass when expressed in transgenic mice (Babij et al. 2003
). In these studies, mice were generated that expressed either wild-type human LRP5 or the G171V mutant, and the mice that expressed the mutant had significantly higher BMD than wild type, despite the fact that levels of transgene expression were similar in both mouse strains. Analysis of the skeletal phenotype in the G171V transgenic mice showed that mineral apposition rate was increased and the rate of osteoblast apoptosis was reduced, whereas eroded surface (reflecting bone resorption) was unaffected. Functional analysis of the HBM-associated mutations of LRP5 has shown that the mutations probably cause activation of
-catenin signaling by inhibiting interactions between LRP5 and Dkk1. An initial study by Boyden et al. (2002)
showed that the G171V mutation did not result in constitutive activation of LRP5 signaling in vitro, but that it impaired Dkk1-mediated inhibition of Wnt-stimulated LRP5 signaling. Another study reached the same conclusion in showing that several HBM-associated mutants (G171V, G171R, A214T, A214V, A242T, T253I, and D111Y) were resistant to Dkk1 inhibition compared with wild-type LRP5, and had a lower affinity for Dkk1 binding (Ai et al. 2005
).
There is evidence that more subtle variations in LRP5 could also underlie variations of BMD in the normal population. Several association studies and family-based studies have been performed in which various common polymorphisms of LRP5 have been related to BMD and/ or osteoporotic fracture. Most of these studies have shown evidence of associations between LRP5 alleles and BMD and, interestingly, these associations have been particularly strong in men (Ferrari et al. 2004
; Koh et al. 2004
; Urano et al. 2004
; van Meurs et al. 2006
). Although many variants have been studied, the most likely functional candidate is an alanine-to-valine amino acid substitution at position 1330 (A1330V). The mechanism by which this variant might affect LRP5 signaling has not been investigated, but evidence of an interaction between the LRP5 A1330V variant and a coding polymorphism of LRP6 (1062V) has also been gained in the Rotterdam study (van Meurs et al. 2006
), where polymorphisms of both genes were found to have an additive effect on fracture susceptibility. Interactions between LRP5 and LRP6 as regulators of BMD have also been observed in preclinical studies where skeletal phenotyping of mice with targeted inactivation of both receptors revealed allele dose-dependent deficits in BMD and limb formation, suggesting functional redundancy between these two genes in bone and limb development (Holmen et al. 2004
).
In conclusion, the data so far indicate that genetic variation in LRP5 and possibly LRP6 plays an important role in regulation of bone mass and osteoporotic fractures in humans. Not only do rare mutations in the LRP5 gene play a major role in regulating BMD, but more subtle polymorphisms also seem to regulate BMD in the normal population. Clearly, other components of this pathwaysuch as the Wnts, Dkks, Frizzled, Kremen, and Sfrpsnow represent prime candidate genes for further study as potential genetic regulators of bone mass.
Transforming growth factor TGF-
1
Members of the TGF-
superfamily of related secreted polypeptides control critical cellular functions during embryonic development and also postnatally. Most attention has focused on TGF-
1, which is encoded by the TGFB1 gene as a regulator of susceptibility to osteoporosis partly because it is particularly abundant in bone and has been shown to have effects on both osteoblast and osteoclast function in vitro (Massagué and Chen 2000
). The receptor system for TGF
molecules consists of two different high-affinity receptors, each of which contains an extracellular ligand-binding domain, a trans-membrane domain, and an intracellular kinase domain. Upon TGF-
binding, the type II receptor (TBR II) homodimer recruits a type I receptor (TBR I) homodimer to form a heterotetrameric receptor complex. Ligand binding to TBR II leads to phosphorylation of a glycine- and serine-rich domain in TBR I, which results in activation of signaling by this receptor. TBR I then phosphorylates Smad2 and Smad3, two so-called regulatory Smads or R-Smads, each of of which forms a complex with a unique common Smad, Smad4. These complexes trans-locate to the nucleus where they interact with other transcription factors and also, in the case of Smad3, with specific DNA sequences to activate a variety of downstream target genes. Inactivation of components of the TGF-
pathway frequently results in embryonic lethality, but TGF-
1-null mice develop generalized inflammation and decreased bone mass with a pronounced reduction of the number of osteoblasts (Geiser et al. 1998
). In addition, Smad3-null mice show osteopenia, which is thought to be due to excessive late-phase osteoblast differentiation into osteocytes and concomitant apoptosis (Borton et al. 2001
).
Polymorphisms of the TGFB1 gene have been extensively studied in relation to osteoporosis. One of the earliest studies was that of Langdahl et al. (1997)
, who identified a polymorphism within intron 4 of TGFB1 that was associated with severe osteoporosis. Subsequent work by the same group evaluated the relationship between several polymorphisms in TGFB1 and osteoporosis in a case-control study and identified an association between a different polymorphism in intron 5 and BMD (Langdahl et al. 2003
). Other research has focused on polymorphisms in the promoter and first exon of TGFB1 in relation to BMD. One of these polymorphisms (a leucine-to-proline substitution at codon 10) has been associated with BMD in some populations (Yamada et al. 1999
, 2001
). Although genetic association studies have yielded somewhat conflicting results, definitive evidence that genetic variation in TGFB1 regulates bone mass in humans comes from the observation that Camurati-Engelmann (CED) disease is caused by mutations in TGFB1 (Janssens et al. 2000a
; Kinoshita et al. 2000
). CED is a rare autosomal dominant genetic disorder characterized by hyperostosis and sclerosis, especially affecting the diaphysis of long bones. Clinically, CED is characterized by radiological osteosclerosis, but affected patients also have increased markers of bone resorption, suggesting that bone resorption is elevated as well as bone formation (Hernandez et al. 1997
). The causal gene for CED was mapped to chromosome 19q13 by linkage analysis (Ghadami et al. 2000
; Janssens et al. 2000b
), and then mutations in TGFB1 were identified as a cause using a positional candidate approach (Janssens et al. 2000a
; Kinoshita et al. 2000
). Most causal mutations lie within a region of the gene that encodes the latency-associated peptide (LAP), although two other mutations have been identified in the signal peptide region, including a duplication of three leucine residues at codon 12 and a tyrosine-to-histidine amino acid change at codon 81. Following transcription and translation, the mature TGF
1 peptide is cleaved from the LAP but is secreted from the cell in an inactive form as a dimer, which is bound noncovalently through disulfide bonds to LAP. These molecules can also form a larger covalent complex with a latent TGF-
-binding protein (LTBP). Activation of TGF
1, which occurs outside the cell, can occur as the result of low pH or proteolytic activity and results in the release of the active TGF
1 dimer from the LAP and other inhibitory molecules so that it can now directly interact with the TGF
1 receptors (Janssens et al. 2003
). Functional studies have shown that CED-causing mutations all act to increase SMAD-dependent TGF
1 signaling, but by different mechanisms. Mutations in the LAP domain appear to inhibit formation of the disulfide bonds that link TGF
1 to the LAP and cause increased release of biologically active TGF
1. However, the other mutations in the signal peptide region prevent TGF
1 secretion and result in intracellular accumulation of TGF-
1 with stimulation of SMAD signaling by an intracrine mechanism that is poorly understood.
Bone morphogenic proteins (BMPs)
BMPs are part of the TGF
superfamily of molecules. Signaling through BMP is regulated extracellularly by antagonists such as noggin and chordin (De Robertis and Kuroda 2004
), and intracellularly by the Smad family of signaling proteins that act downstream from specific BMP receptors (Cao and Chen 2005
). BMPs are strong candidate genes for regulation of bone mass, and naturally occurring mutations in various components of the BMP pathway result in defects of limb formation (Seemann et al. 2005
). Recent studies indicate that polymorphic variation in members of the BMP family may be involved in regulation of bone mass and susceptibility to osteoporosis. A linkage study in Icelanders identified a locus for regulation of BMD on chromosome 20p12 (Styrkarsdottir et al. 2003
), and subsequent positional cloning studies showed that the linkage signal was attributable in part to genetic variation in the BMP2 gene. A coding polymorphism was identified at codon 37 of the BMP2 protein, resulting in a substitution of alanine for a serine residue (Ser37Ala) that was overrepresented in Icelandic patients and in Danish patients with osteoporosis (Styrkarsdottir et al. 2003
). The mechanisms by which this substitution affects BMP signaling to affect bone mass have not yet been investigated. Other workers have reported that common polymorphism, which causes an alanine-to-valine substitution at position 152 in the BMP4 protein (Ala152Val), is associated with low BMD (Ramesh et al. 2005
). These association studies raise the possibility that polymorphic variation in members of the BMP family may be involved in osteoporosis, but further work will be required to confirm these observations and to study the mechanism by which the coding variants so far identified affect bone cell function and modulate bone mass.
Sclerostin
Sclerostin is a cysteine knot-containing protein that shares homology with the DAN family of BMP antagonists (Winkler et al. 2003
; van Bezooijen et al. 2004
). Sclerostin was first identified as a regulator of bone mass by genetic mapping studies in the sclerosing bone dysplasias Sclerosteosis and van Buchem disease. These are rare conditions with a similar phenotype that are inherited in an autosomal recessive manner. Both are characterized by progressive osteosclerosis mainly affecting the skull, mandible, and tubular bones of the extremities (Beighton et al. 1984
). Sclerosteosis differs from van Buchem disease because it is more severe and is associated with hand malformations such as syndactyly, absent or dysplastic nails, and radial deviation of the terminal phalanges (Balemans et al. 2005
). Both diseases were mapped to the same region of chromosome 17q12-21 in the late 1990s (Van Hul et al. 1998
; Balemans et al. 1999
), and positional cloning studies identified recessive loss-of-function mutations in Sclerostin (SOST) gene as the cause of Sclerosteosis in 2001 (Balemans et al. 2001
; Brunkow et al. 2001
). Interestingly, mutation screening of SOST in patients with van Buchem disease showed no mutations, but a homozygous 52-Kb deletion was identified 35 Kb downstream from SOST as the cause of this condition (Balemans et al. 2002b
). Subsequent work has shown that the deleted region contains highly conserved regulatory elements that play a key role in regulating SOST expression in bone (Loots et al. 2005
). Although heterozygous carriers of SOST mutations are asymptomatic, they have significantly increased BMD, indicating that haploinsufficiency of SOST increases bone mass without apparently causing adverse effects (Gardner et al. 2005
).
Polymorphisms of SOST have also been associated with BMD in some population-based studies. In the Rotterdam study (Uitterlinden et al. 2004
), polymorphisms in the SOST promoter and the region deleted in van Buchem disease were both found to be associated with BMD, and evidence was gained to suggest that the effect became more pronounced with increasing age. In another case-control study, SOST polymorphisms were not found to be associated with BMD (Balemans et al. 2002a
), although this was a relatively small study and participants were of a younger age group than in the Rotterdam study.
In view of the structural similarity between SOST and other BMP antagonists, it was initially assumed that SOST regulates bone mass by affecting BMP signaling (Winkler et al. 2003
). In keeping with this view was the observation that Sclerostin binds to both BMP5 and BMP6, competes with BMP6 for binding to the type I BMP receptor, and antagonizes the effects of BMP6 on Smad signaling (Winkler et al. 2003
). Transgenic mice that overexpress SOST are osteopenic and have reduced bone formation, consistent with a model whereby SOST negatively regulates osteoblast function (Winkler et al. 2003
). However, studies by van Bezooijen et al. (2004)
showed that SOST was not a classical BMP antagonist, suggesting that it might interact with other signaling pathways to regulate bone formation. The mechanisms by which SOST inhibits bone formation has been clarified by recent studies, which have shown that SOST directly interacts with LRP5 and LRP6 and antagonizes Wnt signaling both in vitro and in vivo (Li et al. 2005
; Semenov et al. 2005
). However, other workers have proposed that the inhibitory effect of SOST on Wnt signaling is indirect and mediated through BMP production (Winkler et al. 2005
). While further studies will be required to clarify the mechanisms by which SOST affects bone formation, it seems possible that the increased bone mass that results from SOST deficiency is mediated, at least in part, through the LRP5
-catenin pathway.
CBFA1
The CBFA1 gene (also known as Runx2) plays an essential role in regulating osteoblast differentiation, since mice that are deficient in this transcription factor have complete absence of bone (Komori et al. 1997
; Otto et al. 1997
). Moreover, mice with haploinsufficiency of CBFA1 phenocopy the human syndrome of cleidocranial dysplasia (CCD), a skeletal disorder characterized by short stature, hypoplasia or aplasia of the clavicles, patent fontanelles, supernumerary teeth, and other defects in skeletal patterning and growth (Otto et al. 1997
). The human syndrome of CCD is caused by various missense, nonsense, and frameshift mutations of CBFA1 (Lee et al. 1997
; Mundlos et al. 1997
; Quack et al. 1999
). Some of these mutations have been shown to interfere with the DNA-binding activity of CBFA1, whereas others have been found to alter nuclear localization of the protein or to produce a mutant or truncated protein that is biologically inactive. In addition to these rare mutations, various polymorphisms have been identified in CBFA1 and some of these have been associated with bone mass in population-based studies (Vaughan et al. 2002
, 2004
; Doecke et al. 2006
). The best functional candidates lie within the Runx2 promoter (Doecke et al. 2006
) or within polyalanine and polyglutamine repeats in exon 1. The polyalanine and polyglutamine repeats are of interest since they lie within one of the transactivation domains of Runx2. Various polymorphic variations have been identified in this region, including an 18-bp deletion that results in a polyalanine repeat of 11 residues (11 Ala) compared with the more common repeat of 17 residues (17 Ala). Various rare length variants within the polyglutamine repeat have also been identified, resulting in stretches of between 15 and 30 repeats. The strongest association with BMD has been observed with an anonymous polymorphism in the Ala repeat region (Vaughan et al. 2002
, 2004
), although it is thought that this might be due to linkage disequilibrium with polymorphisms in the promoter that have been shown to affect Runx2 transcription in reporter assays (Doecke et al. 2006
). It is currently unclear whether the length variants in the polyalanine and polyglutamine tracts have functional importance, but this is an area of ongoing investigation.
Cathepsin K
Mutations in Cathepsin K are responsible for the rare syndrome of Pycnodysostosis, which is a rare recessive bone dysplasia characterized by osteosclerosis and short stature (Gelb et al. 1996
). Various disease-causing mutations have been identified that result in production of either a truncated or a nonfunctional protein. Studies of more common polymorphic variants of Cathepsin K in relation to BMD regulation in the normal population have so far yielded negative results (Giraudeau et al. 2004
).
TCIRG1
The TCIRG1 gene encodes the APT6i subunit of the osteoclast-specific proton pump, and various inactivating mutations in TCIRGI have been described that are responsible for a subset of patients with recessive osteopetrosis (Frattini et al. 2000
). Recent work indicates that polymorphisms of TCIRG1 might also contribute to regulation of BMD in the normal population; a study by Frattini (Sobacchi et al. 2004
) showed evidence of an association between a polymorphism affecting an AP1-binding site in the TCIRG1 promoter and BMD in perimenopausal women. Functional studies need to be performed to identify the mechanisms that underlie this association, however, and to replicate the finding in other populations.
CLCN7
The CLCN7 gene encodes a chloride channel that is highly expressed in osteoclasts and is essential for acidification of the resorption lacuna. Homozygous inactivation mutations in CLCN7 cause a severe form of recessive osteopetrosis (Kornak et al. 2001
), whereas heterozygous missense mutations of CLCN7 cause autosomal dominant osteopetrosis (Balemans et al. 2005
). Prompted by this observation, Pettersson et al. (2005)
looked for evidence of an association between polymorphisms of CLCN7 and BMD in normal individuals and found that a common polymorphism in exon 15 of CLCN7 that results in a methionine-to-valine amino acid change was associated with BMD in normal women. Further studies will be required to determine if this polymorphism is functionally important and to replicate the observation in other populations.
Vitamin D receptor (VDR)
Vitamin D, through its principal bioactive form 1,25-dihydroxyvitamin D3[1,25-(OH)2D3] plays a crucial role in bone metabolism. The action of 1,25-(OH)2D3 is mediated through a specific hormone-receptor (VDR) that regulates gene expression by forming a heterodimer with retinoic X receptor (RXR) that binds to vitamin D response elements in target genes. Mutations in VDR cause the syndrome of vitamin D-resistant rickets which is a recessive condition characterized by alopecia, hypocalcaemia, hypophosphatemia, and severe rickets, and is resistant to treatment with vitamin D and its active metabolites (Kristjansson et al. 1993
). Targeted inactivation of the VDR gene in mice provides a phenocopy of the human syndrome (Bouillon et al. 2003
), but when these mice are fed a diet rich in calcium and phosphorus to normalize serum calcium and PTH concentrations, they develop normally without bone anomalies. It therefore appears that the bone defects that result from VDR deficiency are due to malabsorption of calcium and phosphate by the intestine rather than absence of 1,25-(OH)2D3 signaling in bone cells. The VDR was the first candidate gene to be studied in relation to osteoporosis, and most attention has focused on polymorphisms situated near the 3' flank of VDR recognized by the restriction enzymes BsmI, ApaI, and TaqI. A recent retrospective meta-analysis of association studies that genotyped the BsmI polymorphism concluded that there was evidence of an association between spine BMD and the BsmI polymorphism, equivalent to ~0.15 Z-score units between the BB genotype and the other genotypes (Thakkinstian et al. 2004
). No association with femoral neck BMD was observed. Another polymorphism affecting exon 2 of VDR has been described that creates an alternative translational start site, resulting in the production of two isoforms of the VDR protein, which differ in length by three amino acids (Gross et al. 1997
). This has been associated with BMD in some, but not other, populations, and functional studies of this polymorphism have yielded inconclusive results (Gross et al. 1998
). A further polymorphism has been identified in the promoter of VDR at a binding site for the transcription factor Cdx-2, which has been associated with BMD in Japanese subjects and appears to be functional (Arai et al. 2001
). This polymorphism has been associated with fracture in other populations, but not with BMD (Fang et al. 2003
). The most comprehensive study of VDR alleles in relation to osteoporosis-related phenotypes was that of Fang et al. (2005)
, who conducted a large-scale study of haplotype tagging SNP of VDR in 6418 participants of the Rotterdam study. As the result of this analysis, the authors identified haplotype alleles in the promoter region and 3' untranslated region (UTR) that were associated with an increased risk of fracture. For a subgroup of individuals who carried risk alleles at both sites, the fracture risk was significantly increased by 48% when compared with control subjects. Functional studies also showed that the promoter haplotype that increased fracture risk was associated with reduced VDR expression in reporter assays, whereas the 3' UTR risk haplotype was associated with increased degradation of VDR mRNA. The data would be consistent with a model whereby the combination of risk haplotypes results in a lower VDR mRNA level due to decreased transcription and increased degradation. Interestingly, the risk alleles for fracture identified in this study were not associated with differences in BMD. In view of this, the mechanism by which these polymorphisms predispose to fracture is unclear, although a possibility, alluded to by the authors, was through an effect on bone geometry. Although this was a large and well-conducted study, the risk estimates were modest, and if correction had been applied for all the combinations of haplotypes tested (and their interactions), the results would not have been significant.
Collagen type I
I
The gene encoding the
I chain of type I collagen (COLIA1) is an important functional candidate for the pathogenesis of osteoporosis, since type I collagen is the major protein of bone and since mutations in this gene cause the syndrome of osteogenesis imperfectaa rare disease characterized by increased bone fragility and reduced BMD (Boyde et al. 1999
). Previous research identified associations between BMD and polymorphisms within the proximal promoter of COLIA1 (Garcia-Giralt et al. 2002
) and within the first intron of COLIA1 (Grant et al. 1996
). Most research has focused on a polymorphism within intron 1, which is situated at a binding site for the transcription factor Sp1. The COLIA1 Sp1-binding site polymorphism has been associated with various osteoporosis-related phenotypes, including bone density (Grant et al. 1996
; Uitterlinden et al. 1998
), fragility fractures (Grant et al. 1996
; Uitterlinden et al. 1998
), post-menopausal bone loss (Harris et al. 2000
; MacDonald et al. 2001
), bone geometry (Qureshi et al. 2001
), bone quality (Mann et al. 2001
), and bone mineralization (Stewart et al. 2005
). Functional analysis has shown that the osteoporosis-associated T allele ("s") of the Sp1 polymorphism is associated with increased DNAprotein binding, increased transcription, and abnormally increased production of the COLIA1 mRNA and protein (Mann et al. 2001
). It is thought that the resulting imbalance between the COLIA1 and COLIA2 chains contributes to impairment of bone strength and reduced bone mass in carriers of the T allele by subtly affecting bone mineralization (Stewart et al. 2005
). Retrospective meta-analyses of published studies have concluded that carriage of the T allele is associated with reduced BMD at the lumbar spine and femoral neck and with vertebral fractures (Efstathiadou et al. 2001
; Mann and Ralston 2003
). Moreover, homozygotes for the T allele of the Sp1 polymorphism were recently found to be associated with BMD and incident vertebral fractures in a large prospective meta-analysis of >20,000 participants of the GENOMOS study (Ralston et al. 2006
). In the GENOMOS study, however, the BMD association was only observed for homozygotes for the T allele, in contrast to previous studies where heterozygotes also showed a reduction in BMD (Mann and Ralston 2003
). The association with vertebral fracture in GENOMOS was similar to previous reports and equivalent to a 40% increase in fracture risk for each copy of the T allele carried. Interestingly, the association between COLIA1 alleles and vertebral fracture reported in GENOMOS and other studies was not fully explained on the basis of reduced bone density, implying that the Sp1 allele also acts as a marker for bone quality. The promoter polymorphisms of COLIA1 are in strong linkage disequilibrium with the Sp1 polymorphism, and two studies have now suggested that an extended haplotype defined by the Sp1 polymorphism and other promoter polymorphism may exert stronger effects on BMD than the individual polymorphisms (Garcia-Giralt et al. 2002
; Stewart et al. 2006
). Evidence has been presented to suggest the promoter polymorphism at position 1663 interacts with the transcription factor NMP4, which plays a role in osteoblast differentiation by interacting with Smads (Garcia-Giralt et al. 2005
). While a functional interaction has been reported between the 1663 and 1997 polymorphisms in regulating COLIA1 transcription (Garcia-Giralt et al. 2005
), it is currently unclear whether the 1997 polymorphism is also at the site of transcription factor binding.
Estrogen receptor 
The estrogen receptor
, encoded by the ESR1 gene, is another important functional candidate for the regulation of bone mass. A large number of investigators have looked for evidence of an association between ESR1 alleles and BMD, mostly focusing on two polymorphisms within intron 1 recognized by the XbaI and PvuII restriction enzymes, and on a TA repeat in the promoter. Studies of ESR1 alleles in relation to BMD have yielded inconsistent results, possibly because most studies have been of small sample size and involved subjects of different ages, menopausal status, and ethnic backgrounds (Ioannidis et al. 2002
). A meta-analysis of published studies showed evidence of an association between the XbaI polymorphisms and both BMD and fractures, with higher BMD values in "XX" homozygotes and a reduced risk of fractures in association with this genotype (Ioannidis et al. 2002
). More recently, a prospective meta-analysis of data from >18,000 subjects in the GENOMOS study confirmed that "XX" homozygotes had a reduced risk of fracture (Ioannidis et al. 2004
). In this study, no association with BMD was observed, indicating that ESR1 influences fracture risk independent of an effect on BMD. The mechanism responsible for this observation remains unclear, but one possibility may be an effect on bone quality, since ESR1 alleles have been associated with ultrasound properties of bone and bone loss (Albagha et al. 2005
). There has been controversy as to whether the associations between BMD and polymorphisms in the 5' region of the ESR1 gene are driven by the intron 1 polymorphisms or the TA repeat polymorphism in the ESR1 promoter. Previous studies have shown that the PvuII polymorphism lies within consensus recognition sites for the AP4 and Myb transcription factors; studies using promoter-reporter assays have shown that the PvuII polymorphism influences Myb-driven transcription in vitro (Herrington et al. 2002
); and other studies have shown that both the XbaI and PvuII polymorphisms also influence reporter gene transcription in vitro (Maruyama et al. 2000
). In this regard, it is of interest that the PvuII and XbaI polymorphisms are located within a region that is 70%80% conserved in the human, mouse, and rat genomes, whereas the TA repeat polymorphism is not conserved to any significant extent across species. These data are consistent with a direct functional effect of the PvuII and XbaI polymorphisms on ESR1 gene transcription, but it remains possible that the intron 1 polymorphisms are simply in linkage disequilibrium (LD) with causal polymorphisms elsewhere in the ESR1 gene.
| Concluding remarks |
|---|
|
|
|---|
| Footnotes |
|---|
E-MAIL stuart.ralston{at}ed.ac.uk; FAX 44-131-651-1085. ![]()
Article is online at http://www.genesdev.org/cgi/doi/10.1101/gad.1449506.
| References |
|---|
|
|
|---|
Albagha, O.M., Pettersson, U., Stewart, A., McGuigan, F.E., MacDonald, H.M., Reid, D.M., Ralston, S.H. 2005. Association of oestrogen receptor
gene polymorphisms with postmenopausal bone loss, bone mass, and quantitative ultrasound properties of bone. J. Med. Genet. 42: 240246.
Altmuller, J., Palmer, L.J., Fischer, G., Scherb, H., Wjst, M. 2001. Genomewide scans of complex human diseases: True linkage is hard to find. Am. J. Hum. Genet. 69: 936950.[CrossRef][Medline]
Andrew, T., Antioniades, L., Scurrah, K.J., MacGregor, A.J., Spector, T.D. 2005. Risk of wrist fracture in women is heritable and is influenced by genes that are largely independent of those influencing BMD. J. Bone Miner. Res. 20: 6774.[CrossRef][Medline]
Arai, H., Miyamoto, K.I., Yoshida, M., Yamamoto, H., Taketani, Y., Morita, K., Kubota, M., Yoshida, S., Ikeda, M., Watabe, F. et al. 2001. The polymorphism in the caudal-related homeodomain protein Cdx-2 binding element in the human vitamin D receptor gene. J. Bone Miner. Res. 16: 12561264.[CrossRef][Medline]
Arden, N.K. and Spector, T.D. 1997. Genetic influences on muscle strength, lean body mass, and bone mineral density: A twin study. J. Bone Miner. Res. 12: 20762081.[CrossRef][Medline]
Arden, N.K., Baker, J., Hogg, C., Baan, K., Spector, T.D. 1996. The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: A study of postmenopausal twins. J. Bone Miner. Res. 11: 530534.[Medline]
Babij, P., Zhao, W., Small, C., Kharode, Y., Yaworsky, P.J., Bouxsein, M.L., Reddy, P.S., Bodine, P.V., Robinson, J.A., Bhat, B. et al. 2003. High bone mass in mice expressing a mutant LRP5 gene. J. Bone Miner. Res. 18: 960974.[CrossRef][Medline]
Bafico, A., Liu, G., Yaniv, A., Gazit, A., Aaronson, S.A. 2001. Novel mechanism of Wnt signalling inhibition mediated by Dickkopf-1 interaction with LRP6/Arrow. Nat. Cell Biol. 3: 683686.[CrossRef][Medline]
Balemans, W., Van Den Ende, J., Freire Paes-Alves, A., Dikkers, F.G., Willems, P.J., Vanhoenacker, F., Almeida-Melo, N., Alves, C.F., Stratakis, C.A., Hill, S.C. et al. 1999. Localization of the gene for sclerosteosis to the van Buchem disease-gene region on chromosome 17q12-q21. Am. J. Hum. Genet. 64: 16611669.[CrossRef][Medline]
Balemans, W., Ebeling, M., Patel, N., Van Hul, E., Olson, P., Dioszegi, M., Lacza, C., Wuyts, W., Van Den Ende, J. Willems, P. et al. 2001. Increased bone density in sclerosteosis is due to the deficiency of a novel secreted protein (SOST). Hum. Mol. Genet. 10: 537543.
Balemans, W., Foernzler, D., Parsons, C., Ebeling, M., Thompson, A., Reid, D.M., Lindpaintner, K., Ralston, S.H., Van Hul, W. 2002a. Lack of association between the SOSTgene and bone mineral density in perimenopausal women: Analysis of five polymorphisms. Bone 31: 515519.[Medline]
Balemans, W., Patel, N., Ebeling, M., Van Hul, E., Wuyts, W., Lacza, C., Dioszegi, M., Dikkers, F.G., Hildering, P., Willems, P.J. et al. 2002b. Identification of a 52 kb deletion downstream of the SOST gene in patients with van Buchem disease. J. Med. Genet. 39: 9197.
Balemans, W., Van Wesenbeeck, L., Van Hul, W. 2005. A clinical and molecular overview of the human osteopetroses. Calcif. Tissue Int. 77: 263274.[CrossRef][Medline]
Beamer, W.G., Shultz, K.L., Donahue, L.R., Churchill, G.A., Sen, S., Wergedal, J.R., Baylink, D.J., Rosen, C.J. 2001. Quantitative trait loci for femoral and lumbar vertebral bone mineral density in C57BL/6J and C3H/HeJ inbred strains of mice. J. Bone Miner. Res. 16: 11951206.[CrossRef][Medline]
Beighton, P., Barnard, A., Hamersma, H., van der Wouden, A. 1984. The syndromic status of sclerosteosis and van Buchem disease. Clin. Genet. 25: 175181.[Medline]
Benes, H., Weinstein, R.S., Zheng, W., Thaden, J.J., Jilka, R.L., Manolagas, S.C., Shmookler Reis, R.J. 2000. Chromosomal mapping of osteopenia-associated quantitative trait loci using closely related mouse strains. J. Bone Miner. Res. 15: 626633.[CrossRef][Medline]
Bilezikian, J.P., Morishima, A., Bell, J., Grumbach, M.M. 1998. Increased bone mass as a result of estrogen therapy in a man with aromatase deficiency. N. Engl. J. Med. 339: 599603.
Borton, A.J., Frederick, J.P., Datto, M.B., Wang, X.F., Weinstein, R.S. 2001. The loss of Smad3 results in a lower rate of bone formation and osteopenia through dysregulation of osteoblast differentiation and apoptosis. J. Bone Miner. Res. 16: 17541764.[CrossRef][Medline]
Bouillon, R., Van Cromphaut, S., Carmeliet, G. 2003. Intestinal calcium absorption: Molecular vitamin D mediated mechanisms. J. Cell. Biochem. 88: 332339.[CrossRef][Medline]
Boyde, A., Travers, R., Glorieux, F.H., Jones, S.J. 1999. The mineralization density of iliac crest bone from children with osteogenesis imperfecta. Calcif. Tissue Int. 64: 185190.[CrossRef][Medline]
Boyden, L.M., Mao, J., Belsky, J., Mitzner, L., Farhi, A., Mitnick, M.A., Wu, D., Insogna, K., Lifton, R.P. 2002. High bone density due to a mutation in LDL-receptor-related protein 5. N. Engl. J. Med. 346: 15131521.
Brunkow, M., Gardner, J., Van Ness, J., Paeper, B., Kovacevich, B., Proll, S., Skonier, J., Zhao, L., Sabo, P., Fu, Y. et al. 2001. Bone dysplasia sclerosteosis results from loss of the sost gene product, a novel cystine knot-containing protein. Am. J. Hum. Genet. 68: 577589.[CrossRef][Medline]
Cao, X. and Chen, D. 2005. The BMP signaling and in vivo bone formation. Gene 357: 18.[CrossRef][Medline]
Cardon, L.R. and Abecasis, G.R. 2003. Using haplotype blocks to map human complex trait loci. Trends Genet. 19: 135140.[CrossRef][Medline]
Christian, J.C., Yu, P.L., Slemenda, C.W., Johnston, C.C. 1989. Heritability of bone mass: A longitudinal study in aging male twins. Am. J. Hum. Genet. 44: 429433.[Medline]
Cong, F., Schweizer, L., Varmus, H. 2004. Wnt signals across the plasma membrane to activate the
-catenin pathway by forming oligomers containing its receptors, Frizzled and LRP. Development 131: 51035115.
Cummings, S.R., Nevitt, M.C., Browner, W.S., Stone, K., Fox, K.M., Ensrud, K.E., Cauley, J., Black, D., Vogt, T.M. 1995. Risk factors for hip fracture in white women. Study of Osteoporotic Fractures Research Group. N. Engl. J. Med. 332: 767773.
Deng, H.W., Chen, W.M., Recker, S., Stegman, M.R., Li, J.L., Davies, K.M., Zhou, Y., Deng, H., Heaney, R., Recker, R.R. 2000. Genetic determination of Colles fracture and differential bone mass in women with and without Colles fracture. J. Bone Miner. Res. 15: 12431252.[CrossRef][Medline]
De Robertis, E.M. and Kuroda, H. 2004. Dorsalventral patterning and neural induction in Xenopus embryos. Annu. Rev. Cell. Dev. Biol. 20: 285308.[CrossRef]