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

A systematic review and cross-database analysis of single nucleotide polymorphisms underlying hip morphology and osteoarthritis reveals shared mechanisms

OBJECTIVE

Understanding the mechanisms of hip disease, such as osteoarthritis (OA), is crucial to advance their treatment. Such hip diseases often involve specific morphological changes. Genetic variations, called single nucleotide polymorphisms (SNPs), influence various hip morphological parameters. This study investigated the biological relevance of SNPs correlated to hip morphology in genome-wide association studies (GWAS). The SNP-associated genes were compared to genes associated with OA in other joints, aiming to see if the same genes play a role in both hip development and the risk of OA in other lower limb joints.

METHODOLOGY

A systematic literature review was conducted to identify SNPs correlated with hip morphology, based on the Population, Intervention, Comparison, Outcome, and Study (PICOS) framework. Afterwards, Gene Ontology (GO) analysis was performed, using EnrichR, on the SNP-associated genes and compared with non-hip OA-associated genes, across different databases.

RESULTS

Reviewing 49 GWAS identified 436 SNPs associated with hip joint morphology, encompassing variance in bone size, structure and shape. Among the SNP-associated genes, SOX9 plays a pivotal role in size, GDF5 impacts bone structure, and BMP7 affects shape. Overall, skeletal system development, regulation of cell differentiation, and chondrocyte differentiation emerged as crucial processes influencing hip morphology. Eighteen percent of GWAS-identified genes related to hip morphology were also associated with non-hip OA.

CONCLUSION

Our findings indicate the existence of multiple shared genetic mechanisms across hip morphology and OA, highlighting the necessity for more extensive research in this area, as in contrast to the hip, the genetic background on knee or foot morphology remains largely understudied.

Bone-cartilage crosstalk informed by aging mouse bone transcriptomics and human osteoarthritis genome-wide association studies

Subchondral bone participates in crosstalk with articular cartilage to maintain joint homeostasis, and disruption of either tissue results in overall joint degeneration. Among the subchondral bone changes observed in osteoarthritis (OA), subchondral bone plate (SBP) thickening has a time-dependent relationship with cartilage degeneration and has recently been shown to be regulated by osteocytes. Here, we evaluate the effect of age on SBP thickness and cartilage degeneration in aging mice. We find that SBP thickness significantly increases by 18-months of age, corresponding temporally with increased cartilage degeneration. To identify factors in subchondral bone that may participate in bone cartilage crosstalk or OA, we leveraged mouse transcriptomic data from one joint tissue compartment - osteocyte-enriched bone - to search for enrichment with human OA in UK Biobank and Arthritis Research UK Osteoarthritis Genetics (arcOGEN) GWAS using the mouse2human (M2H, www.mouse2human.org) strategy. Genes differentially expressed in aging mouse bone are significantly enriched for human OA, showing joint site-specific (knee vs. hip) relationships, exhibit temporal associations with age, and unique gene clusters are implicated in each type of OA. Application of M2H identifies genes with known and unknown functions in osteocytes and OA development that are clinically associated with human OA. Altogether, this work prioritizes genes with a potential role in bone/cartilage crosstalk for further mechanistic study based on their association with human OA in GWAS.

Intergenic Interactions of SBNO1, NFAT5 and GLT8D1 Determine the Susceptibility to Knee Osteoarthritis among Europeans of Russia

