Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel candidate effector genes.

OBJECTIVES: Osteoarthritis is a complex disease with a huge public health burden. Genome-wide association studies (GWAS) have identified hundreds of osteoarthritis-associated sequence variants, but the effector genes underpinning these signals remain largely elusive. Understanding chromosome organisation in three-dimensional (3D) space is essential for identifying long-range contacts between distant genomic features (e.g., between genes and regulatory elements), in a tissue-specific manner. Here, we generate the first whole genome chromosome conformation analysis (Hi-C) map of primary osteoarthritis chondrocytes and identify novel candidate effector genes for the disease. METHODS: Primary chondrocytes collected from 8 patients with knee osteoarthritis underwent Hi-C analysis to link chromosomal structure to genomic sequence. The identified loops were then combined with osteoarthritis GWAS results and epigenomic data from primary knee osteoarthritis chondrocytes to identify variants involved in gene regulation via enhancer-promoter interactions. RESULTS: We identified 345 genetic variants residing within chromatin loop anchors that are associated with 77 osteoarthritis GWAS signals. Ten of these variants reside directly in enhancer regions of 10 newly described active enhancer-promoter loops, identified with multiomics analysis of publicly available chromatin immunoprecipitation sequencing (ChIP-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data from primary knee chondrocyte cells, pointing to two new candidate effector genes SPRY4 and PAPPA (pregnancy-associated plasma protein A) as well as further support for the gene SLC44A2 known to be involved in osteoarthritis. For example, PAPPA is directly associated with the turnover of insulin-like growth factor 1 (IGF-1) proteins, and IGF-1 is an important factor in the repair of damaged chondrocytes. CONCLUSIONS: We have constructed the first Hi-C map of primary human chondrocytes and have made it available as a resource for the scientific community. By integrating 3D genomics with large-scale genetic association and epigenetic data, we identify novel candidate effector genes for osteoarthritis, which enhance our understanding of disease and can serve as putative high-value novel drug targets.

Metagenomics of the Gut Microbiome in Parkinson's Disease: Prodromal Changes.

OBJECTIVE: Prior studies on the gut microbiome in Parkinson's disease (PD) have yielded conflicting results, and few studies have focused on prodromal (premotor) PD or used shotgun metagenomic profiling to assess microbial functional potential. We conducted a nested case-control study within 2 large epidemiological cohorts to examine the role of the gut microbiome in PD. METHODS: We profiled the fecal metagenomes of 420 participants in the Nurses' Health Study and the Health Professionals Follow-up Study with recent onset PD (N = 75), with features of prodromal PD (N = 101), controls with constipation (N = 113), and healthy controls (N = 131) to identify microbial taxonomic and functional features associated with PD and features suggestive of prodromal PD. Omnibus and feature-wise analyses identified bacterial species and pathways associated with prodromal and recently onset PD. RESULTS: We observed depletion of several strict anaerobes associated with reduced inflammation among participants with PD or features of prodromal PD. A microbiome-based classifier had moderate accuracy (area under the curve [AUC] = 0.76 for species and 0.74 for pathways) to discriminate between recently onset PD cases and controls. These taxonomic shifts corresponded with functional shifts indicative of carbohydrate source preference. Similar, but less marked, changes were observed in participants with features of prodromal PD, in both microbial features and functions. INTERPRETATION: PD and features of prodromal PD were associated with similar changes in the gut microbiome. These findings suggest that changes in the microbiome could represent novel biomarkers for the earliest phases of PD. ANN NEUROL 2023;94:486-501.

Transmission of Alzheimer's disease-associated microbiota dysbiosis and its impact on cognitive function: evidence from mice and patients.

Spouses of Alzheimer's disease (AD) patients are at a higher risk of developing incidental dementia. However, the causes and underlying mechanism of this clinical observation remain largely unknown. One possible explanation is linked to microbiota dysbiosis, a condition that has been associated with AD. However, it remains unclear whether gut microbiota dysbiosis can be transmitted from AD individuals to non-AD individuals and contribute to the development of AD pathogenesis and cognitive impairment. We, therefore, set out to perform both animal studies and clinical investigation by co-housing wild-type mice and AD transgenic mice, analyzing microbiota via 16S rRNA gene sequencing, measuring short-chain fatty acid amounts, and employing behavioral test, mass spectrometry, site-mutations and other methods. The present study revealed that co-housing between wild-type mice and AD transgenic mice or administrating feces of AD transgenic mice to wild-type mice resulted in AD-associated gut microbiota dysbiosis, Tau phosphorylation, and cognitive impairment in the wild-type mice. Gavage with Lactobacillus and Bifidobacterium restored these changes in the wild-type mice. The oral and gut microbiota of AD patient partners resembled that of AD patients but differed from healthy controls, indicating the transmission of microbiota. The underlying mechanism of these findings includes that the butyric acid-mediated acetylation of GSK3ß at lysine 15 regulated its phosphorylation at serine 9, consequently impacting Tau phosphorylation. Pending confirmative studies, these results provide insight into a potential link between the transmission of AD-associated microbiota dysbiosis and development of cognitive impairment, which underscore the need for further research in this area.

Interplay between diet and gut microbiome, and circulating concentrations of trimethylamine N-oxide: findings from a longitudinal cohort of US men.

OBJECTIVES: Gut-produced trimethylamine N-oxide (TMAO) is postulated as a possible link between red meat intake and poor cardiometabolic health. We investigated whether gut microbiome could modify associations of dietary precursors with TMAO concentrations and cardiometabolic risk markers among free-living individuals. DESIGN: We collected up to two pairs of faecal samples (n=925) and two blood samples (n=473), 6 months apart, from 307 healthy men in the Men's Lifestyle Validation Study. Diet was assessed repeatedly using food-frequency questionnaires and diet records. We profiled faecal metagenome and metatranscriptome using shotgun sequencing and identified microbial taxonomic and functional features. RESULTS: TMAO concentrations were associated with the overall microbial compositions (permutational analysis of variance (PERMANOVA) test p=0.001). Multivariable taxa-wide association analysis identified 10 bacterial species whose abundance was significantly associated with plasma TMAO concentrations (false discovery rate <0.05). Higher habitual intake of red meat and choline was significantly associated with higher TMAO concentrations among participants who were microbial TMAO-producers (p<0.05), as characterised based on four abundant TMAO-predicting species, but not among other participants (for red meat, P-interaction=0.003; for choline, P-interaction=0.03). Among abundant TMAO-predicting species, Alistipes shahii significantly strengthened the positive association between red meat intake and HbA1c levels (P-interaction=0.01). Secondary analyses revealed that some functional features, including choline trimethylamine-lyase activating enzymes, were associated with TMAO concentrations. CONCLUSION: We identified microbial taxa that were associated with TMAO concentrations and modified the associations of red meat intake with TMAO concentrations and cardiometabolic risk markers. Our data underscore the interplay between diet and gut microbiome in producing potentially bioactive metabolites that may modulate cardiometabolic health.

The oral microbiome in relation to pancreatic cancer risk in African Americans.

BACKGROUND: African Americans have the highest pancreatic cancer incidence of any racial/ethnic group in the United States. The oral microbiome was associated with pancreatic cancer risk in a recent study, but no such studies have been conducted in African Americans. Poor oral health, which can be a cause or effect of microbial populations, was associated with an increased risk of pancreatic cancer in a single study of African Americans. METHODS: We prospectively investigated the oral microbiome in relation to pancreatic cancer risk among 122 African-American pancreatic cancer cases and 354 controls. DNA was extracted from oral wash samples for metagenomic shotgun sequencing. Alpha and beta diversity of the microbial profiles were calculated. Multivariable conditional logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for associations between microbes and pancreatic cancer risk. RESULTS: No associations were observed with alpha or beta diversity, and no individual microbial taxa were differentially abundant between cases and control, after accounting for multiple comparisons. Among never smokers, there were elevated ORs for known oral pathogens: Porphyromonas gingivalis (OR = 1.69, 95% CI: 0.80-3.56), Prevotella intermedia (OR = 1.40, 95% CI: 0.69-2.85), and Tannerella forsythia (OR = 1.36, 95% CI: 0.66-2.77). CONCLUSIONS: Previously reported associations between oral taxa and pancreatic cancer were not present in this African-American population overall.

Dietary lignans, plasma enterolactone levels, and metabolic risk in men: exploring the role of the gut microbiome.

BACKGROUND: The conversion of plant lignans to bioactive enterolignans in the gastrointestinal tract is mediated through microbial processing. The goal of this study was to examine the relationships between lignan intake, plasma enterolactone concentrations, gut microbiome composition, and metabolic risk in free-living male adults. RESULTS: In 303 men participating in the Men's Lifestyle Validation Study (MLVS), lignan intake was assessed using two sets of 7-day diet records, and gut microbiome was profiled through shotgun sequencing of up to 2 pairs of fecal samples (n = 911). A score was calculated to summarize the abundance of bacteria species that were significantly associated with plasma enterolactone levels. Of the 138 filtered species, plasma enterolactone levels were significantly associated with the relative abundances of 18 species at FDR < 0.05 level. Per SD increment of lignan intake was associated with 20.7 nM (SEM: 2.3 nM) higher enterolactone concentrations among participants with a higher species score, whereas the corresponding estimate was 4.0 nM (SEM: 1.7 nM) among participants with a lower species score (P for interaction < 0.001). A total of 12 plasma metabolites were also significantly associated with these enterolactone-predicting species. Of the association between lignan intake and metabolic risk, 19.8% (95%CI: 7.3%-43.6%) was explained by the species score alone, 54.5% (95%CI: 21.8%-83.7%) by both species score and enterolactone levels, and 79.8% (95%CI: 17.7%-98.6%) by further considering the 12 plasma metabolites. CONCLUSION: We identified multiple gut bacteria species that were enriched or depleted at higher plasma levels of enterolactone in men. These species jointly modified the associations of lignan intake with plasma enterolactone levels and explained the majority of association between lignan intake and metabolic risk along with enterolactone levels and certain plasma metabolites.

Multivariable association discovery in population-scale meta-omics studies.

It is challenging to associate features such as human health outcomes, diet, environmental conditions, or other metadata to microbial community measurements, due in part to their quantitative properties. Microbiome multi-omics are typically noisy, sparse (zero-inflated), high-dimensional, extremely non-normal, and often in the form of count or compositional measurements. Here we introduce an optimized combination of novel and established methodology to assess multivariable association of microbial community features with complex metadata in population-scale observational studies. Our approach, MaAsLin 2 (Microbiome Multivariable Associations with Linear Models), uses generalized linear and mixed models to accommodate a wide variety of modern epidemiological studies, including cross-sectional and longitudinal designs, as well as a variety of data types (e.g., counts and relative abundances) with or without covariates and repeated measurements. To construct this method, we conducted a large-scale evaluation of a broad range of scenarios under which straightforward identification of meta-omics associations can be challenging. These simulation studies reveal that MaAsLin 2's linear model preserves statistical power in the presence of repeated measures and multiple covariates, while accounting for the nuances of meta-omics features and controlling false discovery. We also applied MaAsLin 2 to a microbial multi-omics dataset from the Integrative Human Microbiome (HMP2) project which, in addition to reproducing established results, revealed a unique, integrated landscape of inflammatory bowel diseases (IBD) across multiple time points and omics profiles.

Effects of Lifetime Exposures to Environmental Contaminants on the Adult Gut Microbiome.

Emerging experimental evidence indicates that toxicant-induced alterations in gut microbiota composition and activity may affect host homeostasis. However, data from human studies are scarce; to our knowledge, no previous studies have quantified the association of lifetime exposure to environmental chemicals, across multiple time points, with the composition of the adult gut microbiome. Here we studied 124 individuals born in the Faroe Islands in 1986-1987 who were followed approximately every seven years from birth through age 28 years. Organochlorine compounds, including polychlorinated biphenyls (PCBs) and pesticides, perfluoroalkyl substances (PFAS), and mercury (Hg), were measured in cord blood and longitudinally in participants' blood. At age 28, the gut microbiome was assessed using shotgun metagenomic sequencing. Historical contaminant exposures had little direct effect on the adult gut microbiome, while a small number of fastidious anaerobes were weakly linked to recent PFAS/PFOS exposures at age 28. In this cohort, our findings suggest no lasting effects of early life exposures on adult gut microbial composition, but proximal exposures may contribute to gut microbiome alterations. The methods developed and used for this investigation may help in future identification of small but lasting impacts of environmental toxicant exposure on the gut microbiome.

Pathophysiology and Emerging Molecular Therapeutic Targets in Heterotopic Ossification.

The term heterotopic ossification (HO) describes bone formation in tissues where bone is normally not present. Musculoskeletal trauma induces signalling events that in turn trigger cells, probably of mesenchymal origin, to differentiate into bone. The aetiology of HO includes extremely rare but severe, generalised and fatal monogenic forms of the disease; and as a common complex disorder in response to musculoskeletal, neurological or burn trauma. The resulting bone forms through a combination of endochondral and intramembranous ossification, depending on the aetiology, initiating stimulus and affected tissue. Given the heterogeneity of the disease, many cell types and biological pathways have been studied in efforts to find effective therapeutic strategies for the disorder. Cells of mesenchymal, haematopoietic and neuroectodermal lineages have all been implicated in the pathogenesis of HO, and the emerging dominant signalling pathways are thought to occur through the bone morphogenetic proteins (BMP), mammalian target of rapamycin (mTOR), and retinoic acid receptor pathways. Increased understanding of these disease mechanisms has resulted in the emergence of several novel investigational therapeutic avenues, including palovarotene and other retinoic acid receptor agonists and activin A inhibitors that target both canonical and non-canonical signalling downstream of the BMP type 1 receptor. In this article we aim to illustrate the key cellular and molecular mechanisms involved in the pathogenesis of HO and outline recent advances in emerging molecular therapies to treat and prevent HO that have had early success in the monogenic disease and are currently being explored in the common complex forms of HO.

The Kinesin Gene KIF26B Modulates the Severity of Post-Traumatic Heterotopic Ossification.

The formation of pathological bone deposits within soft tissues, termed heterotopic ossification (HO), is common after trauma. However, the severity of HO formation varies substantially between individuals, from relatively isolated small bone islands through to extensive soft tissue replacement by bone giving rise to debilitating symptoms. The aim of this study was to identify novel candidate therapeutic molecular targets for severe HO. We conducted a genome-wide scan in men and women with HO of varying severity following hip replacement for osteoarthritis. HO severity was dichotomized as mild or severe, and association analysis was performed with adjustment for age and sex. We next confirmed expression of the gene encoded by the lead signal in human bone and in primary human mesenchymal stem cells. We then examined the effect of gene knockout in a murine model of osseous trans-differentiation, and finally we explored transcription factor phosphorylation in key pathways perturbed by the gene. Ten independent signals were suggestively associated with HO severity, with KIF26B as the lead. We subsequently confirmed KIF26B expression in human bone and upregulation upon BMP2-induced osteogenic differentiation in primary human mesenchymal stem cells, and also in a rat tendo-Achilles model of post-traumatic HO. CRISPR-Cas9 mediated knockout of Kif26b inhibited BMP2-induced Runx2, Sp7/Osterix, Col1A1, Alp, and Bglap/Osteocalcin expression and mineralized nodule formation in a murine myocyte model of osteogenic trans-differentiation. Finally, KIF26B deficiency inhibited ERK MAP kinase activation during osteogenesis, whilst augmenting p38 and SMAD 1/5/8 phosphorylation. Taken together, these data suggest a role for KIF26B in modulating the severity of post-traumatic HO and provide a potential novel avenue for therapeutic translation.

Obesity and revision surgery, mortality, and patient-reported outcomes after primary knee replacement surgery in the National Joint Registry: A UK cohort study.

BACKGROUND: One in 10 people in the United Kingdom will need a total knee replacement (TKR) during their lifetime. Access to this life-changing operation has recently been restricted based on body mass index (BMI) due to belief that high BMI may lead to poorer outcomes. We investigated the associations between BMI and revision surgery, mortality, and pain/function using what we believe to be the world's largest joint replacement registry. METHODS AND FINDINGS: We analysed 493,710 TKRs in the National Joint Registry (NJR) for England, Wales, Northern Ireland, and the Isle of Man from 2005 to 2016 to investigate 90-day mortality and 10-year cumulative revision. Hospital Episodes Statistics (HES) and Patient Reported Outcome Measures (PROMs) databases were linked to the NJR to investigate change in Oxford Knee Score (OKS) 6 months postoperatively. After adjustment for age, sex, American Society of Anaesthesiologists (ASA) grade, indication for operation, year of primary TKR, and fixation type, patients with high BMI were more likely to undergo revision surgery within 10 years compared to those with "normal" BMI (obese class II hazard ratio (HR) 1.21, 95% CI: 1.10, 1.32 (p < 0.001) and obese class III HR 1.13, 95% CI: 1.02, 1.26 (p = 0.026)). All BMI classes had revision estimates within the recognised 10-year benchmark of 5%. Overweight and obese class I patients had lower mortality than patients with "normal" BMI (HR 0.76, 95% CI: 0.65, 0.90 (p = 0.001) and HR 0.69, 95% CI: 0.58, 0.82 (p < 0.001)). All BMI categories saw absolute increases in OKS after 6 months (range 18-20 points). The relative improvement in OKS was lower in overweight and obese patients than those with "normal" BMI, but the difference was below the minimal detectable change (MDC; 4 points). The main limitations were missing BMI particularly in the early years of data collection and a potential selection bias effect of surgeons selecting the fitter patients with raised BMI for surgery. CONCLUSIONS: Given revision estimates in all BMI groups below the recognised threshold, no evidence of increased mortality, and difference in change in OKS below the MDC, this large national registry shows no evidence of poorer outcomes in patients with high BMI. This study does not support rationing of TKR based on increased BMI.

Linking chondrocyte and synovial transcriptional profile to clinical phenotype in osteoarthritis.

OBJECTIVES: To determine how gene expression profiles in osteoarthritis joint tissues relate to patient phenotypes and whether molecular subtypes can be reproducibly captured by a molecular classification algorithm. METHODS: We analysed RNA sequencing data from cartilage and synovium in 113 osteoarthritis patients, applying unsupervised clustering and Multi-Omics Factor Analysis to characterise transcriptional profiles. We tested the association of the molecularly defined patient subgroups with clinical characteristics from electronic health records. RESULTS: We detected two patient subgroups in low-grade cartilage (showing no/minimal degeneration, cartilage normal/softening only), with differences associated with inflammation, extracellular matrix-related and cell adhesion pathways. The high-inflammation subgroup was associated with female sex (OR 4.12, p=0.0024) and prescription of proton pump inhibitors (OR 4.21, p=0.0040). We identified two independent patient subgroupings in osteoarthritis synovium: one related to inflammation and the other to extracellular matrix and cell adhesion processes. A seven-gene classifier including MMP13, APOD, MMP2, MMP1, CYTL1, IL6 and C15orf48 recapitulated the main axis of molecular heterogeneity in low-grade knee osteoarthritis cartilage (correlation ρ=-0.88, p<10-10) and was reproducible in an independent patient cohort (ρ=-0.85, p<10-10). CONCLUSIONS: These data support the reproducible stratification of osteoarthritis patients by molecular subtype and the exploration of new avenues for tailored treatments.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene.

BACKGROUND: Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology. METHODS: We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees. RESULTS: We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci. CONCLUSIONS: We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

Plant-Based Diet Index and Metabolic Risk in Men: Exploring the Role of the Gut Microbiome.

BACKGROUND: Healthy plant-based diet index (hPDI) is associated with a lower risk of cardiometabolic conditions, but its association as well as interactions with microbiome have not been elucidated. OBJECTIVES: We aimed to investigate the interrelations between hPDI, gut microbiome, and cardiometabolic risk markers. METHODS: hPDI was derived from dietary assessments by a validated FFQ and was examined in relation to metagenomic profiles of 911 fecal samples collected from 303 men aged 71 ± 4 y with an average BMI (in kg/m2) of 25.2 ± 3.6 in the Men's Lifestyle Validation Study. Principal coordinate (PCo) analysis based on Bray-Curtis dissimilarity was conducted, and interactions between hPDI and PCo were examined by using a metabolic risk score composed of blood lipids, BMI, and glycated hemoglobin. RESULTS: After multivariable adjustment, hPDI was significantly associated with the relative abundance of 7 species and 9 pathways. In particular, higher hPDI was significantly associated with a higher relative abundance of Bacteroides cellulosilyticus and Eubacterium eligens, amino acid biosynthesis pathways (l-isoleucine biosynthesis I and III and l-valine biosynthesis), and the pathway of pyruvate fermentation to isobutanol. A favorable association between hPDI and the metabolic risk score was more pronounced among men with a higher PCo characterized by higher abundance of Bacteroides uniformis and lower abundance of Prevotella copri. At the individual species level, a similar interaction was also observed between hPDI and P. copri, as well as with Clostridium clostridioforme or Blautia hydrogenotrophica (all P-interaction < 0.01). CONCLUSION: A greater adherence to a healthy plant-based diet by older men was associated with a microbial profile characterized by a higher abundance of multiple species, including B. cellulosilyticus and E. eligens, as well as pathways in amino acid metabolism and pyruvate fermentation. In addition, inverse associations between healthy plant-based diet and human metabolic risk may partially depend on microbial compositions.

Direct inference of first-year sea ice thickness using broadband acoustic backscattering.

Accurate measurements of sea ice thickness are critical to better understand climate change, to provide situational awareness in ice-covered waters, and to reduce risks for communities that rely on sea ice. Nonetheless, remotely measuring the thickness of sea ice is difficult. The only regularly employed technique that accurately measures the full ice thickness involves drilling a hole through the ice. Other presently used methods are either embedded in or through the ice (e.g., ice mass balance buoys) or calculate thickness from indirect measurements (e.g., ice freeboard from altimetry; ice draft using sonars; total snow and ice thickness using electromagnetic techniques). Acoustic techniques, however, may provide an alternative approach to measure the total ice thickness. Here laboratory-grown sea ice thicknesses, estimated by inverting the time delay between echoes from the water-ice and ice-air interfaces, are compared to those measured using ice cores. A time-domain model capturing the dominant scattering mechanisms is developed to explore the viability of broadband acoustic techniques for measuring sea ice thickness, to compare with experimental measurements, and to investigate optimal frequencies for in situ applications. This approach decouples ice thickness estimates from water column properties and does not preclude ice draft measurements using the same data.

Pathways to understanding the genomic aetiology of osteoarthritis.

Osteoarthritis is a common, complex disease with no curative therapy. In this review, we summarize current knowledge on disease aetiopathogenesis and outline genetics and genomics approaches that are helping catalyse a much-needed improved understanding of the biological underpinning of disease development and progression.

Evaluation of shared genetic aetiology between osteoarthritis and bone mineral density identifies SMAD3 as a novel osteoarthritis risk locus.

Osteoarthritis (OA) is a common complex disease with high public health burden and no curative therapy. High bone mineral density (BMD) is associated with an increased risk of developing OA, suggesting a shared underlying biology. Here, we performed the first systematic overlap analysis of OA and BMD on a genome wide scale. We used summary statistics from the GEFOS consortium for lumbar spine (n = 31,800) and femoral neck (n = 32,961) BMD, and from the arcOGEN consortium for three OA phenotypes (hip, ncases=3,498; knee, ncases=3,266; hip and/or knee, ncases=7,410; ncontrols=11,009). Performing LD score regression we found a significant genetic correlation between the combined OA phenotype (hip and/or knee) and lumbar spine BMD (rg=0.18, P = 2.23 × 10-2), which may be driven by the presence of spinal osteophytes. We identified 143 variants with evidence for cross-phenotype association which we took forward for replication in independent large-scale OA datasets, and subsequent meta-analysis with arcOGEN for a total sample size of up to 23,425 cases and 236,814 controls. We found robustly replicating evidence for association with OA at rs12901071 (OR 1.08 95% CI 1.05-1.11, Pmeta=3.12 × 10-10), an intronic variant in the SMAD3 gene, which is known to play a role in bone remodeling and cartilage maintenance. We were able to confirm expression of SMAD3 in intact and degraded cartilage of the knee and hip. Our findings provide the first systematic evaluation of pleiotropy between OA and BMD, highlight genes with biological relevance to both traits, and establish a robust new OA genetic risk locus at SMAD3.

The interleukin-33 receptor contributes to pulmonary responses to ozone in male mice: role of the microbiome.

BACKGROUND: Interleukin-33 is released in the airways following acute ozone exposure and has the ability to cause airway hyperresponsiveness, a defining feature of asthma. Ozone causes greater airway hyperresponsiveness in male than female mice. Moreover, sex differences in the gut microbiome account for sex differences in this response to ozone. The purpose of this study was to determine whether there were sex differences in the role of interleukin-33 in ozone-induced airway hyperresponsiveness and to examine the role of the microbiome in these events. METHODS: Wildtype mice and mice genetically deficient in ST2, the interleukin-33 receptor, were housed from weaning with either other mice of the same genotype and sex, or with mice of the same sex but opposite genotype. At 15 weeks of age, fecal pellets were harvested for 16S rRNA sequencing and the mice were then exposed to air or ozone. Airway responsiveness was measured and a bronchoalveolar lavage was performed 24 h after exposure. RESULTS: In same-housed mice, ozone-induced airway hyperresponsiveness was greater in male than female wildtype mice. ST2 deficiency reduced ozone-induced airway hyperresponsiveness in male but not female mice and abolished sex differences in the response to ozone. However, sex differences in the role of interleukin-33 were unrelated to type 2 cytokine release: ozone-induced increases in bronchoalveolar lavage interleukin-5 were greater in females than males and ST2 deficiency virtually abolished interleukin-5 in both sexes. Since gut microbiota contribute to sex differences in ozone-induced airway hyperresponsiveness, we examined the role of the microbiome in these ST2-dependent sex differences. To do so, we cohoused wildtype and ST2 deficient mice, a situation that allows for transfer of microbiota among cage-mates. Cohousing altered the gut microbial community structure, as indicated by 16S rRNA gene sequencing of fecal DNA and reversed the effect of ST2 deficiency on pulmonary responses to ozone in male mice. CONCLUSIONS: The data indicate that the interleukin-33 /ST2 pathway contributes to ozone-induced airway hyperresponsiveness in male mice and suggest that the role of interleukin-33 is mediated at the level of the gut microbiome.

Estimating an Individual's Probability of Revision Surgery After Knee Replacement: A Comparison of Modeling Approaches Using a National Data Set.

Tools that provide personalized risk prediction of outcomes after surgical procedures help patients make preference-based decisions among the available treatment options. However, it is unclear which modeling approach provides the most accurate risk estimation. We constructed and compared several parametric and nonparametric models for predicting prosthesis survivorship after knee replacement surgery for osteoarthritis. We used 430,455 patient-procedure episodes between April 2003 and September 2015 from the National Joint Registry for England, Wales, Northern Ireland, and the Isle of Man. The flexible parametric survival and random survival forest models most accurately captured the observed probability of remaining event-free. The concordance index for the flexible parametric model was the highest (0.705, 95% confidence interval (CI): 0.702, 0.707) for total knee replacement and was 0.639 (95% CI: 0.634, 0.643) for unicondylar knee replacement and 0.589 (95% CI: 0.586, 0.592) for patellofemoral replacement. The observed-to-predicted ratios for both the flexible parametric and the random survival forest approaches indicated that models tended to underestimate the risks for most risk groups. Our results show that the flexible parametric model has a better overall performance compared with other tested parametric methods and has better discrimination compared with the random survival forest approach.

A novel variant in GLIS3 is associated with osteoarthritis.

OBJECTIVES: Osteoarthritis (OA) is a complex disease, but its genetic aetiology remains poorly characterised. To identify novel susceptibility loci for OA, we carried out a genome-wide association study (GWAS) in individuals from the largest UK-based OA collections to date. METHODS: We carried out a discovery GWAS in 5414 OA individuals with knee and/or hip total joint replacement (TJR) and 9939 population-based controls. We followed-up prioritised variants in OA subjects from the interim release of the UK Biobank resource (up to 12 658 cases and 50 898 controls) and our lead finding in operated OA subjects from the full release of UK Biobank (17 894 cases and 89 470 controls). We investigated its functional implications in methylation, gene expression and proteomics data in primary chondrocytes from 12 pairs of intact and degraded cartilage samples from patients undergoing TJR. RESULTS: We detect a genome-wide significant association at rs10116772 with TJR (P=3.7×10-8; for allele A: OR (95% CI) 0.97 (0.96 to 0.98)), an intronic variant in GLIS3, which is expressed in cartilage. Variants in strong correlation with rs10116772 have been associated with elevated plasma glucose levels and diabetes. CONCLUSIONS: We identify a novel susceptibility locus for OA that has been previously implicated in diabetes and glycaemic traits.