Cross-talk of subchondral bone osteoblasts and articular cartilage chondrocytes : a new insight in understanding osteoarthritis pathogenesis

Osteoarthritis (OA) is the most common musculoskeletal disorder and represents a major health burden to society. In the course of the pathological development of OA, articular cartilage chondrocytes (ACCs) undergo atypical phenotype changes characterized by the expression of hypertrophic differentiation markers. Also, the adjacent subchondral bone shows signs of abnormal mineral density and enhanced production of bone turnover markers, indicative of osteoblast dysfunction. Collectively these findings indicate that the pathological changes typical of OA, involve alterations of the phenotypic properties of cells in both the subchondral bone and articular cartilage. However, the mechanism(s) by which these changes occur during OA development are not completely understood. The purpose of this project was to address the question of how subchondral bone osteoblasts (SBOs) and ACCs interact with each other with respect to regulation of respective cells’ phenotypic properties and in particular the involvement of mitogen activated protein kinase (MAPK) signalling pathways under normal and OA joint condition. We also endeavoured to test the influence of cross-talk between SBOs and ACCs isolated from normal and OA joint on matrix metalloproteinase (MMP) expression. For this purpose tissues from the knees of OA patients and normal controls were collected to isolate SBOs and ACCs. The cellular cross-talk of SBOs and ACCs were studied by means of both direct and indirect co-culture systems, which made it possible to identify the role of both membrane bound and soluble factors. Histology, immunohistochemistry, qRT-PCR, zymography, ELISA and western blotting were some of the techniques applied to distinguish the changes in the co-cultured vs. non co-cultured cells. The MAPK signalling pathways were probed by using targeted MAPK inhibitors, and their activity monitored by western blot analysis using phospho MAPK specific antibodies. Our co-culture studies demonstrated that OA ACCs enhanced the SBOs differentiation compared to normal ACCs. We demonstrated that OA ACCs induced these phenotypic changes in the SBOs via activating an ERK1/2 signalling pathway. The findings from this study thus provided clear evidence that OA ACCs play an integral role in altering the SBO phenotype. In the second study, we tested the influence of normal SBOs and OA SBOs on ACCs phenotype changes. The results showed that OA SBOs increased the hypertrophic gene expression in co-cultured ACCs compared to normal SBOs, a phenotype which is considered as pathological to the health and integrity of articular cartilage. It was demonstrated that these phenotype changes occurred via de-activation of p38 and activation of ERK1/2 signaling pathways. These findings suggest that the pathological interaction of OA SBOs with ACCs is mediated by cross-talking between ERK1/2 and p38 pathways, resulting in ACCs undergoing hypertrophic differentiation. Subsequent experiments to determine the effect on MMP regulation, of SBOs and ACCs cross-talk, revealed that co-culturing OA SBOs with ACCs significantly enhanced the proteolytic activity and expression of MMP-2 and MMP-9. In turn, co-culture of OA ACCs with SBOs led to abundant MMP-2 expression in SBOs. Furthermore, we showed that the addition of ERK1/2 and JNK inhibitors reversed the elevated MMP-2 and MMP-9 production which otherwise resulted from the interactions of OA SBOs-ACCs. Thus, this study has demonstrated that the altered interactions between OA SBOs-ACCs are capable of triggering the pathological pathways leading to degenerative changes seen in the osteoarthritic joint. In conclusion, the body of work presented in this dissertation has given clear in vitro evidence that the altered bi-directional communication of SBOs and ACCs may play a role in OA development and that this process was mediated by MAPK signalling pathways. Targeting these altered interactions by the use of MAPK inhibitors may provide the scientific rationale for the development of novel therapeutic strategies in the treatment and management of OA.

Optimal fluid injection strategies for in situ mineral leaching in two-dimensions

This paper examines the ground-water flow problem associated with the injection and recovery of certain corrosive fluids into mineral bearing rock. The aim is to dissolve the minerals in situ, and then recover them in solution. In general, it is not possible to recover all the injected fluid, which is of concern economically and environmentally. However, a new strategy is proposed here, that allows all the leaching fluid to be recovered. A mathematical model of the situation is solved approximately using an asymptotic solution, and exactly using a boundary integral approach. Solutions are shown for two-dimensional flow, which is of some practical interest as it is achievable in old mine tunnels, for example.

Vibrational spectroscopic analysis of the mineral crandallite CaAl3(PO4)2(OH)5•(H2O) from the Jenolan Caves, Australia

The mineral crandallite CaAl3(PO4)2(OH)5•(H2O) has been identified in deposits found in the Jenolan Caves, New South Wales, Australia by using a combination of X-ray diffraction and Raman spectroscopic techniques. A comparison is made between the vibrational spectra of crandallite found in the Jenolan Caves and a standard crandallite. Raman and infrared bands are assigned to PO43- and HPO42- stretching and bending modes. The predominant features are the internal vibrations of the PO43 and HPO42- groups. A mechanism for the formation of crandallite is presented and the conditions for the formation are elucidated.

Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data