This study was conducted to examine the associations between genome-wide association studies (GWAS)-important single nucleotide polymorphisms (SNPs) and knee osteoarthritis (KOA) among Europeans of Russia. The present replicative study ("patient-control" design has been used) was carried out on 1000 DNA samples from KOA (n = 500) and KOA-free (n = 500) participants. Ten GWAS-important for KOA SNPs of eight candidate genes (LYPLAL1, GNL3, GLT8D1, SBNO1, WWP2, NFAT5, TGFA, GDF5) were studied. To assess the link between SNPs and KOA susceptibility, logistic regression (to establish independent SNP effects) and MB-MDR (to identify SNP-SNP interactions) were used. As a result of this genetic analysis, the associations of individual SNPs with KOA have not been proven. Eight loci out of ten tested SNPs interacted with each other (within twelve genetic models) and determined susceptibility to KOA. The greatest contribution to the disease development were made by three polymorphisms/genes such as rs6976 (C>T) GLT8D1, rs56116847 (G>A) SBNO1, rs6499244 (T>A) NFAT5 (each was included in 2/3 [8 out 12] KOA-responsible genetic interaction models). A two-locus epistatic interaction of rs56116847 (G >A) SBNO1 × rs6499244 (T>A) NFAT5 determined the maximum percentage (0.86%) of KOA entropy. KOA-associated SNPs are regulatory polymorphisms that affect the expression/splicing level, epigenetic modification of 72 genes in KOA-pathogenetically significant organs such as skeletal muscles, tibial arteries/nerves, thyroid, adipose tissue, etc. These putative KOA-effector genes are mainly involved in the organization/activity of the exoribonuclease complex and antigen processing/presentation pathways. In conclusion, KOA susceptibility among Europeans of Russia is mediated by intergenic interactions (but not the main effects) of GWAS-important SNPs.

Genetics of osteoarthritis

Osteoarthritis genetics has been transformed in the past decade through the application of large-scale genome-wide association scans. So far, over 100 polymorphic DNA variants have been associated with this common and complex disease. These genetic risk variants account for over 20% of osteoarthritis heritability and the vast majority map to non-protein coding regions of the genome where they are presumed to act by regulating the expression of target genes. Statistical fine mapping, in silico analyses of genomics data, and laboratory-based functional studies have enabled the identification of some of these targets, which encode proteins with diverse roles, including extracellular signaling molecules, intracellular enzymes, transcription factors, and cytoskeletal proteins. A large number of the risk variants correlate with epigenetic factors, in particular cartilage DNA methylation changes in cis, implying that epigenetics may be a conduit through which genetic effects on gene expression are mediated. Some of the variants also appear to have been selected as humans adapted to bipedalism, suggesting that a proportion of osteoarthritis genetic susceptibility results from antagonistic pleiotropy, with risk variants having a positive role in joint formation but a negative role in the long-term health of the joint. Although data from an osteoarthritis genetic study has not yet directly led to a novel treatment, some of the osteoarthritis associated genes code for proteins that have available therapeutics. Genetic investigations are therefore revealing fascinating fundamental insights into osteoarthritis and can expose options for translational intervention.

Epigenetic Regulation of Chondrocytes and Subchondral Bone in Osteoarthritis

The aim of this review is to provide an updated review of the epigenetic factors involved in the onset and development of osteoarthritis (OA). OA is a prevalent degenerative joint disease characterized by chronic inflammation, ectopic bone formation within the joint, and physical and proteolytic cartilage degradation which result in chronic pain and loss of mobility. At present, no disease-modifying therapeutics exist for the prevention or treatment of the disease. Research has identified several OA risk factors including mechanical stressors, physical activity, obesity, traumatic joint injury, genetic predisposition, and age. Recently, there has been increased interest in identifying epigenetic factors involved in the pathogenesis of OA. In this review, we detail several of these epigenetic modifications with known functions in the onset and progression of the disease. We also review current therapeutics targeting aberrant epigenetic regulation as potential options for preventive or therapeutic treatment.

Comparison of the Effect of MFAT and MFAT + PRP on Treatment of Hip Osteoarthritis: An Observational, Intention-to-Treat Study at One Year

Hip osteoarthritis (OA) is a major contributor to reduced quality of life and concomitant disability associated with lost working life months. Intra-articular injection of various biological materials has shown promise in alleviating symptoms and potentially slowing down the degenerative process. Here, we compared the effects of treatment of a cohort of 147 patients suffering from grade 1–4 hip OA; with either micro-fragmented adipose tissue (MFAT), or a combination of MFAT with platelet-rich plasma (PRP). We found significant improvements in both the visual analogue score for pain (VAS) and Oxford hip score (OHS) that were similar for both treatments with over 60% having an improvement in the VAS score of 20 points or more. These results suggest a positive role for intra-articular injection of MFAT + PRP as a treatment for hip osteoarthritis which may be important particularly in low body mass index (BMI) patients where the difficulty in obtaining sufficient MFAT for treatment could be offset by using this combination of biologicals.

