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.

Multi-tissue epigenetic analysis of the osteoarthritis susceptibility locus mapping to the plectin gene PLEC

OBJECTIVE

In cartilage, the osteoarthritis (OA) associated single nucleotide polymorphism (SNP) rs11780978 correlates with differential expression of PLEC, and with differential methylation of PLEC CpG dinucleotides, forming eQTLs and mQTLs respectively. This implies that methylation links chondrocyte genotype and phenotype, thus driving the functional effect of this genetic risk signal. PLEC encodes plectin, a cytoskeletal protein that enables tissues to respond to mechanical forces. We sought to assess whether these PLEC functional effects were cartilage specific.

METHOD

Cartilage, fat pad, synovium and peripheral blood were collected from patients undergoing arthroplasty. PLEC CpGs were analysed for mQTLs and allelic expression imbalance (AEI) was performed to test for eQTLs. Plectin was knocked down in a mesenchymal stem cell (MSC) line using CRISPR/Cas9 and cells phenotyped by RNA-sequencing.

RESULTS

mQTLs were discovered in fat pad, synovium and blood. Their effects were however stronger in the joint tissues and of comparable effect between these tissues. We observed AEI in synovium in the same direction as for cartilage and correlations between methylation and PLEC expression. Knocking-down plectin impacted on pathways reported to have a role in OA, including Wnt signalling, glycosaminoglycan biosynthesis and immune regulation.

CONCLUSIONS

Synovium is also a target of the rs11780978 OA association functionally operating on PLEC. In fat pad, mQTLs were identified but these did not correlate with PLEC expression, suggesting the functional effect is not joint-wide. Our study highlights interplay between genetic risk, DNA methylation and gene expression in OA, and reveals clear differences between tissues from the same diseased joint.

Discovery and analysis of methylation quantitative trait loci (mQTLs) mapping to novel osteoarthritis genetic risk signals

OBJECTIVE

Osteoarthritis (OA) is polygenic with over 90 independent genome-wide association loci so far reported. A key next step is the identification of target genes and the molecular mechanisms through which this genetic risk operates. The majority of OA risk-conferring alleles are predicted to act by modulating gene expression. DNA methylation at CpG dinucleotides may be a functional conduit through which this occurs and is detectable by mapping methylation quantitative trait loci, or mQTLs. This approach can therefore provide functional insight into OA risk and will prioritize genes for subsequent investigation. That was our goal, with a focus on the largest set of OA loci yet to be reported.

METHOD

We investigated DNA methylation, genotype and RNA sequencing data derived from the cartilage of patients who had undergone arthroplasty and combined this with in silico analyses of expression quantitative trait loci, epigenomes and chromatin interactions.

RESULTS

We investigated 42 OA risk loci and in ten of these we identified 24 CpGs in which methylation correlated with genotype (false discovery rate (FDR) P-values ranging from 0.049 to 1.73x10^-25). In silico analyses of these mQTLs prioritised genes and regulatory elements at the majority of the ten loci, with COLGALT2 (encoding a collagen galactosyltransferase), COL11A2 (encoding a polypeptide chain of type XI collagen) and WWP2 (encoding a ubiquitin ligase active during chondrogenesis) emerging as particularly compelling target genes.

CONCLUSION

We have highlighted the pivotal role of DNA methylation as a link between genetic risk and OA and prioritized genes for further investigation.

Age-related changes in mesenchymal stem cells identified using a multi-omics approach

Mesenchymal stem cells (MSC) are capable of multipotent differentiation into connective tissues and as such are an attractive source for autologous cell-based treatments for many clinical diseases and injuries. Ageing is associated with various altered cellular phenotypes coupled with a variety of transcriptional, epigenetic and translational changes. Furthermore, the regeneration potential of MSCs is reduced with increasing age and is correlated with changes in cellular functions. This study used a systems biology approach to investigate the transcriptomic (RNASeq), epigenetic (miRNASeq and DNA methylation) and protein alterations in ageing MSCs in order to understand the age-related functional and biological variations, which may affect their applications to regenerative medicine. We identified no change in expression of the cellular senescence markers. Alterations were evident at both the transcriptional and post-transcriptional level in a number of transcription factors. There was enrichment in genes involved in developmental disorders at mRNA and differential methylated loci (DML) level. Alterations in energy metabolism were apparent at the DML and protein level. The microRNA miR-199b-5p, whose expression was reduced in old MSCs, had predicted gene targets involved in energy metabolism and cell survival. Additionally, enrichment of DML and proteins in cell survival was evident. Enrichment in metabolic processes was revealed at the protein level and in genes identified as undergoing alternate splicing. Overall, an altered phenotype in MSC ageing at a number of levels implicated roles for inflamm-ageing and mitochondrial ageing. Identified changes represent novel insights into the ageing process, with implications for stem cell therapies in older patients.

The role of inflammation-related genes in osteoarthritis

