A variant in MCF2L is associated with osteoarthritis.

Osteoarthritis (OA) is a prevalent, heritable degenerative joint disease with a substantial public health impact. We used a 1000-Genomes-Project-based imputation in a genome-wide association scan for osteoarthritis (3177 OA cases and 4894 controls) to detect a previously unidentified risk locus. We discovered a small disease-associated set of variants on chromosome 13. Through large-scale replication, we establish a robust association with SNPs in MCF2L (rs11842874, combined odds ratio [95% confidence interval] 1.17 [1.11-1.23], p = 2.1 × 10(-8)) across a total of 19,041 OA cases and 24,504 controls of European descent. This risk locus represents the third established signal for OA overall. MCF2L regulates a nerve growth factor (NGF), and treatment with a humanized monoclonal antibody against NGF is associated with reduction in pain and improvement in function for knee OA patients.

A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip.

OBJECTIVES: Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS: We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS: We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS: Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

No evidence of an association between mitochondrial DNA variants and osteoarthritis in 7393 cases and 5122 controls.

OBJECTIVES: Osteoarthritis (OA) has a complex aetiology with a strong genetic component. Genome-wide association studies implicate several nuclear genes in the aetiology, but a major component of the heritability has yet to be defined at the molecular level. Initial studies implicate maternally inherited variants of mitochondrial DNA (mtDNA) in subgroups of patients with OA based on gender and specific joint involvement, but these findings have not been replicated. METHODS: The authors studied 138 maternally inherited mtDNA variants genotyped in a two cohort genetic association study across a total of 7393 OA cases from the arcOGEN consortium and 5122 controls genotyped in the Wellcome Trust Case Control consortium 2 study. RESULTS: Following data quality control we examined 48 mtDNA variants that were common in cohort 1 and cohort 2, and found no association with OA. None of the phenotypic subgroups previously associated with mtDNA haplogroups were associated in this study. CONCLUSIONS: We were not able to replicate previously published findings in the largest mtDNA association study to date. The evidence linking OA to mtDNA is not compelling at present.

Evaluation of the genetic overlap between osteoarthritis with body mass index and height using genome-wide association scan data.

OBJECTIVES: Obesity as measured by body mass index (BMI) is one of the major risk factors for osteoarthritis. In addition, genetic overlap has been reported between osteoarthritis and normal adult height variation. We investigated whether this relationship is due to a shared genetic aetiology on a genome-wide scale. METHODS: We compared genetic association summary statistics (effect size, p value) for BMI and height from the GIANT consortium genome-wide association study (GWAS) with genetic association summary statistics from the arcOGEN consortium osteoarthritis GWAS. Significance was evaluated by permutation. Replication of osteoarthritis association of the highlighted signals was investigated in an independent dataset. Phenotypic information of height and BMI was accounted for in a separate analysis using osteoarthritis-free controls. RESULTS: We found significant overlap between osteoarthritis and height (p=3.3×10(-5) for signals with p≤0.05) when the GIANT and arcOGEN GWAS were compared. For signals with p≤0.001 we found 17 shared signals between osteoarthritis and height and four between osteoarthritis and BMI. However, only one of the height or BMI signals that had shown evidence of association with osteoarthritis in the arcOGEN GWAS was also associated with osteoarthritis in the independent dataset: rs12149832, within the FTO gene (combined p=2.3×10(-5)). As expected, this signal was attenuated when we adjusted for BMI. CONCLUSIONS: We found a significant excess of shared signals between both osteoarthritis and height and osteoarthritis and BMI, suggestive of a common genetic aetiology. However, only one signal showed association with osteoarthritis when followed up in a new dataset.

Structure of ubiquitin-fold modifier 1-specific protease UfSP2.

Ubiquitin-fold modifier 1 (Ufm1)-specific protease 2 (UfSP2) is a cysteine protease that is responsible for the release of Ufm1 from Ufm1-conjugated cellular proteins, as well as for the generation of mature Ufm1 from its precursor. The 2.6 Å resolution crystal structure of mouse UfSP2 reveals that it is composed of two domains. The C-terminal catalytic domain is similar to UfSP1 with Cys(294), Asp(418), His(420), Tyr(282), and a regulatory loop participating in catalysis. The novel N-terminal domain shows a unique structure and plays a role in the recognition of its cellular substrate C20orf116 and thus in the recruitment of UfSP2 to the endoplasmic reticulum, where C20orf116 predominantly localizes. Mutagenesis studies were carried out to provide the structural basis for understanding the loss of catalytic activity observed in a recently identified UfSP2 mutation that is associated with an autosomal dominant form of hip dysplasia.

