The IFITM5 mutation c.-14C > T results in an elongated transcript expressed in human bone; And causes varying phenotypic severity of osteogenesis imperfecta type V

Background The genetic mutation resulting in osteogenesis imperfecta (OI) type V was recently characterised as a single point mutation (c.-14C > T) in the 5' untranslated region (UTR) of IFITM5, a gene encoding a transmembrane protein with expression restricted to skeletal tissue. This mutation creates an alternative start codon and has been shown in a eukaryotic cell line to result in a longer variant of IFITM5, but its expression has not previously been demonstrated in bone from a patient with OI type V. Methods Sanger sequencing of the IFITM5 5' UTR was performed in our cohort of subjects with a clinical diagnosis of OI type V. Clinical data was collated from referring clinicians. RNA was extracted from a bone sample from one patient and Sanger sequenced to determine expression of wild-type and mutant IFITM5. Results: All nine subjects with OI type V were heterozygous for the c.-14C > T IFITM5 mutation. Clinically, there was heterogeneity in phenotype, particularly in the manifestation of bone fragility amongst subjects. Both wild-type and mutant IFITM5 mRNA transcripts were present in bone. Conclusions The c.-14C > T IFITM5 mutation does not result in an RNA-null allele but is expressed in bone. Individuals with identical mutations in IFITM5 have highly variable phenotypic expression, even within the same family.

Association study of genes related to bone formation and resorption and the extent of radiographic change in ankylosing spondylitis

Objective: To identify genetic associations with severity of radiographic damage in ankylosing spondylitis (AS). Method: We studied 1537 AS cases of European descent; all fulfilled the modified New York Criteria. Radiographic severity was assessed from digitised lateral radiographs of the cervical and lumbar spine using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). A two-phase genotyping design was used. In phase 1, 498 single nucleotide polymorphisms (SNPs) were genotyped in 688 cases; these were selected to capture >90% of the common haplotypic variation in the exons, exon-intron boundaries, and 5 kb flanking DNA in the 5' and 3' UTR of 74 genes involved in anabolic or catabolic bone pathways. In phase 2, 15 SNPs exhibiting p<0.05 were genotyped in a further cohort of 830 AS cases; results were analysed both separately and in combination with the discovery phase data. Association was tested by contingency tables after separating the samples into 'mild' and 'severe' groups, defined as the bottom and top 40% by mSASSS, adjusted for gender and disease duration. Results: Experiment-wise association was observed with the SNP rs8092336 (combined OR 0.32, p=1.2×10-5), which lies within RANK (receptor activator of NF?B), a gene involved in osteoclastogenesis, and in the interaction between T cells and dendritic cells. Association was also found with the SNP rs1236913 in PTGS1 (prostaglandin-endoperoxide synthase 1, cyclooxygenase 1), giving an OR of 0.53 (p=2.6×10-3). There was no observed association between radiographic severity and HLA-B*27. Conclusions: These findings support roles for bone resorption and prostaglandins pathways in the osteoproliferative changes in AS.

Excessive bone formation in a mouse model of ankylosing spondylitis is associated with decreases in Wnt pathway inhibitors

Introduction: Ankylosing spondylitis (AS) is unique in its pathology where inflammation commences at the entheses before progressing to an osteoproliferative phenotype generating excessive bone formation that can result in joint fusion. The underlying mechanisms of this progression are poorly understood. Recent work has suggested that changes in Wnt signalling, a key bone regulatory pathway, may contribute to joint ankylosis in AS. Using the proteoglycan-induced spondylitis (PGISp) mouse model which displays spondylitis and eventual joint fusion following an initial inflammatory stimulus, we have characterised the structural and molecular changes that underlie disease progression. Methods: PGISp mice were characterised 12 weeks after initiation of inflammation using histology, immunohistochemistry (IHC) and expression profiling. Results: Inflammation initiated at the periphery of the intervertebral discs progressing to disc destruction followed by massively excessive cartilage and bone matrix formation, as demonstrated by toluidine blue staining and IHC for collagen type I and osteocalcin, leading to syndesmophyte formation. Expression levels of DKK1 and SOST, Wnt signalling inhibitors highly expressed in joints, were reduced by 49% and 63% respectively in the spine PGISp compared with control mice (P < 0.05) with SOST inhibition confirmed by IHC. Microarray profiling showed genes involved in inflammation and immune-regulation were altered. Further, a number of genes specifically involved in bone regulation including other members of the Wnt pathway were also dysregulated. Conclusions: This study implicates the Wnt pathway as a likely mediator of the mechanism by which inflammation induces bony ankylosis in spondyloarthritis, raising the potential that therapies targeting this pathway may be effective in preventing this process.

Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

COL1A1 C-propeptide cleavage site mutation causes high bone mass, bone fragility and jaw lesions: A new cause of gnathodiaphyseal dysplasia?