In this review article we examine the role of inflammation-related genes in osteoarthritis (OA) from the perspective of genetics, epigenetics and gene expression. There have been great strides in such genomic analyses of OA in recent years thanks to the study of adequately powered patient cohorts, the detailed analysis of candidate genes, and the application of genome-wide approaches. These have led to some unexpected and therefore exciting discoveries, implicating pathways that would not necessarily have been predicted to have a role in this common arthritis. Inflammatory-related genes sit firmly in the candidate camp based on prior observations that the OA disease process can have an inflammatory component. What is clear from the genetic studies published to date is that there is no compelling evidence that DNA variation in inflammatory genes is an OA risk factor. This conclusion may of course change as ever more powerful association studies are conducted. There is, however, compelling evidence that epigenetic effects involving inflammatory genes are a component of OA and that alteration in the expression of these genes is also highly relevant to the disease process. We may in fact be close to demonstrating, at the genomic level, a clear separation of OA patients into those in whom inflammation is a key driver of the disease and those in whom it is not. This has obvious implications for the design of trials of novel OA interventions and may also guide the intelligent re-purposing of anti-inflammatory therapies.

Identification and analysis of a SMAD3 cis-acting eQTL operating in primary osteoarthritis and in the aneurysms and osteoarthritis syndrome

OBJECTIVE

The TGF-β pathway plays a central role in joint development with polymorphism in TGF-β pathway genes implicated in osteoarthritis susceptibility. One association is to rs12901499, within intron 1 of SMAD3. Since rs12901499 is not in linkage disequilibrium with a non-synonymous polymorphism, it is likely the association is operating by influencing expression of SMAD3.

DESIGN

Using tissues from the joints of primary osteoarthritis patients who had undergone joint replacement we measured the overall expression of SMAD3 by quantitative real-time PCR. We also measured allelic expression of SMAD3 using these tissues and vascular smooth muscle cells from patients with aneurysms and osteoarthritis syndrome, a rare condition featuring early-onset osteoarthritis. We tested the functional effect of SNPs in vitro using luciferase assays and assessed association with osteoarthritis using a large osteoarthritis case-control dataset.

RESULTS

We observed that genotype at rs12901499 did not correlate with overall SMAD3 expression or allelic expression. However, genotype at a 3'UTR SNP, rs8031440, did correlate with SMAD3 expression in cartilage (P = 0.005) which was supported by allelic expression data showing that the G allele correlated with decreased SMAD3 expression in joint tissues and vascular smooth muscle cells. This G allele was underrepresented in osteoarthritis cases vs controls (P = 0.027, odds ratio = 0.921). rs8031440 is in perfect linkage disequilibrium with five other SMAD3 3'UTR SNPs and our luciferase analysis identified rs3743342 and rs12595334 as being functional.

CONCLUSION

SMAD3 is subject to cis-acting regulatory polymorphism in the tissues of relevance to both primary osteoarthritis and the aneurysms-osteoarthritis syndrome.

Identification of the pathogenic pathways in osteoarthritic hip cartilage: commonality and discord between hip and knee OA

OBJECTIVE

To define for the first time the transcriptomes of normal and end-stage osteoarthritis (OA) hip cartilage.

MATERIALS AND METHODS

RNA was isolated from cartilage within 2h of joint replacement surgery. Gene expression was analyzed using Agilent GeneSpring GX 11 following hybridization to Illumina Human HT-12 V3 microarrays. Real-time reverse-transcription polymerase chain reaction (RT-PCR) was used to validate the expression of six genes identified by microarray as differentially expressed. Gene Set Enrichment Analysis (GSEA) and Ingenuity Pathway Analysis (IPA) were used to investigate enriched functions or canonical pathways amongst differentially expressed genes respectively.

RESULTS

In total we identified 998 differentially expressed genes (fold change ≥ ±1.5, P-value ≤ 0.01) between neck of femur fracture (NOF) (n = 10) and OA hip (n = 9) patient cartilage. These differentially expressed genes were enriched within 71 canonical pathways. A comparison between a comparable knee dataset(20) only identified 229 genes similarly differentially expressed although remarkably 34 canonical pathways overlapped between experiments.

CONCLUSIONS

This study is the first to report a comprehensive gene expression analysis of human hip OA cartilage compared to control (NOF) cartilage at the whole-genome level. Our differential gene expression dataset shows excellent correlation with similar defined studies using comparable tissue but reveals discord between hip and knee OA at the individual gene status but with commonality with regards the molecular pathways involved.

Genetic indicators and susceptibility to osteoarthritis

A large number of experiments have been performed to identify genetic loci that influence osteoarthritis (OA) susceptibility, with a particular focus on the primary form of the disease. Unfortunately, the currently reported candidate-gene, genome-wide linkage scans and genome-wide association scans have tended to highlight the heterogeneous nature of OA rather than generate statistically compelling genome signals. Nevertheless, some breakthroughs have been made. For example, genetic susceptibility within genes coding for components of the transforming growth factor β pathway has emerged as a particularly interesting find, while completely novel loci are also being uncovered, such as the signal to a cluster of genes on chromosome 7q22. It also appears that quantitative effects on gene expression, rather than qualitative effects on protein function, are particularly important, and that we need to consider the effects of genetic susceptibility in joint formation as much as we do in joint maintenance. Nevertheless, we are still only at the beginning of our search, and much more sophisticated clinical and laboratory approaches will need to be applied before we get a clear understanding of the genetic indicators that influence OA susceptibility.

Deep sequencing of GDF5 reveals the absence of rare variants at this important osteoarthritis susceptibility locus

OBJECTIVE

