Objectives: To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses... Show moreObjectives: To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients. Methods: RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the genome-wide association studies (GWAS) catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (phi) was calculated for each heterozygous individual with the risk SNP or with the SNP in linkage disequilibrium (LD) with it (r(2) > 0.6). Furthermore, a meta-analysis was performed to generate a meta-phi (null hypothesis median phi = 0.49) and P-value for each SNP. Results: We identified 30 transcript SNPs (28 in cartilage and two in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-phi = 0.54, FDR = 1.7x10(-4)) and ILF3-DT (meta-phi = 0.6, FDR = 1.75x10(-5)). Moreover, 12 drugs were interacting with seven genes displaying AEI, of which seven drugs have been already approved. Conclusions: By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harbouring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further translation towards underlying biological mechanisms. Show less
Almeida, R.C. de; Mahfouz, A.; Mei, H.L.; Houtman, E.; Hollander, W. den; Soul, J.; ... ; Meulenbelt, I. 2021
Objective. To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant... Show moreObjective. To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant characteristics.Methods. This study includes n = 66 primary OA patients (41 knees and 25 hips), who underwent a joint replacement surgery, from which macroscopically unaffected (preserved, n = 56) and lesioned (n = 45) OA articular cartilage were collected [Research Arthritis and Articular Cartilage (RAAK) study]. Unsupervised hierarchical clustering analysis on preserved cartilage transcriptome followed by clinical data integration was performed. Protein-protein interaction (PPI) followed by pathway enrichment analysis were done for genes significant differentially expressed between subgroups with interactions in the PPI network.Results. Analysis of preserved samples (n = 56) resulted in two OA subtypes with n = 41 (cluster A) and n = 15 (cluster B) patients. The transcriptomic profile of cluster B cartilage, relative to cluster A (DE-AB genes) showed among others a pronounced upregulation of multiple genes involved in chemokine pathways. Nevertheless, upon investigating the OA pathophysiology in cluster B patients as reflected by differentially expressed genes between preserved and lesioned OA cartilage (DE-OA-B genes), the chemokine genes were significantly downregulated with OA pathophysiology. Upon integrating radiographic OA data, we showed that the OA phenotype among cluster B patients, relative to cluster A, may be characterized by higher joint space narrowing (JSN) scores and low osteophyte (OP) scores.Conclusion. Based on whole-transcriptome profiling, we identified two robust OA subtypes characterized by unique OA, pathophysiological processes in cartilage as well as a clinical phenotype. We advocate that further characterization, confirmation and clinical data integration is a prerequisite to allow for development of treatments towards personalized care with concurrently more effective treatment response. Show less
Meessen, J.M.T.A.; Saberi-Hosnijeh, F.; Bomer, N.; Hollander, W. den; Bom, J.G. van der; Hilten, J.A. van; ... ; Meulenbelt, I. 2020
Higher body mass index (BMI) is associated with osteoarthritis (OA) in both weight-bearing and non-weight-bearing joints, suggesting a link between OA and poor metabolic health beyond mechanical... Show moreHigher body mass index (BMI) is associated with osteoarthritis (OA) in both weight-bearing and non-weight-bearing joints, suggesting a link between OA and poor metabolic health beyond mechanical loading. This risk may be influenced by systemic factors accompanying BMI. Fluctuations in concentrations of metabolites may mark or even contribute to development of OA. This study explores the association of metabolites with radiographic knee/hip OA prevalence and progression. A H-1-NMR-metabolomics assay was performed on plasma samples of 1564 cases for prevalent OA and 2,125 controls collected from the Rotterdam Study, CHECK, GARP/NORREF and LUMC-arthroplasty cohorts. OA prevalence and 5 to 10 year progression was assessed by means of Kellgren-Lawrence (KL) score and the OARSI-atlas. End-stage knee/hip OA (TJA) was defined as indication for arthroplasty surgery. Controls did not have OA at baseline or follow-up. Principal component analysis of 227 metabolites demonstrated 23 factors, of which 19 remained interpretable after quality-control. Associations of factor scores with OA definitions were investigated with logistic regression. Fatty acids chain length (FALen), which was included in two factors which associated with TJA, was individually associated with both overall OA as well as TJA. Increased Fatty Acid chain Length is associated with OA. Show less
Bosch, M.H.J. van den; Ramos, Y.F.M.; Hollander, W. den; Bomer, N.; Nelissen, R.G.H.H.; Bovee, J.V.M.G.; ... ; Meulenbelt, I. 2019
In dit proefschrift worden de moleculaire mechanismen behandeld die onderliggend zijn aan artrose. Specifiek wordt genoomwijd onderzocht welke genen anders tot expressie komen in aangedaan... Show moreIn dit proefschrift worden de moleculaire mechanismen behandeld die onderliggend zijn aan artrose. Specifiek wordt genoomwijd onderzocht welke genen anders tot expressie komen in aangedaan vergeleken met gezond kraakbeen van artrose patienten. Dit in de context van epigenetische regulatie van gen expressie, specifiek door DNA methylatie in het licht van de lokale genetische context in de vorm van puntmutaties. Show less
Cornelis, F.; Monteagudo, S.; Hollander, W. den; Peeters, T.; Guns, L.A.; Storms, L.; ... ; Lories, R. 2017