The finite element (FE) analysis is an effective method to study the strength and predict the fracture risk of endodontically-treated teeth. This paper presents a rapid method developed to generate a comprehensive tooth FE model using data retrieved from micro-computed tomography (μCT). With this method, the inhomogeneity of material properties of teeth was included into the model without dividing the tooth model into different regions. The material properties of the tooth were assumed to be related to the mineral density. The fracture risk at different tooth portions was assessed for root canal treatments. The micro-CT images of a tooth were processed by a Matlab software programme and the CT numbers were retrieved. The tooth contours were obtained with thresholding segmentation using Amira. The inner and outer surfaces of the tooth were imported into Solidworks and a three-dimensional (3D) tooth model was constructed. An assembly of the tooth model with the periodontal ligament (PDL) layer and surrounding bone was imported into ABAQUS. The material properties of the tooth were calculated from the retrieved CT numbers via ABAQUS user's subroutines. Three root canal geometries (original and two enlargements) were investigated. The proposed method in this study can generate detailed 3D finite element models of a tooth with different root canal enlargements and filling materials, and would be very useful for the assessment of the fracture risk at different tooth portions after root canal treatments.

Nature and genesis of clay minerals of the Rustler Formation in the vicinity of the Waste Isolation Pilot Plant in Souteastern New Mexico

Detailed mineralogical studies of the matrix and fracture-fill materials of a large number of samples from the Rustler Formation have been carried out using x-ray diffraction, high-resolution transmission electron microscopy, electron microprobe analysis, x-ray fluorescence, and atomic absorption spectrophotometry. These analyses indicate the presence of four clay minerals: interstratified chlorite/saponite, illite, chlorite, and serpentine. Corrensite (regularly stratified chlorite/saponite) is the dominant clay mineral in samples from the Culebra dolomite and two shale layers of the lower unnamed member of the Rustler Formation. Within other layers of the Rustler Formation, disordered mixed chlorite/saponite is usually the most abundant clay mineral. Studies of the morphology and composition of clay crystallites suggest that the corrensite was formed by the alteration of detrital dioctahedral smectite in magnesium-rich pore fluids during early diagenesis of the Rustler Formation. This study provides initial estimates of the abundance and nature of the clay minerals in the Culebra dolomite in the vicinity of the Waste Isolation Pilot Plant.

Vibrational spectroscopy of the borate mineral kotoite Mg3(BO3)2

Vibrational spectroscopy has been used to assess the structure of kotoite a borate mineral of magnesium which is isostructural with jimboite. The mineral is orthorhombic with point group: 2/m 2/m 2/m. The mineral has the potential as a new memory insulator material. The mineral has been characterised by a combination of Raman and infrared spectroscopy. The Raman spectrum is dominated by a very intense band at 835 cm−1, assigned to the symmetric stretching mode of tetrahedral boron. Raman bands at 919, 985 and 1015 cm−1 are attributed to the antisymmetric stretching modes of tetrahedral boron. Kotoite is strictly an hydrous borate mineral. An intense Raman band observed at 3559 cm−1 is attributed to the stretching vibration of hydroxyl units, more likely to be associated with the borate mineral hydroxyborate. The lack of observation of water bending modes proves the absence of water in the kotoite structure.

The phosphate mineral arrojadite-(KFe) and its spectroscopic characterization

The arrojadite-(KFe) mineral has been analyzed using a combination of scanning electron microscopy and a combination of Raman and infrared spectroscopy. The origin of the mineral is Rapid Creek sedimentary phosphatic iron formation, northern Yukon. The formula of the mineral was determined as K2.06Na2Ca0.89Na3.23(Fe7.82Mg4.40Mn0.78)Σ13.00Al1.44(PO4)10.85(PO3OH0.23)(OH)2. The complexity of the mineral formula is reflected in the spectroscopy. Raman bands at 975, 991 and 1005 cm−1 with shoulder bands at 951 and 1024 cm−1 are assigned to the View the MathML source ν1 symmetric stretching modes. The Raman bands at 1024, 1066, 1092, 1123, 1148 and 1187 cm−1 are assigned to the View the MathML source ν3 antisymmetric stretching modes. A series of Raman bands observed at 540, 548, 557, 583, 604, 615 and 638 cm−1 are attributed to the ν4 out of plane bending modes of the PO4 and H2PO4 units. The ν2 PO4 and H2PO4 bending modes are observed at 403, 424, 449, 463, 479 and 513 cm−1. Hydroxyl and water stretching bands are readily observed. Vibrational spectroscopy enables new information about the complex phosphate mineral arrojadite-(KFe) to be obtained.