The common single nucleotide polymorphism (SNP) rs143383 in the 5' untranslated region (5'UTR) of growth and differentiation factor 5 (GDF5) is strongly associated with osteoarthritis (OA) and influences GDF5 allelic expression in vitro and in the joint tissues of OA patients. This effect is modulated in cis by another common SNP, also located within the 5'UTR, whilst a common SNP in the 3'UTR influences allelic expression independent of rs143383. DNA variants can be common, rare or extremely rare/unique. To therefore enhance our understanding of the allelic architecture of this very important OA susceptibility locus we sequenced the gene for potentially functional and novel rare variants.

METHOD

Using the Sanger method we sequenced GDF5 in 992 OA patients and 944 controls, with DNA changes identified by sequencing software. We encompassed the protein-coding region of the two GDF5 exons, both untranslated regions and approximately 100 bp of the proximal promoter of the gene.

RESULTS

We detected 13 variants. Six were extremely rare with minor allele frequencies (MAFs) of ≤ 0.0006. One is in a predicted transcription factor binding site in the GDF5 promoter whilst two substitute conserved amino acids. The remaining seven variants were common and are previously known variants, with MAFs ranging from 0.025 to 0.39. There was a complete absence of variants with frequencies in-between the extremely rare (n=6) and the common (n=7).

CONCLUSIONS

This is the first report of the deep sequencing of an OA susceptibility locus. The absence of rare variants informs us that within the regions of the gene that we have sequenced GDF5 does not harbour any novel variants that are able to contribute, at a population level, to the OA association signal mediated by rs143383 nor does it harbour, at a population level, any novel variants that can influence OA susceptibility independent of rs143383.

Osteoarthritis year 2010 in review: genetics

At the 2010 Osteoarthritis Research Society International (OARSI) congress in Brussels I was asked to present on "Genetics" in the "Year in Review" session. This gave me an opportunity to reflect on the talk that I gave on the same topic at the 2007 OARSI congress in Ft Lauderdale, Florida. My 2007 talk was very upbeat and ended with the hostage to fortune statement "The next few years will provide tremendous clarity in our genetic understanding of osteoarthritis".

Insights into the genetic architecture of osteoarthritis from stage 1 of the arcOGEN study

Objectives

The genetic aetiology of osteoarthritis has not yet been elucidated. To enable a well-powered genome-wide association study (GWAS) for osteoarthritis, the authors have formed the arcOGEN Consortium, a UK-wide collaborative effort aiming to scan genome-wide over 7500 osteoarthritis cases in a two-stage genome-wide association scan. Here the authors report the findings of the stage 1 interim analysis.

Methods

The authors have performed a genome-wide association scan for knee and hip osteoarthritis in 3177 cases and 4894 population-based controls from the UK. Replication of promising signals was carried out in silico in five further scans (44 449 individuals), and de novo in 14 534 independent samples, all of European descent.

Results

None of the association signals the authors identified reach genome-wide levels of statistical significance, therefore stressing the need for corroboration in sample sets of a larger size. Application of analytical approaches to examine the allelic architecture of disease to the stage 1 genome-wide association scan data suggests that osteoarthritis is a highly polygenic disease with multiple risk variants conferring small effects.

Conclusions

Identifying loci conferring susceptibility to osteoarthritis will require large-scale sample sizes and well-defined phenotypes to minimise heterogeneity.

Interleukin-1 region meta-analysis with osteoarthritis phenotypes

OBJECTIVE

Several research groups have examined osteoarthritis (OA) association with Interleukin-1 (IL-1) region markers and haplotypes. The results have been suggestive for hand OA, negative for knee OA, and conflicting for hip OA.

DESIGN

Our aim was to address conflicts employing meta-analytical methods on data from 1238 European-descent cases with various OA phenotypes and 1269 European-descent controls from four study centers. We imputed some missing genotype data and reconstructed IL-1 region extended haplotypes. A previously reported 7-marker IL1A-IL1B-IL1RN extended risk haplotype was tested for association with each specific index phenotype.

RESULTS

For hip OA, data from three centers showed heterogeneity of extended-risk-haplotype effect, two panels showing trend toward risk and another showing protection, with overall odds ratio (OR) 1.24 (95% Confidence interval (CI) 0.45-3.41, P 0.67). The heterogeneity fell partly along control ascertainment lines, chiefly between controls ascertained as spouses of arthroplasty patients and controls identified through population radiographic survey. For knee OA, the results showed no heterogeneity and no significant extended-risk-haplotype effect. For hand OA, the results showed little heterogeneity and a modest trend toward positive association (summary OR 1.34, 95% CI 0.83-2.17 P 0.23). Using a Bayesian partition modeling approach, the 7-marker extended haplotypes showed no significant effect on any OA phenotype examined. A 3-single-nucleotide polymorphism (SNP) IL1B-IL1RN haplotype rs1143627-rs16944-rs419598 showed a trend toward hand OA association (posterior probability of association 0.72) with the most prominent feature being protection from a specific haplotype representing a partial mirror image of the extended risk haplotype (OR estimated at 0.46).

CONCLUSIONS

The meta-analysis data do not confirm but only suggest that some hand and hip OA risk could be associated with the IL-1 region, particularly centered in IL1B and possibly also IL1RN.

Association of a nsSNP in ADAMTS14 to some osteoarthritis phenotypes

OBJECTIVE

To investigate the effect in OA (Osteoarthritis) susceptibility of putative damaging changes in ADAM (A Disintegrin And Metalloprotease) and ADAMTS (ADAM with ThromboSpondin motif) proteases.

METHODS

