A follow-up association study of two genetic variants for bone mineral density variation in Caucasians

Wnt16 Promotes Vascular Smooth Muscle Contractile Phenotype and Function via Taz (Wwtr1) Activation in Male LDLR-/- Mice

Wnt16 is expressed in bone and arteries, and maintains bone mass in mice and humans, but its role in cardiovascular physiology is unknown. We show that Wnt16 protein accumulates in murine and human vascular smooth muscle (VSM). WNT16 genotypes that convey risk for bone frailty also convey risk for cardiovascular events in the Dallas Heart Study. Murine Wnt16 deficiency, which causes postnatal bone loss, also reduced systolic blood pressure. Electron microscopy demonstrated abnormal VSM mitochondrial morphology in Wnt16-null mice, with reductions in mitochondrial respiration. Following angiotensin-II (AngII) infusion, thoracic ascending aorta (TAA) dilatation was greater in Wnt16-/- vs Wnt16+/+ mice (LDLR-/- background). Acta2 (vascular smooth muscle alpha actin) deficiency has been shown to impair contractile phenotype and worsen TAA aneurysm with concomitant reductions in blood pressure. Wnt16 deficiency reduced expression of Acta2, SM22 (transgelin), and other contractile genes, and reduced VSM contraction induced by TGFβ. Acta2 and SM22 proteins were reduced in Wnt16-/- VSM as was Ankrd1, a prototypic contractile target of Yap1 and Taz activation via TEA domain (TEAD)-directed transcription. Wnt16-/- VSM exhibited reduced nuclear Taz and Yap1 protein accumulation. SiRNA targeting Wnt16 or Taz, but not Yap1, phenocopied Wnt16 deficiency, and Taz siRNA inhibited contractile gene upregulation by Wnt16. Wnt16 incubation stimulated mitochondrial respiration and contraction (reversed by verteporfin, a Yap/Taz inhibitor). SiRNA targeting Taz inhibitors Ccm2 and Lats1/2 mimicked Wnt16 treatment. Wnt16 stimulated Taz binding to Acta2 chromatin and H3K4me3 methylation. TEAD cognates in the Acta2 promoter conveyed transcriptional responses to Wnt16 and Taz. Wnt16 regulates cardiovascular physiology and VSM contractile phenotype, mediated via Taz signaling.

Multiple Mechanisms Explain Genetic Effects at the CPED1-WNT16 Bone Mineral Density Locus

PURPOSE OF REVIEW

Chromosome region 7q31.31, also known as the CPED1-WNT16 locus, is robustly associated with BMD and fracture risk. The aim of the review is to highlight experimental studies examining the function of genes at the CPED1-WNT16 locus.

RECENT FINDINGS

Genes that reside at the CPED1-WNT16 locus include WNT16, FAM3C, ING3, CPED1, and TSPAN12. Experimental studies in mice strongly support the notion that Wnt16 is necessary for bone mass and strength. In addition, roles for Fam3c and Ing3 in regulating bone morphology in vivo and/or osteoblast differentiation in vitro have been identified. Finally, a role for wnt16 in dually influencing bone and muscle morphogenesis in zebrafish has recently been discovered, which has brought forth new questions related to whether the influence of WNT16 in muscle may conspire with its influence in bone to alter BMD and fracture risk.

Integrating genome-wide association study with regulatory SNP annotations identified novel candidate genes for osteoporosis

Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p_{discovery GWAS} = 1.21 × 10^-25, p_{replication GWAS} = 1.80 × 10^-12), CCDC170 (p_{discovery GWAS} = 1.23 × 10^-11, p_{replication GWAS} = 3.22 × 10^-9), and SOX6 (p_{discovery GWAS} = 4.41 × 10^-15, p_{replication GWAS} = 6.57 × 10^-14). Then, for the 36 common target genes, multiple gene ontology (GO) terms were detected for BMD such as positive regulation of cartilage development (p = 9.27 × 10^-3) and positive regulation of chondrocyte differentiation (p = 9.27 × 10^-3). We explored the potential roles of rSNP in the genetic mechanisms of BMD and identified multiple candidate genes. Our study results support the implication of regulatory genetic variants in the development of OP.

Associations between WNT signaling pathway-related gene polymorphisms and risks of osteoporosis development in Chinese postmenopausal women: a case-control study

BACKGROUND