Single-Cell Transcriptomic Analyses Define Distinct Peripheral B Cell Subsets and Discrete Development Pathways.

Separation of B cells into different subsets has been useful to understand their different functions in various immune scenarios. In some instances, the subsets defined by phenotypic FACS separation are relatively homogeneous and so establishing the functions associated with them is straightforward. Other subsets, such as the "Double negative" (DN, CD19+CD27-IgD-) population, are more complex with reports of differing functionality which could indicate a heterogeneous population. Recent advances in single-cell techniques enable an alternative route to characterize cells based on their transcriptome. To maximize immunological insight, we need to match prior data from phenotype-based studies with the finer granularity of the single-cell transcriptomic signatures. We also need to be able to define meaningful B cell subsets from single cell analyses performed on PBMCs, where the relative paucity of a B cell signature means that defining B cell subsets within the whole is challenging. Here we provide a reference single-cell dataset based on phenotypically sorted B cells and an unbiased procedure to better classify functional B cell subsets in the peripheral blood, particularly useful in establishing a baseline cellular landscape and in extracting significant changes with respect to this baseline from single-cell datasets. We find 10 different clusters of B cells and applied a novel, geometry-inspired, method to RNA velocity estimates in order to evaluate the dynamic transitions between B cell clusters. This indicated the presence of two main developmental branches of memory B cells. A T-independent branch that involves IgM memory cells and two DN subpopulations, culminating in a population thought to be associated with Age related B cells and the extrafollicular response. The other, T-dependent, branch involves a third DN cluster which appears to be a precursor of classical memory cells. In addition, we identify a novel DN4 population, which is IgE rich and closely linked to the classical/precursor memory branch suggesting an IgE specific T-dependent cell population.

Discovery and Toxicological Profiling of Aminopyridines as Orally Bioavailable Selective Inhibitors of PI3-Kinase γ.

Using a novel physiologically relevant in vitro human whole blood neutrophil shape change assay, an aminopyrazine series of selective PI3Kγ inhibitors was identified and prioritized for further optimization. Severe solubility limitations associated with the series leading to low oral bioavailability and poor exposures, especially at higher doses, were overcome by moving to an aminopyridine core. Compound 33, with the optimal balance of on-target activity, selectivity, and pharmacokinetic parameters, progressed into in vivo studies and demonstrated good efficacy (10 mg/kg) in a rat model of airway inflammation. Sufficient exposures were achieved at high doses to support toxicological studies, where unexpected inflammatory cell infiltrates in cardiovascular tissue prevented further compound development.

A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5' UTR of GDF5 with osteoarthritis susceptibility.

We have performed a meta-analysis combining data for more than 11,000 individuals. It provides compelling evidence for a positive association between a functional single-nucleotide polymorphism (SNP) in the 5'-UTR of GDF5 (+104T/C; rs143383) and osteoarthritis (OA) in European and Asian populations. This SNP has recently been reported to be associated with OA in Japanese and Han Chinese populations. Attempts to replicate this association in European samples have been inconclusive, as no association was found in the case-control cohorts from the UK, Spain and Greece when studied individually. However, the pooled data of UK and Spain found an association of the T-allele with an odds ratio (OR) of 1.10. Although the European studies had adequate power to replicate the original findings from the Japanese cohort (OR = 1.79), these results suggest that the role of the GDF5 polymorphism may not be as strong in Europeans. To clarify whether the European studies were hampered by insufficient power, we combined new data from the UK and the Netherlands with the three published studies of Europe and Asia. The results provide strong evidence of a positive association of the GDF5 SNP with knee OA for Europeans as well as for Asians. The combined association for both ethnic groups is highly significant for the allele frequency model (P = 0.0004, OR = 1.21) and the dominant model (P < 0.0001, OR = 1.48). These findings represent the first highly significant evidence for a risk factor for the development of OA which affects two highly diverse ethnic groups.

Superoxide dismutase downregulation in osteoarthritis progression and end-stage disease.

