High microsatellite and SNP genotyping success rate established in a large number of genomic DNA samples extracted from mouth swabs and genotypes

In this article, we present the genomic DNA yield and the microsatellite and single nucleotide polymorphism (SNP) genotyping success rates of genomic DNA extracted from a large number of mouth swab samples. In total, the median yield and quality was determined in 714 individuals and the success rates in 378,480 genotypings of 915 individuals. The median yield of genomic DNA per mouth swab was 4.1 μg (range 0.1-42.2 μg) and was not reduced when mouth swabs were stored for at least 21 months prior to extraction. A maximum of 20 mouth swabs is collected per participant. Mouth swab samples showed in, respectively, 89% for 390 microsatellites and 99% for 24 SNPs a genotyping success rate higher than 75%. A very low success rate of genotyping (0%-10%) was obtained for 3.2% of the 915 mouth swab samples using microsatellite markers. Only 0.005% of the mouth swab samples showed a genotyping success rate lower than 75% (range 58%-71 %) using SNPs. Our results show that mouth swabs can be easily collected, stored by our conditions for months prior to DNA extraction and result in high yield and high-quality DNA appropriate for genotyping with high success rate including whole genome searches using microsatellites or SNPs

Circulating microRNAs highly correlate to expression of cartilage genes potentially reflecting OA susceptibility: towards identification of applicable early OA biomarkers

Objective: To identify and validate circulating micro RNAs (miRNAs) that mark gene expression changes in articular cartilage early in osteoarthritis (OA) pathophysiology process. Methods: Within the ongoing RAAK study, human preserved OA cartilage and plasma (N = 22 paired samples) was collected for RNA sequencing (respectively mRNA and miRNA). Spearman correlation was determined for 114 cartilage genes consistently and significantly differentially expressed early in osteoarthritis and 384 plasma miRNAs. Subsequently, the minimal number of circulating miRNAs serving to discriminate between progressors and non-progressors was assessed by regression analysis and area under receiver operating curves (AUC) was calculated with progression data and plasma miRNA sequencing from the GARP study (N = 71). Results: We identified strong correlations (rho >= |0.7|) among expression levels of 34 unique plasma miRNAs and 21 genes, including 4 genes that correlated with multiple miRNAs. The strongest correlation was between let-7d-5p and EGFLAM (rho = -0.75, P = 6.9 x 10(-5)). Regression analysis of the 34 miRNAs resulted in a set of 7 miRNAs that, when applied to the GARP study, demonstrated clinically relevant predictive value with AUC > 0.8 for OA progression over 2 years and near-clinical value for progression over 5 years- (AUC = 0.8). Conclusions: We show that plasma miRNAs levels reflect gene expression levels in cartilage and can be exploited to represent ongoing pathophysiological processes in articular cartilage. We advocate that identified signature of 7 plasma miRNAs can contribute to direct further studies toward early biomarkers predictive for progression of osteoarthritis over 2 and 5 years.Molecular Epidemiolog

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

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-1b; 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 carti-lage. The determined direct negative effect of warfarin on human explant cultures functionally un-derscores the previously found association between vitamin K deficiency and OA. (c) 2021 The Authors. Published by Elsevier Ltd on behalf of Osteoarthritis Research Society International. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).Molecular Epidemiolog

Genetics and Not Shared Environment Explains Familial Resemblance in Adult Metabolomics Data

Metabolites are small molecules involved in cellular metabolism where they act as reaction substrates or products. The term 'metabolomics' refers to the comprehensive study of these molecules. The concentrations of metabolites in biological tissues are under genetic control, but this is limited by environmental factors such as diet. In adult mono- and dizygotic twin pairs, we estimated the contribution of genetic and shared environmental influences on metabolite levels by structural equation modeling and tested whether the familial resemblance for metabolite levels is mainly explained by genetic or by environmental factors that are shared by family members. Metabolites were measured across three platforms: two based on proton nuclear magnetic resonance techniques and one employing mass spectrometry. These three platforms comprised 237 single metabolic traits of several chemical classes. For the three platforms, metabolites were assessed in 1407, 1037 and 1116 twin pairs, respectively. We carried out power calculations to establish what percentage of shared environmental variance could be detected given these sample sizes. Our study did not find evidence for a systematic contribution of shared environment, defined as the influence of growing up together in the same household, on metabolites assessed in adulthood. Significant heritability was observed for nearly all 237 metabolites; significant contribution of the shared environment was limited to 6 metabolites. The top quartile of the heritability distribution was populated by 5 of the 11 investigated chemical classes. In this quartile, metabolites of the class lipoprotein were significantly overrepresented, whereas metabolites of classes glycerophospholipids and glycerolipids were significantly underrepresented