Gnathodiaphyseal dysplasia (GDD) is a rare autosomal dominant condition characterized by bone fragility, irregular bone mineral density (BMD) and fibro-osseous lesions in the skull and jaw. Mutations in Anoctamin-5 (ANO5) have been identified in some cases. We aimed to identify the causative mutation in a family with features of GDD but no mutation in ANO5, using whole exome capture and massive parallel sequencing (WES). WES of two affected individuals (a mother and son) and the mother's unaffected parents identified a mutation in the C-propeptide cleavage site of COL1A1. Similar mutations have been reported in individuals with osteogenesis imperfecta (OI) and paradoxically increased BMD. C-propeptide cleavage site mutations in COL1A1 may not only cause 'high bone mass OI', but also the clinical features of GDD, specifically irregular sclerotic BMD and fibro-osseous lesions in the skull and jaw. GDD patients negative for ANO5 mutations should be assessed for mutations in type I collagen C-propeptide cleavage sites.

Mimicking breast cancer-induced bone metastasis in vivo: Current transplantation models and advanced humanized strategies

Bone metastasis is a complication that occurs in 80 % of women with advanced breast cancer. Despite the prevalence of bone metastatic disease, the avenues for its clinical management are still restricted to palliative treatment options. In fact, the underlying mechanisms of breast cancer osteotropism have not yet been fully elucidated due to a lack of suitable in vivo models that are able to recapitulate the human disease. In this work, we review the current transplantation-based models to investigate breast cancer-induced bone metastasis and delineate the strengths and limitations of the use of different grafting techniques, tissue sources, and hosts. We further show that humanized xenograft models incorporating human cells or tissue grafts at the primary tumor site or the metastatic site mimic more closely the human disease. Tissue-engineered constructs are emerging as a reproducible alternative to recapitulate functional humanized tissues in these murine models. The development of advanced humanized animal models may provide better platforms to investigate the mutual interactions between human cancer cells and their microenvironment and ultimately improve the translation of preclinical drug trials to the clinic.

Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

The revolution in human monogenic disease mapping

The successful completion of the Human Genome Project (HGP) was an unprecedented scientific advance that has become an invaluable resource in the search for genes that cause monogenic and common (polygenic) diseases. Prior to the HGP, linkage analysis had successfully mapped many disease genes for monogenic disorders; however, the limitations of this approach were particularly evident for identifying causative genes in rare genetic disorders affecting lifespan and/or reproductive fitness, such as skeletal dysplasias. In this review, we illustrate the challenges of mapping disease genes in such conditions through the ultra-rare disorder fibrodysplasia ossificans progressiva (FOP) and we discuss the advances that are being made through current massively parallel (“next generation”) sequencing (MPS) technologies.

The role of IL-17-secreting mast cells in inflammatory joint disease

The proinflammatory cytokine IL-17 has an important role in pathogenesis of several inflammatory diseases. In immune-mediated joint diseases, IL-17 can induce secretion of other proinflammatory cytokines such as IL-1, IL-6 and TNF, as well as matrix metalloproteinase enzymes, leading to inflammation, cartilage breakdown, osteoclastogenesis and bone erosion. In animal models of inflammatory arthritis, mice deficient in IL-17 are less susceptible to development of disease. The list of IL-17-secreting cells is rapidly growing, and mast cells have been suggested to be a dominant source of IL-17 in inflammatory joint disease. However, many other innate sources of IL-17 have been described in both inflammatory and autoinflammatory conditions, raising questions as to the role of mast cells in orchestrating joint inflammation. This article will critically assess the contribution of mast cells and other cell types to IL-17 production in the inflammatory milieu associated with inflammatory arthritis, understanding of which could facilitate targeted therapeutic approaches. © 2013 Macmillan Publishers Limited. All rights reserved.

Genetic control of bone density and turnover: Role of the collagen 1α1, estrogen receptor, and vitamin D receptor genes

Genetic factors are known to influence both the peak bone mass and probably the rate of change in bone density. A range of regulatory and structural genes has been proposed to be involved including collagen 1α1 (COL1A1), the estrogen receptor (ER), and the vitamin D receptor (VDR), but the actual genes involved are uncertain. We therefore studied the role of the COL1A1 and VDR loci in control of bone density by linkage in 45 dizygotic twin pairs and 29 nuclear families comprising 120 individuals. The influences on bone density of polymorphisms of COL1A1, VDR, and ER were studied by association both cross-sectionally and longitudinally in 193 elderly postmenopausal women (average age, 69 years) over a mean follow-up time of 6.3 years. Weak linkage of the COL1A1 locus with bone density was observed in both twins and families (p = 0.02 in both data sets), confirming previous observations of linkage of this locus with bone density. Association between the MscI polymorphism of COL1A1 and rate of lumbar spine bone loss was observed with significant gene-environment interaction related to dietary calcium intake (p = 0.0006). In the lowest tertile of dietary calcium intake, carriers of "s" alleles lost more bone than "SS" homozygotes (p = 0.01), whereas the opposite was observed in the highest dietary calcium intake (p = 0.003). Association also was observed between rate of bone loss at both the femoral neck and the lumbar spine and the TaqI VDR polymorphism (p = 0.03). This association was strongest in those in the lowest tertile of calcium intake, also suggesting the presence of gene-environment interaction involving dietary calcium and VDR, influencing bone turnover. No significant association was observed between the PvuII ER polymorphism alone or in combination with VDR or COL1A1 genotypes, with either bone density or its rate of change. These data support the involvement of COL1A1 in determination of bone density and the interaction of both COL1A1 and VDR with calcium intake in regulation of change of bone density over time.