Non-synonymous single nucleotide polymorphisms (nsSNP) in 18 ADAMTS and 31 ADAM genes were analyzed with two software applications for prediction of functional damage. Four putative damaging nsSNP were found in ADAMTS2, ADAMTS14, ADAMTS16 and ADAM12, respectively. These nsSNPs were analyzed in case-control sample collections with a variety of phenotypes totalling 3217 OA patients and 2214 healthy controls, all of them Caucasians.

RESULTS

No statistically significant differences were found in ADAMTS2, ADAMTS16 and ADAM12 nsSNPs. Conversely, the rare allele of the rs4747096 nsSNP in ADAMTS14 was overrepresented in women requiring joint replacement because of knee OA (O.R.(M-H) (odds ratio. Mantel-Haenszel)=1.41, 95% C.I.=1.1-1.8; P=0.002) and in patients with symptomatic hand OA (O.R.=1.37, 95% C.I.=1.0-1.9; P=0.047). A non significant increase in the frequency of the same allele was also found in patients with hip OA requiring prosthesis (O.R.(M-H)=1.14, 95% C.I.=1.0-1.3; P=0.08). No association was found with other OA phenotypes.

CONCLUSION

Our findings implicate ADAMTS14 in OA, specifically in knee OA requiring joint replacement in women and, possibly, in hand OA. Independent association of ADAMTS14 genetic variation to knee OA in women has been communicated. ADAMTS14 involvement, if confirmed, will open a new area of interest in OA pathogenesis because of its role in the maturation of collagen fibers.

Testing the druggable endothelial differentiation gene 2 knee osteoarthritis genetic factor for replication in a wide range of sample collections

OBJECTIVES

To replicate a previously reported association with osteoarthritis (OA) of the promoter single nucleotide polymorphism (SNP) rs10980705 in the endothelial differentiation gene 2 (EDG2).

METHODS

Five collections of samples, four from Europe and one from China, were studied. They included patients with 3 OA phenotypes: 1501 with knee OA, 1497 with hip OA and 376 with generalised OA. A total of 2521 controls were also studied. Allele and genotype frequencies of the rs10980705 SNP were analysed in each individual sample collection and in pooled data. In addition, a meta-analysis to incorporate results from the original Japanese report was performed.

RESULTS

The association of the rs10980705 SNP with knee OA was not replicated in any of the five sample collections studied or in their combined analysis (odds ratio (OR) 1.10, 95% CI 0.98 to 1.22; p = 0.10). Meta-analysis of all data, including the original Japanese study, did show association with knee OA (OR 1.15, 95% CI 1.06 to 1.26; p = 0.002) but the effect was accounted for by the Japanese data and was less significant than the original report. No association was found with hip OA or with generalised OA.

CONCLUSIONS

The original report of a promising genetic association between a druggable G-protein coupled receptor, EDG2, and knee OA has not been replicated. This lack of replication could be due to a modest effect of the promoter polymorphism that will require even larger studies (the winners curse) although a more pronounced effect in the Asian population vs Europeans cannot be excluded.

Genetic associations in peripheral joint osteoarthritis and spinal degenerative disease: a systematic review

We conducted a systematic review of genetic association studies for osteoarthritis of the peripheral joints (OA) and spinal degenerative disease (SDD). Electronic searches were carried out for any English language article reporting on a gene association study for either OA or SDD published up until the end of 2006. A team of seven reviewers used a standardised template to extract data in duplicate. In all, 90 studies fulfilled our inclusion criteria, reporting a total of 94 significant associations from 83 different genes. We found relatively few instances in which a specific gene-disease association had been analysed by more than one study, and there were 14 cases in which significant associations were replicated in independent studies (at joints associated with the AGC1, ASPN, COL9A2, COL9A3, COL11A2, ESR1, FZRB, HFE, IL1A, IL1RN, PTGS2 and VDR genes).

METHOD

logical and reporting problems were widespread, including failure to report full results, missing population details, multiple testing, and over-reliance on subgroup analysis. In summary, the complex phenotypes of OA and SDD may have made it difficult for researchers to focus their efforts. The field is dominated by isolated analyses of disparate potential associations, a problem that is amplified by the frequent analysis of different polymorphisms within individual genes. Flaws in study methodology and interpretation undoubtedly increase the risk of publication bias. Closer adherence to published recommendations (in particular those produced by HuGENet) will help to ensure that future studies are well-designed and build on current understanding, rather than simply adding to the growing bank of potential associations.

Allelic expression analysis suggests that cis-acting polymorphism of FRZB expression does not contribute to osteoarthritis susceptibility

OBJECTIVE

The FRZB gene codes for secreted frizzled-related protein 3 (sFRP3), a soluble antagonist of wnt signalling that is required for maintaining cartilage integrity. Two common single nucleotide polymorphisms (SNPs) that code for the substitution of conserved arginine residues have previously been identified in FRZB and found to be associated with osteoarthritis (OA). Functional studies revealed that the arginine substitutions lead to a loss-of-function of sFRP3 activity. We set out to assess whether cis-acting polymorphism in the regulation of FRZB expression exists, as this may be an additional mechanism through which sFRP3 activity could be modulated.

METHODS

RNA was extracted from the articular cartilage chondrocytes of 25 individuals who had undergone joint replacement for OA and who were heterozygous for one of the two FRZB SNPs. Allelic output was measured by single base extension in to the SNPs and deviations from the expected 1:1 pattern were assessed using the Mann-Whitney U test.