The WNT signaling pathway is involved in the regulation of bone homeostasis, and the effect of WNT signaling pathway-related gene (WNT16 and LRP5) polymorphisms on osteoporosis risk among Chinese postmenopausal women is still unknown. Hence, we performed a case-control study to assess the association of WNT signaling pathway-related gene polymorphisms and osteoporosis risk.

METHODS

A total of 1026 women (515 osteoporosis patients and 511 controls) of postmenopausal age who were randomly sampled from Xi'an 630 Hospital (Shaanxi Province, China) were involved in this study. Seven genetic polymorphisms in WNT16 (rs3779381, rs3801387, rs917727 and rs7776725) and LRP5 (rs2291467, rs11228240 and rs12272917) were selected and genotyped using the Agena MassARRAY iPLEX system. The association of the genetic polymorphisms and osteoporosis risk was assessed by odds ratios and 95% confidence intervals. The multifactor dimensionality reduction (MDR) method was conducted to analyze single nucleotide polymorphism (SNP)-SNP interaction.

RESULTS

We found that LRP5 polymorphisms (rs2291467, rs11228240 and rs12272917) were significantly associated with a decreased risk of osteoporosis in homozygote, recessive and additive models (p < 0.05). Stratification analysis showed that LRP5 polymorphisms (rs2291467, rs11228240 and rs12272917) significantly decreased the osteoporosis risk in the subgroup of body mass index (BMI) ≤ 24 (p < 0.05) and that individuals carrying a heterozygote genotype of WNT16 polymorphisms (rs3779381, rs3801387, rs917727 and rs7776725) had a higher osteoporosis risk in the subgroup of BMI > 24 (p < 0.05). Two haplotypes (haplotype 1: rs3779381, rs3801387, rs917727 and rs7776725; haplotype 2: rs2291467 and rs11228240) were observed, yet only T_rs2291467T_rs11228240 and C_rs2291467C_rs11228240 had a strong association with a decreased risk of osteoporosis (p < 0.05). Additionally, MDR analysis revealed that LRP5 rs2291467 was the best model in single-locus MDR analysis. A seven-locus model including rs3779381-AG, rs7776725-TC, rs3801387-GA and rs917727-TC in WNT16 and rs11228240-CC, rs12272917-TC and rs2291467-CC in LRP5 was the best model in multiple-loci MDR analysis (p < 0.001). These two best models were the most significantly associated with osteoporosis risk.

CONCLUSIONS

Our findings suggested that WNT16 and LRP5 genetic polymorphisms are associated with osteoporosis risk among Chinese postmenopausal women.

FAM3C: an emerging biomarker and potential therapeutic target for cancer

FAM3C is a member of the FAM3 family. Recently, overexpression of FAM3C has been reported in numerous types of cancer, including breast and colon cancer. Increasing evidence suggests that elevated FAM3C and its altered subcellular localization are closely associated with tumor formation, invasion, metastasis and poor survival. Moreover, FAM3C has been found to be the regulator of various proteins that associate with cancer, including Ras, STAT3, TGF-β and LIFR. This review summarizes the current knowledge regarding FAM3C, including its structure, expression patterns, regulation, physiological roles and regulatory functions in various malignancies. These findings highlight the importance of FAM3C in cancer development and provide evidence that FAM3C is a novel biomarker and potential therapeutic target for various cancers.

Genetic variants in the FAM3C gene are associated with lipid traits in Chinese children

BACKGROUND

Previous studies have related FAM3C gene with childhood bone health, and the regulation of lipid metabolism in hepatocytes. The present case-control study aimed to analyze the association of FAM3C genetic variants with overweight/obesity and lipid traits among Chinese children.

METHODS

Two genetic variants (rs7776725 and rs7793554) within the FAM3C gene were genotyped in 3305 Chinese children aged 6-18 years.

RESULTS

In the whole study population, the T-allele of rs7776725 and A-allele of rs7793554 within the FAM3C gene were associated with 40.2% (95% CI: 11.6-76.1%; P = 0.004) and 29.1% (6.9-56.0%; P = 0.008) increased risk of dyslipidemia, higher triglyceride (P = 0.014 and P = 0.001) and lower HDL-C (P = 0.015 and P = 0.003). In addition, we found that rs7776725 interacted with sex on dyslipidemia (P_{for interaction} = 0.004), and sex-stratified analyses showed that it was significantly associated with dyslipidemia only in girls (P = 8.78 × 10^-5). The variant also showed nominally significant interactions with sex on total cholesterol and LDL-C (P_{for interaction} = 0.012 and 0.008).