BACKGROUND: Oxidative stress is proposed as an important factor in osteoarthritis (OA). OBJECTIVE: To investigate the expression of the three superoxide dismutase (SOD) antioxidant enzymes in OA. METHODS: SOD expression was determined by real-time PCR and immunohistochemistry using human femoral head cartilage. SOD2 expression in Dunkin-Hartley guinea pig knee articular cartilage was determined by immunohistochemistry. The DNA methylation status of the SOD2 promoter was determined using bisulphite sequencing. RNA interference was used to determine the consequence of SOD2 depletion on the levels of reactive oxygen species (ROS) using MitoSOX and collagenases, matrix metalloproteinase 1 (MMP-1) and MMP-13, gene expression. RESULTS: All three SOD were abundantly expressed in human cartilage but were markedly downregulated in end-stage OA cartilage, especially SOD2. In the Dunkin-Hartley guinea pig spontaneous OA model, SOD2 expression was decreased in the medial tibial condyle cartilage before, and after, the development of OA-like lesions. The SOD2 promoter had significant DNA methylation alterations in OA cartilage. Depletion of SOD2 in chondrocytes increased ROS but decreased collagenase expression. CONCLUSION: This is the first comprehensive expression profile of all SOD genes in cartilage and, importantly, using an animal model, it has been shown that a reduction in SOD2 is associated with the earliest stages of OA. A decrease in SOD2 was found to be associated with an increase in ROS but a reduction of collagenase gene expression, demonstrating the complexities of ROS function.

Oxidative stress modulates theophylline effects on steroid responsiveness.

Oxidative stress is a central factor in many chronic inflammatory diseases such as severe asthma and chronic obstructive pulmonary disease (COPD). Oxidative stress reduces the anti-inflammatory corticosteroid action and may therefore contribute to the relative corticosteroid insensitivity seen in these diseases. Low concentrations of theophylline can restore the anti-inflammatory action of corticosteroids in oxidant exposed cells, however the mechanism remains unknown. Here, we demonstrate that a low concentration of theophylline restores corticosteroid repression of pro-inflammatory mediator release and histone acetylation in oxidant exposed cells. Global gene expression analysis shows that theophylline regulates distinct pathways in naïve and oxidant exposed cells and reverses oxidant mediated modulated of pathways. Furthermore, quantitative chemoproteomics revealed that theophylline has few high affinity targets in naive cells but an elevated affinity in oxidant stressed cells. In conclusion, oxidative stress alters theophylline binding profile and gene expression which may result in restoration of corticosteroid function.

Col2a1 lineage tracing reveals that the meniscus of the knee joint has a complex cellular origin.

The knee joint consists of multiple interacting tissues that are prone to injury- and disease-related degeneration. Although much is known about the structure and function of the knee's constituent tissues, relatively little is known about their cellular origin and the mechanisms governing their segregation. To investigate the origin and segregation of knee tissues in vivo we performed lineage tracing using a Col2a1-Cre/R26R mouse model system and compared the data obtained with actual Col2a1 expression. These studies demonstrated that at E13.5 the interzone at the presumptive joint site forms when cells within the Col2a1-expressing anlagen cease expression of Col2a1 and not through cellular invasion into the anlagen. Later in development these interzone cells form the cruciate ligament and inner medial meniscus of the knee. At E14.5, after interzone formation, cells that had never expressed Col2a1 appeared in the joint and formed the lateral meniscus. Furthermore, cells with a Col2a1-positive expression history combined with the negative cells to form the medial meniscus. The invading cells started to express Col2a1 1 week after birth, resulting in all cells within the meniscus synthesizing collagen II. These findings support a model of knee development in which cells present in the original anlagen combine with invading cells in the formation of this complex joint.

The function and interrelationship between GDF5 and ERG-010 during chondrogenesis in vitro.