RESULTS

Differential allelic expression was observed in six of the 25 individuals. However, the average fold difference in allelic expression in the six was only 1.19.

CONCLUSIONS

The presence of a small degree of differential allelic expression in a low proportion (24%) of the individuals studied suggests that polymorphism in FRZB cis-acting regulatory elements can be discounted as a major factor that could influence the development of OA.

Genetic association analysis of LRCH1 as an osteoarthritis susceptibility locus

OBJECTIVE

A genetic association with knee osteoarthritis (OA) of a single nucleotide polymorphism (SNP) in intron 1 of the LRCH1 gene was recently reported in a UK Caucasian case-control sample and confirmed in a Newfoundland Caucasian sample. Our objective was to assess whether the SNP was associated with OA in our large UK Caucasian sample.

METHODS

The SNP was genotyped in 1521 cases that had undergone elective joint replacement of the hip (1098 cases), of the knee (340 cases) or of the hip and knee (83 cases) due to end-stage primary OA. The SNP was also genotyped in 736 controls of similar ages in the cases.

RESULTS

There was no significant difference (all P-values >0.05) in genotype or allele frequencies between our cases and our controls. There was also no significant difference when the cases were stratified by sex, by joint replaced or by sex combined with joint replaced.

CONCLUSION

Our data on 2257 individuals implies that the LRCH1 intron 1 SNP is not a risk factor for OA aetiology of the knee or of the hip in our UK Caucasian sample.

The CALM1 core promoter polymorphism is not associated with hip osteoarthritis in a United Kingdom Caucasian population

OBJECTIVE

A convincing genetic association with hip osteoarthritis (OA) of a functional single nucleotide polymorphism (SNP) in the core promoter of the calmodulin 1 gene CALM1 was recently reported in a Japanese population. The T-allele of the SNP encoded OA susceptibility and this was mediated by a reduced expression of CALM1. Our objective was to assess whether the SNP was also associated with hip OA in UK Caucasians.

METHODS

The SNP was genotyped in 920 cases that had undergone elective joint replacement of the hip due to end-stage primary OA and in 752 age-matched controls.

RESULTS

Our study had greater than 97% power to observe an effect comparable to that seen in the Japanese study. However, there was no significant difference (P< or =0.05) in genotype or allele frequencies between our cases and our controls. There was also no significant difference when the cases were stratified by sex.

CONCLUSION

Our data on a cohort of 1672 individuals implies that the CALM1 core promoter polymorphism is not a risk factor for OA etiology in Caucasians. Our study does not call in to question the veracity of the Japanese report. Instead it highlights the heterogeneous nature of OA genetic susceptibility.

Genetic association analysis of the IGFBP7, ADAMTS3, and IL8 genes as the potential osteoarthritis susceptibility that maps to chromosome 4q

OBJECTIVES

To determine by genetic association analysis whether the 4q functional candidate genes IGFBP7, ADAMTS3, and IL8 might encode for susceptibility to osteoarthritis (OA).

METHODS

Expression of IGFBP7, ADAMTS3, and IL8 in adult OA articular cartilage chondrocytes was demonstrated by reverse transcription-polymerase chain reaction. The genes were screened for common polymorphic DNA variants by direct sequencing of exons, intron-exon boundaries, and the 5' and 3' untranslated regions. The variants were genotyped in the female probands from the 146 families which each contained two or more sisters who had undergone total hip replacement (THR) and in 375 female controls matched for age. Variants showing evidence for association were subsequently genotyped in 244 female-THR patients with OA. Allele frequencies between the probands (or patients) and the controls were compared by chi(2) analysis.

RESULTS

Fourteen common variants were identified in the three genes. An ADAMTS3 single nucleotide polymorphism was associated in the probands (p = 0.015) and an ADAMTS3 insertion/deletion approached significance (p = 0.059). However, neither variant was associated in the additional 244 patients with hip OA, with p values of 0.12 and 0.19, respectively.

CONCLUSIONS

The analysis implies that the chromosome 4q female hip OA susceptibility is not coded for by polymorphism within the functional candidates IGFBP7, ADAMTS3, or IL8.

Mutations in the region encoding the von Willebrand factor A domain of matrilin-3 are associated with multiple epiphyseal dysplasia

Multiple epiphyseal dysplasia (MED) is a relatively mild and clinically variable osteochondrodysplasia, primarily characterized by delayed and irregular ossification of the epiphyses and early-onset osteoarthritis. Mutations in the genes encoding cartilage oligomeric matrix protein (COMP) and type IX collagen (COL9A2 and COL9A3) have previously been shown to cause different forms of MED (refs. 4-13). These dominant forms of MED (EDM1-3) are caused by mutations in the genes encoding structural proteins of the cartilage extracellular matrix (ECM); these proteins interact with high affinity in vitro. A recessive form of MED (EDM4) has also been reported; it is caused by a mutation in the diastrophic dysplasia sulfate transporter gene (SLC26A). A genomewide screen of family with autosomal-dominant MED not linked to the EDM1-3 genes provides significant genetic evidence for a MED locus on the short arm of chromosome 2 (2p24-p23), and a search for candidate genes identified MATN3 (ref. 18), encoding matrilin-3, within the critical region. Matrilin-3 is an oligomeric protein that is present in the cartilage ECM. We have identified two different missense mutations in the exon encoding the von Willebrand factor A (vWFA) domain of matrilin-3 in two unrelated families with MED (EDM5). These are the first mutations to be identified in any of the genes encoding the matrilin family of proteins and confirm a role for matrilin-3 in the development and homeostasis of cartilage and bone.

