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 20.11 standard deviations [SD] per C allele, P = 6.2610 29 ). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly.Arg), also had genome-wide significant association with forearm BMD (20.14 SD per C allele, P = 2.3610 212 , and 20.16 SD per G allele, P = 1.2610 215 , 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.3610 29 ), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9610 26 and rs2707466: OR = 1.22, P = 7.2610 26 ). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16 2/2 mice had 27% (P,0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%–61% (6.5610 213 ,P,5.9610 24 ) 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.

/pdf/wnt16-influences-bone-mineral-density-cortical-bone-494zq7zu5p.pdf

WNT16 Influences Bone Mineral Density, Cortical Bone Thickness, Bone Strength, and Osteoporotic Fracture Risk

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-11 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-31 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-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-16) and rs7801723 (P = 8.9×10-28), also mapping to C7orf58 (r2 = 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.

/pdf/meta-analysis-of-genome-wide-scans-for-total-body-bmd-in-x12arfypf1.pdf

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

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. 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. 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.

/pdf/proteomic-analysis-of-human-osteoarthritis-synovial-fluid-3xi498rgrb.pdf

Proteomic analysis of human osteoarthritis synovial fluid

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. 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. 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.

/pdf/genomic-structure-and-marker-derived-gene-networks-for-1t1v1ukz11.pdf

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

Fam3c modulates osteogenic differentiation by down-regulating Runx2

Summary
We tested whether two genetic variants were associated with BMD at multiple clinically relevant skeletal sites in Caucasians. We found that variant rs7776725 is consistently associated with hip, spine, wrist and whole-body BMD, which highlights the potential importance of this variant or linked variants for osteoporosis.

Y.-P. Zheng

Papers

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