CONCLUSION

We found that FAM3C genetic variants were associated with dyslipidemia and lipid traits among Chinese children. In addition, we found significant gene-by-sex interactions. Our findings provided evidence supporting the role of FAM3C gene in regulating lipid metabolism in humans.

IMPACT

FAM3C genetic variants were associated with dyslipidemia and lipid traits among Chinese children. In addition, we found significant gene-by-sex interactions. FAM3C/rs7776725 was associated with dyslipidemia and lipid traits only in girls. Our findings provided evidence supporting the role of FAM3C gene in regulating lipid metabolism in humans.

Fam3c modulates osteogenic differentiation by down-regulating Runx2

Fam3c, a cytokine-like protein, is a member of the Fam3 family (family with sequence similarity 3) and has been implicated to play a crucial role in Epithelial-to- mesenchymal transition (EMT) and subsequent metastasis during cancer progression. A few independent genome-wide association studies on different population cohorts predicted the gene locus of Fam3c to be associated with bone mineral density and fractures. In this study, we examined the role of Fam3c during osteoblast differentiation. Fam3c was found to be expressed during osteogenic differentiation of both primary bone marrow stromal cells and MC3T3-E1 pre-osteoblasts. In differentiating osteoblasts, knockdown of Fam3c increased alkaline phosphatase expression and activity whereas overexpression of Fam3c reduced it. Furthermore, overexpression of Fam3c caused reduction of Runx2 expression at both mRNA and protein levels. Fam3c was localized in the cytoplasm and it was not secreted outside the cell during osteoblast differentiation and therefore, may function intracellularly. Furthermore, Fam3c and TGF-β1 were found to regulate each other reciprocally. Our findings therefore suggest a functional role of Fam3c in the regulation of osteoblast differentiation.

Fam3c modulates osteogenic cell differentiation and affects bone volume and cortical bone mineral density

Fam3c, a cytokine-like growth factor, has been suggested to have a role in epithelial-to-mesenchymal transition (EMT), tumor growth and metastasis. A single-nucleotide polymorphism affecting bone mineral density has been found in the first intron of the Fam3c gene in a study analyzing an Asian population cohort. Other independent studies on different population cohorts have found the fam3c locus to be associated with bone mineral density and fractures. In order to investigate the role of Fam3c in bone biology, we have generated a Fam3c knock-out (KO) mouse strain. The Fam3c KO mice were found to have normal appearance, behavior and fertility, but small changes in bone morphology and content were also observed. Micro-CT analysis of tibiae of the female mice revealed decreased number of trabeculae. In male mice the changes in the bone phenotype were smaller, but hematological changes were observed. Furthermore, there was a negative correlation between body weight and tibial trabecular and cortical bone volume in the male KO mice. There was a small increase in cortical bone mineral density, but in the lateral direction of tibiae the breaking strength was reduced. Fam3c KO bone marrow cells showed accelerated osteogenic differentiation and mineralization in vitro. The reduced number of bone trabeculae in Fam3c KO mice and the stimulated osteogenic differentiation indicate a role for Fam3c in osteoblast differentiation and bone homeostasis.

Genomic structure and marker-derived gene networks for growth and meat quality traits of Brazilian Nelore beef cattle

BACKGROUND

Nelore is the major beef cattle breed in Brazil with more than 130 million heads. Genome-wide association studies (GWAS) are often used to associate markers and genomic regions to growth and meat quality traits that can be used to assist selection programs. An alternative methodology to traditional GWAS that involves the construction of gene network interactions, derived from results of several GWAS is the AWM (Association Weight Matrices)/PCIT (Partial Correlation and Information Theory). With the aim of evaluating the genetic architecture of Brazilian Nelore cattle, we used high-density SNP genotyping data (~770,000 SNP) from 780 Nelore animals comprising 34 half-sibling families derived from highly disseminated and unrelated sires from across Brazil. The AWM/PCIT methodology was employed to evaluate the genes that participate in a series of eight phenotypes related to growth and meat quality obtained from this Nelore sample.

RESULTS