The Genetic Epidemiology of Joint Shape and the Development of Osteoarthritis

Congruent, low-friction relative movement between the articulating elements of a synovial joint is an essential pre-requisite for sustained, efficient, function. Where disorders of joint formation or maintenance exist, mechanical overloading and osteoarthritis (OA) follow. The heritable component of OA accounts for ~ 50% of susceptible risk. Although almost 100 genetic risk loci for OA have now been identified, and the epidemiological relationship between joint development, joint shape and osteoarthritis is well established, we still have only a limited understanding of the contribution that genetic variation makes to joint shape and how this modulates OA risk. In this article, a brief overview of synovial joint development and its genetic regulation is followed by a review of current knowledge on the genetic epidemiology of established joint shape disorders and common shape variation. A summary of current genetic epidemiology of OA is also given, together with current evidence on the genetic overlap between shape variation and OA. Finally, the established genetic risk loci for both joint shape and osteoarthritis are discussed.

Separating positional noise from neutral alignment in multicomponent statistical shape models

Given sufficient training samples, statistical shape models can provide detailed population representations for use in anthropological and computational genetic studies, injury biomechanics, musculoskeletal disease models or implant design optimization. While the technique has become extremely popular for the description of isolated anatomical structures, it suffers from positional interference when applied to coupled or articulated input data. In the present manuscript we describe and validate a novel approach to extract positional noise from such coupled data. The technique was first validated and then implemented in a multicomponent model of the lower limb. The impact of noise on the model itself as well as on the description of sexual dimorphism was evaluated. The novelty of our methodology lies in the fact that no rigid transformations are calculated or imposed on the data by means of idealized joint definitions and by extension the models obtained from them.

Statistical Shape Modeling of Skeletal Anatomy for Sex Discrimination: Their Training Size, Sexual Dimorphism, and Asymmetry

Purpose: Statistical shape modeling provides a powerful tool for describing and analyzing human anatomy. By linearly combining the variance of the shape of a population of a given anatomical entity, statistical shape models (SSMs) identify its main modes of variation and may approximate the total variance of that population to a selected threshold, while reducing its dimensionality. Even though SSMs have been used for over two decades, they lack in characterization of their goodness of prediction, in particular when defining whether these models are actually representative for a given population.

Methods: The current paper presents, to the authors' knowledge, the most extent lower limb anatomy shape model considering the pelvis, femur, patella, tibia, fibula, talus, and calcaneum to date. The present study includes the segmented training shapes (n = 542) obtained from 271 lower limb CT scans. The different models were evaluated in terms of accuracy, compactness, generalizability as well as specificity.

Results: The size of training samples needed in each model so that it can be considered population covering was estimated to approximate around 200 samples, based on the generalizability properties of the different models. Simultaneously differences in gender and patterns in left-right asymmetry were identified and characterized. Size was found to be the most pronounced sexual discriminator whereas intra-individual variations in asymmetry were most pronounced at the insertion site of muscles.

Conclusion: For models aimed at population covering descriptive studies, the number of training samples required should amount a sizeable 200 samples. The geometric morphometric method for sex discrimination scored excellent, however, it did not largely outperformed traditional methods based on discrete measures.

Unpicking observational relationships between hip shape and osteoarthritis: hype or hope?

PURPOSE OF REVIEW

To review recent findings concerning the observational relationship between hip shape and hip osteoarthritis (HOA) and their shared genetic influences, and the potential for clinical application.

RECENT FINDINGS