Joint formation begins with the establishment of an interzone within the cartilaginous anlagen of the future skeleton. Both GDF5 and ERG are proposed as regulators of chondrocyte differentiation during and post interzone formation. The aim of this study was to examine the relationship between Gdf5 and Erg expression and downstream effects on chondrocyte gene expression. Erg expression was identified in mouse knee joints at E13.5. Expression analyses were performed using micromass cultures of murine C3H10T1/2 mesenchymal cells undergoing induced chondrogenesis in the presence and absence of GDF5 and ERG. At E13.5, Erg expression was found to surround epiphyseal chondrocytes and span the interzone up to the intermediate zone. Erg splice forms were expressed in micromass cultures, and their expression profile was altered by the addition of recombinant GDF5 depending on the stage of differentiation. Overexpression of Erg-010 resulted in a downregulation of Col2a1 and Col10a1. Microarray analysis following Erg-010 overexpression identified two potential downstream targets, Ube2b and Osr2, which were also differentially regulated by GDF5. Erg regulation by GDF5 in induced mesenchymal cells in vitro is dependent on the stage of chondrogenesis, and its expression in vivo demarcates chondrocytes that are not destined to be consumed by endochondral ossification. Functionally, Erg expression causes downregulation of Col2a1 and Col10a1 expression and this effect is potentially mediated by Osr2 and/or Ube2b. Combined, these data suggest a possible pathway linking GDF5, ERG and downstream factors in the processes of chondrocyte differentiation during articular joint formation.

Identification of a mutation in the ubiquitin-fold modifier 1-specific peptidase 2 gene, UFSP2, in an extended South African family with Beukes hip dysplasia.

BACKGROUND: Beukes hip dysplasia (BHD) is an autosomal dominant disorder of variable penetrance that was originally identified in a large South African family of European origin. BHD is characterised by bilateral dysmorphism of the proximal femur, which results in severe degenerative osteoarthropathy. Previous studies mapped the disorder to a 3.34 Mb region on chromosome 4q35. OBJECTIVE: To fine-map the BHD locus and identify the disease-causing mutation by direct sequencing. RESULTS: The linked BHD allele was refined to 1.33 Mb, reducing the number of candidate genes from 25 to 16. Analysis of protein coding and invariant splice-site sequences in three distantly related individuals identified a single-candidate disease-causing variant c.868T>C within exon 8 of the ubiquitin-fold modifier 1 (Ufm1)-specific peptidase 2 gene, UFSP2. The presence of this unique mutation was confirmed in all 17 affected members of the BHD family who were genotyped. The mutation segregated with the BHD phenotype in the extended family with a two-point (single marker) LOD score of 10.4 (θ=0.0 and 80% penetrance). The mutation predicts the substitution of a highly conserved amino acid, p.Tyr290His, in the encoded protein. In vitro functional assays performed using purified recombinant wild-type and mutant UFSP2 protein demonstrated that the BHD mutation abolishes UFSP2-mediated C-terminal cleavage of its substrate, Ufm1. CONCLUSION: We report a unique UFSP2 mutation that segregates with the BHD phenotype. The predicted amino acid substitution inactivates UFSP2 proteolytic function, thus implicating the ubiquitin-fold modifier 1 cascade in this form of severe hip osteoarthropathy. The facile polymerase chain reaction-based assay we describe could be used to confirm the diagnosis of BHD, or for presymptomatic testing of members of the extended BHD family.

Investigation of association between hip osteoarthritis susceptibility loci and radiographic proximal femur shape.

OBJECTIVE: To test whether previously reported hip morphology or osteoarthritis (OA) susceptibility loci are associated with proximal femur shape as represented by statistical shape model (SSM) modes and as univariate or multivariate quantitative traits. METHODS: We used pelvic radiographs and genotype data from 929 subjects with unilateral hip OA who had been recruited previously for the Arthritis Research UK Osteoarthritis Genetics Consortium genome-wide association study. We built 3 SSMs capturing the shape variation of the OA-unaffected proximal femur in the entire mixed-sex cohort and for male/female-stratified cohorts. We selected 41 candidate single-nucleotide polymorphisms (SNPs) previously reported as being associated with hip morphology (for replication analysis) or OA (for discovery analysis) and for which genotype data were available. We performed 2 types of analysis for genotype-phenotype associations between these SNPs and the modes of the SSMs: 1) a univariate analysis using individual SSM modes and 2) a multivariate analysis using combinations of SSM modes. RESULTS: The univariate analysis identified association between rs4836732 (within the ASTN2 gene) and mode 5 of the female SSM (P = 0.0016) and between rs6976 (within the GLT8D1 gene) and mode 7 of the mixed-sex SSM (P = 0.0003). The multivariate analysis identified association between rs5009270 (near the IFRD1 gene) and a combination of modes 3, 4, and 9 of the mixed-sex SSM (P = 0.0004). Evidence of associations remained significant following adjustment for multiple testing. All 3 SNPs had previously been associated with hip OA. CONCLUSION: These de novo findings suggest that rs4836732, rs6976, and rs5009270 may contribute to hip OA susceptibility by altering proximal femur shape.