Genetic epidemiology of primary osteoarthritis

Primary osteoarthritis (OA) is a late onset disease that fits most accurately into the oligogenic, multifactorial class of genetic diseases. Twin pair and family risk studies have highlighted a surprisingly large genetic component to OA and have prompted the search for predisposing genes. These searches have taken three forms: (1) parametric linkage analysis of rare families in which OA segregates as a Mendelian trait, (2) model-free linkage analysis of affected sibling pairs, and (3) association analysis of known candidate genes. Within the past year linkage analysis studies have highlighted that chromosomes 2, 4, 6, 7, 11, 16, and the X may each harbor an OA susceptibility gene. Chromosomes 2, 4, and 16 were identified in multiple genome scans and are therefore the most likely to encode susceptibility. Association analysis of candidates suggests that the syntenic genes for type II collagen and the vitamin D receptor (12q12--q13.1) may also encode for OA susceptibility.

Eukaryotic translation initiation factor 5 (eIF5) acts as a classical GTPase-activator protein

GTP hydrolysis occurs at several specific stages during the initiation, elongation, and termination stages of mRNA translation. However, it is unclear how GTP hydrolysis occurs; it has previously been suggested to involve a GTPase active center in the ribosome, although proof for this is lacking. Alternatively, it could involve the translation factors themselves, e.g., be similar to the situation for small G in which the GTPase active site involves arginine residues contributed by a further protein termed a GTPase-activator protein (GAP). During translation initiation in eukaryotes, initiation factor eIF5 is required for hydrolysis of GTP bound to eIF2 (the protein which brings the initiator Met-tRNA(i) to the 40S subunit). Here we show that eIF5 displays the hallmarks of a classical GAP (e.g., RasGAP). Firstly, its interaction with eIF2 is enhanced by AlF(4)(-). Secondly, eIF5 possesses a conserved arginine (Arg15) which, like the "arginine fingers" of classical GAPs, is flanked by hydrophobic residues. Mutation of Arg15 to methionine abolishes the ability of eIF5 either to stimulate GTP hydrolysis or to support mRNA translation in vitro. Mutation studies suggest that a second conserved arginine (Arg48) also contributes to the GTPase active site of the eIF2.eIF5 complex. Our data thus show that eIF5 behaves as a classical GAP and that GTP hydrolysis during translation involves proteins extrinsic to the ribosome. Indeed, inspection of their sequences suggests that other translation factors may also act as GAPs.

Linkage analysis of chromosome 2q in osteoarthritis

BACKGROUND

In independent linkage studies chromosome 2q11-q24 and chromosome 2q23-35 have previously been implicated as regions potentially harbouring susceptibility loci for osteoarthritis (OA).

OBJECTIVE

To test chromosome 2q for linkage to idiopathic osteoarthritis.

METHODS

Using a cohort of 481 OA families that each contained at least one affected sibling pair with severe end-stage disease (ascertained by hip or knee joint replacement surgery), we conducted a linkage analysis of chromosome 2q using 16 polymorphic microsatellite markers at an average spacing of one marker every 8.5 cM.

RESULTS

Our results provide suggestive evidence for a locus at 2q31 with a maximum multipoint logarithm of the odds score (MLS) of 1.22 which increased to 2.19 in those families concordant for hip-only disease (n = 311). This suggestive linkage was greater in male-hip families (MLS = 1.57, n = 71) than in female-hip families (MLS = 0.71, n = 132).

CONCLUSIONS

Chromosome 2q is likely to contain at least one susceptibility locus for OA.

Linkage analysis of candidate genes as susceptibility loci for osteoarthritis-suggestive linkage of COL9A1 to female hip osteoarthritis

OBJECTIVE

To examine 11 candidate genes as susceptibility loci for osteoarthritis (OA).

METHODS

A total of 481 families have been ascertained in which at least two siblings have had joint replacement surgery of the hip, or knee, or hip and knee for idiopathic OA. Each candidate gene was targeted using one or more intragenic or closely linked microsatellite marker. The linkage data were analysed unstratified and following stratification by sex and by joint replaced (hip or knee).

RESULTS

The analyses revealed suggestive linkage of the type IX collagen gene COL9A1 (6q12-q13) to a subset of 132 families that contained affected females who were concordant for hip OA (female-hip) with a P-value of 0.00053 and logarithm of the odds (LOD) score of 2.33 [corrected P-value of 0. 0016, corrected LOD score of 1.85].

CONCLUSIONS

COL9A1 may therefore be a susceptibility locus for female hip OA. In addition, there was weak evidence of linkage to HLA/COL11A2 (6p21.3) in female hip OA with a corrected P-value of 0.016.

Association analysis of the vitamin D receptor gene, the type I collagen gene COL1A1, and the estrogen receptor gene in idiopathic osteoarthritis

OBJECTIVE

