PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Bone. Author manuscript; available in PMC Aug 1, 2011.
Published in final edited form as:
PMCID: PMC2902654
NIHMSID: NIHMS203738
A Common Variant in Fibroblast Growth Factor Binding Protein 1 (FGFBP1) is Associated with Bone Mineral Density and Influences Gene Expression In Vitro
Nicole Hoppman,1 John C. McLenithan,1 Daniel J. McBride,1 Haiqing Shen,1 Jan Bruder,2 Richard L. Bauer,2 John R. Shaffer,3 Jie Liu,1 Elizabeth A. Streeten,1 Alan R. Shuldiner,1,4 Candace M. Kammerer,3 and Braxton D. Mitchell1
1Division of Endocrinology, Diabetes and Nutrition, University of Maryland School of Medicine, Baltimore, MD 21201, USA
2 Department of Medicine, University of Texas Health Science Center, San Antonio, TX 78207, USA
3Department of Human Genetics, University of Pittsburgh Graduate School of Public Health, Pittsburgh, PA 15261, USA.
4Geriatric Research, Education, and Clinical Center, Baltimore Veterans Administration Medical Center, Baltimore, MD 21201, USA
Corresponding Author: Nicole Hoppman 971 Hilton 200 1st Street SW Rochester, MN 55905 (507) 538-8863; Fax: (507) 284-0670 ; hoppmanchaney.nicole/at/mayo.edu
Nicole Hoppman: hoppmanchaney.nicole/at/mayo.edu John McLenithan: jmcle001/at/umaryland.edu Daniel McBride: dmcbride/at/UDel.Edu Haiqing Shen: hshen/at/medicine.umaryland.edu Jan Bruder: bruder/at/uthscsa.edu Richard Bauer: bauer/at/uthscsa.edu John Shaffer: jrs51/at/pitt.edu Jie Liu: jliu/at/medicine.umaryland.edu Elizabeth Streeten: estreete/at/medicine.umaryland.edu Alan Shuldiner: ashuldin/at/medicine.umaryland.edu Candace Kammerer: cmk3/at/pitt.edu Braxton Mitchell: bmitchel/at/medicine.umaryland.edu
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 invitro 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.
Keywords: FGFBP1, polymorphisms, bone mineral density, osteoporosis, testosterone
Osteoporosis, a disease characterized by low bone mass and structural deterioration of bone, is a major threat to public health. It is estimated that approximately 10 million Americans are affected while another 34 million are at risk due to low BMD. Fractures resulting from osteoporosis are very costly in terms of public health dollars as well as quality of life. Osteoporosis is estimated to be responsible for more than 1.5 million fractures annually, of which 300,000 are fractures of the hip and 250,000 are fractures of the forearm. The estimated cost for osteoporotic hip fractures in the US was $18 billion dollars in 2002, and the cost is rising [1].
Because BMD is such an important predictor of future fracture, efforts are underway to understand the genetic and environmental factors that influence BMD. Twin and family studies suggest that BMD has a strong genetic component [2-5]. Several chromosomal regions potentially harboring quantitative trait loci (QTL) that influence BMD have been identified through linkage analysis, though few have been replicated across multiple studies. The paucity of replicated linkages is not surprising, since locus-specific effects on BMD are likely to be small as bone homeostasis is a complex process that is influenced by many factors. Nevertheless, determining the genetic causes and mechanism of this disease may make it possible to improve diagnosis, treatment, and offer better screening and prevention to at-risk individuals.
We previously performed a genome-wide linkage scan in 664 Mexican Americans comprising 29 large pedigrees and detected very strong evidence for linkage of forearm (radius midpoint) BMD to an 18 megabase region on chromosome 4p with the peak LOD near marker D4S2639 (LOD = 4.33, genome-wide p = 0.006) [6]. This region includes over 30 known genes and several hypothetical proteins.
An attractive candidate gene in this region is fibroblast growth factor binding protein 1 (FGFBP1). FGFBP1, also known as heparin-binding growth factor binding protein 17 (HBP17), is located at chromosome 4p15.32, which is directly under our linkage peak in Mexican Americans. The gene is approximately 2.8 kilobases long, and its 234 amino acid protein product is secreted from the cell. FGFBP1 has been shown to bind various fibroblast growth factors (FGFs), such as FGF-1 and FGF-2 [7], as well as FGFs 7, 10, and 22 [8]. Additional studies suggest that FGFBP1 has the ability to mobilize FGF-2 from the extracellular matrix and enhance the biological activities of both FGF-1 and FGF-2 [9, 10]. The ability of FGFBP1 to bind and modify the properties of FGF-2 is particularly relevant to the study of BMD for many reasons. FGF-2 is essential in regulating the differentiation, proliferation, and apoptosis of human osteoblasts [11]. FGF-2 binds cell surface receptors (FGFRs) on osteoblasts and activates mitogen-activated protein (MAP) kinase signal transduction pathways involving various MAP kinases [12-14]. Each of these pathways ultimately activates transcription of genes involved in osteoblast differentiation and function, such as type I collagen and alkaline phosphatase. Since FGFBP1 is secreted it has the potential to act on FGF-2 upstream of this pathway before it binds its receptor on osteoblast cell membranes.
In this follow-up to our linkage study, we report identification of several SNPs associated with BMD in our Mexican American population and replicate the association in an independent population. We then demonstrate by luciferase reporter gene assay that the three common haplotypes show marked differences in FGFBP1 expression with a promoter polymorphism responsible for abolishing this gene's response to testosterone in vitro. Therefore, variation in FGFBP1 may be important in expression of this gene and possibly subsequent risk of osteoporosis.
Mexican American Subjects
The San Antonio Family Osteoporosis Study (SAFOS) was initiated in 1997 to identify the genetic determinants of BMD in Mexican Americans. The baseline phase of the SAFOS included 897 individuals from 34 families (872 of whom had BMD measured at the hip and/or spine and are included in the present analyses) who were invited to participate because of their concurrent participation in a follow-up examination of the San Antonio Family Heart Study (SAFHS). Eligibility criteria for probands were that they be 40-60 years of age and have large families in the San Antonio area. First, second, and third degree relatives of each proband and the proband's spouse were invited to participate irrespective of medical history. Details of sampling and recruitment procedures were previously described [15].
Examinations were conducted at the General Clinical Research Center located at the Audie Murphy Veteran's Administration Hospital on the campus of the Health Science Center in San Antonio. Body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared. BMD was assessed using a Hologic Model 1500W dual energy X-ray absorptiometer (DXA, Hologic, Inc., Bedford, MA). Relevant to this report, we measured BMD at the forearm (radius one-third), lumbar spine (L1-L4), and hip (trochanter and femoral neck). Areal BMD (g/cm2) was determined by the manufacturer's software by dividing the bone mineral content (BMC, g) by the projected area of the region scanned (cm2). The short-term in vivo precision of the BMD was determined on 27 subjects who were examined twice on the same day. The precision of the lumbar spine was 0.009 g/cm2 (CV% = 1.0%) and precision of the total hip was 0.007 g/cm2 (CV% = 0.87%). The precision of the manufacturer's spine phantom was 0.0017 g/cm2 (CV% = 0.17%).
Amish Replication Subjects
The Amish Family Osteoporosis Study (AFOS) was also initiated in 1997 with the goal of identifying the genetic determinants of osteoporosis. Details of ascertainment, phenotyping and clinical characteristics of participants have been reported previously [16-19]. Ascertainment for this study was based on families with a proband believed to be at risk for osteoporosis by virtue of fracture history or prior bone mineral density measurements. The AFOS also comprised large extended families and included adults aged 20 years and older.
Exams were conducted at the Amish Research Clinic in Strasburg, PA. Height was measured by using a stadiometer and weight was recorded without shoes. BMI was calculated as weight in kilograms divided by height in meters squared. BMD was measured by DXA using a Hologic Model 4500 W (Hologic Inc., Bedford, MA) at the same sites as for the Mexican American study.
The protocols and procedures for the San Antonio Family Osteoporosis Study and Amish Family Osteoporosis Study were approved by the Institutional Review Boards at the respective institutions and all subjects gave written informed consent.
Sequencing and Genotyping
FGFBP1 is 2.8 kilobases in length and contains two exons. Primers were designed using Primer3 software (available at http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) to amplify approximately 1500 nucleotides of the promoter, both exons including 5′ and 3′ untranslated regions (UTR), and splice junctions. We sequenced PCR products from 40 Mexican Americans. To enhance our chances of detecting a highly penetrant rare allele, we selected 20 subjects with high BMD and 20 with low BMD from the families showing the greatest evidence of linkage to 4p. There was an equal number of males and females in the sequencing subset. The probability of not observing a variant in 40 individuals (80 chromosomes) whose frequency in the full set is 4% or higher was 3.8%. Both strands were sequenced from all samples on an ABI 3700 DNA sequencer and analyzed with Sequence Analysis 3.2 software (Applied Biosystems). We used Sequencher (GeneCodes) to align individual sequences and identify polymorphisms.
All SNPs were genotyped using SNPlex Genotyping System (Applied Biosystems) according to the manufacturer's protocol [20]. SNPs were analyzed using GeneMapper version 4.0 software (Applied Biosystems). Genotyping error rates ranged from 0-2% in Mexican Americans (n = 50 replicated samples) and from 0-1.5% in Amish (n = 180 replicated samples).
Association Analyses
All genotypes were checked for Mendelian errors using Pedcheck software [21] and errors were resolved or removed prior to analysis. Linkage disequilibrium (LD) among SNPs was assessed using the program Haploview [22].
Because both the Mexican American and Amish studies included large extended families, we utilized variance component methods to test for associations between each SNP genotype and BMD measures at different sites because these methods account for the relatedness among family members. Briefly, we used SOLAR to estimate the effects of genotype on BMD while incorporating the effects of sex, age, age squared, age by sex interactions, and BMI, as well as the residual correlations in BMD existing between related individuals. We initially tested for significance using additive genetic models, that is, linear allele dosage effects. When statistically significant genotype effects were observed, we then evaluated dominant and recessive genetic models that compared trait values between carriers and noncarriers of each allele. Only a recessive model was used in the Amish replication study based on findings in the Mexican American study population. Statistical significance was assessed by likelihood ratio tests in which we compared the likelihood of the data given the pedigree structure for a model in which genotype effects were estimated to that of a model in which genotype effects were constrained to be zero. Twice the difference between the logarithm of the likelihoods of the two models is distributed asymptotically as a χ2 statistic with degrees of freedom equal to the difference in the numbers of independent factors in the models being compared. In the initial Mexican American sample, we used a Bonferroni correction to account for the number of `independent' SNPs tested by excluding SNPs in high LD (r2 > 0.80) with another tested SNP. We did not adjust for the multiple BMD phenotypes tested because of the high correlations in BMD across sites.
Vector Preparation and Cloning
There were three common haplotypes at FGFBP1, which comprised 92% of haplotypes in Mexican Americans and 99% of haplotypes in the Amish. A 3.12 kb region of the gene (promoter, exon 1 (untranslated), intron 1, and first 7 bases of exon 2 prior to the translation start site) was amplified from DNA samples isolated from individuals homozygous for each haplotype. PCR was performed with Advantage 2 Polymerase Mix (BD Biosciences), according to the manufacturer's protocol. The amplified fragments and pGL3-Basic Reporter vector were digested with XhoI (New England Biolabs) and PCR products were cloned into the linearized vector using BD In-Fusion PCR Cloning Kit (BD Biosciences). Constructs were transformed into E. coli Fusion-Blue chemically competent cells (BD Biosciences) and isolated using QIAprep Spin Miniprep Kit (Qiagen). Constructs with single SNPs were created from the HapA construct using the QuikChange Site-Directed Mutagenesis Kit (Stratagene). All recombinant plasmid vectors were sequenced to confirm the sequence and orientation of inserts. Correct constructs were transformed into chemically competent E. coli cells and plasmids were isolated using HiSpeed Plasmid Maxi Kit (Qiagen).
Cell Culture and Transfection
Human bronchial epithelial cells (ATCC designation HBE135-E6E7), which were shown to express FGFBP1 at high levels by quantitative RT-PCR (data not shown), were obtained from American Type Culture Collection (ATCC). Osteoblast cell lines were not used because they did not express FGFBP1 at sufficient levels for in vitro studies to be performed (data not shown). Cells were subcultured and plated onto 12- or 24-well culture dishes (Costar) and grown in keratinocyte serum-free media (GIBCO) supplemented with 5 ng/mL human recombinant EGF (GIBCO), 0.05 mg/ml bovine pituitary extract (GIBCO), 0.005 mg/ml insulin (Sigma), and 500 ng/ml hydrocortisone (Sigma). Cells were incubated at 37°C with 95% air and 5% CO 2 and were transfected at 50-70% confluence. For each construct, at least two wells were transfected using FuGENE HD Reagent (Roche) according to manufacturer's protocol using a 6:2 ratio (6 uL reagent: 2 ug DNA). 0.1 ug of pRL-CMV was cotransfected per well in order to normalize for transfection efficiency. Cells were incubated for 2 days, after which media was removed and cells were lysed for reporter gene assays as described below. Dexamethasone (100 nM) (Sigma), 3,3′,5-Triiodo-L-thyronine sodium salt (thyroid hormone/T3, 50 nM) (Sigma), testosterone (100 nM) (Sigma), and 9-cis-retinoic acid (100 nM) (Sigma) were added to cells 24 hours post-transfection at the concentrations indicated; cells were lysed and assayed for luciferase activity 24 hours after addition.
Luciferase Activity
The Dual-Luciferase Reporter Assay System (Promega) was used to quantitate luciferase activity in transfected cells. Firefly luciferase activity was determined by mixing 20 uL of cell lysate with 100 uL of Luciferase Assay Reagent II. Relative light units were measured with a TD-20/20 Luminometer (Turner Designs). Renilla luciferase activity was then measured by the addition of 100 uL of Stop & Glo Reagent and subsequent measurement on the TD-20/20. Each well was read three times, and the average of the readings was used in statistical analyses. Relative luciferase activity was corrected for Renilla luciferase activity, expression levels among wells were standardized within experiments relative to HapA, and at least 3 experiments were combined for statistical analysis. The Student's T-test was used to test for expression differences between luciferase reporter gene constructs.
Quantitative RT-PCR (qRT-PCR)
Total RNA was extracted from HBE cells stimulated with testosterone or not for 24 hours (100nM) using RNeasy Mini Kits (Qiagen,). Total RNA was reverse transcribed using the Roche Transcriptor cDNA kit (Roche Applied Science). cDNA was quantified using the LightCycler 480, LightCycler 480 Probes Master kit (Roche Applied Science) and Taqman probe/primer sets for FGFBP1 and cyclophilin (Applied Biosystems). Normalized gene expression of FGFBP1 was expressed as a ratio of FGFBP1/cyclophilin mRNA.
SNP Discovery and Association with BMD in Mexican Americans
Characteristics of the 872 Mexican American (initial) and 972 Amish (replication) samples are given in Table 1. Compared to the Amish sample, the Mexican American sample was on average younger (42-43 versus 50-51 yrs) and heavier (mean BMI 29.5 and 31.1 kg/m2 in men and women, respectively versus. 26.4 and 28.0 kg/m2. Consistent with their younger age and higher BMI, Mexican Americans also had higher BMD at the spine and femoral neck.
Table 1
Table 1
Characteristics of Mexican American and Amish study subjects (mean ± SD).
Sequencing of FGFBP1 in 40 Mexican American subjects revealed 9 polymorphisms (Figure 1). All of the variants detected by sequencing were present in the HapMap database with the lone exception of rs2246196, which lies over 1kb upstream of the FGFBP1 promoter and is monomorphic in the Amish cohort. We genotyped 7 of the 9 SNPs in the full set of SAFOS participants; the two remaining SNPs, rs2246196 and rs732244, failed genotyping but were in near-perfect correlation (r2 > 0.90) with other genotyped SNPs in the 40 sequencing samples (rs2246196 in LD with rs2531174 and rs2072313; and rs732244 in LD with rs732245 and rs12503796). In the full SAFOS sample there was extensive linkage disequilibrium among the 7 genotyped SNPs, indicating that FGFBP1 comprised essentially one large haplotype block (Figure 2A).
Fig.1
Fig.1
Polymorphisms identified in Mexican Americans and the Amish. Exons are numbered below the gene schematic; black regions correspond to UTR and white regions correspond to coding sequence. Direction of transcription is indicated by arrow. Abbreviations: (more ...)
Fig. 2
Fig. 2
Linkage disequilibrium (D' top, r2 bottom) pattern for all SNPs genotyped in Mexican Americans (a) and the Amish (b). For both D' and r2, values range from 0 (no LD) to 100 (complete LD). Schematic of the FGFBP1 gene is shown above to indicate approximate (more ...)
Two of the seven genotyped SNPs (rs2531174 and rs2072313) had minor allele frequencies ≤ 0.02 and were not analyzed further (Figure 1). Associations of the remaining five SNPs with BMD are shown in Figure 3. In the full set of Mexican Americans, four of the five SNPs were significantly associated with BMD at one or more sites with statistically significant p-values for the recessive model ranging from 0.004 to 0.042. The two most strongly associated SNPs, rs2245964 and rs16892645, were in high LD (r2 = 0.86) and subjects homozygous for the minor allele (frequency = 0.18-0.20; n = 32 subjects for rs16892645) had higher BMD than subjects with one or two copies of the major allele. Two additional SNPs were also in high LD with each other (rs12503796 and rs732245), and so we considered our threshold p-value for statistical significance in the initial Mexican American cohort to be 0.05/3, or 0.017.
Fig. 3
Fig. 3
P-values reflecting association of FGFBP1 SNPs with BMD in Mexican Americans and Amish. P-values indicated correspond to a recessive genetic model. Covariates include age, age squared, sex, BMI (and in Mexican Americans diabetes status and physical activity (more ...)
Replication Study in Old Order Amish
We next sought to replicate the associations in an independent sample consisting of 972 Amish individuals. The LD patterns of the five common FGFBP1 SNPs in the Amish closely resembled those estimated in Mexican Americans (Figure 2B) as well as the HapMap CEPH population (data not shown). As shown in Figure 3, both SNPs were associated with trochanter BMD, again with subjects homozygous for the minor allele having higher BMD than subjects with one or two copies of the major allele (p=0.02 for each; 2-tailed p-value with adjustments for age, sex, and BMI). Modest associations (p < 0.05) were also observed in the Amish between the other three FGFBP1 SNPs and BMD at other sites.
Figure 4 demonstrates the association between genotype and BMD at the trochanter in both populations for one representative SNP (rs16892645), adjusted for age, sex, BMI, and other covariates, demonstrating a similar effect in both populations. This particular SNP accounts for 1.2% and 0.3% of the total residual variation in trochanter BMD after accounting for the effects of age, sex, and BMI (and other covariates) in Mexican Americans and Amish respectively.
Fig. 4
Fig. 4
Mean BMD (± standard error) at the trochanter according to rs16892645 genotype in Mexican Americans and the Amish. BMD is adjusted for age, sex, BMI, diabetes status, and metabolic equivalents in Mexican Americans, and for age, sex, and BMI in (more ...)
Because BMD is known to differ between men and women, we also assessed whether these SNPs were differentially associated with BMD in men and women in both the Mexican American and Amish cohorts. Sex-specific analyses in both populations revealed significant associations with BMD in males while very weak or even lacking in females. In Mexican Americans, for example, the associations with rs16892645 remained highly statistically significant in men despite a markedly reduced sample size (p: 0.037-0.001; n = 359), but not in women (p > 0.05 for all; n = 528) (Figure 3). However, only 15 males and 17 females in Mexican Americans, and only 3 males and 4 females in the Amish, had the high BMD genotype. A formal genotype × sex interaction was carried out in Mexican Americans and revealed no strong evidence supporting a stronger genotype effect on BMD in men than in women (interaction p-value > 0.05). There was insufficient power to detect a genotype × sex interaction involving rs16892645 in the Amish due to the low frequency of this SNP (MAF = 0.06) and small number of subjects with the high BMD-associated genotype.
In Vitro Functional Studies
In order to determine if one or more of the associated SNPs in FGFBP1 influence gene expression we created luciferase reporter gene constructs which included 1.6 kb of the promoter, the entire 5′ UTR, and entire intron of FGFBP1, totaling 3.12 kb in length. Constructs included the three common haplotypes formed by SNPs rs12503796, rs2245964, and rs16892645 in this region. These haplotypes, hereafter referred to as HapA (GCG), HapB (TCG), and HapC (GGA), occurred in frequencies of 0.31, 0.48, and 0.18 in the population. Because of the high LD existing between the UTR and intronic SNPs (rs2245964 and rs16892645), only the CG and GA haplotypes were observed at these two SNPs (Figure 5a). We observed significantly decreased luciferase expression with HapB as compared to HapA (p = 6.2 × 10−6, while HapC showed only slightly decreased expression (p = 0.017) (Figure 5b). Expression of HapB is approximately 25% lower than HapA. To confirm this finding, a construct containing only the promoter SNP (designated PRO), created by site-directed mutagenesis of HapA, also showed significantly decreased expression compared to HapA, confirming the effect of the promoter SNP on expression (p = 6.51 × 10−5, data not shown). Site-directed mutagenesis was also used to separate rs2245964 and rs16892645 in order to determine whether one or both influence expression when isolated; however we observed no significant difference in expression from HapA for either SNP (data not shown).
Fig. 5
Fig. 5
A. Diagram of native haplotype expression vector constructs. Differences from HapA are shown in red. B. Luciferase activity values (adjusted for Renilla luciferase activity and normalized to HapA) for native haplotypes (n = 11 wells). Bars represent standard (more ...)
Table 2 lists several transcription factors that are predicted to bind to the sequence containing the promoter polymorphism according to three different prediction software programs—Patch [23], Promo [24, 25], and TFsearch [26]. It is possible that one or more of these transcription factors may directly or indirectly influence bone homeostasis. To determine whether any of these factors does indeed bind to one of these polymorphic loci, we added several nuclear hormone receptor agonists to the HBE cells transfected with HapA, HapB, and HapC and measured luciferase expression. Three agonists, dexamethasone, thyroid hormone and retinoic acid, did not influence expression of the luciferase reporter gene. However, when testosterone was added to cultures, we observed an increase in luciferase expression for both HapA (p = 0.00087) and HapC (p = 0.0098), but not HapB (p = 0.752), implying that the promoter polymorphism (rs12503796) found in HapB may inhibit response to testosterone. To confirm this result, we performed the experiment with the PRO construct and observed the same lack of response to testosterone (p = 0.59, Figure 5c).
Table 2
Table 2
Potential transcription factor binding sites for sequences containing polymorphisms. Those relevant to bone homeostasis and therefore tested in vitro are shown in bold. No relevant sites were predicted for the two other polymorphisms tested (rs2245964 (more ...)
We next measured endogenous gene expression of FGFBP1 in HBE cells by qRT-PCR to verify that the regulation of the endogenous FGFBP1 gene by testosterone was similar to the observed regulation of FGFBP1 promoter/ luciferase constructs. FGFBP1 gene expression (as a ratio of cyclophilin mRNA) increased 1.5 fold over basal expression after addition of testosterone (0.93 ± 0.07 vs 1.43 ± 0.2, respectively, p < 0.05, n=12). Thus, testosterone regulation of the endogenous FGFBP1 gene was of similar magnitude and direction as the promoter/luciferase constructs (data not shown).
Despite the high heritability of BMD, efforts to identify individual genes that explain a large proportion of the variation in this trait have been largely unsuccessful. Numerous candidate genes have been proposed, some of which have been associated with BMD in multiple populations [27, 28], and recent genome-wide association studies have identified a small number of additional loci that appear to be robustly associated with BMD [29-34]. However, very few association studies have resulted in the identification of functional variants with demonstrable effects on gene expression. We have performed an association study in two very different populations, a genetically heterogenous population with a high degree of admixture and a genetically homogenous population of European ancestry, and followed up with functional studies of associated haplotypes. The fact that the FGFBP1 associations demonstrated in our study were observed across two populations, taken together with our in vitro studies showing a dramatic effect on gene expression, make this gene a very convincing candidate for bone homeostasis.
The mechanisms for differences in FGFBP1 expression are not clear; however, we provide evidence that the promoter polymorphism (rs12503796) may be involved. Rs12503796 is located in the promoter approximately 280 bases upstream of the transcription start site. It is not located in a known transcription factor binding site; however, it is a predicted binding site for several transcription factors that may be involved in bone metabolism. In this regard we did, in fact, observe enhanced in vitro luciferase expression when transfected cells were co-cultured with testosterone, although the testosterone-enhancing effect was obliterated in constructs having the rs12503796 T allele, suggesting that this allele may destroy a binding site for the androgen receptor. An attractive interpretation of this hypothesis is that loss of a testosterone-enhancing effect on BMD might be expected to have a greater effect on BMD in men than in women, which is also consistent with our observed association results. However, the “T” allele is predicted to bind the androgen receptor, while the “G” allele is not, so if the prediction is correct, other mechanisms to account for the differences in haplotype expression must be in play. For example, there may be other transcription factors not tested, either alone or through interaction with the androgen receptor, that mediate this effect. In fact, the promoter polymorphism is also a predicted binding site for the general transcription factor TFII-I and the rs12503796 polymorphism might alter TFII-I binding, possibly explaining why basal expression levels of the constructs were different in the absence of hormones. Finally, it is also possible that FGFBP1 expression differences may have nothing to do with the rs12503796 promoter polymorphism because FGFBP1 expression was decreased in HapC, despite the presence of the G allele at this promoter site. It is also possible that this effect is not true as this construct did not have significantly different expression from HapA in the presence of testosterone, which is likely more representative of true physiological conditions in which testosterone is always present at some level.
Despite the parallel associations observed in these data of FGFBP1 SNPs with BMD and of FGFBP1 haplotypes with gene expression, there are some limitations and inconsistencies in our study. First, Mexican Americans are an admixed population and our analysis of FGFBP1 SNPS did not account for population substructure. However, replication in a non-admixed population, the Amish, mitigates this limitation. Second, the expression experiments may not have been performed in the most physiologically relevant tissue. FGFBP1 expression in osteoblasts was measured; expression was detected but at insufficient amounts to carry out in vitro studies (data not shown). However, because FGFBP1 is secreted, it is possible that high expression in tissues other than bone could certainly influence bone mineral density. Third, the allele most strongly associated with higher BMD (e.g. rs16892645 A allele) is not by itself associated with gene expression. Fourth, the T allele at the promoter polymorphism (rs12503796) appears to block the increased expression of this gene in the presence of testosterone. How this potential decrease in activity and/or lack of stimulation leads to increased BMD is unclear since FGFBP1 is thought to enhance FGF-2 binding to its receptors, which should lead to osteoblast stimulation. Therefore, one might expect decreased expression or function of FGFBP1 to lead to decreased BMD. It is possible that in bone, FGFBP1 acts to sequester FGF-2 instead of enhancing binding. Indeed, some data suggests that while FGFBP1 stimulates FGF-2 at low concentrations, it can also inhibit FGF-2 at higher concentrations [7, 35]. Likewise, FGF-2 has also been found to have dual effects on osteoblasts. FGF-2 induces proliferation and differentiation in immature osteoblasts but has also been found to block mineralization and induce apoptosis in differentiated osteoblasts [13]. FGF-2 is also thought to influence osteoclast maturation and activity [36-40]; hence we would expect in this scenario that decreased expression of FGFBP1 would lead to increased BMD, as we observed, via decreased osteoclast activity. These hypotheses require further testing in order to determine the mechanism of this effect. Additionally, the cell type used in our in vitro experiments may not be representative of the pathway in relevant tissues in vitro.
The possibility that the effects of FGFBP1 SNPs may manifest differently between men and women is intriguing, but must be considered speculative. While the intronic rs16892645 variant was more strongly associated in Mexican American men than women, the estimated effect size in men was not statistically greater than that observed in women in a formal sex by SNP interaction analysis. Moreover, we had little power to test for sex differences in the association of this SNP with BMD in the Amish because of the low frequency of the rs16892645 variant in this population. Nevertheless, coupled with the observation that expression was significantly enhanced in the presence of testosterone, the possibility of a sex interaction makes an interesting hypothesis. Future studies of FGFBP1 should consider the possibility of differential effects between men and women.
It seems unlikely that variation in FGFBP1 is responsible for our previously observed linkage peak due to the modest, albeit replicated, associations observed with BMD. Although linkage analysis can detect rare genetic effects, our prior linkage was exclusively with forearm BMD, while the current analyses indicate association across multiple BMD sites, particularly at the hip. It is unclear why our previous linkage study detected little evidence for linkage to hip BMD; likewise, despite modest association of rs2245964 and rs16892645 with forearm BMD, a subsequent (not shown) linkage analysis indicated essentially no change in the evidence for linkage to forearm BMD after accounting for the associated SNPs from this study. Other genes in this region must therefore be considered as candidates for their role in driving the forearm BMD linkage result.
In summary, we have identified several polymorphisms in FGFBP1 that are associated with BMD at several sites in both Mexican Americans and the Old Order Amish. The minor alleles of these polymorphisms tend to be associated with higher BMD in both populations and there is a suggestion that effects may be stronger in males than in females. A promoter polymorphism (rs12503796) was shown to influence gene expression in vitro, possibly due to alteration of a transcription factor binding site as demonstrated by abolishment of response to testosterone in vitro. However, while these data provide evidence that FGFBP1 is involved in bone remodeling, the precise mechanism governing the regulation of gene expression is not clear. This mechanism likely involves an interaction between sequence variation and hormonal or other influences on transcription. Further elucidation of the factors regulating expression of this gene and the downstream effects of altered gene expression of skeletal phenotypes seems warranted.
Acknowledgements
We are deeply grateful for the cooperation of the families participating in the SAFOS and AFOS. The San Antonio Family Osteoporosis Study is supported by research grants RO1-AR43351, RO1- HL69313, and PO1-HL45522 awarded by the National Institutes of Health. Support for the Frederic C. Bartter General Clinical Research Center was made available by Clinical National Institutes of Health Grant MO1-RR-01346. The Amish Family Osteoporosis Study is supported by research grants R01-AR46838 awarded by the National Institutes of Health, University of Maryland General Clinical Research Center Grant M01 RR165001, General Clinical Research Centers Program, National Center for Research Resources (NCRR), and the Geriatrics Research and Education Clinical Center, Baltimore Veterans Administration Medical Center. We would like to thank Faith Rowland for help with sequencing and Jessica Pray and Samuel Woodle for help with the functional analyses.
Footnotes
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
Conflict of Interest: All authors have no conflicts of interest.
1. Bone Health and Osteoporosis: A Report of the Surgeon General. U.S. Department of Health and Human Services, Office of the Surgeon General; Rockville, MD: 2004.
2. Arden NK, Spector TD. Genetic influences on muscle strength, lean body mass, and bone mineral density: a twin study. J Bone Miner Res. 1997;12:2076–81. [PubMed]
3. Dequeker J, Nijs J, Verstraeten A, Geusens P, Gevers G. Genetic determinants of bone mineral content at the spine and radius: a twin study. Bone. 1987;8:207–9. [PubMed]
4. Pocock NA, Eisman JA, Hopper JL, Yeates MG, Sambrook PN, Eberl S. Genetic determinants of bone mass in adults. A twin study. J Clin Invest. 1987;80:706–10. [PMC free article] [PubMed]
5. Smith DM, Nance WE, Kang KW, Christian JC, Johnston CC., Jr Genetic factors in determining bone mass. J Clin Invest. 1973;52:2800–8. [PMC free article] [PubMed]
6. Kammerer CM, Schneider JL, Cole SA, Hixson JE, Samollow PB, O'Connell JR, Perez R, Dyer TD, Almasy L, Blangero J, Bauer RL, Mitchell BD. Quantitative trait loci on chromosomes 2p, 4p, and 13q influence bone mineral density of the forearm and hip in Mexican Americans. J Bone Miner Res. 2003;18:2245–52. [PubMed]
7. Wu DQ, Kan MK, Sato GH, Okamoto T, Sato JD. Characterization and molecular cloning of a putative binding protein for heparin-binding growth factors. J Biol Chem. 1991;266:16778–85. [PubMed]
8. Beer HD, Bittner M, Niklaus G, Munding C, Max N, Goppelt A, Werner S. The fibroblast growth factor binding protein is a novel interaction partner of FGF-7, FGF-10 and FGF-22 and regulates FGF activity: implications for epithelial repair. Oncogene. 2005;24:5269–77. [PubMed]
9. Aigner A, Butscheid M, Kunkel P, Krause E, Lamszus K, Wellstein A, Czubayko F. An FGF-binding protein FGF-BP exerts its biological function by parallel paracrine stimulation of tumor cell and endothelial cell proliferation through FGF-2 release. Int J Cancer. 2001;92:510–7. [PubMed]
10. Tassi E, Al-Attar A, Aigner A, Swift MR, McDonnell K, Karavanov A, Wellstein A. Enhancement of fibroblast growth factor FGF activity by an FGF-binding protein. J Biol Chem. 2001;276:40247–53. [PubMed]
11. Marie PJ, Debiais F, Hay E. Regulation of human cranial osteoblast phenotype by FGF-2, FGFR-2 and BMP-2 signaling. Histol Histopathol. 2002;17:877–85. [PubMed]
12. Debiais F, Lemonnier J, Hay E, Delannoy P, Caverzasio J, Marie PJ. Fibroblast growth factor-2 FGF-2 increases N-cadherin expression through protein kinase C and Src-kinase pathways in human calvaria osteoblasts. J Cell Biochem. 2001;81:68–81. [PubMed]
13. Mansukhani A, Bellosta P, Sahni M, Basilico C. Signaling by fibroblast growth factors FGF and fibroblast growth factor receptor 2 FGFR2-activating mutations blocks mineralization and induces apoptosis in osteoblasts. J Cell Biol. 2000;149:1297–308. [PMC free article] [PubMed]
14. Newberry EP, Willis D, Latifi T, Boudreaux JM, Towler DA. Fibroblast growth factor receptor signaling activates the human interstitial collagenase promoter via the bipartite Ets-AP1 element. Mol Endocrinol. 1997;11:1129–44. [PubMed]
15. Mitchell BD, Kammerer CM, Schneider JL, Perez R, Bauer RL. Genetic and environmental determinants of bone mineral density in Mexican Americans: results from the San Antonio Family Osteoporosis Study. Bone. 2003;33:839–46. [PubMed]
16. Brown LB, Streeten EA, Shapiro JR, McBride D, Shuldiner AR, Peyser PA, Mitchell BD. Genetic and environmental influences on bone mineral density in pre- and post-menopausal women. Osteoporos Int. 2005;16:1849–56. [PubMed]
17. Brown LB, Streeten EA, Shuldiner AR, Almasy LA, Peyser PA, Mitchell BD. Assessment of sex-specific genetic and environmental effects on bone mineral density. Genet Epidemiol. 2004;27:153–61. [PubMed]
18. Streeten EA, McBride DJ, Lodge AL, Pollin TI, Stinchcomb DG, Agarwala R, Schaffer AA, Shapiro JR, Shuldiner AR, Mitchell BD. Reduced incidence of hip fracture in the Old Order Amish. J Bone Miner Res. 2004;19:308–13. [PubMed]
19. Streeten EA, McBride DJ, Pollin TI, Ryan K, Shapiro J, Ott S, Mitchell BD, Shuldiner AR, O'Connell JR. Quantitative trait loci for BMD identified by autosome-wide linkage scan to chromosomes 7q and 21q in men from the Amish Family Osteoporosis Study. J Bone Miner Res. 2006;21:1433–42. [PubMed]
20. De la Vega FM, Lazaruk KD, Rhodes MD, Wenz MH. Assessment of two flexible and compatible SNP genotyping platforms: TaqMan SNP Genotyping Assays and the SNPlex Genotyping System. Mutat Res. 2005;573:111–35. [PubMed]
21. O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998;63:259–66. [PubMed]
22. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5. [PubMed]
23. PATCH: Pattern Search for Transcription Factor Binding Sites. Available at: http://www.gene-regulation.com/pub/programs.html. Last accessed January 2007.
24. Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, Messeguer X. Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003;31:3651–3. [PMC free article] [PubMed]
25. Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM. PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002;18:333–4. [PubMed]
26. Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel OV, Ignatieva EV, Ananko EA, Podkolodnaya OA, Kolpakov FA, Podkolodny NL, Kolchanov NA. Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998;26:362–7. [PMC free article] [PubMed]
27. Ralston SH. Genetic determinants of osteoporosis. Curr Opin Rheumatol. 2005;17:475–9. [PubMed]
28. Liu YJ, Shen H, Xiao P, Xiong DH, Li LH, Recker RR, Deng HW. Molecular genetic studies of gene identification for osteoporosis: a 2004 update. J Bone Miner Res. 2006;21:1511–35. [PMC free article] [PubMed]
29. Cho YS, Go MJ, Kim YJ, Heo JY, Oh JH, Ban HJ, Yoon D, Lee MH, Kim DJ, Park M, Cha SH, Kim JW, Han BG, Min H, Ahn Y, Park MS, Han HR, Jang HY, Cho EY, Lee JE, Cho NH, Shin C, Park T, Park JW, Lee JK, Cardon L, Clarke G, McCarthy MI, Lee JY, Lee JK, Oh B, Kim HL. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet. 2009;41:527–34. [PubMed]
30. Kiel DP, Demissie S, Dupuis J, Lunetta KL, Murabito JM, Karasik D. Genome-wide association with bone mass and geometry in the Framingham Heart Study. BMC Med Genet. 2007;8(Suppl 1):S14. [PMC free article] [PubMed]
31. Richards JB, Rivadeneira F, Inouye M, Pastinen TM, Soranzo N, Wilson SG, Andrew T, Falchi M, Gwilliam R, Ahmadi KR, Valdes AM, Arp P, Whittaker P, Verlaan DJ, Jhamai M, Kumanduri V, Moorhouse M, van Meurs JB, Hofman A, Pols HA, Hart D, Zhai G, Kato BS, Mullin BH, Zhang F, Deloukas P, Uitterlinden AG, Spector TD. Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study. Lancet. 2008;371:1505–12. [PMC free article] [PubMed]
32. Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, Ingvarsson T, Jonsdottir T, Saemundsdottir J, Center JR, Nguyen TV, Bagger Y, Gulcher JR, Eisman JA, Christiansen C, Sigurdsson G, Kong A, Thorsteinsdottir U, Stefansson K. Multiple genetic loci for bone mineral density and fractures. N Engl J Med. 2008;358:2355–65. [PubMed]
33. Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, Ingvarsson T, Jonsdottir T, Saemundsdottir J, Snorradottir S, Center JR, Nguyen TV, Alexandersen P, Gulcher JR, Eisman JA, Christiansen C, Sigurdsson G, Kong A, Thorsteinsdottir U, Stefansson K. New sequence variants associated with bone mineral density. Nat Genet. 2009;41:15–7. [PubMed]
34. Xiong DH, Liu XG, Guo YF, Tan LJ, Wang L, Sha BY, Tang ZH, Pan F, Yang TL, Chen XD, Lei SF, Yerges LM, Zhu XZ, Wheeler VW, Patrick AL, Bunker CH, Guo Y, Yan H, Pei YF, Zhang YP, Levy S, Papasian CJ, Xiao P, Lundberg YW, Recker RR, Liu YZ, Liu YJ, Zmuda JM, Deng HW. Genome-wide association and follow-up replication studies identified ADAMTS18 and TGFBR3 as bone mass candidate genes in different ethnic groups. Am J Hum Genet. 2009;84:388–98. [PubMed]
35. Chen JH, Wang XC, Sato JD. Bifunctional effects of heparin-binding protein HBp17 on DNA synthesis in cells. Cell Biol Int. 2001;25:567–70. [PubMed]
36. Chikazu D, Hakeda Y, Ogata N, Nemoto K, Itabashi A, Takato T, Kumegawa M, Nakamura K, Kawaguchi H. Fibroblast growth factor FGF-2 directly stimulates mature osteoclast function through activation of FGF receptor 1 and p42/p44 MAP kinase. J Biol Chem. 2000;275:31444–50. [PubMed]
37. Kawaguchi H, Chikazu D, Nakamura K, Kumegawa M, Hakeda Y. Direct and indirect actions of fibroblast growth factor 2 on osteoclastic bone resorption in cultures. J Bone Miner Res. 2000;15:466–73. [PubMed]
38. Nakano K, Okada Y, Saito K, Tanaka Y. Induction of RANKL expression and osteoclast maturation by the binding of fibroblast growth factor 2 to heparan sulfate proteoglycan on rheumatoid synovial fibroblasts. Arthritis Rheum. 2004;50:2450–8. [PubMed]
39. Shimoaka T, Ogasawara T, Yonamine A, Chikazu D, Kawano H, Nakamura K, Itoh N, Kawaguchi H. Regulation of osteoblast, chondrocyte, and osteoclast functions by fibroblast growth factor FGF-18 in comparison with FGF-2 and FGF-10. J Biol Chem. 2002;277:7493–500. [PubMed]
40. Zuo J, Jiang J, Dolce C, Holliday LS. Effects of basic fibroblast growth factor on osteoclasts and osteoclast-like cells. Biochem Biophys Res Commun. 2004;318:162–7. [PubMed]