Recent observational studies have strengthened the evidence that specific shape deformities, such as cam and acetabular dysplasia, are related to HOA. Statistical shape modelling has emerged as a method to measure hip shape holistically, with the added advantage that this can be applied to dual X-ray absorptiometry scan images. This has led to several additional aspects of hip shape variation being identified, such as a wider femoral neck and larger lesser trochanter, in association with HOA. Furthermore, this method has formed the basis of genetic studies identifying novel genetic influences on hip shape, several of which are shared with known genetic risk factors for HOA.

SUMMARY

Shared genetic influences of hip shape and HOA raise the possibility that hip shape plays a casual role in the development of HOA, justifying preventive approaches aiming to combat these adverse consequences.

Collagen glycosylation

Despite the ubiquity of collagens in the animal kingdom, little is known about the biology of the disaccharide Glc(α1-2)Gal(β1-O) bound to hydroxylysine across collagens from sponges to mammals. The extent of collagen glycosylation varies by the types of collagen, with basement membrane collagen type IV being more glycosylated than fibrillar collagens. Beyond true collagens, proteins including collagen domains such as the complement protein 1Q and the hormone adiponectin also feature glycosylated hydroxylysine. Collagen glycosylation is initiated in the endoplasmic reticulum by the galactosyltransferases COLGALT1 and COLGALT2. Mutations in the COLGALT1 gene cause cerebral small vessel abnormality and porencephaly, which are common in collagen type IV deficiency. Beyond the strongly conserved Glc(α1-2)Gal(β1-O) glycan, additional forms of collagen glycosylation have been described in the deep-sea worm Riftia pachyptila and in the giant virus Mimivirus, thereby suggesting that further forms of collagen glycosylation are likely to be identified in the future.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Unfolding of hidden white blood cell count phenotypes for gene discovery using latent class mixed modeling

Resting-state white blood cell (WBC) count is a marker of inflammation and immune system health. There is evidence that WBC count is not fixed over time and there is heterogeneity in WBC trajectory that is associated with morbidity and mortality. Latent class mixed modeling (LCMM) is a method that can identify unobserved heterogeneity in longitudinal data and attempts to classify individuals into groups based on a linear model of repeated measurements. We applied LCMM to repeated WBC count measures derived from electronic medical records of participants of the National Human Genetics Research Institute (NHRGI) electronic MEdical Record and GEnomics (eMERGE) network study, revealing two WBC count trajectory phenotypes. Advancing these phenotypes to GWAS, we found genetic associations between trajectory class membership and regions on chromosome 1p34.3 and chromosome 11q13.4. The chromosome 1 region contains CSF3R, which encodes the granulocyte colony-stimulating factor receptor. This protein is a major factor in neutrophil stimulation and proliferation. The association on chromosome 11 contain genes RNF169 and XRRA1; both involved in the regulation of double-strand break DNA repair.

Significant morphological change in osteoarthritic hips identified over 6-12 months using statistical shape modelling

OBJECTIVE

Predicting who will develop osteoarthritis, assessing how rapidly their disease will progress and monitoring early responses to treatment are key to the development of therapeutic agents able to treat this crippling disease and to their future clinical use. Statistical Shape Modelling (SSM) enables quantification of variations in multiple geometric measures describing the whole hip joint to be considered in concert. This prospective study evaluates the responsiveness of SSM to changes in hip-shape within 1 year.

METHODS

Sixty-two people, mean age 67.1 yrs, were recruited. Dual-energy X-ray Absorptiometry images were taken at three timepoints (baseline, 6 and 12 months). Based on Kellgren-Lawrence grading (KLG) of their baseline images, subjects were classified into control/doubtful OA: KLG < 1 in both hips; moderate OA: KLG = 2; and severe OA: KLG ≥ 3 in their most severe hip. Morphology was quantified using SSM and changes in shape were assessed using generalised estimating equations. Standardized response means (SRMs) were calculated for the first and second 6 month periods, then the full 12 months.

RESULTS

Disease severity ranged from KLG0-KLG4 in the 124 hips assessed at baseline. Three SSM modes (Modes 1, 3 and 4) were associated with OA severity. Across the whole cohort, SRM magnitudes ranged from 0.16 to 0.63. The greatest subgroup SRM (magnitude 0.91) was observed over 12 months in those subjects with moderate OA (KLG2).