Evidence has accumulated supporting a role for genes in the etiology of osteoarthritis (OA). Several candidates have been targeted as potential susceptibility loci including genes that are involved in the regulation of bone density. Genetic association analysis has suggested a role for the vitamin D receptor gene (VDR) and the estrogen receptor gene (ER) in susceptibility. Such findings must be tested in additional independent cohorts. We tested for association of these 2 genes, plus a third gene implicated in bone density, COL1A1, with idiopathic OA.

METHODS

A case-control cohort of 371 affected probands and 369 unaffected spouses was used. Association was tested using 4 intragenic single nucleotide polymorphisms (SNP), one each for the VDR and COL1A1 genes, and 2 for the ER gene. The VDR and ER SNP are the same SNP that have been associated with OA. All 4 SNP affect restriction enzyme sites and were genotyped using polymerase chain reaction and enzyme digestion. Allele and genotype distributions for each SNP were compared between cases and controls and analyzed using Fisher's exact test.

RESULTS

There was no evidence of association of the VDR or the ER gene SNP to OA. There was weak evidence of association of the COL1A1 SNP in female cases (p = 0.017), reflected by a difference in the distribution of genotypes at this SNP between female cases and controls (p = 0.027). However, when corrected for multiple testing, these results were not significant.

CONCLUSION

If the VDR, ER, or COL1A1 genes do encode predisposition to OA then the 4 SNP tested are not associated with major susceptibility alleles at these 3 loci.

Osteoarthritis-susceptibility locus on chromosome 11q, detected by linkage

We present a two-stage genomewide scan for osteoarthritis-susceptibility loci, using 481 families that each contain at least one affected sibling pair. The first stage, with 272 microsatellite markers and 297 families, involved a sparse map covering 23 chromosomes at intervals of approximately 15 cM. Sixteen markers that showed evidence of linkage at nominal P</=.05 were then taken through to the second stage, with an additional 184 families. This second stage confirmed evidence of linkage for markers on chromosome 11q. Additional markers from this region were then typed to create a denser map. We obtained a maximum single-point LOD score, at D11S901, of 2.40 (P=.0004) and a maximum multipoint-LOD score of 3.15, between markers D11S1358 and D11S35. A subset of 196 of the 481 families, comprising affected female sibling pairs, generated a corrected LOD score of 2.54 (P=.0003) for marker D11S901, with evidence for linkage extending 12 cM proximal to this marker. When we stratified for affected male sibling pairs there was no evidence of linkage to chromosome 11. Our data suggest that a female-specific susceptibility gene for idiopathic osteoarthritis is located on chromosome 11q.

Brachydactyly type B: linkage to chromosome 9q22 and evidence for genetic heterogeneity

Brachydactyly type B (BDB), an autosomal dominant disorder, is the most severe of the brachydactylies and is characterized by hypoplasia or absence of the terminal portions of the index to little fingers, usually with absence of the nails. The thumbs may be of normal length but are often flattened and occasionally are bifid. The feet are similarly but less severely affected. We have performed a genomewide linkage analysis of three families with BDB, two English and one Portugese. The two English families show linkage to the same region on chromosome 9 (combined multipoint maximum LOD score 8.69 with marker D9S257). The 16-cM disease interval is defined by recombinations with markers D9S1680 and D9S1786. These two families share an identical disease haplotype over 18 markers, inclusive of D9S278-D9S280. This provides strong evidence that the English families have the same ancestral mutation, which reduces the disease interval to <12.7 cM between markers D9S257 and D9S1851 in chromosome band 9q22. In the Portuguese family, we excluded linkage to this region, a result indicating that BDB is genetically heterogeneous. Reflecting this, there were atypical clinical features in this family, with shortening of the thumbs and absence or hypoplasia of the nails of the thumb and hallux. These results enable a refined classification of BDB and identify a novel locus for digit morphogenesis in 9q22.

European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA

We have used Y-chromosomal polymorphisms to trace paternal lineages in Polynesians by use of samples previously typed for mtDNA variants. A genealogical approach utilizing hierarchical analysis of eight rare-event biallelic polymorphisms, seven microsatellite loci, and internal structural analysis of the hypervariable minisatellite, MSY1, has been used to define three major paternal-lineage clusters in Polynesians. Two of these clusters, both defined by novel MSY1 modular structures and representing 55% of the Polynesians studied, are also found in coastal Papua New Guinea. Reduced Polynesian diversity, relative to that in Melanesians, is illustrated by the presence of several examples of identical MSY1 codes and microsatellite haplotypes within these lineage clusters in Polynesians. The complete lack of Y chromosomes having the M4 base substitution in Polynesians, despite their prevalence (64%) in Melanesians, may also be a result of the multiple bottleneck events during the colonization of this region of the world. The origin of the M4 mutation has been dated by use of two independent methods based on microsatellite-haplotype and minisatellite-code diversity. Because of the wide confidence limits on the mutation rates of these loci, the M4 mutation cannot be conclusively dated relative to the colonization of Polynesia, 3,000 years ago. The other major lineage cluster found in Polynesians, defined by a base substitution at the 92R7 locus, represents 27% of the Polynesians studied and, most probably, originates in Europe. This is the first Y-chromosomal evidence of major European admixture with indigenous Polynesian populations and contrasts sharply with the picture given by mtDNA evidence.

Critical ovarian hyperstimulation syndrome in a coasted in-vitro fertilization patient

