1
|
Bloks NGC, Dicks A, Harissa Z, Nelissen RGHH, Hajmousa G, Ramos YFM, de Almeida RC, Guilak F, Meulenbelt I. Hyper-physiologic mechanical cues, as an osteoarthritis disease-relevant environmental perturbation, cause a critical shift in set points of methylation at transcriptionally active CpG sites in neo-cartilage organoids. Clin Epigenetics 2024; 16:64. [PMID: 38730337 PMCID: PMC11087253 DOI: 10.1186/s13148-024-01676-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Accepted: 04/30/2024] [Indexed: 05/12/2024] Open
Abstract
BACKGROUND Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (Guilak in Best Pract Res Clin Rheumatol 25:815-823, 2011). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state. RESULTS We showed that hyper-physiological loading evokes consistent changes in CpGs associated with expression changes (ML-tCpGs) in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology (OA-tCpGs), we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity. CONCLUSION Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCpGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.
Collapse
Affiliation(s)
- Niek G C Bloks
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Einthovenweg 20, 2333 ZC, Leiden, The Netherlands
| | - Amanda Dicks
- Washington University, Saint Louis, MO, USA
- Shriners Hospitals for Children, Saint Louis, MO, USA
| | - Zainab Harissa
- Washington University, Saint Louis, MO, USA
- Shriners Hospitals for Children, Saint Louis, MO, USA
| | - Rob G H H Nelissen
- Department Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Ghazaleh Hajmousa
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Einthovenweg 20, 2333 ZC, Leiden, The Netherlands
| | - Yolande F M Ramos
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Einthovenweg 20, 2333 ZC, Leiden, The Netherlands
| | - Rodrigo Coutinho de Almeida
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Einthovenweg 20, 2333 ZC, Leiden, The Netherlands
| | - Farshid Guilak
- Washington University, Saint Louis, MO, USA
- Shriners Hospitals for Children, Saint Louis, MO, USA
| | - Ingrid Meulenbelt
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Einthovenweg 20, 2333 ZC, Leiden, The Netherlands.
| |
Collapse
|
2
|
Korthagen NM, Houtman E, Boone I, Coutinho de Almeida R, Sivasubramaniyan K, Mahdad R, Nelissen RGHH, Ramos YFM, Tessari MA, Meulenbelt I. Thyroid hormone induces ossification and terminal maturation in a preserved OA cartilage biomimetic model. Arthritis Res Ther 2024; 26:91. [PMID: 38664820 PMCID: PMC11044551 DOI: 10.1186/s13075-024-03326-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2023] [Accepted: 04/21/2024] [Indexed: 04/29/2024] Open
Abstract
OBJECTIVE To characterize aspects of triiodothyronine (T3) induced chondrocyte terminal maturation within the molecular osteoarthritis pathophysiology using the previously established T3 human ex vivo osteochondral explant model. DESIGNS RNA-sequencing was performed on explant cartilage obtained from OA patients (n = 8), that was cultured ex vivo with or without T3 (10 ng/ml), and main findings were validated using RT-qPCR in an independent sample set (n = 22). Enrichment analysis was used for functional clustering and comparisons with available OA patient RNA-sequencing and GWAS datasets were used to establish relevance for OA pathophysiology by linking to OA patient genomic profiles. RESULTS Besides the upregulation of known hypertrophic genes EPAS1 and ANKH, T3 treatment resulted in differential expression of 247 genes with main pathways linked to extracellular matrix and ossification. CCDC80, CDON, ANKH and ATOH8 were among the genes found to consistently mark early, ongoing and terminal maturational OA processes in patients. Furthermore, among the 37 OA risk genes that were significantly affected in cartilage by T3 were COL12A1, TNC, SPARC and PAPPA. CONCLUSIONS RNA-sequencing results show that metabolic activation and recuperation of growth plate morphology are induced by T3 in OA chondrocytes, indicating terminal maturation is accelerated. The molecular mechanisms involved in hypertrophy were linked to all stages of OA pathophysiology and will be used to validate disease models for drug testing.
Collapse
Affiliation(s)
- N M Korthagen
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - E Houtman
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - I Boone
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - R Coutinho de Almeida
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - K Sivasubramaniyan
- Galapagos BV, Willem Einthovenstraat 13, Oegstgeest, 2342 BH, The Netherlands
| | - R Mahdad
- Alrijne hospital, Simon Smitweg 1, Leiderdorp, 2353 GA, The Netherlands
| | - R G H H Nelissen
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - Y F M Ramos
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands
| | - M A Tessari
- Galapagos BV, Willem Einthovenstraat 13, Oegstgeest, 2342 BH, The Netherlands
| | - I Meulenbelt
- Department Biomedical Data Sciences, Section of Molecular Epidemiology, LUMC, Einthovenweg 20, Postzone S05-P, 2333 ZC, Leiden, The Netherlands.
| |
Collapse
|
3
|
Govindaraj K, Meteling M, van Rooij J, Becker M, van Wijnen AJ, van den Beucken JJJP, Ramos YFM, van Meurs J, Post JN, Leijten J. Osmolarity-Induced Altered Intracellular Molecular Crowding Drives Osteoarthritis Pathology. Adv Sci (Weinh) 2024; 11:e2306722. [PMID: 38213111 PMCID: PMC10953583 DOI: 10.1002/advs.202306722] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 11/08/2023] [Indexed: 01/13/2024]
Abstract
Osteoarthritis (OA) is a multifactorial degenerative joint disease of which the underlying mechanisms are yet to be fully understood. At the molecular level, multiple factors including altered signaling pathways, epigenetics, metabolic imbalance, extracellular matrix degradation, production of matrix metalloproteinases, and inflammatory cytokines, are known to play a detrimental role in OA. However, these factors do not initiate OA, but are mediators or consequences of the disease, while many other factors causing the etiology of OA are still unknown. Here, it is revealed that microenvironmental osmolarity can induce and reverse osteoarthritis-related behavior of chondrocytes via altered intracellular molecular crowding, which represents a previously unknown mechanism underlying OA pathophysiology. Decreased intracellular crowding is associated with increased sensitivity to proinflammatory triggers and decreased responsiveness to anabolic stimuli. OA-induced lowered intracellular molecular crowding could be renormalized via exposure to higher extracellular osmolarity such as those found in healthy joints, which reverse OA chondrocyte's sensitivity to catabolic stimuli as well as its glycolytic metabolism.
Collapse
Affiliation(s)
- Kannan Govindaraj
- Department of Developmental BioengineeringFaculty of Science and Technology, Technical Medical CentreUniversity of TwenteDrienerlolaan 5Enschede7522NBThe Netherlands
| | - Marieke Meteling
- Department of Developmental BioengineeringFaculty of Science and Technology, Technical Medical CentreUniversity of TwenteDrienerlolaan 5Enschede7522NBThe Netherlands
| | - Jeroen van Rooij
- Department of Internal MedicineErasmus MCDr. Molewaterplein 40Rotterdam3015GDThe Netherlands
| | - Malin Becker
- Department of Developmental BioengineeringFaculty of Science and Technology, Technical Medical CentreUniversity of TwenteDrienerlolaan 5Enschede7522NBThe Netherlands
| | | | | | - Yolande F. M. Ramos
- Department of Biomedical Data SciencesSection Molecular EpidemiologyLUMCEinthovenweg 20Leiden2333 ZCThe Netherlands
| | - Joyce van Meurs
- Department of Internal MedicineErasmus MCDr. Molewaterplein 40Rotterdam3015GDThe Netherlands
- Department of Orthopedics & Sports MedicineErasmus MCDr. Molewaterplein 40Rotterdam3015GDThe Netherlands
| | - Janine N. Post
- Department of Developmental BioengineeringFaculty of Science and Technology, Technical Medical CentreUniversity of TwenteDrienerlolaan 5Enschede7522NBThe Netherlands
| | - Jeroen Leijten
- Department of Developmental BioengineeringFaculty of Science and Technology, Technical Medical CentreUniversity of TwenteDrienerlolaan 5Enschede7522NBThe Netherlands
| |
Collapse
|
4
|
Ramos YFM, Rice SJ, Ali SA, Pastrello C, Jurisica I, Rai MF, Collins KH, Lang A, Maerz T, Geurts J, Ruiz-Romero C, June RK, Thomas Appleton C, Rockel JS, Kapoor M. Evolution and advancements in genomics and epigenomics in OA research: How far we have come. Osteoarthritis Cartilage 2024:S1063-4584(24)00054-2. [PMID: 38428513 DOI: 10.1016/j.joca.2024.02.656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 01/29/2024] [Accepted: 02/25/2024] [Indexed: 03/03/2024]
Abstract
OBJECTIVE Osteoarthritis (OA) is the most prevalent musculoskeletal disease affecting articulating joint tissues, resulting in local and systemic changes that contribute to increased pain and reduced function. Diverse technological advancements have culminated in the advent of high throughput "omic" technologies, enabling identification of comprehensive changes in molecular mediators associated with the disease. Amongst these technologies, genomics and epigenomics - including methylomics and miRNomics, have emerged as important tools to aid our biological understanding of disease. DESIGN In this narrative review, we selected articles discussing advancements and applications of these technologies to OA biology and pathology. We discuss how genomics, deoxyribonucleic acid (DNA) methylomics, and miRNomics have uncovered disease-related molecular markers in the local and systemic tissues or fluids of OA patients. RESULTS Genomics investigations into the genetic links of OA, including using genome-wide association studies, have evolved to identify 100+ genetic susceptibility markers of OA. Epigenomic investigations of gene methylation status have identified the importance of methylation to OA-related catabolic gene expression. Furthermore, miRNomic studies have identified key microRNA signatures in various tissues and fluids related to OA disease. CONCLUSIONS Sharing of standardized, well-annotated omic datasets in curated repositories will be key to enhancing statistical power to detect smaller and targetable changes in the biological signatures underlying OA pathogenesis. Additionally, continued technological developments and analysis methods, including using computational molecular and regulatory networks, are likely to facilitate improved detection of disease-relevant targets, in-turn, supporting precision medicine approaches and new treatment strategies for OA.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Dept. Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Sarah J Rice
- Biosciences Institute, International Centre for Life, Newcastle University, Newcastle upon Tyne, United Kingdom
| | - Shabana Amanda Ali
- Henry Ford Health + Michigan State University Health Sciences, Detroit, MI, USA
| | - Chiara Pastrello
- Osteoarthritis Research Program, Division of Orthopedic Surgery, Schroeder Arthritis Institute, UHN, Toronto, Ontario, Canada
| | - Igor Jurisica
- Osteoarthritis Research Program, Division of Orthopedic Surgery, Schroeder Arthritis Institute, UHN, Toronto, Ontario, Canada; Departments of Medical Biophysics and Computer Science, University of Toronto, Toronto, Ontario, Canada
| | - Muhammad Farooq Rai
- Department of Biological Sciences, Center for Biotechnology, College of Medicine & Health Sciences, Khalifa University, Abu Dhabi, United Arab Emirates
| | - Kelsey H Collins
- Department of Orthopaedic Surgery, University of California San Francisco, San Francisco, CA, USA
| | - Annemarie Lang
- Departments of Orthopaedic Surgery and Bioengineering, University of Pennsylvania, Philadelphia, PA, USA
| | - Tristan Maerz
- Department of Orthopaedic Surgery, University of Michigan, Ann Arbor, MI, USA
| | - Jeroen Geurts
- Rheumatology, Department of Musculoskeletal Medicine, Lausanne University Hospital, Lausanne, Switzerland
| | - Cristina Ruiz-Romero
- Grupo de Investigación de Reumatología (GIR), Unidad de Proteómica, INIBIC -Hospital Universitario A Coruña, SERGAS, A Coruña, Spain
| | - Ronald K June
- Department of Mechanical & Industrial Engineering, Montana State University, Bozeman, MT, USA
| | - C Thomas Appleton
- Department of Medicine, University of Western Ontario, London, Ontario, Canada
| | - Jason S Rockel
- Osteoarthritis Research Program, Division of Orthopedic Surgery, Schroeder Arthritis Institute, UHN, Toronto, Ontario, Canada
| | - Mohit Kapoor
- Osteoarthritis Research Program, Division of Orthopedic Surgery, Schroeder Arthritis Institute, UHN, Toronto, Ontario, Canada.
| |
Collapse
|
5
|
Timmermans RGM, Blom AB, Nelissen RGHH, Broekhuis D, van der Kraan PM, Meulenbelt I, van den Bosch MHJ, Ramos YFM. Mechanical stress and inflammation have opposite effects on Wnt signaling in human chondrocytes. J Orthop Res 2024; 42:286-295. [PMID: 37525432 DOI: 10.1002/jor.25673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Revised: 07/12/2023] [Accepted: 07/26/2023] [Indexed: 08/02/2023]
Abstract
Dysregulation of Wingless and Int-1 (Wnt) signaling has been strongly associated with development and progression of osteoarthritis (OA). Here, we set out to investigate the independent effects of either mechanical stress (MS) or inflammation on Wnt signaling in human neocartilage pellets, and to relate this Wnt signaling to OA pathophysiology. OA synovium-conditioned media (OAS-CM) was collected after incubating synovium from human end-stage OA joints for 24 h in medium. Cytokine levels in the OAS-CM were determined with a multiplex immunoassay (Luminex). Human neocartilage pellets were exposed to 20% MS, 2% OAS-CM or 1 ng/mL Interleukin-1β (IL-1β). Effects on expression levels of Wnt signaling members were determined by reverse transcription-quantitative polymerase chain reaction. Additionally, the expression of these members in articular cartilage from human OA joints was analyzed in association with joint space narrowing (JSN) and osteophyte scores. Protein levels of IL-1β, IL-6, IL-8, IL-10, tumor necrosis factor α, and granulocyte-macrophage colony-stimulating factor positively correlated with each other. MS increased noncanonical WNT5A and FOS expression. In contrast, these genes were downregulated upon stimulation with OAS-CM or IL-1β. Furthermore, Wnt inhibitors DKK1 and FRZB decreased in response to OAS-CM or IL-1β exposure. Finally, expression of WNT5A in OA articular cartilage was associated with increased JSN scores, but not osteophyte scores. Our results demonstrate that MS and inflammatory stimuli have opposite effects on canonical and noncanonical Wnt signaling in human neocartilage. Considering the extent to which MS and inflammation contribute to OA in individual patients, we hypothesize that targeting specific Wnt pathways offers a more effective, individualized approach.
Collapse
Affiliation(s)
- Ritchie G M Timmermans
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Arjen B Blom
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Demiën Broekhuis
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Peter M van der Kraan
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
6
|
Todtenhaupt P, Franken LA, Groene SG, van Hoolwerff M, van der Meeren LE, van Klink JMM, Roest AAW, de Bruin C, Ramos YFM, Haak MC, Lopriore E, Heijmans BT, van Pel M. A robust and standardized method to isolate and expand mesenchymal stromal cells from human umbilical cord. Cytotherapy 2023; 25:1057-1068. [PMID: 37516948 DOI: 10.1016/j.jcyt.2023.07.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 06/22/2023] [Accepted: 07/14/2023] [Indexed: 07/31/2023]
Abstract
BACKGROUND AIMS Human umbilical cord-derived mesenchymal stromal cells (hUC-MSCs) are increasingly used in research and therapy. To obtain hUC-MSCs, a diversity of isolation and expansion methods are applied. Here, we report on a robust and standardized method for hUC-MSC isolation and expansion. METHODS Using 90 hUC donors, we compared and optimized critical variables during each phase of the multi-step procedure involving UC collection, processing, MSC isolation, expansion and characterization. Furthermore, we assessed the effect of donor-to-donor variability regarding UC morphology and donor attributes on hUC-MSC characteristics. RESULTS We demonstrated robustness of our method across 90 UC donors at each step of the procedure. With our method, UCs can be collected up to 6 h after birth, and UC-processing can be initiated up to 48 h after collection without impacting on hUC-MSC characteristics. The removal of blood vessels before explant cultures improved hUC-MSC purity. Expansion in Minimum essential medium α supplemented with human platelet lysate increased reproducibility of the expansion rate and MSC characteristics as compared with Dulbecco's Modified Eagle's Medium supplemented with fetal bovine serum. The isolated hUC-MSCs showed a purity of ∼98.9%, a viability of >97% and a high proliferative capacity. Trilineage differentiation capacity of hUC-MSCs was reduced as compared with bone marrow-derived MSCs. Functional assays indicated that the hUC-MSCs were able to inhibit T-cell proliferation demonstrating their immune-modulatory capacity. CONCLUSIONS We present a robust and standardized method to isolate and expand hUC-MSCs, minimizing technical variability and thereby lay a foundation to advance reliability and comparability of results obtained from different donors and different studies.
Collapse
Affiliation(s)
- Pia Todtenhaupt
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands; Neonatology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Laura A Franken
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Sophie G Groene
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands; Neonatology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Marcella van Hoolwerff
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Lotte E van der Meeren
- Department of Pathology, Leiden University Medical Center, Leiden, The Netherlands; Department of Pathology, Erasmus Medical Center, Leiden, The Netherlands
| | - Jeanine M M van Klink
- Neonatology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Arno A W Roest
- Pediatric Cardiology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Christiaan de Bruin
- Pediatric Endocrinology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Monique C Haak
- Fetal Medicine, Department of Obstetrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Enrico Lopriore
- Neonatology, Willem-Alexander Children's Hospital, Department of Pediatrics, Leiden University Medical Center, Leiden, The Netherlands
| | - Bastiaan T Heijmans
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Melissa van Pel
- NecstGen, Leiden, The Netherlands; Department of Internal Medicine, Leiden University Medical Center, Leiden, The Netherlands.
| |
Collapse
|
7
|
Timmermans RGM, Blom AB, Bloks NGC, Nelissen RGHH, van der Linden EHMJ, van der Kraan PM, Meulenbelt I, Ramos YFM, van den Bosch MHJ. CCN4/WISP1 Promotes Migration of Human Primary Osteoarthritic Chondrocytes. Cartilage 2023; 14:67-75. [PMID: 36546648 PMCID: PMC10076902 DOI: 10.1177/19476035221144747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
OBJECTIVES Previously, we have shown the involvement of cellular communication network factor 4/Wnt-activated protein Wnt-1-induced signaling protein 1 (CCN4/WISP1) in osteoarthritic (OA) cartilage and its detrimental effects on cartilage. Here, we investigated characteristics of CCN4 in chondrocyte biology by exploring correlations of CCN4 with genes expressed in human OA cartilage with functional follow-up. DESIGN Spearman correlation analysis was performed for genes correlating with CCN4 using our previously established RNA sequencing dataset of human preserved OA cartilage of the RAAK study, followed by a pathway enrichment analysis for genes with ρ ≥|0.6.| Chondrocyte migration in the absence or presence of CCN4 was determined in a scratch assay, measuring scratch size using a live cell imager for up to 36 h. Changes in expression levels of 12 genes, correlating with CCN4 and involved in migratory processes, were determined with reverse transcription-quantitative polymerase chain reaction (RT-qPCR). RESULTS Correlation of CCN4 with ρ ≥|0.6| was found for 58 genes in preserved human OA cartilage. Pathway analysis revealed "neural crest cell migration" as most significant enriched pathway, containing among others CORO1C, SEMA3C, and SMO. Addition of CCN4 to primary chondrocytes significantly enhance chondrocyte migration as demonstrated by reduced scratch size over the course of 36 h, but at the timepoints measured no effect was observed on mRNA expression of the 12 genes. CONCLUSION CCN4 increases cell migration of human primary OA chondrocytes. Since WISP1 expression is known to be increased in OA cartilage, this may serve to direct chondrocytes toward cartilage defects and orchestrate repair.
Collapse
Affiliation(s)
- Ritchie G M Timmermans
- Experimental Rheumatology, Radboud university medical center, Nijmegen, The Netherlands
- Section Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
| | - Arjen B Blom
- Experimental Rheumatology, Radboud university medical center, Nijmegen, The Netherlands
| | - Niek G C Bloks
- Section Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Centre, Leiden, The Netherlands
| | | | - Peter M van der Kraan
- Experimental Rheumatology, Radboud university medical center, Nijmegen, The Netherlands
| | - Ingrid Meulenbelt
- Section Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
| | - Yolande F M Ramos
- Section Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
| | | |
Collapse
|
8
|
van Hoolwerff M, Tuerlings M, Wijnen IJL, Suchiman HED, Cats D, Mei H, Nelissen RGHH, van der Linden-van der Zwaag HMJ, Ramos YFM, Coutinho de Almeida R, Meulenbelt I. Identification and functional characterization of imbalanced osteoarthritis-associated fibronectin splice variants. Rheumatology (Oxford) 2023; 62:894-904. [PMID: 35532170 PMCID: PMC9891405 DOI: 10.1093/rheumatology/keac272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 03/31/2022] [Accepted: 04/24/2022] [Indexed: 02/04/2023] Open
Abstract
OBJECTIVE To identify FN1 transcripts associated with OA pathophysiology and investigate the downstream effects of modulating FN1 expression and relative transcript ratio. METHODS FN1 transcriptomic data was obtained from our previously assessed RNA-seq dataset of lesioned and preserved OA cartilage samples from the Research osteoArthritis Articular Cartilage (RAAK) study. Differential transcript expression analysis was performed on all 27 FN1 transcripts annotated in the Ensembl database. Human primary chondrocytes were transduced with lentiviral particles containing short hairpin RNA (shRNA) targeting full-length FN1 transcripts or non-targeting shRNA. Subsequently, matrix deposition was induced in our 3D in vitro neo-cartilage model. Effects of changes in the FN1 transcript ratio on sulphated glycosaminoglycan (sGAG) deposition were investigated by Alcian blue staining and dimethylmethylene blue assay. Moreover, gene expression levels of 17 cartilage-relevant markers were determined by reverse transcription quantitative polymerase chain reaction. RESULTS We identified 16 FN1 transcripts differentially expressed between lesioned and preserved cartilage. FN1-208, encoding migration-stimulating factor, was the most significantly differentially expressed protein coding transcript. Downregulation of full-length FN1 and a concomitant increased FN1-208 ratio resulted in decreased sGAG deposition as well as decreased ACAN and COL2A1 and increased ADAMTS-5, ITGB1 and ITGB5 gene expression levels. CONCLUSION We show that full-length FN1 downregulation and concomitant relative FN1-208 upregulation was unbeneficial for deposition of cartilage matrix, likely due to decreased availability of the classical RGD (Arg-Gly-Asp) integrin-binding site of fibronectin.
Collapse
Affiliation(s)
| | - Margo Tuerlings
- Department of Biomedical Data Sciences, Section Molecular Epidemiology
| | - Imke J L Wijnen
- Department of Biomedical Data Sciences, Section Molecular Epidemiology
| | - H Eka D Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology
| | | | | | - Rob G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology
| | | | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology
| |
Collapse
|
9
|
Tuerlings M, Janssen GMC, Boone I, van Hoolwerff M, Rodriguez Ruiz A, Houtman E, Suchiman HED, van der Wal RJP, Nelissen RGHH, Coutinho de Almeida R, van Veelen PA, Ramos YFM, Meulenbelt I. WWP2 confers risk to osteoarthritis by affecting cartilage matrix deposition via hypoxia associated genes. Osteoarthritis Cartilage 2023; 31:39-48. [PMID: 36208715 DOI: 10.1016/j.joca.2022.09.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 09/12/2022] [Accepted: 09/28/2022] [Indexed: 11/06/2022]
Abstract
OBJECTIVE To explore the co-expression network of the osteoarthritis (OA) risk gene WWP2 in articular cartilage and study cartilage characteristics when mimicking the effect of OA risk allele rs1052429-A on WWP2 expression in a human 3D in vitro model of cartilage. METHOD Co-expression behavior of WWP2 with genes expressed in lesioned OA articular cartilage (N = 35 samples) was explored. By applying lentiviral particle mediated WWP2 upregulation in 3D in vitro pellet cultures of human primary chondrocytes (N = 8 donors) the effects of upregulation on cartilage matrix deposition was evaluated. Finally, we transfected primary chondrocytes with miR-140 mimics to evaluate whether miR-140 and WWP2 are involved in similar pathways. RESULTS Upon performing Spearman correlations in lesioned OA cartilage, 98 highly correlating genes (|ρ| > 0.7) were identified. Among these genes, we identified GJA1, GDF10, STC2, WDR1, and WNK4. Subsequent upregulation of WWP2 on 3D chondrocyte pellet cultures resulted in a decreased expression of COL2A1 and ACAN and an increase in EPAS1 expression. Additionally, we observed a decreased expression of GDF10, STC2, and GJA1. Proteomics analysis identified 42 proteins being differentially expressed with WWP2 upregulation, which were enriched for ubiquitin conjugating enzyme activity. Finally, upregulation of miR-140 in 2D chondrocytes resulted in significant upregulation of WWP2 and WDR1. CONCLUSIONS Mimicking the effect of OA risk allele rs1052429-A on WWP2 expression initiates detrimental processes in the cartilage shown by a response in hypoxia associated genes EPAS1, GDF10, and GJA1 and a decrease in anabolic markers, COL2A1 and ACAN.
Collapse
Affiliation(s)
- M Tuerlings
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - G M C Janssen
- Center for Proteomics and Metabolomics, Leiden University Medical Center, Leiden, the Netherlands.
| | - I Boone
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - M van Hoolwerff
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - A Rodriguez Ruiz
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - E Houtman
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - H E D Suchiman
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - R J P van der Wal
- Dept. Orthopaedics, Leiden University Medical Center, Leiden, the Netherlands.
| | - R G H H Nelissen
- Dept. Orthopaedics, Leiden University Medical Center, Leiden, the Netherlands.
| | - R Coutinho de Almeida
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - P A van Veelen
- Center for Proteomics and Metabolomics, Leiden University Medical Center, Leiden, the Netherlands.
| | - Y F M Ramos
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - I Meulenbelt
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| |
Collapse
|
10
|
Ramos YFM, Tertel T, Shaw G, Staubach S, de Almeida RC, Suchiman E, Kuipers TB, Mei H, Barry F, Murphy M, Giebel B, Meulenbelt I. Characterizing the secretome of licensed hiPSC-derived MSCs. Stem Cell Res Ther 2022; 13:434. [PMID: 36056373 PMCID: PMC9438242 DOI: 10.1186/s13287-022-03117-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 08/04/2022] [Indexed: 11/23/2022] Open
Abstract
Although mesenchymal stromal cells (MSCs) from primary tissues have been successfully applied in the clinic, their expansion capabilities are limited and results are variable. MSCs derived from human-induced pluripotent stem cells (hiMSCs) are expected to overcome these limitations and serve as a reproducible and sustainable cell source. We have explored characteristics and therapeutic potential of hiMSCs in comparison to hBMSCs. RNA sequencing confirmed high resemblance, with average Pearson correlation of 0.88 and Jaccard similarity index of 0.99, and similar to hBMSCs the hiMSCs released extracellular vesicles with in vitro immunomodulatory properties. Potency assay with TNFα and IFNγ demonstrated an increase in well-known immunomodulatory genes such as IDO1, CXCL8/IL8, and HLA-DRA which was also highlighted by enhanced secretion in the media. Notably, expression of 125 genes increased more than 1000-fold. These genes were predicted to be regulated by NFΚB signaling, known to play a central role in immune response. Altogether, our data qualify hiMSCs as a promising source for cell therapy and/or cell-based therapeutic products. Additionally, the herewith generated database will add to our understanding of the mode of action of regenerative cell-based therapies and could be used to identify relevant potency markers.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, LUMC Postzone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands.
| | - Tobias Tertel
- Institute for Transfusion Medicine, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
| | - Georgina Shaw
- National University of Ireland Galway, Galway, Ireland
| | - Simon Staubach
- Institute for Transfusion Medicine, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
| | - Rodrigo Coutinho de Almeida
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, LUMC Postzone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands
| | - Eka Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, LUMC Postzone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands
| | | | - Hailiang Mei
- LUMC, Sequencing Analysis Support Core, Leiden, The Netherlands
| | - Frank Barry
- National University of Ireland Galway, Galway, Ireland
| | - Mary Murphy
- National University of Ireland Galway, Galway, Ireland
| | - Bernd Giebel
- Institute for Transfusion Medicine, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, LUMC Postzone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands
| |
Collapse
|
11
|
Rodríguez Ruiz A, van Hoolwerff M, Sprangers S, Suchiman E, Schoenmaker T, Dibbets-Schneider P, Bloem JL, Nelissen RGHH, Freund C, Mummery C, Everts V, de Vries TJ, Ramos YFM, Meulenbelt I. Mutation in the CCAL1 locus accounts for bidirectional process of human subchondral bone turnover and cartilage mineralization. Rheumatology (Oxford) 2022; 62:360-372. [PMID: 35412619 PMCID: PMC9788812 DOI: 10.1093/rheumatology/keac232] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Revised: 03/11/2022] [Accepted: 03/25/2022] [Indexed: 01/06/2023] Open
Abstract
OBJECTIVES To study the mechanism by which the readthrough mutation in TNFRSF11B, encoding osteoprotegerin (OPG) with additional 19 amino acids at its C-terminus (OPG-XL), causes the characteristic bidirectional phenotype of subchondral bone turnover accompanied by cartilage mineralization in chondrocalcinosis patients. METHODS OPG-XL was studied by human induced pluripotent stem cells expressing OPG-XL and two isogenic CRISPR/Cas9-corrected controls in cartilage and bone organoids. Osteoclastogenesis was studied with monocytes from OPG-XL carriers and matched healthy controls followed by gene expression characterization. Dual energy X-ray absorptiometry scans and MRI analyses were used to characterize the phenotype of carriers and non-carriers of the mutation. RESULTS Human OPG-XL carriers relative to sex- and age-matched controls showed, after an initial delay, large active osteoclasts with high number of nuclei. By employing hiPSCs expressing OPG-XL and isogenic CRISPR/Cas9-corrected controls to established cartilage and bone organoids, we demonstrated that expression of OPG-XL resulted in excessive fibrosis in cartilage and high mineralization in bone accompanied by marked downregulation of MGP, encoding matrix Gla protein, and upregulation of DIO2, encoding type 2 deiodinase, gene expression, respectively. CONCLUSIONS The readthrough mutation at CCAL1 locus in TNFRSF11B identifies an unknown role for OPG-XL in subchondral bone turnover and cartilage mineralization in humans via DIO2 and MGP functions. Previously, OPG-XL was shown to affect binding between RANKL and heparan sulphate (HS) resulting in loss of immobilized OPG-XL. Therefore, effects may be triggered by deficiency in the immobilization of OPG-XL Since the characteristic bidirectional pathophysiology of articular cartilage calcification accompanied by low subchondral bone mineralization is also a hallmark of OA pathophysiology, our results are likely extrapolated to common arthropathies.
Collapse
Affiliation(s)
| | | | | | - Eka Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden
| | - Ton Schoenmaker
- Department of Oral Cell Biology,Department of Periodontology, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and Vrije Universiteit , Amsterdam
| | | | | | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | | | | | | | - Teun J de Vries
- Department of Oral Cell Biology,Department of Periodontology, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and Vrije Universiteit , Amsterdam
| | - Yolande F M Ramos
- Correspondence to: Department of Molecular Epidemiology, Leiden University Medical Center, LUMC Postzone S-05-P, P.O. Box 9600, 2300 RC Leiden, The Netherlands. E-mail:
| | | |
Collapse
|
12
|
Houtman E, Tuerlings M, Suchiman HED, Lakenberg N, Cornelis FMF, Mei H, Broekhuis D, Nelissen RGHH, Coutinho de Almeida R, Ramos YFM, Lories RJ, Cruz LJ, Meulenbelt I. Inhibiting thyroid activation in aged human explants prevents mechanical induced detrimental signalling by mitigating metabolic processes. Rheumatology (Oxford) 2022; 62:457-466. [PMID: 35383365 PMCID: PMC9788824 DOI: 10.1093/rheumatology/keac202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 03/24/2022] [Accepted: 03/24/2022] [Indexed: 01/01/2023] Open
Abstract
OBJECTIVES To investigate whether the deiodinase inhibitor iopanoic acid (IOP) has chondroprotective properties, a mechanical stress induced model of human aged explants was used to test both repeated dosing and slow release of IOP. METHODS Human osteochondral explants subjected to injurious mechanical stress (65%MS) were treated with IOP or IOP encapsulated in poly lactic-co-glycolic acid-polyethylene glycol nanoparticles (NP-IOP). Changes to cartilage integrity and signalling were determined by Mankin scoring of histology, sulphated glycosaminoglycan (sGAG) release and expression levels of catabolic, anabolic and hypertrophic markers. Subsequently, on a subgroup of samples, RNA sequencing was performed on 65%MS (n = 14) and 65%MS+IOP (n = 7) treated cartilage to identify IOP's mode of action. RESULTS Damage from injurious mechanical stress was confirmed by increased cartilage surface damage in the Mankin score, increased sGAG release, and consistent upregulation of catabolic markers and downregulation of anabolic markers. IOP and, though less effective, NP-IOP treatment, reduced MMP13 and increased COL2A1 expression. In line with this, IOP and NP-IOP reduced cartilage surface damage induced by 65%MS, while only IOP reduced sGAG release from explants subjected to 65%MS. Lastly, differential expression analysis identified 12 genes in IOP's mode of action to be mainly involved in reducing metabolic processes (INSIG1, DHCR7, FADS1 and ACAT2) and proliferation and differentiation (CTGF, BMP5 and FOXM1). CONCLUSION Treatment with the deiodinase inhibitor IOP reduced detrimental changes of injurious mechanical stress. In addition, we identified that its mode of action was likely on metabolic processes, cell proliferation and differentiation.
Collapse
Affiliation(s)
- Evelyn Houtman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Margo Tuerlings
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - H Eka D Suchiman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Nico Lakenberg
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Frederique M F Cornelis
- Department of Development and Regeneration, Skeletal Biology and Engineering Research Centre, Laboratory of Tissue Homeostasis and Disease, KU Leuven, Leuven, Belgium
| | | | - Demiën Broekhuis
- Department of Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Rodrigo Coutinho de Almeida
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Rik J Lories
- Department of Development and Regeneration, Skeletal Biology and Engineering Research Centre, Laboratory of Tissue Homeostasis and Disease, KU Leuven, Leuven, Belgium,Division of Rheumatology, University Hospitals Leuven, Leuven, Belgium
| | - Luis J Cruz
- Translational Nanobiomaterials and Imaging (TNI) Group, Department of Radiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Ingrid Meulenbelt
- Correspondence to: Ingrid Meulenbelt, Molecular Epidemiology, Department of Biomedical Data Sciences Postzone J-11-R, Albinusdreef 2, 2333 ZA Leiden, The Netherlands. E-mail:
| |
Collapse
|
13
|
van Hoolwerff M, Rodríguez Ruiz A, Bouma M, Suchiman HED, Koning RI, Jost CR, Mulder AA, Freund C, Guilak F, Ramos YFM, Meulenbelt I. High-impact FN1 mutation decreases chondrogenic potential and affects cartilage deposition via decreased binding to collagen type II. Sci Adv 2021; 7:eabg8583. [PMID: 34739320 PMCID: PMC8570604 DOI: 10.1126/sciadv.abg8583] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Accepted: 09/17/2021] [Indexed: 06/13/2023]
Abstract
Osteoarthritis is the most prevalent joint disease worldwide, yet progress in development of effective disease-modifying treatments is slow because of lack of insight into the underlying disease pathways. Therefore, we aimed to identify the causal pathogenic mutation in an early-onset osteoarthritis family, followed by functional studies in human induced pluripotent stem cells (hiPSCs) in an in vitro organoid cartilage model. We demonstrated that the identified causal missense mutation in the gelatin-binding domain of the extracellular matrix protein fibronectin resulted in significant decreased binding capacity to collagen type II. Further analyses of formed hiPSC-derived neo-cartilage tissue highlighted that mutated fibronectin affected chondrogenic capacity and propensity to a procatabolic osteoarthritic state. Together, we demonstrate that binding of fibronectin to collagen type II is crucial for fibronectin downstream gene expression of chondrocytes. We advocate that effective treatment development should focus on restoring or maintaining proper binding between fibronectin and collagen type II.
Collapse
Affiliation(s)
- Marcella van Hoolwerff
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Alejandro Rodríguez Ruiz
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Marga Bouma
- LUMC hiPSC Hotel, Leiden University Medical Center, Leiden, Netherlands
- Department of Anatomy and Embryology, Leiden University Medical Center, Leiden, Netherlands
| | - H. Eka D. Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Roman I. Koning
- Section Electron Microscopy, Department of Cell and Chemical Biology, Leiden University Medical Center, Leiden, Netherlands
| | - Carolina R. Jost
- Section Electron Microscopy, Department of Cell and Chemical Biology, Leiden University Medical Center, Leiden, Netherlands
| | - Aat A. Mulder
- Section Electron Microscopy, Department of Cell and Chemical Biology, Leiden University Medical Center, Leiden, Netherlands
| | - Christian Freund
- LUMC hiPSC Hotel, Leiden University Medical Center, Leiden, Netherlands
- Department of Anatomy and Embryology, Leiden University Medical Center, Leiden, Netherlands
| | - Farshid Guilak
- Department of Orthopedic Surgery, Washington University and Shriners Hospitals for Children, St. Louis, MO, USA
| | - Yolande F. M. Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| |
Collapse
|
14
|
Tuerlings M, van Hoolwerff M, van Bokkum JM, Suchiman HED, Lakenberg N, Broekhuis D, Nelissen RGHH, Ramos YFM, Mei H, Cats D, Coutinho de Almeida R, Meulenbelt I. Long non-coding RNA expression profiling of subchondral bone reveals AC005165.1 modifying FRZB expression during osteoarthritis. Rheumatology (Oxford) 2021; 61:3023-3032. [PMID: 34730803 PMCID: PMC9258540 DOI: 10.1093/rheumatology/keab826] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Revised: 10/29/2021] [Indexed: 12/21/2022] Open
Abstract
Objective To gain insight in the expression profile of long non-coding RNAs (lncRNAs) in OA subchondral bone. Methods RNA sequencing data of macroscopically preserved and lesioned OA subchondral bone of patients that underwent joint replacement surgery due to OA (N = 22 pairs; 5 hips, 17 knees, Research osteoArthrits Articular Tissue (RAAK study) was run through an in-house pipeline to detect expression of lncRNAs. Differential expression analysis between preserved and lesioned bone was performed. Spearman correlations were calculated between differentially expressed lncRNAs and differentially expressed mRNAs identified previously in the same samples. Primary osteogenic cells were transfected with locked nucleic acid (LNA) GapmeRs targeting AC005165.1 lncRNA, to functionally investigate its potential mRNA targets. Results In total, 2816 lncRNAs were well-expressed in subchondral bone and we identified 233 lncRNAs exclusively expressed in knee and 307 lncRNAs exclusively in hip. Differential expression analysis, using all samples (N = 22 pairs; 5 hips, 17 knees), resulted in 21 differentially expressed lncRNAs [false discovery rate (FDR) < 0.05, fold change (FC) range 1.19–7.39], including long intergenic non-protein coding RNA (LINC) 1411 (LINC01411, FC = 7.39, FDR = 2.20 × 10−8), AC005165.1 (FC = 0.44, FDR = 2.37 × 10−6) and empty spiracles homeobox 2 opposite strand RNA (EMX2OS, FC = 0.41, FDR = 7.64 × 10−3). Among the differentially expressed lncRNAs, five were also differentially expressed in articular cartilage, including AC005165.1, showing similar direction of effect. Downregulation of AC005165.1 in primary osteogenic cells resulted in consistent downregulation of highly correlated frizzled related protein (FRZB). Conclusion The current study identified a novel lncRNA, AC005165.1, being dysregulated in OA articular cartilage and subchondral bone. Downregulation of AC005165.1 caused a decreased expression of OA risk gene FRZB, an important member of the wnt pathway, suggesting that AC005165.1 could be an attractive potential therapeutic target with effects in articular cartilage and subchondral bone.
Collapse
Affiliation(s)
- Margo Tuerlings
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Marcella van Hoolwerff
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Jessica M van Bokkum
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - H Eka D Suchiman
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Nico Lakenberg
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Demiën Broekhuis
- Dept. Orthopaedics Leiden, University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Dept. Orthopaedics Leiden, University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Hailiang Mei
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | - Davy Cats
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| | | | - Ingrid Meulenbelt
- Dept. of Biomedical Data Sciences, Leiden, The Netherlands, Leiden University Medical Center
| |
Collapse
|
15
|
Houtman E, Coutinho de Almeida R, Tuerlings M, Suchiman HED, Broekhuis D, Nelissen RGHH, Ramos YFM, van Meurs JBJ, Meulenbelt I. Characterization of dynamic changes in Matrix Gla Protein (MGP) gene expression as function of genetic risk alleles, osteoarthritis relevant stimuli, and the vitamin K inhibitor warfarin. Osteoarthritis Cartilage 2021; 29:1193-1202. [PMID: 33984465 DOI: 10.1016/j.joca.2021.05.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 04/28/2021] [Accepted: 05/05/2021] [Indexed: 02/02/2023]
Abstract
OBJECTIVE We here aimed to characterize changes of Matrix Gla Protein (MGP) expression in relation to its recently identified OA risk allele rs1800801-T in OA cartilage, subchondral bone and human ex vivo osteochondral explants subjected to OA related stimuli. Given that MGP function depends on vitamin K bioavailability, we studied the effect of frequently prescribed vitamin K antagonist warfarin. METHODS Differential (allelic) mRNA expression of MGP was analyzed using RNA-sequencing data of human OA cartilage and subchondral bone. Human osteochondral explants were used to study exposures to interleukin one beta (IL-1β; inflammation), triiodothyronine (T3; Hypertrophy), warfarin, or 65% mechanical stress (65%MS) as function of rs1800801 genotypes. RESULTS We confirmed that the MGP risk allele rs1800801-T was associated with lower expression and that MGP was significantly upregulated in lesioned as compared to preserved OA tissues, mainly in risk allele carriers, in both cartilage and subchondral bone. Moreover, MGP expression was downregulated in response to OA like triggers in cartilage and subchondral bone and this effect might be reduced in carriers of the rs1800801-T risk allele. Finally, warfarin treatment in cartilage increased COL10A1 and reduced SOX9 and MMP3 expression and in subchondral bone reduced COL1A1 and POSTN expression. DISCUSSION & CONCLUSIONS Our data highlights that the genetic risk allele lowers MGP expression and upon OA relevant triggers may hamper adequate dynamic changes in MGP expression, mainly in cartilage. The determined direct negative effect of warfarin on human explant cultures functionally underscores the previously found association between vitamin K deficiency and OA.
Collapse
Affiliation(s)
- E Houtman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands
| | - R Coutinho de Almeida
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands
| | - M Tuerlings
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands
| | - H E D Suchiman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands
| | - D Broekhuis
- Department of Orthopaedics, Leiden University Medical Center, Leiden, the Netherlands
| | - R G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden, the Netherlands
| | - Y F M Ramos
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands
| | - J B J van Meurs
- Department of Internal Medicine, Erasmus MC, University Medical Center, Rotterdam, the Netherlands
| | - I Meulenbelt
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, the Netherlands.
| |
Collapse
|
16
|
Ramos YFM, Meulenbelt I, Meijer J. Clock genes for joint health: if we could turn back time. Rheumatology (Oxford) 2021; 61:3-5. [PMID: 34260695 DOI: 10.1093/rheumatology/keab550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Affiliation(s)
- Yolande F M Ramos
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Ingrid Meulenbelt
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Joke Meijer
- Dept. Cellular and Chemical Biology, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
17
|
Coutinho de Almeida R, Mahfouz A, Mei H, Houtman E, den Hollander W, Soul J, Suchiman E, Lakenberg N, Meessen J, Huetink K, Nelissen RGHH, Ramos YFM, Reinders M, Meulenbelt I. Identification and characterization of two consistent osteoarthritis subtypes by transcriptome and clinical data integration. Rheumatology (Oxford) 2021; 60:1166-1175. [PMID: 32885253 PMCID: PMC7937023 DOI: 10.1093/rheumatology/keaa391] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 06/04/2020] [Indexed: 01/13/2023] Open
Abstract
OBJECTIVE 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.
Collapse
Affiliation(s)
- Rodrigo Coutinho de Almeida
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Ahmed Mahfouz
- Delft Bioinformatics Lab, Delft University of Technology, Delft, The Netherlands.,Leiden Computational Biology Center, Leiden, The Netherlands
| | - Hailiang Mei
- Sequence Analysis Support Core, Leiden University Medical Center, Leiden, The Netherlands
| | - Evelyn Houtman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Wouter den Hollander
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Jamie Soul
- Skeletal Research Group, Institute of Genetic Medicine, Newcastle University, Central Parkway, Newcastle upon Tyne, UK
| | - Eka Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Nico Lakenberg
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Jennifer Meessen
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands.,Department Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Kasper Huetink
- Department Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Marcel Reinders
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands.,Delft Bioinformatics Lab, Delft University of Technology, Delft, The Netherlands.,Leiden Computational Biology Center, Leiden, The Netherlands
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
18
|
Ruiz AR, Tuerlings M, Das A, de Almeida RC, Eka Suchiman H, Nelissen RGHH, Ramos YFM, Meulenbelt I. The role of TNFRSF11B in development of osteoarthritic cartilage. Rheumatology (Oxford) 2021; 61:856-864. [PMID: 33989379 PMCID: PMC8824428 DOI: 10.1093/rheumatology/keab440] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2021] [Revised: 05/12/2021] [Indexed: 11/13/2022] Open
Abstract
OBJECTIVES Osteoarthritis (OA) is a complex genetic disease with different risk factors contributing to its development. One of the genes, TNFRSF11B, previously identified with gain-of-function mutation in a family with early-onset OA with chondrocalcinosis, is among the highest upregulated genes in lesioned OA cartilage (RAAK-study). Here, we determined the role of TNFRSF11B overexpression in development of OA. METHODS Human primary articular chondrocytes (9 donors RAAK study) were transduced using lentiviral particles with or without TNFRSF11B. Cells were cultured for 1 week in a 3D in-vitro chondrogenic model . TNFRSF11B overexpression was confirmed by RT-qPCR, immunohistochemistry and ELISA. Effects of TNFRSF11B overexpression on cartilage matrix deposition, matrix mineralization, and genes highly correlated to TNFRSF11B in RNA-sequencing dataset (r>|0.75|) were determined by RT-qPCR. Additionally, glycosaminoglycans and collagen deposition were visualized with Alcian blue staining and immunohistochemistry (COL1 and COL2). RESULTS Overexpression of TNFRSF11B resulted in strong upregulation of MMP13, COL2A1 and COL1A1. Likewise, mineralization and osteoblast characteristic markers RUNX2, ASPN and OGN showed a consistent increase. Among 30 genes highly correlated to TNFRSF11B, expression of only 8 changed significantly, with BMP6 showing highest increase (9-fold) while expression of RANK and RANKL remained unchanged indicating previously unknown downstream pathways of TNFRSF11B in cartilage. CONCLUSION Results of our 3D in vitro chondrogenesis model indicate that upregulation of TNFRSF11B in lesioned OA cartilage may act as a direct driving factor for chondrocyte to osteoblast transition observed in OA pathophysiology. This transition does not appear to act via the OPG/RANK/RANKL triad common in bone remodeling. ETHICS APPROVAL AND CONSENT TO PARTICIPATE The Medical Ethics Committee of the LUMC gave approval for the RAAK study (P08.239). Written informed consent was obtained from all donors.
Collapse
Affiliation(s)
- Alejandro Rodríguez Ruiz
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Margo Tuerlings
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Ankita Das
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Rodrigo Coutinho de Almeida
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - H Eka Suchiman
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Rob G H H Nelissen
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Yolande F M Ramos
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| | - Ingrid Meulenbelt
- Dept. of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands; Dept. Orthopaedics, LUMC
| |
Collapse
|
19
|
Affiliation(s)
- Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.
| | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands
| | - S Rubina Baglio
- Department of Pathology, Cancer Center Amsterdam, Amsterdam University Medical Centers, Amsterdam, the Netherlands
| | - D Michiel Pegtel
- Department of Pathology, Cancer Center Amsterdam, Amsterdam University Medical Centers, Amsterdam, the Netherlands.
| |
Collapse
|
20
|
Houtman E, van Hoolwerff M, Lakenberg N, Suchiman EHD, van der Linden-van der Zwaag E, Nelissen RGHH, Ramos YFM, Meulenbelt I. Human Osteochondral Explants: Reliable Biomimetic Models to Investigate Disease Mechanisms and Develop Personalized Treatments for Osteoarthritis. Rheumatol Ther 2021; 8:499-515. [PMID: 33608843 PMCID: PMC7991015 DOI: 10.1007/s40744-021-00287-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Accepted: 01/30/2021] [Indexed: 02/07/2023] Open
Abstract
Introduction Likely due to ignored heterogeneity in disease pathophysiology, osteoarthritis (OA) has become the most common disabling joint disease, without effective disease-modifying treatment causing a large social and economic burden. In this study we set out to explore responses of aged human osteochondral explants upon different OA-related perturbing triggers (inflammation, hypertrophy and mechanical stress) for future tailored biomimetic human models. Methods Human osteochondral explants were treated with IL-1β (10 ng/ml) or triiodothyronine (T3; 10 nM) or received 65% strains of mechanical stress (65% MS). Changes in chondrocyte signalling were determined by expression levels of nine genes involved in catabolism, anabolism and hypertrophy. Breakdown of cartilage was measured by sulphated glycosaminoglycans (sGAGs) release, scoring histological changes (Mankin score) and mechanical properties of cartilage. Results All three perturbations (IL-1β, T3 and 65% MS) resulted in upregulation of the catabolic genes MMP13 and EPAS1. IL-1β abolished COL2A1 and ACAN gene expression and increased cartilage degeneration, reflected by increased Mankin scores and sGAGs released. Treatment with T3 resulted in a high and significant upregulation of the hypertrophic markers COL1A1, COL10A1 and ALPL. However, 65% MS increased sGAG release and detrimentally altered mechanical properties of cartilage. Conclusion We present consistent and specific output on three different triggers of OA. Perturbation with the pro-inflammatory IL-1β mainly induced catabolic chondrocyte signalling and cartilage breakdown, while T3 initiated expression of hypertrophic and mineralization markers. Mechanical stress at a strain of 65% induced catabolic chondrocyte signalling and changed cartilage matrix integrity. The major strength of our ex vivo models was that they considered aged, preserved, human cartilage of a heterogeneous OA patient population. As a result, the explants may reflect a reliable biomimetic model prone to OA onset allowing for development of different treatment modalities. Supplementary Information The online version contains supplementary material available at 10.1007/s40744-021-00287-y.
Collapse
Affiliation(s)
- Evelyn Houtman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Marcella van Hoolwerff
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Nico Lakenberg
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Eka H D Suchiman
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Rob G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Ingrid Meulenbelt
- Molecular Epidemiology, Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands.
| |
Collapse
|
21
|
Ramos YFM, Mobasheri A. MicroRNAs and Regulation of Autophagy in Chondrocytes. Methods Mol Biol 2021; 2245:179-194. [PMID: 33315203 DOI: 10.1007/978-1-0716-1119-7_13] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Chondrocytes are the main cells responsible for the maintenance of cartilage homeostasis and integrity. During development, extracellular matrix (ECM) macromolecules are produced and deposited by chondrocyte precursors. Autophagy, a highly dynamic process aimed at degradation of dysfunctional or pathogenic proteins, organelles, and intracellular microbes that can damage tissues, is one of the key processes required for sustained cartilage homeostasis. In different cell types it has been shown that, among others, autophagy is regulated by epigenetic mechanisms such as small noncoding RNAs (miRNAs, ~22 base pairs). Increasing evidence suggests that miRNAs are also involved in the regulation of autophagy in chondrocytes. Based on our previous research of gene and miRNA expression in articular cartilage, in this chapter we provide a summary of the tools models to direct in vitro and in vivo studies aimed at gaining a better understanding of the regulatory roles of miRNAs in chondrocyte autophagy.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands.
| | - Ali Mobasheri
- Department of Regenerative Medicine, State Research Institute Centre for Innovative Medicine, Vilnius, Lithuania.,Research Unit of Medical Imaging, Physics and Technology, Faculty of Medicine, University of Oulu, Oulu, Finland.,Centre for Sport, Exercise and Osteoarthritis Research Versus Arthritis, Queen's Medical Centre, Nottingham, UK.,Departments of Orthopedics, Rheumatology and Clinical Immunology, University Medical Centre Utrecht, Utrecht, Netherlands
| |
Collapse
|
22
|
Pelsma ICM, Claessen KMJA, Slagboom PE, van Heemst D, Pereira AM, Kroon HM, Ramos YFM, Kloppenburg M, Biermasz NR, Meulenbelt IM. Variants of FOXO3 and RPA3 genes affecting IGF-1 levels alter the risk of development of primary osteoarthritis. Eur J Endocrinol 2021; 184:29-39. [PMID: 33112260 DOI: 10.1530/eje-20-0904] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Accepted: 10/05/2020] [Indexed: 11/08/2022]
Abstract
INTRODUCTION Pathologically high growth hormone (GH) and insulin-like growth factor-1 (IGF-1) levels in patients with acromegaly are associated with arthropathy. Several studies highlight the potential role of the GH/IGF-1 axis in primary osteoarthritis (OA). We aimed to disentangle the role of IGF-1 levels in primary OA pathogenesis. METHODS Patients from the Genetics osteoARthritis and Progression (GARP) Study with familial, generalized, symptomatic OA (n = 337, mean age: 59.8 ± 7.4 years, 82% female) were compared to Leiden Longevity Study (LLS) controls (n = 456, mean age: 59.8 ± 6.8 years, 51% female). Subjects were clinically and radiographically assessed, serum IGF-1 levels were measured, and 10 quantitative trait loci (QTL) in the FOXO3, IGFBP3/TNS3, RPA3, SPOCK2 genes, previously related to serum IGF-1 levels, were genotyped. Linear or binary logistic generalized estimating equation models were performed. RESULTS Serum IGF-1 levels were increased in OA patients, with male patients exhibiting the strongest effect (males OR = 1.10 (1.04-1.17), P=0.002 vs females OR = 1.04 (1.01-1.07), P = 0.02). Independent of the increased IGF-1 levels, male carriers of the minor allele of FOXO3 QTL rs4946936 had a lower risk to develop hip OA (OR = 0.41 (0.18-0.90), P = 0.026). Additionally, independent of IGF-1 levels, female carriers of the minor alleles of RPA3 QTL rs11769597 had a higher risk to develop knee OA (OR = 1.90 (1.20-2.99), P = 0.006). CONCLUSION Patients with primary OA had significantly higher IGF-1 levels compared to controls. Moreover, SNPs in the FOXO3 and RPA3 genes were associated with an altered risk of OA. Therefore, altered IGF-1 levels affect the development of OA, and are potentially the result of the pathophysiological OA process.
Collapse
Affiliation(s)
- I C M Pelsma
- Division of Endocrinology, Department of Medicine and Center for Endocrine Tumors Leiden
| | - K M J A Claessen
- Division of Endocrinology, Department of Medicine and Center for Endocrine Tumors Leiden
| | - P E Slagboom
- Department of Biomedical Data Science, Section Molecular Epidemiology
| | - D van Heemst
- Department of Geriatrics and Gerontology, Leiden University Medical Center, Leiden, the Netherlands
| | - A M Pereira
- Division of Endocrinology, Department of Medicine and Center for Endocrine Tumors Leiden
| | - H M Kroon
- Department of Radiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Y F M Ramos
- Department of Biomedical Data Science, Section Molecular Epidemiology
| | - M Kloppenburg
- Department of Rheumatology, Leiden University Medical Center, Leiden, the Netherlands
- Department of Epidemiology, Leiden University Medical Center, Leiden, the Netherlands
| | - N R Biermasz
- Division of Endocrinology, Department of Medicine and Center for Endocrine Tumors Leiden
| | - I M Meulenbelt
- Department of Biomedical Data Science, Section Molecular Epidemiology
| |
Collapse
|
23
|
van Hoolwerff M, Metselaar PI, Tuerlings M, Suchiman HED, Lakenberg N, Ramos YFM, Cats D, Nelissen RGHH, Broekhuis D, Mei H, de Almeida RC, Meulenbelt I. Elucidating Epigenetic Regulation by Identifying Functional cis-Acting Long Noncoding RNAs and Their Targets in Osteoarthritic Articular Cartilage. Arthritis Rheumatol 2020; 72:1845-1854. [PMID: 32840049 PMCID: PMC7702083 DOI: 10.1002/art.41396] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 06/02/2020] [Indexed: 12/23/2022]
Abstract
Objective To identify robustly differentially expressed long noncoding RNAs (lncRNAs) with osteoarthritis (OA) pathophysiology in cartilage and to explore potential target messenger RNA (mRNA) by establishing coexpression networks, followed by functional validation. Methods RNA sequencing was performed on macroscopically lesioned and preserved OA cartilage from patients who underwent joint replacement surgery due to OA (n = 98). Differential expression analysis was performed on lncRNAs that were annotated in GENCODE and Ensembl databases. To identify potential interactions, correlations were calculated between the identified differentially expressed lncRNAs and the previously reported differentially expressed protein‐coding genes in the same samples. Modulation of chondrocyte lncRNA expression was achieved using locked nucleic acid GapmeRs. Results By applying our in‐house pipeline, we identified 5,053 lncRNAs that were robustly expressed, of which 191 were significantly differentially expressed (according to false discovery rate) between lesioned and preserved OA cartilage. Upon integrating mRNA sequencing data, we showed that intergenic and antisense differentially expressed lncRNAs demonstrate high, positive correlations with their respective flanking sense genes. To functionally validate this observation, we selected P3H2‐AS1, which was down‐regulated in primary chondrocytes, resulting in the down‐regulation of P3H2 gene expression levels. As such, we can confirm that P3H2‐AS1 regulates its sense gene P3H2. Conclusion By applying an improved detection strategy, robustly differentially expressed lncRNAs in OA cartilage were detected. Integration of these lncRNAs with differential mRNA expression levels in the same samples provided insight into their regulatory networks. Our data indicate that intergenic and antisense lncRNAs play an important role in regulating the pathophysiology of OA.
Collapse
Affiliation(s)
| | | | | | | | - Nico Lakenberg
- Leiden University Medical Center, Leiden, The Netherlands
| | | | - Davy Cats
- Leiden University Medical Center, Leiden, The Netherlands
| | | | | | - Hailiang Mei
- Leiden University Medical Center, Leiden, The Netherlands
| | | | | |
Collapse
|
24
|
Ross AK, Coutinho de Almeida R, Ramos YFM, Li J, Meulenbelt I, Guilak F. The miRNA-mRNA interactome of murine induced pluripotent stem cell-derived chondrocytes in response to inflammatory cytokines. FASEB J 2020; 34:11546-11561. [PMID: 32767602 DOI: 10.1096/fj.202000889r] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 06/03/2020] [Accepted: 06/05/2020] [Indexed: 12/24/2022]
Abstract
Osteoarthritis (OA) is a degenerative joint disease, and inflammation within an arthritic joint plays a critical role in disease progression. Pro-inflammatory cytokines, specifically IL-1 and TNF-α, induce aberrant expression of catabolic and degradative enzymes and inflammatory cytokines in OA and result in a challenging environment for cartilage repair and regeneration. MicroRNAs (miRNAS) are small noncoding RNAs and are important regulatory molecules that act by binding to target messenger RNAs (mRNAs) to reduce protein synthesis and have been implicated in many diseases, including OA. The goal of this study was to understand the mechanisms of miRNA regulation of the transcriptome of tissue-engineered cartilage in response to IL-1β and TNF-α using an in vitro murine induced pluripotent stem cell (miPSC) model system. We performed miRNA and mRNA sequencing to determine the temporal and dynamic responses of genes to specific inflammatory cytokines as well as miRNAs that are differentially expressed (DE) in response to both cytokines or exclusively to IL-1β or TNF-α. Through integration of mRNA and miRNA sequencing data, we created networks of miRNA-mRNA interactions which may be controlling the response to inflammatory cytokines. Within the networks, hub miRNAs, miR-29b-3p, miR-17-5p, and miR-20a-5p, were identified. As validation of these findings, we found that delivery of miR-17-5p and miR-20a-5p mimics significantly decreased degradative enzyme activity levels while also decreasing expression of inflammation-related genes in cytokine-treated cells. This study utilized an integrative approach to determine the miRNA interactome controlling the response to inflammatory cytokines and novel mediators of inflammation-driven degradation in tissue-engineered cartilage.
Collapse
Affiliation(s)
- Alison K Ross
- Department of Orthopaedic Surgery, Washington University, St. Louis, MO, USA.,Department of Biomedical Engineering, Washington University, St. Louis, MO, USA.,Shriners Hospitals for Children, St. Louis, MO, USA.,Center of Regenerative Medicine, Washington University, St. Louis, MO, USA
| | - Rodrigo Coutinho de Almeida
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Jiehan Li
- Department of Orthopaedic Surgery, Washington University, St. Louis, MO, USA.,Department of Biomedical Engineering, Washington University, St. Louis, MO, USA.,Shriners Hospitals for Children, St. Louis, MO, USA.,Center of Regenerative Medicine, Washington University, St. Louis, MO, USA
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Farshid Guilak
- Department of Orthopaedic Surgery, Washington University, St. Louis, MO, USA.,Department of Biomedical Engineering, Washington University, St. Louis, MO, USA.,Shriners Hospitals for Children, St. Louis, MO, USA.,Center of Regenerative Medicine, Washington University, St. Louis, MO, USA
| |
Collapse
|
25
|
van den Bosch MHJ, Ramos YFM, den Hollander W, Bomer N, Nelissen RGHH, Bovée JVMG, van den Berg WB, van Lent PLEM, Blom AB, van der Kraan PM, Meulenbelt I. Increased WISP1 expression in human osteoarthritic articular cartilage is epigenetically regulated and decreases cartilage matrix production. Rheumatology (Oxford) 2020; 58:1065-1074. [PMID: 30649473 DOI: 10.1093/rheumatology/key426] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2018] [Revised: 11/21/2018] [Indexed: 01/08/2023] Open
Abstract
OBJECTIVES Previously, we have shown the involvement of Wnt-activated protein Wnt-1-induced signaling protein 1 (WISP1) in the development of OA in mice. Here, we aimed to characterize the relation between WISP1 expression and human OA and its regulatory epigenetic determinants. METHODS Preserved and lesioned articular cartilage from end-stage OA patients and non-OA-diagnosed individuals was collected. WISP1 expression was determined using immunohistochemistry and damage was classified using Mankin scoring. RNA expression and DNA methylation were assessed in silico from genome-wide datasets (microarray analysis and RNA sequencing, and 450 k-methylationarrays, respectively). Effects of WISP1 were tested in pellet cultures of primary human chondrocytes. RESULTS WISP1 expression in cartilage of OA patients was increased compared with non-OA-diagnosed controls and, within OA patients, WISP1 was even higher in lesioned compared with preserved regions, with expression strongly correlating with Mankin score. In early symptomatic OA patients with disease progression, higher synovial WISP1 expression was observed as compared with non-progressors. Notably, increased WISP1 expression was inversely correlated with methylation levels of a positional CpG-dinucleotide (cg10191240), with lesioned areas showing strong hypomethylation for this CpG as compared with preserved cartilage. Additionally, we observed that methylation levels were allele-dependent for an intronic single-nucleotide polymorphism nearby cg10191240. Finally, addition of recombinant WISP1 to pellets of primary chondrocytes strongly inhibited deposition of extracellular matrix as reflected by decreased pellet circumference, proteoglycan content and decreased expression of matrix components. CONCLUSION Increased WISP1 expression is found in lesioned human articular cartilage, and appears epigenetically regulated via DNA methylation. In vitro assays suggest that increased WISP1 is detrimental for cartilage integrity.
Collapse
Affiliation(s)
| | - Yolande F M Ramos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Wouter den Hollander
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Nils Bomer
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Judith V M G Bovée
- Department of Pathology, Leiden University Medical Center, Leiden, The Netherlands
| | - Wim B van den Berg
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Peter L E M van Lent
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Arjen B Blom
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Peter M van der Kraan
- Experimental Rheumatology, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
26
|
den Hollander W, Pulyakhina I, Boer C, Bomer N, van der Breggen R, Arindrarto W, Couthino de Almeida R, Lakenberg N, Sentner T, Laros JFJ, ‘t Hoen PAC, Slagboom EPE, Nelissen RGHH, van Meurs J, Ramos YFM, Meulenbelt I. Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage. Arthritis Rheumatol 2019; 71:561-570. [PMID: 30298554 PMCID: PMC6593438 DOI: 10.1002/art.40748] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Accepted: 10/02/2018] [Indexed: 01/10/2023]
Abstract
OBJECTIVE Multiple single-nucleotide polymorphisms (SNPs) conferring susceptibility to osteoarthritis (OA) mark imbalanced expression of positional genes in articular cartilage, reflected by unequally expressed alleles among heterozygotes (allelic imbalance [AI]). We undertook this study to explore the articular cartilage transcriptome from OA patients for AI events to identify putative disease-driving genetic variation. METHODS AI was assessed in 42 preserved and 5 lesioned OA cartilage samples (from the Research Arthritis and Articular Cartilage study) for which RNA sequencing data were available. The count fraction of the alternative alleles among the alternative and reference alleles together (φ) was determined for heterozygous individuals. A meta-analysis was performed to generate a meta-φ and P value for each SNP with a false discovery rate (FDR) correction for multiple comparisons. To further validate AI events, we explored them as a function of multiple additional OA features. RESULTS We observed a total of 2,070 SNPs that consistently marked AI of 1,031 unique genes in articular cartilage. Of these genes, 49 were found to be significantly differentially expressed (fold change <0.5 or >2, FDR <0.05) between preserved and paired lesioned cartilage, and 18 had previously been reported to confer susceptibility to OA and/or related phenotypes. Moreover, we identified notable highly significant AI SNPs in the CRLF1, WWP2, and RPS3 genes that were related to multiple OA features. CONCLUSION We present a framework and resulting data set for researchers in the OA research field to probe for disease-relevant genetic variation that affects gene expression in pivotal disease-affected tissue. This likely includes putative novel compelling OA risk genes such as CRLF1, WWP2, and RPS3.
Collapse
Affiliation(s)
| | - Irina Pulyakhina
- Radboud University Medical Center Nijmegen, The Netherlands, and Wellcome Trust Centre for Human GeneticsOxfordUK
| | - Cindy Boer
- Erasmus Medical CenterRotterdamThe Netherlands
| | - Nils Bomer
- Leiden University Medical CenterLeidenThe Netherlands
| | | | | | | | | | - Thom Sentner
- Leiden University Medical CenterLeidenThe Netherlands
| | | | | | | | | | | | | | | |
Collapse
|
27
|
Coutinho de Almeida R, Ramos YFM, Meulenbelt I. Involvement of epigenetics in osteoarthritis. Best Pract Res Clin Rheumatol 2019; 31:634-648. [PMID: 30509410 DOI: 10.1016/j.berh.2018.03.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2018] [Accepted: 03/02/2018] [Indexed: 02/06/2023]
Abstract
Osteoarthritis (OA) is the most prevalent chronic age-related arthritic disease that mainly affects the diarthrodial joints. Nevertheless, there is no treatment currently available that can effectively reduce symptoms or slow down or stop disease progression. The lack of disease-modifying therapies could be explained by the complex pathogenesis of OA, which is still not completely understood. Intertwined epigenetic mechanisms such as DNA methylation, histone modifications, and noncoding RNAs (ncRNAs) have been indicated as important cellular tools to maintain tissue homeostasis upon environmental challenges. The current review illustrates that dysfunctional epigenetic control mechanisms in the articular cartilage likely play an important role in driving OA pathophysiology.
Collapse
Affiliation(s)
- Rodrigo Coutinho de Almeida
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Post-zone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands
| | - Yolande F M Ramos
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Post-zone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands
| | - Ingrid Meulenbelt
- Dept. Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Post-zone S-05-P, P.O. Box 9600, 2300 RC, Leiden, The Netherlands.
| |
Collapse
|
28
|
Coutinho de Almeida R, Ramos YFM, Mahfouz A, den Hollander W, Lakenberg N, Houtman E, van Hoolwerff M, Suchiman HED, Rodríguez Ruiz A, Slagboom PE, Mei H, Kiełbasa SM, Nelissen RGHH, Reinders M, Meulenbelt I. RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage. Ann Rheum Dis 2018; 78:270-277. [PMID: 30504444 PMCID: PMC6352405 DOI: 10.1136/annrheumdis-2018-213882] [Citation(s) in RCA: 109] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2018] [Revised: 11/09/2018] [Accepted: 11/12/2018] [Indexed: 01/24/2023]
Abstract
Objective To uncover the microRNA (miRNA) interactome of the osteoarthritis (OA) pathophysiological process in the cartilage. Methods We performed RNA sequencing in 130 samples (n=35 and n=30 pairs for messenger RNA (mRNA) and miRNA, respectively) on macroscopically preserved and lesioned OA cartilage from the same patient and performed differential expression (DE) analysis of miRNA and mRNAs. To build an OA-specific miRNA interactome, a prioritisation scheme was applied based on inverse Pearson’s correlations and inverse DE of miRNAs and mRNAs. Subsequently, these were filtered by those present in predicted (TargetScan/microT-CDS) and/or experimentally validated (miRTarBase/TarBase) public databases. Pathway enrichment analysis was applied to elucidate OA-related pathways likely mediated by miRNA regulatory mechanisms. Results We found 142 miRNAs and 2387 mRNAs to be differentially expressed between lesioned and preserved OA articular cartilage. After applying prioritisation towards likely miRNA-mRNA targets, a regulatory network of 62 miRNAs targeting 238 mRNAs was created. Subsequent pathway enrichment analysis of these mRNAs (or genes) elucidated that genes within the ‘nervous system development’ are likely mediated by miRNA regulatory mechanisms (familywise error=8.4×10−5). Herein NTF3 encodes neurotrophin-3, which controls survival and differentiation of neurons and which is closely related to the nerve growth factor. Conclusions By an integrated approach of miRNA and mRNA sequencing data of OA cartilage, an OA miRNA interactome and related pathways were elucidated. Our functional data demonstrated interacting levels at which miRNA affects expression of genes in the cartilage and exemplified the complexity of functionally validating a network of genes that may be targeted by multiple miRNAs.
Collapse
Affiliation(s)
- Rodrigo Coutinho de Almeida
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Yolande F M Ramos
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Ahmed Mahfouz
- The Delft Bioinformatics Lab, Delft University of Technology, Delft, The Netherlands.,Leiden Computational Biology Center, Leiden University Medical Center, Leiden, The Netherlands
| | - Wouter den Hollander
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Nico Lakenberg
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Evelyn Houtman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Marcella van Hoolwerff
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - H Eka D Suchiman
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Alejandro Rodríguez Ruiz
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - P Eline Slagboom
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Hailiang Mei
- Sequence Analysis Support Core, Leiden University Medical Center, Leiden, The Netherlands
| | - Szymon M Kiełbasa
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands.,Sequence Analysis Support Core, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Marcel Reinders
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands.,The Delft Bioinformatics Lab, Delft University of Technology, Delft, The Netherlands
| | - Ingrid Meulenbelt
- Department of Biomedical Data Sciences, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
29
|
den Hollander W, Boer CG, Hart DJ, Yau MS, Ramos YFM, Metrustry S, Broer L, Deelen J, Cupples LA, Rivadeneira F, Kloppenburg M, Peters M, Spector TD, Hofman A, Slagboom PE, Nelissen RGHH, Uitterlinden AG, Felson DT, Valdes AM, Meulenbelt I, van Meurs JJB. Genome-wide association and functional studies identify a role for matrix Gla protein in osteoarthritis of the hand. Ann Rheum Dis 2017; 76:2046-2053. [PMID: 28855172 PMCID: PMC5788019 DOI: 10.1136/annrheumdis-2017-211214] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2017] [Revised: 07/20/2017] [Accepted: 07/31/2017] [Indexed: 12/26/2022]
Abstract
OBJECTIVE Osteoarthritis (OA) is the most common form of arthritis and the leading cause of disability in the elderly. Of all the joints, genetic predisposition is strongest for OA of the hand; however, only few genetic risk loci for hand OA have been identified. Our aim was to identify novel genes associated with hand OA and examine the underlying mechanism. METHODS We performed a genome-wide association study of a quantitative measure of hand OA in 12 784 individuals (discovery: 8743, replication: 4011). Genome-wide significant signals were followed up by analysing gene and allele-specific expression in a RNA sequencing dataset (n=96) of human articular cartilage. RESULTS We found two significantly associated loci in the discovery set: at chr12 (p=3.5 × 10-10) near the matrix Gla protein (MGP) gene and at chr12 (p=6.1×10-9) near the CCDC91 gene. The DNA variant near the MGP gene was validated in three additional studies, which resulted in a highly significant association between the MGP variant and hand OA (rs4764133, Betameta=0.83, Pmeta=1.8*10-15). This variant is high linkage disequilibrium with a coding variant in MGP, a vitamin K-dependent inhibitor of cartilage calcification. Using RNA sequencing data from human primary cartilage tissue (n=96), we observed that the MGP RNA expression of the hand OA risk allele was significantly lowercompared with the MGP RNA expression of the reference allele (40.7%, p<5*10-16). CONCLUSIONS Our results indicate that the association between the MGP variant and increased risk for hand OA is caused by a lower expression of MGP, which may increase the burden of hand OA by decreased inhibition of cartilage calcification.
Collapse
Affiliation(s)
- Wouter den Hollander
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Cindy G Boer
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Deborah J Hart
- Department of Twin Research and Genetic Epidemiology, King’s College London, London, UK
| | - Michelle S Yau
- Institute for Aging Research, Hebrew SeniorLife, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
- Clinical Epidemiology Research and Training Unit, Boston University School of Medicine, Boston, Massachusetts, USA
| | - Yolande F M Ramos
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Sarah Metrustry
- Department of Twin Research and Genetic Epidemiology, King’s College London, London, UK
| | - Linda Broer
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Joris Deelen
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- Max Planck Institute for Biology of Ageing, Cologne, Germany
| | - L Adrienne Cupples
- Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, USA
| | - Fernando Rivadeneira
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Margreet Kloppenburg
- Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands
| | - Marjolein Peters
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Tim D Spector
- Department of Twin Research and Genetic Epidemiology, King’s College London, London, UK
| | - Albert Hofman
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
| | - P Eline Slagboom
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | - André G Uitterlinden
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - David T Felson
- Arthritis Research UK Epidemiology Unit, University of Manchester, Manchester, UK
| | - Ana M Valdes
- School of Medicine, University of Nottingham, Nottingham, UK
| | - Ingrid Meulenbelt
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Joyce J B van Meurs
- Department of Internal Medicine, Genetic Laboratory, Erasmus Medical Center, Rotterdam, The Netherlands
| |
Collapse
|
30
|
Steinberg J, Ritchie GRS, Roumeliotis TI, Jayasuriya RL, Clark MJ, Brooks RA, Binch ALA, Shah KM, Coyle R, Pardo M, Le Maitre CL, Ramos YFM, Nelissen RGHH, Meulenbelt I, McCaskie AW, Choudhary JS, Wilkinson JM, Zeggini E. Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis. Sci Rep 2017; 7:8935. [PMID: 28827734 PMCID: PMC5566454 DOI: 10.1038/s41598-017-09335-6] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2017] [Accepted: 07/17/2017] [Indexed: 11/09/2022] Open
Abstract
Osteoarthritis (OA) is a common disease characterized by cartilage degeneration and joint remodeling. The underlying molecular changes underpinning disease progression are incompletely understood. We investigated genes and pathways that mark OA progression in isolated primary chondrocytes taken from paired intact versus degraded articular cartilage samples across 38 patients undergoing joint replacement surgery (discovery cohort: 12 knee OA, replication cohorts: 17 knee OA, 9 hip OA patients). We combined genome-wide DNA methylation, RNA sequencing, and quantitative proteomics data. We identified 49 genes differentially regulated between intact and degraded cartilage in at least two -omics levels, 16 of which have not previously been implicated in OA progression. Integrated pathway analysis implicated the involvement of extracellular matrix degradation, collagen catabolism and angiogenesis in disease progression. Using independent replication datasets, we showed that the direction of change is consistent for over 90% of differentially expressed genes and differentially methylated CpG probes. AQP1, COL1A1 and CLEC3B were significantly differentially regulated across all three -omics levels, confirming their differential expression in human disease. Through integration of genome-wide methylation, gene and protein expression data in human primary chondrocytes, we identified consistent molecular players in OA progression that replicated across independent datasets and that have translational potential.
Collapse
Affiliation(s)
- Julia Steinberg
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.,Cancer Research Division, Cancer Council NSW, Sydney, NSW, 2011, Australia
| | - Graham R S Ritchie
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.,European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK.,Usher Institute of Population Health Sciences & Informatics, University of Edinburgh, Edinburgh, EH16 4UX, UK.,MRC Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Theodoros I Roumeliotis
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Raveen L Jayasuriya
- Department of Oncology and Metabolism, University of Sheffield, Beech Hill Road, Sheffield, S10 2RX, UK
| | - Matthew J Clark
- Department of Oncology and Metabolism, University of Sheffield, Beech Hill Road, Sheffield, S10 2RX, UK
| | - Roger A Brooks
- Division of Trauma & Orthopaedic Surgery, University of Cambridge, Box 180, Addenbrooke's Hospital, Hills Road, Cambridge, CB2 0QQ, UK
| | - Abbie L A Binch
- Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, S1 1WB, UK
| | - Karan M Shah
- Department of Oncology and Metabolism, University of Sheffield, Beech Hill Road, Sheffield, S10 2RX, UK
| | - Rachael Coyle
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Mercedes Pardo
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Christine L Le Maitre
- Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, S1 1WB, UK
| | - Yolande F M Ramos
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, 2300RC, The Netherlands
| | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, 2300RC, The Netherlands
| | - Ingrid Meulenbelt
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, 2300RC, The Netherlands
| | - Andrew W McCaskie
- Division of Trauma & Orthopaedic Surgery, University of Cambridge, Box 180, Addenbrooke's Hospital, Hills Road, Cambridge, CB2 0QQ, UK
| | - Jyoti S Choudhary
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - J Mark Wilkinson
- Department of Oncology and Metabolism, University of Sheffield, Beech Hill Road, Sheffield, S10 2RX, UK.
| | - Eleftheria Zeggini
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.
| |
Collapse
|
31
|
Castaño-Betancourt MC, Evans DS, Ramos YFM, Boer CG, Metrustry S, Liu Y, den Hollander W, van Rooij J, Kraus VB, Yau MS, Mitchell BD, Muir K, Hofman A, Doherty M, Doherty S, Zhang W, Kraaij R, Rivadeneira F, Barrett-Connor E, Maciewicz RA, Arden N, Nelissen RGHH, Kloppenburg M, Jordan JM, Nevitt MC, Slagboom EP, Hart DJ, Lafeber F, Styrkarsdottir U, Zeggini E, Evangelou E, Spector TD, Uitterlinden AG, Lane NE, Meulenbelt I, Valdes AM, van Meurs JBJ. Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis. PLoS Genet 2016; 12:e1006260. [PMID: 27701424 PMCID: PMC5049763 DOI: 10.1371/journal.pgen.1006260] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2016] [Accepted: 07/26/2016] [Indexed: 12/31/2022] Open
Abstract
Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10-8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.
Collapse
Affiliation(s)
| | - Dan S. Evans
- California Pacific Medical Center Research Institute, San Francisco, California, United States of America
| | - Yolande F. M. Ramos
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology. Leiden University Medical Center, Leiden, The Netherlands
| | - Cindy G. Boer
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Sarah Metrustry
- Department of Twins Research and Genetic Epidemiology Unit, King’s College London, London, United Kingdom
| | - Youfang Liu
- Thurston Arthritis Research Center, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Wouter den Hollander
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology. Leiden University Medical Center, Leiden, The Netherlands
| | - Jeroen van Rooij
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Virginia B. Kraus
- Duke Molecular Physiology Institute and Division of Rheumatology. Duke University School of Medicine, Durham, North Carolina, United States of America
| | - Michelle S. Yau
- Departments of Medicine and Epidemiology & Public Health, University of Maryland School of Medicine, Baltimore, Maryland, United States of America
| | - Braxton D. Mitchell
- Departments of Medicine and Epidemiology & Public Health, University of Maryland School of Medicine, Baltimore, Maryland, United States of America
- Geriatrics Research and Education Clinical Center, Baltimore Veterans Administration Medical Center, Baltimore, Maryland, United States of America
| | - Kenneth Muir
- Health Sciences Research Institute, University of Warwick, Warwick, United Kingdom
| | - Albert Hofman
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Harvard T.H. School of Public Health, Boston, Massachusetts, United States of America
| | - Michael Doherty
- School of Medicine, University of Nottingham, Nottingham, United Kingdom
| | - Sally Doherty
- School of Medicine, University of Nottingham, Nottingham, United Kingdom
| | - Weiya Zhang
- School of Medicine, University of Nottingham, Nottingham, United Kingdom
| | - Robert Kraaij
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Fernando Rivadeneira
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Elizabeth Barrett-Connor
- Epidemiology Division, Family Medicine and Public Health Department, University of California, San Diego, La Jolla, California
| | - Rose A. Maciewicz
- Respiratory, Inflammation, Autoimmunity Innovative Medicines, AstraZeneca AB, Mölndal, Sweden
| | - Nigel Arden
- Nuffield Department of Orthopaedics, Rheumatology and musculoskeletal sciences, University of Oxford, United Kingdom
| | - Rob G. H. H. Nelissen
- Department of Orthopaedics, Leiden University Medical Center, Leiden The Netherlands
| | - Margreet Kloppenburg
- Department of Rheumatology and Department of Clinical Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Joanne M. Jordan
- Thurston Arthritis Research Center, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Michael C. Nevitt
- University of California at San Francisco, San Francisco, California
| | - Eline P. Slagboom
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology. Leiden University Medical Center, Leiden, The Netherlands
| | - Deborah J. Hart
- Department of Twins Research and Genetic Epidemiology Unit, King’s College London, London, United Kingdom
| | - Floris Lafeber
- University Medical Center Utrecht, Utrecht, The Netherlands
| | | | | | - Evangelos Evangelou
- Department of Hygiene & Epidemiology, University of Ioannina School of Medicine, Ioannina, Greece
- Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, United Kingdom
| | - Tim D. Spector
- Department of Twins Research and Genetic Epidemiology Unit, King’s College London, London, United Kingdom
| | - Andre G. Uitterlinden
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Nancy E. Lane
- University of California at San Francisco, San Francisco, California
- School of Medicine, University of California, Davis, Sacramento, California
| | - Ingrid Meulenbelt
- Department of Medical Statistics and Bioinformatics, Section Molecular Epidemiology. Leiden University Medical Center, Leiden, The Netherlands
| | - Ana M. Valdes
- School of Medicine, University of Nottingham, Nottingham, United Kingdom
| | | |
Collapse
|
32
|
Bomer N, den Hollander W, Suchiman H, Houtman E, Slieker RC, Heijmans BT, Slagboom PE, Nelissen RGHH, Ramos YFM, Meulenbelt I. Neo-cartilage engineered from primary chondrocytes is epigenetically similar to autologous cartilage, in contrast to using mesenchymal stem cells. Osteoarthritis Cartilage 2016; 24:1423-30. [PMID: 26995110 DOI: 10.1016/j.joca.2016.03.009] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/22/2015] [Revised: 02/16/2016] [Accepted: 03/10/2016] [Indexed: 02/02/2023]
Abstract
OBJECTIVES To compare the epigenetic landscape of 3D cell models of human primary articular chondrocytes (hPACs) and human bone-marrow derived mesenchymal stem cells (hBMSCs) and their respective autologous articular cartilage. DESIGN Using Illumina Infinium HumanMethylation450 BeadChip arrays, the DNA methylation landscape of the different cell sources and autologous cartilage was determined. Pathway enrichment was analyzed using DAVID. RESULTS Principal Component Analysis (PCA) of methylation data revealed separate clustering of hBMSC samples. Between hBMSCs and autologous cartilage 86,881 cytosine-phosphate-guanine dinucleotides (CpGs) (20.2%), comprising 3,034 differentially methylated regions (DMRs; Δβ > 0.1; with the same direction of effect), were significantly differentially methylated. In contrast, between hPACs and autologous cartilage only 5,706 CpGs (1.33%) were differentially methylated. Of interest was the finding of the transcriptionally active, hyper-methylation of a Cartilage Intermediate Layer Protein (CILP) annotated DMR (Δβ = 0.16) in PAC-cartilage, corresponding to a profound decrease in CILP expression after in vitro culturing of hPACs as compared to autologous cartilage. CONCLUSIONS In vitro engineered neo-cartilage tissue from primary chondrocytes, hPACs, exhibits a DNA methylation landscape that is almost identical (99% similarity) to autologous cartilage, in contrast to neo-cartilage engineered from bone marrow-derived mesenchymal stem cells (MSCs). Although hBMSCs are widely used for cartilage engineering purposes the effects of these vast differences on cartilage regeneration and long term consequences of implantation, are not known. The use of hBMSCs or hPACs for future cartilage tissue regeneration purposes should therefore be investigated in more depth in future endeavors to better understand the consequences of the differential methylome on neo-cartilage.
Collapse
Affiliation(s)
- N Bomer
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands; IDEAL, LUMC, Leiden, The Netherlands
| | - W den Hollander
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - H Suchiman
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - E Houtman
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - R C Slieker
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands; IDEAL, LUMC, Leiden, The Netherlands
| | - B T Heijmans
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands; IDEAL, LUMC, Leiden, The Netherlands
| | - P E Slagboom
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands; IDEAL, LUMC, Leiden, The Netherlands
| | | | - Y F M Ramos
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - I Meulenbelt
- Dept. of Molecular Epidemiology, LUMC, Leiden, The Netherlands.
| |
Collapse
|
33
|
Bomer N, Cornelis FMF, Ramos YFM, den Hollander W, Lakenberg N, van der Breggen R, Storms L, Slagboom PE, Lories RJU, Meulenbelt I. Aberrant Calreticulin Expression in Articular Cartilage of Dio2 Deficient Mice. PLoS One 2016; 11:e0154999. [PMID: 27163789 PMCID: PMC4862667 DOI: 10.1371/journal.pone.0154999] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2016] [Accepted: 04/22/2016] [Indexed: 11/19/2022] Open
Abstract
OBJECTIVE To identify intrinsic differences in cartilage gene expression profiles between wild-type- and Dio2-/--mice, as a mechanism to investigate factors that contribute to prolonged healthy tissue homeostasis. METHODS Previously generated microarray-data (Illumina MouseWG-6 v2) of knee cartilage of wild-type and Dio2 -/- -mice were re-analyzed to identify differential expressed genes independent of mechanical loading conditions by forced treadmill-running. RT-qPCR and western blot analyses of overexpression and knockdown of Calr in mouse chondro-progenitor cells (ATDC5) were applied to assess the direct effect of differential Calr expression on cartilage deposition. RESULTS Differential expression analyses of articular cartilage of Dio2-/- (N = 9) and wild-type-mice (N = 11) while applying a cutoff threshold (P < 0.05 (FDR) and FC > |1,5|) resulted in 1 probe located in Calreticulin (Calr) that was found significantly downregulated in Dio2-/- mice (FC = -1.731; P = 0.044). Furthermore, overexpression of Calr during early chondrogenesis in ATDC5 cells leads to decreased proteoglycan deposition and corresponding lower Aggrecan expression, whereas knocking down Calr expression does not lead to histological differences of matrix composition. CONCLUSION We here demonstrate that the beneficial homeostatic state of articular cartilage in Dio2-/- mice is accompanied with significant lower expression of Calr. Functional analyses further showed that upregulation of Calr expression could act as an initiator of cartilage destruction. The consistent association between Calr and Dio2 expression suggests that enhanced expression of these genes facilitate detrimental effects on cartilage integrity.
Collapse
Affiliation(s)
- Nils Bomer
- Department of Molecular Epidemiology, LUMC, Leiden, Netherlands
| | - Frederique M. F. Cornelis
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium
| | | | | | - Nico Lakenberg
- Department of Molecular Epidemiology, LUMC, Leiden, Netherlands
| | | | - Lies Storms
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium
| | | | - Rik J. U. Lories
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium
- Division of Rheumatology, University Hospitals Leuven, Leuven, Belgium
| | | |
Collapse
|
34
|
Abstract
Osteoarthritis (OA) is the most common age-related arthritic disorder and is characterized by aberrant extracellular matrix (ECM) content and surface disruptions that range from fibrillation, clefting and delamination, leading to articular surface erosion. Worldwide, over 20% of the population is affected with OA and 80% of these patients have limitations in movement, whereas 25% experience inhibition in major daily activities of life. OA is the most common disabling arthritic disease; nevertheless, no disease-modifying treatment is available except for the expensive total joint replacement surgery at end-stage disease. Lack of insight into the underlying pathophysiological mechanisms of OA has considerably contributed to the inability of the scientific community to develop disease-modifying drugs. To overcome this critical barrier, focus should be on translation of identified robust gene deviations towards the underlying biological mechanisms.
Collapse
Affiliation(s)
- Nils Bomer
- Dept. Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Centre, LUMC Post-zone S-05-P, P.O. Box 9600, 2300 RC Leiden, The Netherlands.
| | - Wouter den Hollander
- Dept. Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Centre, LUMC Post-zone S-05-P, P.O. Box 9600, 2300 RC Leiden, The Netherlands.
| | - Yolande F M Ramos
- Dept. Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Centre, LUMC Post-zone S-05-P, P.O. Box 9600, 2300 RC Leiden, The Netherlands.
| | - Ingrid Meulenbelt
- Dept. Medical Statistics and Bioinformatics, Section Molecular Epidemiology, Leiden University Medical Centre, LUMC Post-zone S-05-P, P.O. Box 9600, 2300 RC Leiden, The Netherlands.
| |
Collapse
|
35
|
Bomer N, Cornelis FMF, Ramos YFM, den Hollander W, Storms L, van der Breggen R, Lakenberg N, Slagboom PE, Meulenbelt I, Lories RJL. The effect of forced exercise on knee joints in Dio2(-/-) mice: type II iodothyronine deiodinase-deficient mice are less prone to develop OA-like cartilage damage upon excessive mechanical stress. Ann Rheum Dis 2016; 75:571-7. [PMID: 25550340 DOI: 10.1136/annrheumdis-2014-206608] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2014] [Accepted: 12/07/2014] [Indexed: 01/10/2023]
Abstract
OBJECTIVE To further explore deiodinase iodothyronine type 2 (DIO2) as a therapeutic target in osteoarthritis (OA) by studying the effects of forced mechanical loading on in vivo joint cartilage tissue homeostasis and the modulating effect herein of Dio2 deficiency. METHODS Wild-type and C57BL/6-Dio2(-/-) -mice were subjected to a forced running regime for 1 h per day for 3 weeks. Severity of OA was assessed by histological scoring for cartilage damage and synovitis. Genome-wide gene expression was determined in knee cartilage by microarray analysis (Illumina MouseWG-6 v2). STRING-db analyses were applied to determine enrichment for specific pathways and to visualise protein-protein interactions. RESULTS In total, 158 probes representing 147 unique genes showed significantly differential expression with a fold-change ≥1.5 upon forced exercise. Among these are genes known for their association with OA (eg, Mef2c, Egfr, Ctgf, Prg4 and Ctnnb1), supporting the use of forced running as an OA model in mice. Dio2-deficient mice showed significantly less cartilage damage and signs of synovitis. Gene expression response upon exercise between wild-type and knockout mice was significantly different for 29 genes. CONCLUSIONS Mice subjected to a running regime have significant increased cartilage damage and synovitis scores. Lack of Dio2 protected against cartilage damage in this model and was reflected in a specific gene expression profile, and either mark a favourable effect in the Dio2 knockout (eg, Gnas) or an unfavourable effect in wild-type cartilage homeostasis (eg, Hmbg2 and Calr). These data further support DIO2 activity as a therapeutic target in OA.
Collapse
MESH Headings
- Animals
- Cartilage, Articular/metabolism
- Cartilage, Articular/pathology
- Gene Expression Profiling
- Iodide Peroxidase/genetics
- Knee Joint/metabolism
- Knee Joint/pathology
- Mice
- Mice, Inbred C57BL
- Mice, Knockout
- Oligonucleotide Array Sequence Analysis
- Osteoarthritis, Knee/genetics
- Osteoarthritis, Knee/metabolism
- Osteoarthritis, Knee/pathology
- Physical Conditioning, Animal
- RNA, Messenger/metabolism
- Real-Time Polymerase Chain Reaction
- Stress, Mechanical
- Iodothyronine Deiodinase Type II
Collapse
Affiliation(s)
- Nils Bomer
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands Integrated Research of Developmental Determinants of Ageing and Longevity (IDEAL), Leiden, Netherlands
| | - Frederique M F Cornelis
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium
| | | | | | - Lies Storms
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium
| | | | - Nico Lakenberg
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - P Eline Slagboom
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands Integrated Research of Developmental Determinants of Ageing and Longevity (IDEAL), Leiden, Netherlands The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Rik J L Lories
- Laboratory of Tissue Homeostasis and Disease, Skeletal Biology and Engineering Research Centre, KU Leuven, Leuven, Belgium Division of Rheumatology, University Hospitals Leuven, Leuven, Belgium
| |
Collapse
|
36
|
Peters MJ, Ramos YFM, den Hollander W, Schiphof D, Hofman A, Uitterlinden AG, Oei EHG, Slagboom PE, Kloppenburg M, Bloem JL, Bierma-Zeinstra SMA, Meulenbelt I, van Meurs JBJ. Associations between joint effusion in the knee and gene expression levels in the circulation: a meta-analysis. F1000Res 2016; 5:109. [PMID: 27134727 PMCID: PMC4837985 DOI: 10.12688/f1000research.7763.1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 01/19/2016] [Indexed: 01/16/2023] Open
Abstract
Objective: To identify molecular biomarkers for early knee osteoarthritis (OA), we examined whether joint effusion in the knee associated with different gene expression levels in the circulation. Materials and Methods: Joint effusion grades measured with magnetic resonance (MR) imaging and gene expression levels in blood were determined in women of the Rotterdam Study (N=135) and GARP (N=98). Associations were examined using linear regression analyses, adjusted for age, fasting status, RNA quality, technical batch effects, blood cell counts, and BMI. To investigate enriched pathways and protein-protein interactions, we used the DAVID and STRING webtools. Results: In a meta-analysis, we identified 257 probes mapping to 189 unique genes in blood that were nominally significantly associated with joint effusion grades in the knee. Several compelling genes were identified such as
C1orf38 and
NFATC1. Significantly enriched biological pathways were: response to stress, gene expression, negative regulation of intracellular signal transduction, and antigen processing and presentation of exogenous pathways. Conclusion: Meta-analyses and subsequent enriched biological pathways resulted in interesting candidate genes associated with joint effusion that require further characterization. Associations were not transcriptome-wide significant most likely due to limited power. Additional studies are required to replicate our findings in more samples, which will greatly help in understanding the pathophysiology of OA and its relation to inflammation, and may result in biomarkers urgently needed to diagnose OA at an early stage.
Collapse
Affiliation(s)
| | - Yolande F M Ramos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Wouter den Hollander
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Dieuwke Schiphof
- Department of General Practice, Erasmus MC, Rotterdam, Netherlands
| | - Albert Hofman
- Department of Epidemiology, Erasmus MC, Rotterdam, Netherlands
| | - André G Uitterlinden
- Department of Internal Medicine, Erasmus MC, Rotterdam, Netherlands; Department of Epidemiology, Erasmus MC, Rotterdam, Netherlands
| | - Edwin H G Oei
- Department of Radiology, Erasmus MC University Medical Center, Rotterdam, Netherlands
| | - P Eline Slagboom
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | - Margreet Kloppenburg
- Department of Clinical Epidemiology and Rheumatology, Leiden University Medical Center, Leiden, Netherlands
| | - Johan L Bloem
- Department of Radiology, Leiden University Medical Center, Leiden, Netherlands
| | - Sita M A Bierma-Zeinstra
- Department of General Practice, Erasmus MC, Rotterdam, Netherlands; Department of Orthopedics, Erasmus MC, Rotterdam, Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, Netherlands
| | | |
Collapse
|
37
|
Bomer N, den Hollander W, Ramos YFM, Bos SD, van der Breggen R, Lakenberg N, Pepers BA, van Eeden AE, Darvishan A, Tobi EW, Duijnisveld BJ, van den Akker EB, Heijmans BT, van Roon-Mom WMC, Verbeek FJ, van Osch GJVM, Nelissen RGHH, Slagboom PE, Meulenbelt I. Underlying molecular mechanisms of DIO2 susceptibility in symptomatic osteoarthritis. Ann Rheum Dis 2015; 74:1571-9. [PMID: 24695009 PMCID: PMC4516000 DOI: 10.1136/annrheumdis-2013-204739] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2013] [Revised: 03/14/2014] [Accepted: 03/15/2014] [Indexed: 11/30/2022]
Abstract
OBJECTIVES To investigate how the genetic susceptibility gene DIO2 confers risk to osteoarthritis (OA) onset in humans and to explore whether counteracting the deleterious effect could contribute to novel therapeutic approaches. METHODS Epigenetically regulated expression of DIO2 was explored by assessing methylation of positional CpG-dinucleotides and the respective DIO2 expression in OA-affected and macroscopically preserved articular cartilage from end-stage OA patients. In a human in vitro chondrogenesis model, we measured the effects when thyroid signalling during culturing was either enhanced (excess T3 or lentiviral induced DIO2 overexpression) or decreased (iopanoic acid). RESULTS OA-related changes in methylation at a specific CpG dinucleotide upstream of DIO2 caused significant upregulation of its expression (β=4.96; p=0.0016). This effect was enhanced and appeared driven specifically by DIO2 rs225014 risk allele carriers (β=5.58, p=0.0006). During in vitro chondrogenesis, DIO2 overexpression resulted in a significant reduced capacity of chondrocytes to deposit extracellular matrix (ECM) components, concurrent with significant induction of ECM degrading enzymes (ADAMTS5, MMP13) and markers of mineralisation (ALPL, COL1A1). Given their concurrent and significant upregulation of expression, this process is likely mediated via HIF-2α/RUNX2 signalling. In contrast, we showed that inhibiting deiodinases during in vitro chondrogenesis contributed to prolonged cartilage homeostasis as reflected by significant increased deposition of ECM components and attenuated upregulation of matrix degrading enzymes. CONCLUSIONS Our findings show how genetic variation at DIO2 could confer risk to OA and raised the possibility that counteracting thyroid signalling may be a novel therapeutic approach.
Collapse
Affiliation(s)
- Nils Bomer
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- IDEAL, The Netherlands
| | | | | | - Steffan D Bos
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- Genomics Initiative, sponsored by the NCHA, Leiden, The Netherlands
| | | | - Nico Lakenberg
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
| | - Barry A Pepers
- Department of Human Genetics, LUMC, Leiden, The Netherlands
| | | | - Arash Darvishan
- Department of Imaging & BioInformatics, LIACS, Leiden, The Netherlands
| | - Elmar W Tobi
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- IDEAL, The Netherlands
| | | | - Erik B van den Akker
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- The Delft Bioinformatics Lab, Delft University of Technology, Delft, The Netherlands
| | - Bastiaan T Heijmans
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- Genomics Initiative, sponsored by the NCHA, Leiden, The Netherlands
| | | | - Fons J Verbeek
- Department of Imaging & BioInformatics, LIACS, Leiden, The Netherlands
| | - Gerjo J V M van Osch
- Department of Orthopaedics, Erasmus MC, Rotterdam, The Netherlands
- Deptartment of Otorhinolaryngology, Erasmus MC, Rotterdam, The Netherlands
| | | | - P Eline Slagboom
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- IDEAL, The Netherlands
- Genomics Initiative, sponsored by the NCHA, Leiden, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands
- Genomics Initiative, sponsored by the NCHA, Leiden, The Netherlands
| |
Collapse
|
38
|
den Hollander W, Ramos YFM, Bomer N, Elzinga S, van der Breggen R, Lakenberg N, de Dijcker WJ, Suchiman HED, Duijnisveld BJ, Houwing-Duistermaat JJ, Slagboom PE, Bos SD, Nelissen RGHH, Meulenbelt I. Transcriptional Associations of Osteoarthritis-Mediated Loss of Epigenetic Control in Articular Cartilage. Arthritis Rheumatol 2015; 67:2108-16. [DOI: 10.1002/art.39162] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2014] [Accepted: 04/14/2015] [Indexed: 12/31/2022]
Affiliation(s)
| | | | - Nils Bomer
- Leiden University Medical Center; Leiden The Netherlands
| | - Stefan Elzinga
- Leiden University Medical Center; Leiden The Netherlands
| | | | - Nico Lakenberg
- Leiden University Medical Center; Leiden The Netherlands
| | | | | | | | | | - P. Eline Slagboom
- Leiden University Medical Center, Leiden, The Netherlands, and The Netherlands Genomics Initiative-sponsored Netherlands Consortium for Healthy Aging; Rotterdam The Netherlands
| | - Steffan D. Bos
- Leiden University Medical Center, Leiden, The Netherlands, and The Netherlands Genomics Initiative-sponsored Netherlands Consortium for Healthy Aging; Rotterdam The Netherlands
| | | | - Ingrid Meulenbelt
- Leiden University Medical Center, Leiden, The Netherlands, and The Netherlands Genomics Initiative-sponsored Netherlands Consortium for Healthy Aging; Rotterdam The Netherlands
| |
Collapse
|
39
|
|
40
|
Evangelou E, Kerkhof HJ, Styrkarsdottir U, Ntzani EE, Bos SD, Esko T, Evans DS, Metrustry S, Panoutsopoulou K, Ramos YFM, Thorleifsson G, Tsilidis KK, Arden N, Aslam N, Bellamy N, Birrell F, Blanco FJ, Carr A, Chapman K, Day-Williams AG, Deloukas P, Doherty M, Engström G, Helgadottir HT, Hofman A, Ingvarsson T, Jonsson H, Keis A, Keurentjes JC, Kloppenburg M, Lind PA, McCaskie A, Martin NG, Milani L, Montgomery GW, Nelissen RGHH, Nevitt MC, Nilsson PM, Ollier WER, Parimi N, Rai A, Ralston SH, Reed MR, Riancho JA, Rivadeneira F, Rodriguez-Fontenla C, Southam L, Thorsteinsdottir U, Tsezou A, Wallis GA, Wilkinson JM, Gonzalez A, Lane NE, Lohmander LS, Loughlin J, Metspalu A, Uitterlinden AG, Jonsdottir I, Stefansson K, Slagboom PE, Zeggini E, Meulenbelt I, Ioannidis JPA, Spector TD, van Meurs JBJ, Valdes AM. A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip. Ann Rheum Dis 2014; 73:2130-6. [PMID: 23989986 PMCID: PMC4251181 DOI: 10.1136/annrheumdis-2012-203114] [Citation(s) in RCA: 94] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Accepted: 07/26/2013] [Indexed: 12/15/2022]
Abstract
OBJECTIVES Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS Novel genetic loci for hip OA were found in this meta-analysis of GWAS.
Collapse
Affiliation(s)
- Evangelos Evangelou
- Department of Hygiene and Epidemiology, University of Ioannina Medical School, Ioannina, Greece
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - Hanneke J Kerkhof
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | | | - Evangelia E Ntzani
- Department of Hygiene and Epidemiology, University of Ioannina Medical School, Ioannina, Greece
| | - Steffan D Bos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- Netherlands Consortium for Healthy Ageing, The Netherlands
| | - Tonu Esko
- Estonian Genome Center, University of Tartu, Tartu, Estonia
- Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia
| | - Daniel S Evans
- California Pacific Medical Center Research Institute, San Francisco, USA
| | - Sarah Metrustry
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | | | - Yolande F M Ramos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Konstantinos K Tsilidis
- Department of Hygiene and Epidemiology, University of Ioannina Medical School, Ioannina, Greece
| | | | - Nigel Arden
- NIHR Biomedical Research Unit and ARUK Centre of excellence for Sport, Exercise and Osteoarthritis, University of Oxford, Oxford, UK
- MRC Epidemiology Resource Centre, University of Southampton, Southampton, UK
| | - Nadim Aslam
- Worcestershire Royal Hospital, Worcestershire Acute Hospitals NHS Trust, Worcester, UK
| | - Nicholas Bellamy
- Centre of National Research on Disability and Rehabilitation Medicine, The University of Queensland, Brisbane, Australia
| | - Fraser Birrell
- Musculoskeletal Research Group, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, UK
- Wansbeck General Hospital, Northumbria Healthcare NHS Foundation Trust, Ashington, UK
| | - Francisco J Blanco
- Rheumatology Division, Instituto de Investigación Biomédica-Hospital Universitario A Coruña, A Corunna, Spain
| | - Andrew Carr
- NIHR Biomedical Research Unit and ARUK Centre of excellence for Sport, Exercise and Osteoarthritis, University of Oxford, Oxford, UK
| | - Kay Chapman
- NIHR Biomedical Research Unit and ARUK Centre of excellence for Sport, Exercise and Osteoarthritis, University of Oxford, Oxford, UK
| | | | - Panos Deloukas
- Department of Human Genetics, Wellcome Trust Sanger Institute, Hinxton, UK
| | - Michael Doherty
- Department of Academic Rheumatology, University of Nottingham, Nottingham, UK
| | - Gunnar Engström
- Department of Clinical Sciences Malmo, Lund University, Malmo, Sweden
| | | | - Albert Hofman
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Thorvaldur Ingvarsson
- Department of Orthopedic Surgery, Akureyri Hospital, Akureyri, Iceland
- School of Health Sciences, University of Akureyri, Akureyri, Iceland
| | - Helgi Jonsson
- Department of Medicine, The National University Hospital of Iceland, Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Reykjavik, Iceland
| | - Aime Keis
- Department of Public Health, University of Tartu, Tartu, Estonia
- Orthopedic Surgeons, Elva Hospital, Elva, Estonia
| | | | - Margreet Kloppenburg
- Department of Rheumatology, Leiden University Medical Center, Leiden, The Netherlands
| | - Penelope A Lind
- Department of Quantitative Genetics, QIMR Berghofer Medical Research Institute, Brisbane, Australia
| | - Andrew McCaskie
- Musculoskeletal Research Group, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, UK
| | - Nicholas G Martin
- Department of Genetic Epidemiology, QIMR Berghofer Medical Research Institute, Brisbane, Australia
| | - Lili Milani
- Estonian Genome Center, University of Tartu, Tartu, Estonia
| | - Grant W Montgomery
- Department of Molecular Epidemiology, QIMR Berghofer Medical Research Institute, Brisbane, Australia
| | - Rob G H H Nelissen
- Department of Orthopedics, Leiden University Medical Center, Leiden, The Netherlands
| | - Michael C Nevitt
- Department of Epidemiology and Biostatistics, University of California, San Francisco, San Francisco, California, USA
| | - Peter M Nilsson
- Department of Clinical Sciences Malmo, Lund University, Malmo, Sweden
| | - William ER Ollier
- Centre for Integrated Genomic Medical Research, University of Manchester, Manchester, UK
| | - Neeta Parimi
- California Pacific Medical Center Research Institute, San Francisco, USA
| | - Ashok Rai
- Worcestershire Acute Hospitals NHS Trust, Worcester, UK
| | - Stuart H Ralston
- Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, UK
| | - Mike R Reed
- Wansbeck General Hospital, Northumbria Healthcare NHS Foundation Trust, Ashington, UK
| | - Jose A Riancho
- Department of Internal Medicine, Hospital U.M. Valdecilla-IFIMAV, University of Cantabria, Santander, Spain
| | - Fernando Rivadeneira
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Cristina Rodriguez-Fontenla
- Laboratorio Investigacion 10 and Rheumatology Unit, Instituto de Investigacion Sanitaria—Hospital Clinico Universitario de Santiago, Santiago de Compostela, Spain
| | - Lorraine Southam
- Department of Human Genetics, Wellcome Trust Sanger Institute, Hinxton, UK
| | - Unnur Thorsteinsdottir
- Department of Population Genetics, deCODE Genetics, Reykjavik, Iceland
- Department of Medicine, The National University Hospital of Iceland, Reykjavik, Iceland
| | - Aspasia Tsezou
- Department of Biology, University of Thessaly, Medical School, Larissa, Greece
| | - Gillian A Wallis
- Wellcome Trust Centre for Cell-Matrix Research, University of Manchester, Manchester, UK
| | - J Mark Wilkinson
- Department of Human Metabolism, University of Sheffield, Sheffield, UK
| | - Antonio Gonzalez
- Laboratorio Investigacion 10 and Rheumatology Unit, Instituto de Investigacion Sanitaria—Hospital Clinico Universitario de Santiago, Santiago de Compostela, Spain
| | - Nancy E Lane
- Department of Medicine, University of California at Davis, Sacramento, USA
| | - L Stefan Lohmander
- Research Unit for Musculoskeletal Function and Physiotherapy, and Department of Orthopedics and Traumatology, University of Southern Denmark, Odense, Denmark
- Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| | - John Loughlin
- Musculoskeletal Research Group, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, UK
| | - Andres Metspalu
- Estonian Genome Center, University of Tartu, Tartu, Estonia
- Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia
| | - Andre G Uitterlinden
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
- Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Ingileif Jonsdottir
- Department of Population Genetics, deCODE Genetics, Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Reykjavik, Iceland
| | - Kari Stefansson
- Department of Population Genetics, deCODE Genetics, Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Reykjavik, Iceland
| | - P Eline Slagboom
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- Netherlands Consortium for Healthy Ageing, The Netherlands
| | - Eleftheria Zeggini
- Department of Human Genetics, Wellcome Trust Sanger Institute, Hinxton, UK
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- Netherlands Consortium for Healthy Ageing, The Netherlands
| | - John PA Ioannidis
- Department of Hygiene and Epidemiology, University of Ioannina Medical School, Ioannina, Greece
- Stanford Prevention Research Center, Stanford University School of Medicine, Stanford, USA
| | - Tim D Spector
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - Joyce B J van Meurs
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Ana M Valdes
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
- Department of Academic Rheumatology, University of Nottingham, Nottingham, UK
| |
Collapse
|
41
|
den Hollander W, Ramos YFM, Bos SD, Bomer N, van der Breggen R, Lakenberg N, de Dijcker WJ, Duijnisveld BJ, Slagboom PE, Nelissen RGHH, Meulenbelt I. Knee and hip articular cartilage have distinct epigenomic landscapes: implications for future cartilage regeneration approaches. Ann Rheum Dis 2014; 73:2208-12. [DOI: 10.1136/annrheumdis-2014-205980] [Citation(s) in RCA: 88] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
|
42
|
Ramos YFM, den Hollander W, Bovée JVMG, Bomer N, van der Breggen R, Lakenberg N, Keurentjes JC, Goeman JJ, Slagboom PE, Nelissen RGHH, Bos SD, Meulenbelt I. Genes involved in the osteoarthritis process identified through genome wide expression analysis in articular cartilage; the RAAK study. PLoS One 2014; 9:e103056. [PMID: 25054223 PMCID: PMC4108379 DOI: 10.1371/journal.pone.0103056] [Citation(s) in RCA: 120] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Accepted: 06/27/2014] [Indexed: 11/19/2022] Open
Abstract
Objective Identify gene expression profiles associated with OA processes in articular cartilage and determine pathways changing during the disease process. Methods Genome wide gene expression was determined in paired samples of OA affected and preserved cartilage of the same joint using microarray analysis for 33 patients of the RAAK study. Results were replicated in independent samples by RT-qPCR and immunohistochemistry. Profiles were analyzed with the online analysis tools DAVID and STRING to identify enrichment for specific pathways and protein-protein interactions. Results Among the 1717 genes that were significantly differently expressed between OA affected and preserved cartilage we found significant enrichment for genes involved in skeletal development (e.g. TNFRSF11B and FRZB). Also several inflammatory genes such as CD55, PTGES and TNFAIP6, previously identified in within-joint analyses as well as in analyses comparing preserved cartilage from OA affected joints versus healthy cartilage were among the top genes. Of note was the high up-regulation of NGF in OA cartilage. RT-qPCR confirmed differential expression for 18 out of 19 genes with expression changes of 2-fold or higher, and immunohistochemistry of selected genes showed a concordant change in protein expression. Most of these changes associated with OA severity (Mankin score) but were independent of joint-site or sex. Conclusion We provide further insights into the ongoing OA pathophysiological processes in cartilage, in particular into differences in macroscopically intact cartilage compared to OA affected cartilage, which seem relatively consistent and independent of sex or joint. We advocate that development of treatment could benefit by focusing on these similarities in gene expression changes and/or pathways.
Collapse
Affiliation(s)
- Yolande F. M. Ramos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
- * E-mail:
| | - Wouter den Hollander
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Nils Bomer
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Ruud van der Breggen
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Nico Lakenberg
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Jelle J. Goeman
- Department of Biostatistics and Bioinformatics, Leiden University Medical Center, Leiden, The Netherlands
| | - P. Eline Slagboom
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Rob G. H. H. Nelissen
- Department of Orthopeadics, Leiden University Medical Center, Leiden, The Netherlands
| | - Steffan D. Bos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
- The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| |
Collapse
|
43
|
Ramos YFM, Metrustry S, Arden N, Bay-Jensen AC, Beekman M, de Craen AJM, Cupples LA, Esko T, Evangelou E, Felson DT, Hart DJ, Ioannidis JPA, Karsdal M, Kloppenburg M, Lafeber F, Metspalu A, Panoutsopoulou K, Slagboom PE, Spector TD, van Spil EWE, Uitterlinden AG, Zhu Y, Valdes AM, van Meurs JBJ, Meulenbelt I. Meta-analysis identifies loci affecting levels of the potential osteoarthritis biomarkers sCOMP and uCTX-II with genome wide significance. J Med Genet 2014; 51:596-604. [PMID: 25057126 DOI: 10.1136/jmedgenet-2014-102478] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
BACKGROUND Research for the use of biomarkers in osteoarthritis (OA) is promising, however, adequate discrimination between patients and controls may be hampered due to innate differences. We set out to identify loci influencing levels of serum cartilage oligomeric protein (sCOMP) and urinary C-telopeptide of type II collagen (uCTX-II). METHODS Meta-analysis of genome-wide association studies was applied to standardised residuals of sCOMP (N=3316) and uCTX-II (N=4654) levels available in 6 and 7 studies, respectively, from TreatOA. Effects were estimated using a fixed-effects model. Six promising signals were followed up by de novo genotyping in the Cohort Hip and Cohort Knee study (N = 964). Subsequently, their role in OA susceptibility was investigated in large-scale genome-wide association studies meta-analyses for OA. Differential expression of annotated genes was assessed in cartilage. RESULTS Genome-wide significant association with sCOMP levels was found for a SNP within MRC1 (rs691461, p = 1.7 × 10(-12)) and a SNP within CSMD1 associated with variation in uCTX-II levels with borderline genome-wide significance (rs1983474, p = 8.5 × 10(-8)). Indication for association with sCOMP levels was also found for a locus close to the COMP gene itself (rs10038, p = 7.1 × 10(-6)). The latter SNP was subsequently found to be associated with hip OA whereas COMP expression appeared responsive to the OA pathophysiology in cartilage. CONCLUSIONS We have identified genetic loci affecting either uCTX-II or sCOMP levels. The genome wide significant association of MRC1 with sCOMP levels was found likely to act independent of OA subtypes. Increased sensitivity of biomarkers with OA may be accomplished by taking genetic variation into account.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands
| | - Sarah Metrustry
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - Nigel Arden
- NIHR Musculoskeletal Biomedical Research Unit, University of Oxford, Oxford, UK Arthritis Research UK, Sport, Exercise and Osteoarthritis Centre of Excellence, London, UK
| | | | - Marian Beekman
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands
| | - Anton J M de Craen
- Department of Gerontology and Geriatrics, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands
| | - L Adrienne Cupples
- Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, USA The Framingham Heart Study, Framingham, Massachusetts, USA
| | - Tõnu Esko
- Institute of Molecular and Cell Biology and Estonian Genome Center, University of Tartu, Tartu, Estonia Department of Endocrinology, Children's Hospital Boston, Boston, Massachusetts, USA Broad Institute, Cambridge, Massachusetts, USA
| | - Evangelos Evangelou
- Department of Hygiene & Epidemiology, University of Ioannina School of Medicine, Ioannina, Greece Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - David T Felson
- Clinical Epidemiology Unit, Boston University School of Medicine, Boston, Massachusetts, USA
| | - Deborah J Hart
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - John P A Ioannidis
- Department of Hygiene & Epidemiology, University of Ioannina School of Medicine, Ioannina, Greece Department of Medicine, Stanford Prevention Research Center, Stanford, USA Department of Health Research and Policy, Stanford University School of Medicine, Stanford, USA Department of Statistics, Stanford University School of Humanities and Sciences, Stanford, USA
| | - Morten Karsdal
- Department of Rheumatology, Nordic Bioscience, Herlev, Denmark
| | - Margreet Kloppenburg
- Department of Rheumatology & Clinical Epidemiology, LUMC, Leiden, The Netherlands
| | - Floris Lafeber
- Department of Rheumatology & Clinical Immunology, UMC Utrecht, Utrecht, The Netherlands
| | - Andres Metspalu
- Institute of Molecular and Cell Biology and Estonian Genome Center, University of Tartu, Tartu, Estonia
| | | | - P Eline Slagboom
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands
| | - Tim D Spector
- Department of Twin Research & Genetic Epidemiology, King's College London, London, UK
| | - Erwin W E van Spil
- Department of Rheumatology & Clinical Immunology, UMC Utrecht, Utrecht, The Netherlands
| | - Andre G Uitterlinden
- The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Yanyan Zhu
- Global Analytical Science, Takeda Pharmaceuticals International Co., Cambridge, Massachusetts, USA
| | | | | | - Ana M Valdes
- Academic Rheumatology, University of Nottingham, Nottingham, UK
| | - Joyce B J van Meurs
- The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, LUMC, Leiden, The Netherlands The Netherlands Genomics Initiative-Sponsored Netherlands Consortium for Healthy Aging, Leiden and Rotterdam, The Netherlands
| |
Collapse
|
44
|
Ramos YFM, Bos SD, van der Breggen R, Kloppenburg M, Ye K, Lameijer EWEMW, Nelissen RGHH, Slagboom PE, Meulenbelt I. A gain of function mutation inTNFRSF11Bencoding osteoprotegerin causes osteoarthritis with chondrocalcinosis. Ann Rheum Dis 2014; 74:1756-62. [DOI: 10.1136/annrheumdis-2013-205149] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2013] [Accepted: 03/23/2014] [Indexed: 01/16/2023]
|
45
|
Ramos YFM, Bos SD, Lakenberg N, Böhringer S, den Hollander WJ, Kloppenburg M, Slagboom PE, Meulenbelt I. Genes expressed in blood link osteoarthritis with apoptotic pathways. Ann Rheum Dis 2013; 73:1844-53. [PMID: 23864235 DOI: 10.1136/annrheumdis-2013-203405] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
OBJECTIVE To identify novel gene expression networks in blood of osteoarthritis patients compared to controls. METHODS A comprehensive exploration of gene expression in peripheral blood was performed by microarray analysis for a subset of osteoarthritis patients from the Genetics osteoARthritis and Progression (GARP) study in comparison with sex and age-matched healthy controls. To identify pathways, we performed gene enrichment analyses (database for annotation, visualisation and integrated discovery and search tool for the retrieval of interacting genes). Quantitative PCR analysis in overlapping and in additional osteoarthritis samples was performed for prioritised genes to validate and replicate findings. Classification of cases and controls was explored by applying statistical models. RESULTS 741 probes representing 694 unique genes were differentially expressed between cases and controls, including 86 genes expressed with at least a 1.5-fold difference. ATF4, GPR18 and H3F3B were among the top genes identified (p<4.5 × 10(-8)). We found that in the blood of osteoarthritis patients the apoptosis pathway, including the well-known gene CASP3, was significantly enriched among the differentially expressed genes. Our findings were validated in independent samples and when using a small subset of the validated genes, we could accurately distinguish patients from controls (area under the curve 98%). CONCLUSIONS In the current study, we have identified specific gene expression networks, in the easily accessible tissue blood, which associated consistently with osteoarthritis among GARP study cases. Our data further hint at the relevance of apoptosis as an aetiological factor in osteoarthritis onset, thereby qualifying expression profiling of blood as a useful tool to understand the underlying molecular mechanisms of osteoarthritis.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Steffan D Bos
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Nico Lakenberg
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Stefan Böhringer
- Department of Biostatistics and Bioinformatics, Leiden University Medical Center, Leiden, The Netherlands
| | - Wouter J den Hollander
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Margreet Kloppenburg
- Departments of Clinical Epidemiology and Rheumatology, Leiden University Medical Center, Leiden, The Netherlands
| | - P Eline Slagboom
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| | - Ingrid Meulenbelt
- Department of Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands The Netherlands Genomics Initiative, sponsored by the NCHA, Leiden-Rotterdam, The Netherlands
| |
Collapse
|
46
|
Bos SD, Bovée JVMG, Duijnisveld BJ, Raine EVA, van Dalen WJ, Ramos YFM, van der Breggen R, Nelissen RGHH, Slagboom PE, Loughlin J, Meulenbelt I. Increased type II deiodinase protein in OA-affected cartilage and allelic imbalance of OA risk polymorphism rs225014 at DIO2 in human OA joint tissues. Ann Rheum Dis 2012; 71:1254-8. [PMID: 22492780 DOI: 10.1136/annrheumdis-2011-200981] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
OBJECTIVE Genetic variation at the type II deiodinase (D2) gene (DIO2) was previously identified as osteoarthritis (OA) risk factor. To investigate mechanisms possibly underlying this association, we assessed D2 protein in healthy and OA-affected cartilage and investigated allelic balance of the OA risk polymorphism rs225014 at DIO2 in human OA joints. METHODS Immunohistochemical staining of healthy and OA-affected cartilage was performed for D2. We then assessed allelic balance of DIO2 mRNA within OA-affected cartilage both at and away from the lesion, ligaments and subchondral bone. Allelic balance was measured by the amount of alleles 'C' and 'T' of the intragenic OA risk polymorphism rs225014 in heterozygous carriers. RESULTS A markedly higher amount of D2 positive cells and staining intensity was observed in OA cartilage. A significant, 1.3-fold higher presence was observed for the OA-associated rs225014 'C' allele relative to the 'T' allele of DIO2, which was significant in 28 of 31 donors. CONCLUSION In OA cartilage, D2 protein presence is increased. The allelic imbalance of the DIO2 mRNA transcript, with the OA risk allele 'C' of rs225014 more abundant than the wild-type 'T' allele in heterozygote carriers provides a possible mechanism by which genetic variation at DIO2 confers OA risk.
Collapse
Affiliation(s)
- Steffan D Bos
- LUMC, Molecular Epidemiology, Einthovenweg 20, Leiden 2333 ZC, The Netherlands.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
47
|
Ramos YFM, Hestand MS, Verlaan M, Krabbendam E, Ariyurek Y, van Galen M, van Dam H, van Ommen GJB, den Dunnen JT, Zantema A, 't Hoen PAC. Genome-wide assessment of differential roles for p300 and CBP in transcription regulation. Nucleic Acids Res 2010; 38:5396-408. [PMID: 20435671 PMCID: PMC2938195 DOI: 10.1093/nar/gkq184] [Citation(s) in RCA: 110] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2009] [Revised: 02/26/2010] [Accepted: 03/04/2010] [Indexed: 12/20/2022] Open
Abstract
Despite high levels of homology, transcription coactivators p300 and CREB binding protein (CBP) are both indispensable during embryogenesis. They are largely known to regulate the same genes. To identify genes preferentially regulated by p300 or CBP, we performed an extensive genome-wide survey using the ChIP-seq on cell-cycle synchronized cells. We found that 57% of the tags were within genes or proximal promoters, with an overall preference for binding to transcription start and end sites. The heterogeneous binding patterns possibly reflect the divergent roles of CBP and p300 in transcriptional regulation. Most of the 16 103 genes were bound by both CBP and p300. However, after stimulation 89 and 1944 genes were preferentially bound by CBP or p300, respectively. Target genes were found to be primarily involved in the regulation of metabolic and developmental processes, and transcription, with CBP showing a stronger preference than p300 for genes active in negative regulation of transcription. Analysis of transcription factor binding sites suggest that CBP and p300 have many partners in common, but AP-1 and Serum Response Factor (SRF) appear to be more prominent in CBP-specific sequences, whereas AP-2 and SP1 are enriched in p300-specific targets. Taken together, our findings further elucidate the distinct roles of coactivators p300 and CBP in transcriptional regulation.
Collapse
Affiliation(s)
- Yolande F M Ramos
- Department of Molecular Cell Biology, Leiden University Medical Center, Postzone S4-0P, PO Box 9600, 2300 RC Leiden, The Netherlands.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
48
|
de Graaf P, Little NA, Ramos YFM, Meulmeester E, Letteboer SJF, Jochemsen AG. Hdmx protein stability is regulated by the ubiquitin ligase activity of Mdm2. J Biol Chem 2003; 278:38315-24. [PMID: 12874296 DOI: 10.1074/jbc.m213034200] [Citation(s) in RCA: 114] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The stability of the p53 tumor suppressor protein is critically regulated by the Hdm2 and Hdmx proteins. Hdm2 protein levels are auto-regulated by the self-ubiquitination activity of Hdm2 and on the transcriptional level by p53-activated transcription of the hdm2 gene. Little is known about the regulation of Hdmx expression levels, apart from the observation that the Mdmx protein can be cleaved by caspase-3 in a p53-inducible manner. In the functional analysis of two mutant Hdmx proteins, products of two alternatively spliced mRNAs, it was found that Hdmx proteins are targets for ubiquitination by Mdm2. The stability of the Hdmx protein is partly dependent on the presence of its internal acidic domain. Mdm2 appears only to require an intact RING domain to be able to ubiquitinate Hdmx and target it for proteasomal degradation. These findings highlight the intricate functional relationships between p53, Mdm2, and Hdmx.
Collapse
Affiliation(s)
- Petra de Graaf
- Leiden University Medical Center, Department of Molecular and Cell Biology and Center for Biomedical Genetics, P. O. Box 9503, 2300 RA Leiden, The Netherlands
| | | | | | | | | | | |
Collapse
|