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

Meta-analysis of genome-wide studies identifies WNT16 and ESR1 SNPs associated with bone mineral density in premenopausal women

Previous genome-wide association studies (GWAS) have identified common variants in genes associated with variation in bone mineral density (BMD), although most have been carried out in combined samples of older women and men. Meta-analyses of these results have identified numerous single-nucleotide polymorphisms (SNPs) of modest effect at genome-wide significance levels in genes involved in both bone formation and resorption, as well as other pathways. We performed a meta-analysis restricted to premenopausal white women from four cohorts (n = 4061 women, aged 20 to 45 years) to identify genes influencing peak bone mass at the lumbar spine and femoral neck. After imputation, age- and weight-adjusted bone-mineral density (BMD) values were tested for association with each SNP. Association of an SNP in the WNT16 gene (rs3801387; p = 1.7 × 10(-9) ) and multiple SNPs in the ESR1/C6orf97 region (rs4870044; p = 1.3 × 10(-8) ) achieved genome-wide significance levels for lumbar spine BMD. These SNPs, along with others demonstrating suggestive evidence of association, were then tested for association in seven replication cohorts that included premenopausal women of European, Hispanic-American, and African-American descent (combined n = 5597 for femoral neck; n = 4744 for lumbar spine). When the data from the discovery and replication cohorts were analyzed jointly, the evidence was more significant (WNT16 joint p = 1.3 × 10(-11) ; ESR1/C6orf97 joint p = 1.4 × 10(-10) ). Multiple independent association signals were observed with spine BMD at the ESR1 region after conditioning on the primary signal. Analyses of femoral neck BMD also supported association with SNPs in WNT16 and ESR1/C6orf97 (p < 1 × 10(-5) ). Our results confirm that several of the genes contributing to BMD variation across a broad age range in both sexes have effects of similar magnitude on BMD of the spine in premenopausal women. These data support the hypothesis that variants in these genes of known skeletal function also affect BMD during the premenopausal period.

Rare coding variants in ALPL are associated with low serum alkaline phosphatase and low bone mineral density

Alkaline phosphatase (ALP) plays an essential role in the regulation of tissue mineralization, and its activity is highly heritable. Guided by genetic associations discovered in a murine model, we hypothesized a role for rare coding variants in determining serum ALP level and bone mineral density (BMD) in humans. We sequenced the coding regions of the ALP gene (ALPL) in men with low and normal serum ALP activity levels. Single-nucleotide ALPL variants, including 19 rare nonsynonymous variants (minor allele frequency <1%), were much more frequent among the low ALP group (33.8%) than the normal group (1.4%, p = 1 × 10(-11)). Within the low ALP group, men with a rare, nonsynonymous variant had 11.2% lower mean serum ALP (p = 3.9 × 10(-4)), 6.7% lower BMD (p = 0.03), and 11.1% higher serum phosphate (p = 0.002) than those without. In contrast, common nonsynonymous variants had no association with serum ALP, phosphate, or BMD. Multiple rare ALPL coding variants are present in the general population, and nonsynonymous coding variants may be responsible for heritable differences in mineralization and thus BMD.

Molecular genetic studies of gene identification for osteoporosis: the 2009 update

Osteoporosis is a complex human disease that results in increased susceptibility to fragility fractures. It can be phenotypically characterized using several traits, including bone mineral density, bone size, bone strength, and bone turnover markers. The identification of gene variants that contribute to osteoporosis phenotypes, or responses to therapy, can eventually help individualize the prognosis, treatment, and prevention of fractures and their adverse outcomes. Our previously published reviews have comprehensively summarized the progress of molecular genetic studies of gene identification for osteoporosis and have covered the data available to the end of September 2007. This review represents our continuing efforts to summarize the important and representative findings published between October 2007 and November 2009. The topics covered include genetic association and linkage studies in humans, transgenic and knockout mouse models, as well as gene-expression microarray and proteomics studies. Major results are tabulated for comparison and ease of reference. Comments are made on the notable findings and representative studies for their potential influence and implications on our present understanding of the genetics of osteoporosis.

PedHunter 2.0 and its usage to characterize the founder structure of the Old Order Amish of Lancaster County

BACKGROUND

Because they are a closed founder population, the Old Order Amish (OOA) of Lancaster County have been the subject of many medical genetics studies. We constructed four versions of Anabaptist Genealogy Database (AGDB) using three sources of genealogies and multiple updates. In addition, we developed PedHunter, a suite of query software that can solve pedigree-related problems automatically and systematically.

METHODS