CONCLUSIONS

We have demonstrated that SSM can capture changes in hip shape over 6 and 12 months across the entire hip joint providing a sensitive measure of hip OA progression.

The importance of hip shape in predicting hip osteoarthritis

Purpose of review

This narrative review summarizes the evidence relating hip shape and risk of osteoarthritis at the hip, with a focus on the most recent body of work.

Recent findings

Hip OA is a prevalent and potentially disabling condition with few effective non-surgical treatment options. Risk factors for hip OA appear to differ somewhat from those at other sites. Variations in hip morphology, whether assessed through standard geometric measures or statistical modeling methods, seem to increase hip OA risk and may provide a novel approach to interventions to reduce or prevent OA. Such variations have also led to focused surgical interventions to "correct" abnormal shape, although comparisons with non-surgical management are lacking.

Summary

There remains a lack of understanding regarding the optimal management, whether surgical, non-surgical, or a combination, for FAI syndrome. Even less is known regarding other potential morphologic variations that may contribute to OA risk. Additionally, many individuals who have shape variations that would seem to increase their risk will never develop hip OA. Questions remain regarding key risk factors for hip OA development, which individuals should be targeted for therapies, whether directed at symptoms, function, or prevention, and which therapies should be studied and offered. Trials are underway to help address some of these questions.

DXA-derived hip shape is related to osteoarthritis: findings from in the MrOS cohort

OBJECTIVE

Statistical shape modelling (SSM) of radiographs has been used to explore relationships between altered joint shape and hip osteoarthritis (OA). We aimed to apply SSM to Dual-energy X-ray Absorptiometry (DXA) hip scans, and examine associations between resultant hip shape modes (HSMs), radiographic hip OA (RHOA), and hip pain, in a large population based cohort.

METHOD

SSM was performed on baseline hip DXA scans from the Osteoporotic Fractures in Men (MrOS) Study. Associations between the top ten HSMs, and prevalent RHOA from pelvic radiographs obtained 4.6 years later, were analysed in 4100 participants. RHOA was defined as Croft score ≥2. Hip pain was based on pain on walking, hip pain on examination, and Western Ontario and McMaster Universities Arthritis Index (WOMAC).

RESULTS

The five HSMs associated with RHOA showed features of either pincer- or cam-type deformities. HSM 1 (increased pincer-type deformity) was positively associated with RHOA [1.23 (1.09, 1.39)] [odds ratio (OR) and 95% CI]. HSM 8 (reduced pincer-type deformity) was inversely associated with RHOA [0.79 (0.70, 0.89)]. HSM 10 (increased cam-type deformity) was positively associated with RHOA [1.21 (1.07, 1.37)]. HSM 3 and HSM 4 (reduced cam-type deformity) were inversely associated with RHOA [0.73 (0.65, 0.83) and 0.82 (0.73, 0.93), respectively]. HSM 3 was inversely related to pain on examination [0.84 (0.76, 0.92)] and walking [0.88, (0.81, 0.95)], and to WOMAC score [0.87 (0.80, 0.93)].

CONCLUSIONS

DXA-derived measures of hip shape are associated with RHOA, and to a lesser extent hip pain, possibly reflecting their role in the pathogenesis of hip OA.

Statistical shape modelling of hip and lumbar spine morphology and their relationship in the MRC National Survey of Health and Development