Our results indicate a lack of structuring between the individuals studied since principal component analyses were not able to differentiate families by its sires or by its ancestral lineages. The application of the AWM/PCIT methodology revealed a trio of transcription factors (comprising VDR, LHX9 and ZEB1) which in combination connected 66 genes through 359 edges and whose biological functions were inspected, some revealing to participate in biological growth processes in literature searches.

CONCLUSIONS

The diversity of the Nelore sample studied is not high enough to differentiate among families neither by sires nor by using the available ancestral lineage information. The gene networks constructed from the AWM/PCIT methodology were a useful alternative in characterizing genes and gene networks that were allegedly influential in growth and meat quality traits in Nelore cattle.

Applying personal genetic data to injury risk assessment in athletes

Recent studies have identified genetic markers associated with risk for certain sports-related injuries and performance-related conditions, with the hope that these markers could be used by individual athletes to personalize their training and diet regimens. We found that we could greatly expand the knowledge base of sports genetic information by using published data originally found in health and disease studies. For example, the results from large genome-wide association studies for low bone mineral density in elderly women can be re-purposed for low bone mineral density in young endurance athletes. In total, we found 124 single-nucleotide polymorphisms associated with: anterior cruciate ligament tear, Achilles tendon injury, low bone mineral density and stress fracture, osteoarthritis, vitamin/mineral deficiencies, and sickle cell trait. Of these single nucleotide polymorphisms, 91% have not previously been used in sports genetics. We conducted a pilot program on fourteen triathletes using this expanded knowledge base of genetic variants associated with sports injury. These athletes were genotyped and educated about how their individual genetic make-up affected their personal risk profile during an hour-long personal consultation. Overall, participants were favorable of the program, found it informative, and most acted upon their genetic results. This pilot program shows that recent genetic research provides valuable information to help reduce sports injuries and to optimize nutrition. There are many genetic studies for health and disease that can be mined to provide useful information to athletes about their individual risk for relevant injuries.

Is GSN significant for hip BMD in female Caucasians?

Low bone mineral density (BMD) is a risk factor for osteoporosis. Osteoporosis is more prevalent in females than in males. So far, the pathophysiological mechanisms underlying osteoporosis are unclear. Peripheral blood monocytes (PBMs) are precursors of bone-resorbing osteoclasts. This study aims to identify PBM-expressed proteins (genes) influencing hip BMD in humans. We utilized three independent study cohorts (N=34, 29, 40), including premenopausal Caucasians with discordant hip BMD. We studied PBM proteome-wide protein expression profiles in cohort 1 and identified 57 differentially expressed proteins (DEPs) between low vs. high BMD subjects. One protein gelsolin (GSN), after validation by Western blotting, was subject to follow-up. We compared GSN mRNA level in PBM between low vs. high BMD subjects in cohorts 2 and 3. We genotyped SNPs across GSN in 2286 unrelated Caucasians (cohort 4) and 1627 Chinese (cohort 5) and tested their association with hip BMD in females and males, respectively. We discovered and validated that GSN protein expression level in PBM was down-regulated 3.0-fold in low vs. high BMD subjects (P<0.05). Down-regulation of GSN in PBM in low BMD subjects was also observed at mRNA level in both cohort 2 and cohort 3. We identified that SNP rs767770 was significantly associated with hip BMD in female Caucasians (P=0.0003) only. Integrating analyses of the datasets at DNA, RNA, and protein levels from female Caucasians substantiated that GSN is highly significant for hip BMD (P=0.0001). We conclude that GSN is a significant gene influencing hip BMD in female Caucasians.

Proteomic analysis of human osteoarthritis synovial fluid

Background

Osteoarthritis is a chronic musculoskeletal disorder characterized mainly by progressive degradation of the hyaline cartilage. Patients with osteoarthritis often postpone seeking medical help, which results in the diagnosis being made at an advanced stage of cartilage destruction. Sustained efforts are needed to identify specific markers that might help in early diagnosis, monitoring disease progression and in improving therapeutic outcomes. We employed a multipronged proteomic approach, which included multiple fractionation strategies followed by high resolution mass spectrometry analysis to explore the proteome of synovial fluid obtained from osteoarthritis patients. In addition to the total proteome, we also enriched glycoproteins from synovial fluid using lectin affinity chromatography.

Results