We report an instance of critical ovarian hyperstimulation syndrome in a highly responsive in-vitro fertilization patient despite the preventive measure of a 4 day 'coast' interval during which no gonadotrophins were administered while gonadotrophin-releasing hormone agonist therapy continued until serum oestradiol concentrations fell below 3000 pg/ml.

Mutations in fibroblast growth-factor receptor 3 in sporadic cases of achondroplasia occur exclusively on the paternally derived chromosome

More than 97% of achondroplasia cases are caused by one of two mutations (G1138A and G1138C) in the fibroblast growth factor receptor 3 (FGFR3) gene, which results in a specific amino acid substitution, G380R. Sporadic cases of achondroplasia have been associated with advanced paternal age, suggesting that these mutations occur preferentially during spermatogenesis. We have determined the parental origin of the achondroplasia mutation in 40 sporadic cases. Three distinct 1-bp polymorphisms were identified in the FGFR3 gene, within close proximity to the achondroplasia mutation site. Ninety-nine families, each with a sporadic case of achondroplasia in a child, were analyzed in this study. In this population, the achondroplasia mutation occurred on the paternal chromosome in all 40 cases in which parental origin was unambiguous. This observation is consistent with the clinical observation of advanced paternal age resulting in new cases of achondroplasia and suggests that factors influencing DNA replication or repair during spermatogenesis, but not during oogenesis, may predispose to the occurrence of the G1138 FGFR3 mutations.

Genetic influences in end-stage osteoarthritis. Sibling risks of hip and knee replacement for idiopathic osteoarthritis

From a prospective, cross-sectional survey of 402 patients who had a total hip (THR) or a total knee (TKR) replacement for idiopathic osteoarthritis (OA) at a major centre, we determined the prevalence of these replacements for idiopathic OA in their 1171 siblings and 376 spouses. Using spouses as controls, the relative risk of THR in siblings was 1.86 (95% CI 0.93 to 3.69). The relative risk for TKR in siblings v spouses was 4.8 (95 % CI 0.64 to 36.4) whereas the risk for the combined outcome measure of THR or TKR was 2.32 (95% CI 1.22 to 4.43) when siblings and spouses over 64 years of age were compared. Using a threshold liability model (Falconer), the heritability of end-stage OA of the hip was estimated at 27%. The increased risks of joint replacement for severe, idiopathic OA which we found in siblings suggest that genetic influences are important in end-stage OA of the hip and knee.

GENETIC INFLUENCES IN END-STAGE OSTEOARTHRITIS

From a prospective, cross-sectional survey of 402 patients who had a total hip (THR) or a total knee (TKR) replacement for idiopathic osteoarthritis (OA) at a major centre, we determined the prevalence of these replacements for idiopathic OA in their 1171 siblings and 376 spouses. Using spouses as controls, the relative risk of THR in siblings was 1.86 (95% CI 0.93 to 3.69). The relative risk for TKR in siblings v spouses was 4.8 (95% CI 0.64 to 36.4) whereas the risk for the combined outcome measure of THR or TKR was 2.32 (95% CI 1.22 to 4.43) when siblings and spouses over 64 years of age were compared. Using a threshold liability model (Falconer), the heritability of end-stage OA of the hip was estimated at 27%.

The increased risks of joint replacement for severe, idiopathic OA which we found in siblings suggest that genetic influences are important in end-stage OA of the hip and knee.

Women's health services: restructuring for Medicaid managed care

New Jersey health care providers face the need to change dramatically the way health care is delivered as it enters a new era of managed care. This year, more than 24% of New Jersey's total population is enrolled in commercial managed care plans (New Jersey Department of Insurance, 1996). In addition, the state's Medicaid agency took steps to improve the delivery of health services to recipients by initiating implementation activities to transition from the traditional Medicaid program to a managed care model. Eighty-two percent of New Jersey's Aid to Families with Dependent Children (AFDC) and related populations have already been enrolled in managed care. The state plans to expand enrollment in managed care to the remaining 400,000 Medicaid beneficiaries. Communities with high Medicaid populations are challenged with the need to move through the managed care evolution at an accelerated rate.

Changes in protein and mRNA levels of growth factor/growth factor receptors in rat livers after administration of phenobarbitone or methylclofenapate

The effects of phenobarbitone and methylclofenapate were studied on the expression of growth factor and growth factor receptors in livers of male Wistar rats. The major findings were: (1) a significant reduction in epidermal growth factor receptor (EGFR) protein observed with both treatments, and (2) levels of EGFR transcripts were only slightly decreased with both compounds. The reduction in the receptor level therefore does not occur via regulation of transcription. Mannose-6-phosphate receptors (M6PR, also called insulin-like growth factor II receptor) and M6PR transcripts remained unchanged in both experimental groups. Hepatocyte growth factor receptor (HGFR) transcripts were also unchanged in both experimental groups. Transcript levels of transforming growth factor-beta 1 (TGF-beta 1) were lower in both treatment groups compared with the control; the reduction was significant in the methylclofenapate group. This may have relevance to the finding by others that nafenopin, another peroxisome proliferator, suppresses rat hepatocyte apoptosis. Another finding of general interest was that the three "housekeeping genes", namely albumin, actin and glyceraldehyde-3-phosphate dehydrogenase, were influenced by both treatments thus limiting their use as controls for gel loading. The adaptation of a growth regulatory mechanism via EGFR and its ligands may provide conditions such that cells with aberrant growth control have a selective growth advantage over normal cells thus promoting tumorigenesis.