We report on how we have used new features in PedHunter to quantify the number and expected genetic contribution of founders to the OOA. The queries and utility of PedHunter programs are illustrated by examples using AGDB in this paper. For example, we calculated the number of founders expected to be contributing genetic material to the present-day living OOA and estimated the mean relative founder representation for each founder. New features in PedHunter also include pedigree trimming and pedigree renumbering, which should prove useful for studying large pedigrees.

RESULTS

With PedHunter version 2.0 querying AGDB version 4.0, we identified 34,160 presumed living OOA individuals and connected them into a 14-generation pedigree descending from 554 founders (332 females and 222 males) after trimming. From the analysis of cumulative mean relative founder representation, 128 founders (78 females and 50 males) accounted for over 95% of the mean relative founder contribution among living OOA descendants.

DISCUSSION/CONCLUSIONS

The OOA are a closed founder population in which a modest number of founders account for the genetic variation present in the current OOA population. Improvements to the PedHunter software will be useful in future studies of both the OOA and other populations with large and computerized genealogies.

Functional and association analysis of frizzled 1 (FZD1) promoter haplotypes with femoral neck geometry

Frizzleds are receptors for Wnt signaling and are involved in skeletal morphogenesis. Little is known about the transcriptional regulation of frizzleds in bone cells. In the current study, we determined if two common and potentially functional genetic variants (rs2232157, rs2232158) in the frizzled-1 (FZD1) promoter region and their haplotypes influence FZD1 promoter activity in human osteoblast-like cells. We also determined if these variants are associated with femoral neck bone mineral density (BMD) and geometry in 1319 African ancestry men aged > or =40 years. Real-time quantitative PCR and western blot analysis demonstrated FZD1 mRNA and protein expression in the human osteoblast-like cell lines, MG63 and SaOS-2. Promoter activity was next assessed by transient expression of haplotype specific FZD1 promoter reporter plasmids in these cells. In comparison to the common GG haplotype, promoter activity was 3-fold higher for the TC haplotype in both cell lines (p<0.05). We previously demonstrated that rs2232158 is associated with differential FZD1 promoter activity and Egr1 binding and thus focused further functional analyses on the rs2232157 G-to-T polymorphism. Electrophoretic mobility shift assay demonstrated that distinct nuclear protein complexes were associated with rs2232157 in an allele specific manner. Bioinformatics analysis predicted that the G to T transversion creates an E2F1 binding site, further supporting the functional significance of rs2232157 in FZD1 promoter regulation. Individual SNPs and haplotypes were not associated with femoral neck BMD. The TC haplotype was associated with larger subperiosteal width and greater CSMI (p<0.05). These results suggest that FZD1 expression is regulated in a haplotype-dependent manner in osteoblasts and that these same haplotypes may be associated with biomechanical indices of bone strength.

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.

Quantifying the material and structural determinants of bone strength

The ability of a bone to resist fracture depends on the amount of bone present, the spatial distribution of the bone mass as cortical and trabecular bone and the intrinsic properties of the bone material. Whereas low areal bone mineral density (aBMD) predicts fractures, its sensitivity and specificity is low, as over 50% of fractures occur in persons without osteoporosis by BMD testing and most women with osteoporosis do not sustain a fracture. New non-invasive imaging techniques, including three-dimensional (3D) assessments of bone density and geometry, microarchitecture and integrated measurements of bone strength such as finite element analysis (FEA), provide estimates of bone strength that can be used to increase the sensitivity and specificity of fracture risk assessment. Initial observations have shown that these techniques provide information that will improve our understanding of the pathophysiology of skeletal fragility and suggest that these techniques are likely to have a role in the clinical management of individuals at risk for fracture.

Differentially expressed genes strongly correlated with femur strength in rats

The region of chromosome 1q33-q54 harbors quantitative trait loci (QTL) for femur strength in COPxDA and F344xLEW F2 rats. The purpose of this study is to identify the genes within this QTL region that contribute to the variation in femur strength. Microarray analysis was performed using RNA extracted from femurs of COP, DA, F344 and LEW rats. Genes differentially expressed in the 1q33-q54 region among these rat strains were then ranked based on the strength of correlation with femur strength in F2 animals derived from these rats. A total of 214 genes in this QTL region were differentially expressed among all rat strains, and 81 genes were found to be strongly correlated (r(2)>0.50) with femur strength. Of these, 12 candidate genes were prioritized for further validation, and 8 of these genes (Ifit3, Ppp2r5b, Irf7, Mpeg1, Bloc1s2, Pycard, Sec23ip, and Hps6) were confirmed by quantitative PCR (qPCR). Ingenuity Pathway Analysis suggested that these genes were involved in interferon alpha, nuclear factor-kappa B (NFkB), extracellular signal-related kinase (ERK), hepatocyte nuclear factor 4 alpha (HNF4A) and tumor necrosis factor (TNF) pathways.