Human osteochondral explants: reliable biomimetic models to investigate disease mechanisms and develop personalized treatments for osteoarthritis

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 beta (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 beta, T3 and 65% MS) resulted in upregulation of the catabolic genes MMP13 and EPAS1. IL-1 beta 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 beta 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.Orthopaedics, Trauma Surgery and Rehabilitatio

Capturing essential physiological aspects of interacting cartilage and bone tissue with osteoarthritis pathophysiology: a human osteochondral unit-on-a-chip model

Given the multi-tissue aspects of osteoarthritis (OA) pathophysiology, translation of OA susceptibility genes towards underlying biological mechanism and eventually drug target discovery requires appropriate human in vitro OA models that incorporate both functional bone and cartilage tissue units. Therefore, a microfluidic chip is developed with an integrated fibrous polycaprolactone matrix in which neo-bone and cartilage are produced, that could serve as a tailored human in vitro disease model of the osteochondral unit of joints. The model enables to evaluate OA-related environmental perturbations to (individual) tissue units by controlling environmental cues, for example by adding biochemical agents. After establishing the co-culture in the system, a layer of cartilaginous matrix is deposited in the chondrogenic compartment, while a bone-like matrix is deposited between the fibers, indicated by both histology and gene expression levels of collagen type 2 and osteopontin, respectively. As proof-of-principle, the bone and cartilaginous tissue are exposed to active thyroid hormone, creating an OA disease model. This results in increased expression levels of hypertrophy markers integrin-binding sialoprotein and alkaline phosphatase in both cartilage and bone, as expected. Altogether, this model could contribute to enhanced translation from OA risk genes towards novel OA therapies.Molecular Epidemiolog

Elucidating mechano-pathology of osteoarthritis: transcriptome-wide differences in mechanically stressed aged human cartilage explants

Background Failing of intrinsic chondrocyte repair after mechanical stress is known as one of the most important initiators of osteoarthritis. Nonetheless, insight into these early mechano-pathophysiological processes in age-related human articular cartilage is still lacking. Such insights are needed to advance clinical development. To highlight important molecular processes of osteoarthritis mechano-pathology, the transcriptome-wide changes following injurious mechanical stress on human aged osteochondral explants were characterized. Methods Following mechanical stress at a strain of 65% (65%MS) on human osteochondral explants (n(65%MS) = 14 versus n(control) = 14), RNA sequencing was performed. Differential expression analysis between control and 65%MS was performed to determine mechanical stress-specific changes. Enrichment for pathways and protein-protein interactions was analyzed with Enrichr and STRING. Results We identified 156 genes significantly differentially expressed between control and 65%MS human osteochondral explants. Of note, IGFBP5 (FC = 6.01; FDR = 7.81 x 10(-3)) and MMP13 (FC = 5.19; FDR = 4.84 x 10(-2)) were the highest upregulated genes, while IGFBP6 (FC = 0.19; FDR = 3.07 x 10(-4)) was the most downregulated gene. Protein-protein interactions were significantly higher than expected by chance (P = 1.44 x 10(-15) with connections between 116 out of 156 genes). Pathway analysis showed, among others, enrichment for cellular senescence, insulin-like growth factor (IGF) I and II binding, and focal adhesion. Conclusions Our results faithfully represent transcriptomic wide consequences of mechanical stress in human aged articular cartilage with MMP13, IGF binding proteins, and cellular senescence as the most notable results. Acquired knowledge on the as such identified initial, osteoarthritis-related, detrimental responses of chondrocytes may eventually contribute to the development of effective disease-modifying osteoarthritis treatments.Molecular Epidemiolog

Identification and characterization of two consistent osteoarthritis subtypes by transcriptome and clinical data integration

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.Molecular Technology and Informatics for Personalised Medicine and Healt

Characterizing the secretome of licensed hiPSC-derived MSCs

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 alpha and IFN gamma 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 NFKB 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.Molecular Epidemiolog

Inhibiting thyroid activation in aged human explants prevents mechanical induced detrimental signalling by mitigating metabolic processes

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.Molecular Epidemiolog