Annexin VIII is differentially expressed by chondrocytes in the mammalian growth plate during endochondral ossification and in osteoarthritic cartilage.

Endochondral ossification is the developmental process that leads to the formation and coordinated longitudinal growth of the majority of the vertebrate skeleton. Central to this process is chondrocyte differentiation occurring in the growth plate that lies at the junction between the epiphyseal cartilage and the bone. To identify novel factors involved in this differentiation process, suppression subtractive hybridization was performed to amplify preferentially cDNAs uniquely expressed in fetal bovine growth plate chondrocytes as opposed to epiphyseal chondrocytes. The subtracted product was used to screen a fetal bovine chondrocyte cDNA library. One of the cDNA clones identified encoded the bovine orthologue of annexin VIII, a protein not previously described in the growth plate. Northern and Western blotting confirmed that annexin VIII was expressed by growth plate chondrocytes and not by epiphyseal chondrocytes. Immunohistochemistry of the fetal bovine growth plate identified a gradient of increasing annexin VIII protein from the proliferative to the hypertrophic zone. Immunofluorescence localized annexin VIII largely to the chondrocyte cell membrane. In a preliminary study, we examined the distribution of annexin VIII in normal and osteoarthritic (OA) articular cartilage. In OA cartilage, the protein was located in a subset of mid- to deep zone chondrocytes and in the matrix surrounding these cells; no annexin VIII was detected in normal articular cartilage. Thus annexin VIII is a marker for chondrocyte differentiation during normal endochondral ossification and may act as a marker for cells undergoing inappropriate differentiation in OA.

An enhancer complex confers both high-level and cell-specific expression of the human type X collagen gene.

Type X collagen expression is restricted to hypertrophic chondrocytes in the endochondral growth plate. Transient transfection of reporter constructs containing the human collagen X promoter into primary growth plate chondrocytes identified a cis-acting positive regulatory DNA element(s) that has cell-specific enhancer properties and binds a nuclear protein expressed specifically in growth plate chondrocytes. Functional disruption of this region results in a significant reduction in the activation of reporter gene transcription. The identified enhancer is a major element controlling both high-level and cell-specific expression of type X collagen gene.

Assessment of osteoarthritis candidate genes in a meta-analysis of nine genome-wide association studies.

OBJECTIVE: To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA. METHODS: A total of 199 candidate genes for association with OA were identified using Human Genome Epidemiology (HuGE) Navigator. All of their single-nucleotide polymorphisms (SNPs) with an allele frequency of >5% were assessed by fixed-effects meta-analysis of 9 genome-wide association studies (GWAS) that included 5,636 patients with knee OA and 16,972 control subjects and 4,349 patients with hip OA and 17,836 control subjects of European ancestry. An additional 5,921 individuals were genotyped for significantly associated SNPs in the meta-analysis. After correction for the number of independent tests, P values less than 1.58 × 10(-5) were considered significant. RESULTS: SNPs at only 2 of the 199 candidate genes (COL11A1 and VEGF) were associated with OA in the meta-analysis. Two SNPs in COL11A1 showed association with hip OA in the combined analysis: rs4907986 (P = 1.29 × 10(-5) , odds ratio [OR] 1.12, 95% confidence interval [95% CI] 1.06-1.17) and rs1241164 (P = 1.47 × 10(-5) , OR 0.82, 95% CI 0.74-0.89). The sex-stratified analysis also showed association of COL11A1 SNP rs4908291 in women (P = 1.29 × 10(-5) , OR 0.87, 95% CI 0.82-0.92); this SNP showed linkage disequilibrium with rs4907986. A single SNP of VEGF, rs833058, showed association with hip OA in men (P = 1.35 × 10(-5) , OR 0.85, 95% CI 0.79-0.91). After additional samples were genotyped, association at one of the COL11A1 signals was reinforced, whereas association at VEGF was slightly weakened. CONCLUSION: Two candidate genes, COL11A1 and VEGF, were significantly associated with OA in this focused meta-analysis. The remaining candidate genes were not associated.