We identified 677 proteins from synovial fluid of patients with osteoarthritis of which 545 proteins have not been previously reported. These novel proteins included ADAM-like decysin 1 (ADAMDEC1), alanyl (membrane) aminopeptidase (ANPEP), CD84, fibulin 1 (FBLN1), matrix remodelling associated 5 (MXRA5), secreted phosphoprotein 2 (SPP2) and spondin 2 (SPON2). We identified 300 proteins using lectin affinity chromatography, including the glycoproteins afamin (AFM), attractin (ATRN), fibrillin 1 (FBN1), transferrin (TF), tissue inhibitor of metalloproteinase 1 (TIMP1) and vasorin (VSN). Gene ontology analysis confirmed that a majority of the identified proteins were extracellular and are mostly involved in cell communication and signaling. We also confirmed the expression of ANPEP, dickkopf WNT signaling pathway inhibitor 3 (DKK3) and osteoglycin (OGN) by multiple reaction monitoring (MRM) analysis of osteoarthritis synovial fluid samples.

Conclusions

We present an in-depth analysis of the synovial fluid proteome from patients with osteoarthritis. We believe that the catalog of proteins generated in this study will further enhance our knowledge regarding the pathophysiology of osteoarthritis and should assist in identifying better biomarkers for early diagnosis.

Association between the SPRY1 gene polymorphism and obesity-related traits and osteoporosis in Korean women

BACKGROUND

Emerging evidence has revealed a close relationship between obesity and osteoporosis. It was reported recently that conditional knockout of the Spry1 gene in mice adipocytes causes an increase in body fat and a decrease in bone mass, and that these phenotypes are rescued by Spry1 overexpression in adipose tissue. In this study, we investigated whether genetic variation in the human SPRY1 gene is associated with obesity-related phenotypes and/or osteoporosis in humans.

METHODS

We performed a candidate gene association analysis between the four single nucleotide polymorphisms (SNPs) and 14 imputed SNPs in the SPRY1 gene and obesity-related traits and osteoporosis in a Korean women cohort (3013 subjects).

RESULTS

All four SPRY1 gene SNPs were significantly associated with either obesity-related traits or osteoporosis. The TGCC haplotype in the SRPY1 gene showed simultaneous association with an increased risk for obesity-related traits, percentage body fat (p=0.0087) and percentage abdominal fat (p=0.047), and osteoporosis (odds ratio=1.50; p=0.025) in the recessive genetic model.

CONCLUSIONS

Our results support a previous finding in conditional Spry1 gene knockout mice and suggest that the SPRY1 gene is an important genetic factor for determining the risk of both obesity and osteoporosis in humans.

Pourcain B, Xiao Y, Timpson NJ, Smith GD, Reid IR, Ring SM, Sambrook PN, Karlsson M, Dennison EM, Kemp JP, Danoy P, Sayers A, Wilson SG, Nethander M, McCloskey E, Vandenput L, Eastell R, Liu J, Spector T, Mitchell BD, Streeten EA, Brommage R, Pettersson-Kymmer U, Brown MA, Ohlsson C, Richards JB, Lorentzon M. WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)<P<5.9 × 10(-4)) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture.

Pourcain B, Hofman A, Jaddoe VWV, Smith GD, Lorentzon M, Gautvik KM, Uitterlinden AG, Brommage R, Ohlsson C, Tobias JH, Rivadeneira F. Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1×10⁻¹¹ observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ±500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6×10⁻³¹ and an effect size explaining between 0.6%–1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42×10⁻¹⁰) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9×10⁻¹⁶) and rs7801723 (P = 8.9×10⁻²⁸), also mapping to C7orf58 (r² = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

Genetic investigations on bone mineral density (BMD) variation in children allow the identification of factors determining peak bone mass and their influence on developing osteoporosis later in life. We ran a genome-wide association study (GWAS) for total body BMD based on 2,660 children of different ethnic backgrounds, followed by replication in an additional 12,066 individuals comprising children, young adults, and elderly populations. Our GWAS meta-analysis identified two independent signals in the 7q31.31 locus, arising from SNPs in the vicinity of WNT16, FAM3C, and C7orf58. These variants were also associated with skull BMD, a skeletal trait with much less environmental influence for which one of the signals displayed age-specific effects. Integration of functional studies in a Wnt16 knockout mouse model and gene expression profiles in human bone tissue provided additional evidence that WNT16 and C7orf58 underlie the described associations. All together our findings demonstrate the relevance of these factors for bone biology, the attainment of peak bone mass, and their likely impact on bone fragility later in life.