The anatomical shape of bones and joints is important for their proper function but quantifying this, and detecting pathological variations, is difficult to do. Numerical descriptions would also enable correlations between joint shapes to be explored. Statistical shape modelling (SSM) is a method of image analysis employing pattern recognition statistics to describe and quantify such shapes from images; it uses principal components analysis to generate modes of variation describing each image in terms of a set of numerical scores after removing global size variation. We used SSM to quantify the shapes of the hip and the lumbar spine in dual-energy x-ray absorptiometry (DXA) images from 1511 individuals in the MRC National Survey of Health and Development at ages 60-64 years. We compared shapes of both joints in men and women and hypothesised that hip and spine shape would be strongly correlated. We also investigated associations with height, weight, body mass index (BMI) and local (hip or lumber spine) bone mineral density. In the hip, all except one of the first 10 modes differed between men and women. Men had a wider femoral neck, smaller neck-shaft angle, increased presence of osteophytes and a loss of the femoral head/neck curvature compared with women. Women presented with a flattening of the femoral head and greater acetabular coverage of the femoral head. Greater weight was associated with a shorter, wider femoral neck and larger greater and lesser trochanters. Taller height was accompanied by a flattening of the curve between superior head and neck and a larger lesser trochanter. Four of the first eight modes describing lumbar spine shape differed between men and women. Women tended to have a more lordotic spine than men with relatively smaller but caudally increasing anterior-posterior (a-p) vertebral diameters. Men were more likely to have a straighter spine with larger vertebral a-p diameters relative to vertebral height than women, increasing cranially. A weak correlation was found between body weight and a-p vertebral diameter. No correlations were found between shape modes and height in men, whereas in women there was a weak positive correlation between height and evenness of spinal curvature. Linear relationships between hip and spine shapes were weak and inconsistent in both sexes, thereby offering little support for our hypothesis. In conclusion, men and women entering their seventh decade have small but statistically significant differences in the shapes of their hips and their spines. Associations with height, weight, BMI and BMD are small and correspond to subtle variations whose anatomical significance is not yet clear. Correlations between hip and spine shapes are small.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis

OBJECTIVE

Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA.

METHODS

A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage.

RESULTS

We report suggestive evidence (p<5×10-6) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the LRCH1 gene (rs754106, p=1.49×10-7, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between STT3B and GADL1 (rs6766414, p=3.13×10-6, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). STT3B was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10-4).

CONCLUSIONS

Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

Analysis of the three-dimensional anatomical variance of the distal radius using 3D shape models

BACKGROUND

Various medical fields rely on detailed anatomical knowledge of the distal radius. Current studies are limited to two-dimensional analysis and biased by varying measurement locations. The aims were to 1) generate 3D shape models of the distal radius and investigate variations in the 3D shape, 2) generate and assess morphometrics in standardized cut planes, and 3) test the model's classification accuracy.

METHODS

The local radiographic database was screened for CT-scans of intact radii. 1) The data sets were segmented and 3D surface models generated. Statistical 3D shape models were computed (overall, gender and side separate) and the 3D shape variation assessed by evaluating the number of modes. 2) Anatomical landmarks were assigned and used to define three standardized cross-sectional cut planes perpendicular to the main axis. Cut planes were generated for the mean shape models and each individual radius. For each cut plane, the following morphometric parameters were calculated and compared: maximum width and depth, perimeter and area. 3) The overall shape model was utilized to evaluate the predictive value (leave one out cross validation) for gender and side identification within the study population.

RESULTS

Eighty-six radii (45 left, 44% female, 40 ± 18 years) were included. 1) Overall, side and gender specific statistical 3D models were successfully generated. The first mode explained 37% of the overall variance. Left radii had a higher shape variance (number of modes: 20 female / 23 male) compared to right radii (number of modes: 6 female / 6 male). 2) Standardized cut planes could be defined using anatomical landmarks. All morphometric parameters decreased from distal to proximal. Male radii were larger than female radii with no significant side difference. 3) The overall shape model had a combined median classification probability for side and gender of 80%.

CONCLUSIONS

Statistical 3D shape models of the distal radius can be generated using clinical CT-data sets. These models can be used to assess overall bone variance, define and analyze standardized cut-planes, and identify the gender of an unknown sample. These data highlight the potential of shape models to assess the 3D anatomy and anatomical variance of human bones.

Osteoarthritis year in review 2016: genetics, genomics and epigenetics