Infrared and Raman spectroscopic characterization of the carbonate mineral huanghoite - and in comparison with selected rare earth carbonates

Raman spectroscopy complimented with infrared spectroscopy has been used to study the rare earth based mineral huanghoite with possible formula given as BaCe(CO3)2F and compared with the Raman spectra of a series of selected natural halogenated carbonates from different origins including bastnasite, parisite and northupite. The Raman spectrum of huanghoite displays three bands are at 1072, 1084 and 1091 cm−1 attributed to the symmetric stretching vibration. The observation of three symmetric stretching vibrations is very unusual. The position of symmetric stretching vibration varies with mineral composition. Infrared spectroscopy of huanghoite show bands at 1319, 1382, 1422 and 1470 cm−1. No Raman bands of huanghoite were observed in these positions. Raman spectra of bastnasite, parisite and northupite show a single band at 1433, 1420 and 1554 cm−1 assigned to the ν3 (CO3)2− antisymmetric stretching mode. The observation of additional Raman bands for the ν3 modes for some halogenated carbonates is significant in that it shows distortion of the carbonate anion in the mineral structure. Four Raman bands for huanghoite are observed at 687, 704, 718 and 730 cm−1and assigned to the (CO3)2− ν2 bending modes. Raman bands are observed for huanghoite at around 627 cm−1 and are assigned to the (CO3)2− ν4 bending modes. Raman bands are observed for the carbonate ν4 in phase bending modes at 722 cm−1 for bastnasite, 736 and 684 cm−1 for parisite, 714 cm−1 for northupite. Raman bands for huanghoite observed at 3259, 3484 and 3589 cm−1 are attributed to water stretching bands. Multiple bands are observed in the OH stretching region for bastnasite and parisite indicating the presence of water and OH units in their mineral structure. Vibrational spectroscopy enables new information on the structure of huanghoite to be assessed.

A new approach to spatial data interpolation using higher-order statistics

Interpolation techniques for spatial data have been applied frequently in various fields of geosciences. Although most conventional interpolation methods assume that it is sufficient to use first- and second-order statistics to characterize random fields, researchers have now realized that these methods cannot always provide reliable interpolation results, since geological and environmental phenomena tend to be very complex, presenting non-Gaussian distribution and/or non-linear inter-variable relationship. This paper proposes a new approach to the interpolation of spatial data, which can be applied with great flexibility. Suitable cross-variable higher-order spatial statistics are developed to measure the spatial relationship between the random variable at an unsampled location and those in its neighbourhood. Given the computed cross-variable higher-order spatial statistics, the conditional probability density function (CPDF) is approximated via polynomial expansions, which is then utilized to determine the interpolated value at the unsampled location as an expectation. In addition, the uncertainty associated with the interpolation is quantified by constructing prediction intervals of interpolated values. The proposed method is applied to a mineral deposit dataset, and the results demonstrate that it outperforms kriging methods in uncertainty quantification. The introduction of the cross-variable higher-order spatial statistics noticeably improves the quality of the interpolation since it enriches the information that can be extracted from the observed data, and this benefit is substantial when working with data that are sparse or have non-trivial dependence structures.

The global burden attributable to low bone mineral density

Introduction The Global Burden of Disease Study 2010 estimated the worldwide health burden of 291 diseases and injuries and 67 risk factors by calculating disability-adjusted life years (DALYs). Osteoporosis was not considered as a disease, and bone mineral density (BMD) was analysed as a risk factor for fractures, which formed part of the health burden due to falls. Objectives To calculate (1) the global distribution of BMD, (2) its population attributable fraction (PAF) for fractures and subsequently for falls, and (3) the number of DALYs due to BMD. Methods A systematic review was performed seeking population-based studies in which BMD was measured by dual-energy X-ray absorptiometry at the femoral neck in people aged 50 years and over. Age- and sex-specific mean ± SD BMD values (g/cm2) were extracted from eligible studies. Comparative risk assessment methodology was used to calculate PAFs of BMD for fractures. The theoretical minimum risk exposure distribution was estimated as the age- and sex-specific 90th centile from the Third National Health and Nutrition Examination Survey (NHANES III). Relative risks of fractures were obtained from a previous meta-analysis. Hospital data were used to calculate the fraction of the health burden of falls that was due to fractures. Results Global deaths and DALYs attributable to low BMD increased from 103 000 and 3 125 000 in 1990 to 188 000 and 5 216 000 in 2010, respectively. The percentage of low BMD in the total global burden almost doubled from 1990 (0.12%) to 2010 (0.21%). Around one-third of falls-related deaths were attributable to low BMD. Conclusions Low BMD is responsible for a growing global health burden, only partially representative of the real burden of osteoporosis.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

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.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.