Accurate bone segmentation in 2D radiographs using fully automatic shape model matching based on regression-voting.

Recent work has shown that using Random Forests (RFs) to vote for the optimal position of model feature points leads to robust and accurate shape model matching. This paper applies RF regression-voting as part of a fully automatic shape model matching (FASMM) system to three different radiograph segmentation problems: the proximal femur, the bones of the knee joint and the joints of the hand. We investigate why this approach works so well and demonstrate that the performance comes from a combination of three properties: (i) The integration of votes from multiple regions around the model point. (ii) The combination of multiple independent votes from each tree. (iii) The use of a coarse to fine strategy. We show that each property can improve performance, and that the best performance comes from using all three. We demonstrate that FASMM based on RF regression-voting generalises well across application areas, achieving state of the art performance in each of the three segmentation problems. This FASMM system provides an accurate and time-efficient way for the segmentation of bony structures in radiographs.

The effect of genome-wide association scan quality control on imputation outcome for common variants.

Imputation is an extremely valuable tool in conducting and synthesising genome-wide association studies (GWASs). Directly typed SNP quality control (QC) is thought to affect imputation quality. It is, therefore, common practise to use quality-controlled (QCed) data as an input for imputing genotypes. This study aims to determine the effect of commonly applied QC steps on imputation outcomes. We performed several iterations of imputing SNPs across chromosome 22 in a dataset consisting of 3177 samples with Illumina 610 k (Illumina, San Diego, CA, USA) GWAS data, applying different QC steps each time. The imputed genotypes were compared with the directly typed genotypes. In addition, we investigated the correlation between alternatively QCed data. We also applied a series of post-imputation QC steps balancing elimination of poorly imputed SNPs and information loss. We found that the difference between the unQCed data and the fully QCed data on imputation outcome was minimal. Our study shows that imputation of common variants is generally very accurate and robust to GWAS QC, which is not a major factor affecting imputation outcome. A minority of common-frequency SNPs with particular properties cannot be accurately imputed regardless of QC stringency. These findings may not generalise to the imputation of low frequency and rare variants.

Large-scale analysis of association between GDF5 and FRZB variants and osteoarthritis of the hip, knee, and hand.

OBJECTIVE: GDF5 and FRZB have been proposed as genetic loci conferring susceptibility to osteoarthritis (OA); however, the results of several studies investigating the association of OA with the rs143383 polymorphism of the GDF5 gene or the rs7775 and rs288326 polymorphisms of the FRZB gene have been conflicting or inconclusive. To examine these associations, we performed a large-scale meta-analysis of individual-level data. METHODS: Fourteen teams contributed data on polymorphisms and knee, hip, and hand OA. For rs143383, the total number of cases and controls, respectively, was 5,789 and 7,850 for hip OA, 5,085 and 8,135 for knee OA, and 4,040 and 4,792 for hand OA. For rs7775, the respective sample sizes were 4,352 and 10,843 for hip OA, 3,545 and 6,085 for knee OA, and 4,010 and 5,151 for hand OA, and for rs288326, they were 4,346 and 8,034 for hip OA, 3,595 and 6,106 for knee OA, and 3,982 and 5,152 for hand OA. For each individual study, sex-specific odds ratios (ORs) were calculated for each OA phenotype that had been investigated. The ORs for each phenotype were synthesized using both fixed-effects and random-effects models for allele-based effects, and also for haplotype effects for FRZB. RESULTS: A significant random-effects summary OR for knee OA was demonstrated for rs143383 (1.15 [95% confidence interval 1.09-1.22]) (P=9.4x10(-7)), with no significant between-study heterogeneity. Estimates of effect sizes for hip and hand OA were similar, but a large between-study heterogeneity was observed, and statistical significance was borderline (for OA of the hip [P=0.016]) or absent (for OA of the hand [P=0.19]). Analyses for FRZB polymorphisms and haplotypes did not reveal any statistically significant signals, except for a borderline association of rs288326 with hip OA (P=0.019). CONCLUSION: Evidence of an association between the GDF5 rs143383 polymorphism and OA is substantially strong, but the genetic effects are consistent across different populations only for knee OA. Findings of this collaborative analysis do not support the notion that FRZB rs7775 or rs288326 has any sizable genetic effect on OA phenotypes.