The purpose of this narrative review is to provide an overview of last year's publications in the field of genetics, genomics and epigenetics in the osteoarthritis (OA) field. Major themes arising from a Pubmed search on (epi)genetics in OA were identified. In addition, general developments in the fast evolving field of (epi)genetics are reviewed and relevance for the OA field is summarized. In the last 5 years, a number of genome-wide association studies have identified a modest number of genetic loci associated to OA. Continued functional research into these DNA variants is showing putative biological mechanisms underlying these associations. Over the last year, no additional large genome-wide association studies were published, but there clearly remains much to be discovered in the OA genetic field. A lot of research has been done into the epigenetics of OA over the last year. Several genome-wide screens examining the methylome of osteoarthritic cartilage were done. Pathway analysis confirmed deregulation of developmental and extracellular pathways in OA cartilage. Over the last year many microRNAs (miRNAs) have been identified that potentially play important roles in cartilage homeostasis and/or OA process. Continued research will learn whether these identified miRNAs are truly causal and can be used in clinical applications. Many of the epigenetic findings need further confirmation, but they highlight potential novel pathways involved in cartilage biology and OA.

The complex genetics of gait speed: genome-wide meta-analysis approach

Emerging evidence suggests that the basis for variation in late-life mobility is attributable, in part, to genetic factors, which may become increasingly important with age. Our objective was to systematically assess the contribution of genetic variation to gait speed in older individuals. We conducted a meta-analysis of gait speed GWASs in 31,478 older adults from 17 cohorts of the CHARGE consortium, and validated our results in 2,588 older adults from 4 independent studies. We followed our initial discoveries with network and eQTL analysis of candidate signals in tissues. The meta-analysis resulted in a list of 536 suggestive genome wide significant SNPs in or near 69 genes. Further interrogation with Pathway Analysis placed gait speed as a polygenic complex trait in five major networks. Subsequent eQTL analysis revealed several SNPs significantly associated with the expression of PRSS16, WDSUB1 and PTPRT, which in addition to the meta-analysis and pathway suggested that genetic effects on gait speed may occur through synaptic function and neuronal development pathways. No genome-wide significant signals for gait speed were identified from this moderately large sample of older adults, suggesting that more refined physical function phenotypes will be needed to identify the genetic basis of gait speed in aging.

Identification of DNA methylation changes associated with disease progression in subchondral bone with site-matched cartilage in knee osteoarthritis

Subchondral bone plays a key role in the development of osteoarthritis, however, epigenetics of subchondral bone has not been extensively studied. In this study, we examined the genome-wide DNA methylation profiles of subchondral bone from three regions on tibial plateau representing disease progression using HumanMethylation450 BeadChip to identify progression associated DNA methylation alterations. Significant differential methylated probes (DMPs) and differential methylated genes (DMGs) were identified in the intermediate and late stages and during the transition from intermediate to late stage of OA in the subchondral bone. Over half of the DMPs were hyper-methylated. Genes associated with OA and bone remodeling were identified. DMGs were enriched in morphogenesis and development of skeletal system, and HOX transcription factors. Comparison of DMGs identified in subchondral bone and site-matched cartilage indicated that DNA methylation changes occurred earlier in subchondral bone and identified different methylation patterns at the late stage of OA. However, shared DMPs, DMGs and common pathways that implicated the tissue reparation were also identified. Methylation is one key mechanism to regulate the crosstalk between cartilage and subchondral bone.

Analysis of genetics and DNA methylation in osteoarthritis: What have we learnt about the disease?

Osteoarthritis (OA) is a chronic musculoskeletal disease characterised by the destruction of articular cartilage, synovial inflammation and bone remodelling. Disease aetiology is complex and highly heritable, with genetic variation estimated to contribute to 50% of OA occurrence. Epigenetic alterations, including DNA methylation changes, have also been implicated in OA pathophysiology. This review examines what genetic and DNA methylation studies have taught us about the genes and pathways involved in OA pathology. The influence of DNA methylation on the molecular mechanisms underlying OA genetic risk and the consequence of this interaction on disease susceptibility and penetrance are also discussed.