Inflammation-driven bone formation in a mouse model of ankylosing spondylitis: Sequential not parallel processes

Background Ankylosing spondylitis (AS) is an immune-mediated arthritis particularly targeting the spine and pelvis and is characterised by inflammation, osteoproliferation and frequently ankylosis. Current treatments that predominately target inflammatory pathways have disappointing efficacy in slowing disease progression. Thus, a better understanding of the causal association and pathological progression from inflammation to bone formation, particularly whether inflammation directly initiates osteoproliferation, is required. Methods The proteoglycan-induced spondylitis (PGISp) mouse model of AS was used to histopathologically map the progressive axial disease events, assess molecular changes during disease progression and define disease progression using unbiased clustering of semi-quantitative histology. PGISp mice were followed over a 24-week time course. Spinal disease was assessed using a novel semi-quantitative histological scoring system that independently evaluated the breadth of pathological features associated with PGISp axial disease, including inflammation, joint destruction and excessive tissue formation (osteoproliferation). Matrix components were identified using immunohistochemistry. Results Disease initiated with inflammation at the periphery of the intervertebral disc (IVD) adjacent to the longitudinal ligament, reminiscent of enthesitis, and was associated with upregulated tumor necrosis factor and metalloproteinases. After a lag phase, established inflammation was temporospatially associated with destruction of IVDs, cartilage and bone. At later time points, advanced disease was characterised by substantially reduced inflammation, excessive tissue formation and ectopic chondrocyte expansion. These distinct features differentiated affected mice into early, intermediate and advanced disease stages. Excessive tissue formation was observed in vertebral joints only if the IVD was destroyed as a consequence of the early inflammation. Ectopic excessive tissue was predominantly chondroidal with chondrocyte-like cells embedded within collagen type II- and X-rich matrix. This corresponded with upregulation of mRNA for cartilage markers Col2a1, sox9 and Comp. Osteophytes, though infrequent, were more prevalent in later disease. Conclusions The inflammation-driven IVD destruction was shown to be a prerequisite for axial disease progression to osteoproliferation in the PGISp mouse. Osteoproliferation led to vertebral body deformity and fusion but was never seen concurrent with persistent inflammation, suggesting a sequential process. The findings support that early intervention with anti-inflammatory therapies will be needed to limit destructive processes and consequently prevent progression of AS.

Mutations in known monogenic high bone mass loci only explain a small proportion of high bone mass cases

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM ( approximately prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion ( approximately 3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.

The role of 25-hydroxyvitamin D deficiency in promoting insulin resistance and inflammation in patients with chronic kidney disease: A randomised controlled trial

BACKGROUND Approximately 50% of patients with stage 3 Chronic Kidney Disease are 25-hydroxyvitamin D insufficient, and this prevalence increases with falling glomerular filtration rate. Vitamin D is now recognised as having pleiotropic roles beyond bone and mineral homeostasis, with the vitamin D receptor and metabolising machinery identified in multiple tissues. Worryingly, recent observational data has highlighted an association between hypovitaminosis D and increased cardiovascular mortality, possibly mediated via vitamin D effects on insulin resistance and inflammation. The main hypothesis of this study is that oral Vitamin D supplementation will ameliorate insulin resistance in patients with Chronic Kidney Disease stage 3 when compared to placebo. Secondary hypotheses will test whether this is associated with decreased inflammation and bone/adipocyte-endocrine dysregulation. METHODS/DESIGN This study is a single-centre, double-blinded, randomised, placebo-controlled trial. Inclusion criteria include; estimated glomerular filtration rate 30-59 ml/min/1.73 m(2); aged >or=18 on entry to study; and serum 25-hydroxyvitamin D levels <75 nmol/L. Patients will be randomised 1:1 to receive either oral cholecalciferol 2000IU/day or placebo for 6 months. The primary outcome will be an improvement in insulin sensitivity, measured by hyperinsulinaemic euglycaemic clamp. Secondary outcome measures will include serum parathyroid hormone, cytokines (Interleukin-1beta, Interleukin-6, Tumour Necrosis Factor alpha), adiponectin (total and High Molecular Weight), osteocalcin (carboxylated and under-carboxylated), peripheral blood mononuclear cell Nuclear Factor Kappa-B p65 binding activity, brachial artery reactivity, aortic pulse wave velocity and waveform analysis, and indirect calorimetry. All outcome measures will be performed at baseline and end of study. DISCUSSION To date, no randomised controlled trial has been performed in pre-dialysis CKD patients to study the correlation between vitamin D status with supplementation, insulin resistance and markers of adverse cardiovascular risk. We remain hopeful that cholecalciferol may be a safe intervention, with health benefits beyond those related to bone-mineral homeostasis. TRIAL REGISTRATION Australian and New Zealand Clinical Trials Registry ACTRN12609000246280.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.