Depressive Symptoms and Plasma Markers of Alzheimer's Disease and Neurodegeneration: A Coordinated Meta-Analysis of 8 Cohort Studies.

BACKGROUND: Depressive symptoms are associated with an increased risk of Alzheimer's disease (AD). There has been a recent emergence in plasma biomarkers for AD pathophysiology, such as amyloid-beta (Aß) and phosphorylated tau (p-tau), as well as for axonal damage (neurofilament light, NfL) and astrocytic activation (glial fibrillary acidic protein, GFAP). Hypothesizing that depressive symptoms may occur along the AD process, we investigated associations between plasma biomarkers of AD with depressive symptoms in individuals without dementia. METHODS: A two-stage meta-analysis was performed on 2 clinic-based and 6 population-based cohorts (N = 7210) as part of the Netherlands Consortium of Dementia Cohorts. Plasma markers (Aß42/40, p-tau181, NfL, and GFAP) were measured using Single Molecular Array (Simoa; Quanterix) assays. Depressive symptoms were measured with validated questionnaires. We estimated the cross-sectional association of each standardized plasma marker (determinants) with standardized depressive symptoms (outcome) using linear regressions, correcting for age, sex, education, and APOE ε4 allele presence, as well as subgrouping by sex and APOE ε4 allele. Effect estimates were entered into a random-effects meta-analysis. RESULTS: Mean age of participants was 71 years. The prevalence of clinically relevant depressive symptoms ranged from 1% to 22%. None of the plasma markers were associated with depressive symptoms in the meta-analyses. However, NfL was associated with depressive symptoms only in APOE ε4 carriers (ß 0.11; 95% CI: 0.05-0.17). CONCLUSIONS: Late-life depressive symptoms did not show an association to plasma biomarkers of AD pathology. However, in APOE ε4 allele carriers, a more profound role of neurodegeneration was suggested with depressive symptoms.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis.

OBJECTIVES: To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients. METHODS: RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the genome-wide association studies (GWAS) catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP. RESULTS: We identified 30 transcript SNPs (28 in cartilage and two in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7×10-4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75×10-5). Moreover, 12 drugs were interacting with seven genes displaying AEI, of which seven drugs have been already approved. CONCLUSIONS: By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harbouring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further translation towards underlying biological mechanisms.

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.

Long non-coding RNA expression profiling of subchondral bone reveals AC005165.1 modifying FRZB expression during osteoarthritis.

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.

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.

A data-driven methodology reveals novel myofiber clusters in older human muscles.

Skeletal muscles control posture, mobility and strength, and influence whole-body metabolism. Muscles are built of different types of myofibers, each having specific metabolic, molecular, and contractile properties. Fiber classification is, therefore, regarded the key for understanding muscle biology, (patho-) physiology. The expression of three myosin heavy chain (MyHC) isoforms, MyHC-1, MyHC-2A, and MyHC-2X, marks myofibers in humans. Typically, myofiber classification is performed by an eye-based histological analysis. This classical approach is insufficient to capture complex fiber classes, expressing more than one MyHC-isoform. We, therefore, developed a methodological procedure for high-throughput characterization of myofibers on the basis of multiple isoforms. The mean fluorescence intensity of the three most abundant MyHC isoforms was measured per myofiber in muscle biopsies of 56 healthy elderly adults, and myofiber classes were identified using computational biology tools. Unsupervised clustering revealed the existence of six distinct myofiber clusters. A comparison with the visual assessment of myofibers using the same images showed that some of these myofiber clusters could not be detected or were frequently misclassified. The presence of these six clusters was reinforced by RNA expressions levels of sarcomeric genes. In addition, one of the clusters, expressing all three MyHC isoforms, correlated with histological measures of muscle health. To conclude, this methodological procedure enables deep characterization of the complex muscle heterogeneity. This study opens opportunities to further investigate myofiber composition in comparative studies.

RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage.

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.

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.

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.

Underlying molecular mechanisms of DIO2 susceptibility in symptomatic osteoarthritis.

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.

Genes expressed in blood link osteoarthritis with apoptotic pathways.

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.

Genome-wide association study (GWAS)-identified disease risk alleles do not compromise human longevity.

A set of currently known alleles increasing the risk for coronary artery disease, cancer, and type 2 diabetes as identified by genome-wide association studies was tested for compatibility with human longevity. Here, we show that nonagenarian siblings from long-lived families and singletons older than 85 y of age from the general population carry the same number of disease risk alleles as young controls. Longevity in this study population is not compromised by the cumulative effect of this set of risk alleles for common disease.

Large replication study and meta-analyses of DVWA as an osteoarthritis susceptibility locus in European and Asian populations.

Recently, through a genome wide association study in Japanese knee osteoarthritis (OA) cases, a previously unknown gene, DVWA, was identified. The non-synonymous single nucleotide polymorphism (SNP) rs7639618 was subsequently found to be consistent and most significantly associated in Japanese and Han Chinese knee OA studies and functional relevant. Here, the association of the DVWA polymorphisms (rs7639618, rs11718863 and rs9864422) was genotyped in 1120 knee OA cases, 1482 hip OA cases and 2147 controls, all of white European descent from the Netherlands, the UK, Spain and Greece. Random effect DerSimonian and Laird meta-analyses were performed to assess the association in the different strata. To assess a more global effect, the original Japanese and Chinese data were included with the European. The meta-analyses provided evidence for global association of rs7639618 with knee OA with an odds ratio (OR) of 1.29, 95% confidence interval (CI) of 1.15-1.45 and a P-value of 2.70 x 10(-5). This effect, however, showed moderate heterogeneity, and rs7639618 was not independently associated with knee OA in Europeans, with an OR of 1.16, 95% CI of 0.99-1.35 and a P-value of 0.063. Furthermore, no association was observed with hip OA in Europeans, with a P-value of 0.851. Our results suggest that there may be global relevance for the DVWA SNP rs7639618 among knee OA cases, however, the apparent lower effect size in combination with the higher risk allele frequency in the European samples highlights again the ethnic differences in effects of discovered OA susceptibility genes.

Interleukin-1 gene cluster variants with innate cytokine production profiles and osteoarthritis in subjects from the Genetics, Osteoarthritis and Progression Study.

OBJECTIVE: To assess whether genetic variation in the interleukin-1 (IL-1) gene cluster contributes to familial osteoarthritis (OA) by influencing innate ex vivo production of IL-1beta or IL-1 receptor antagonist (IL-1Ra). METHODS: Innate ex vivo IL-1beta and IL-1Ra production upon lipopolysaccharide (LPS) stimulation of whole blood cells was measured in subjects from the Genetics, Osteoarthritis and Progression (GARP) Study, which includes sibling pairs in which at least one sibling has symptomatic OA at multiple sites. Radiographic OA (ROA) was assessed by Kellgren/Lawrence score. Subjects from the GARP Study and controls from the Rotterdam Study were genotyped for 7 single-nucleotide polymorphisms (SNPs) encompassing the IL-1 gene cluster on chromosome 2q13. Linkage disequilibrium analysis and genotype and haplotype association analysis were performed to assess the relationship between the IL-1 gene cluster SNPs, innate ex vivo cytokine production, and OA. RESULTS: Among subjects in the GARP Study, the haplotype variable-number tandem repeat in intron 2/T+8006C/T+11100C 2/2/1 of the IL1RN gene was significantly associated with reduced innate ex vivo bioavailability of IL-1beta upon LPS stimulation (P = 0.026) and with ROA at the highest number of joint locations. CONCLUSION: These results show that genetic variation at the IL-1 gene cluster is associated with lower IL-1beta bioavailability and with OA at a large number of joint locations. The data further indicate that, among subjects with OA affecting the highest number of joints, the innate immune system may be activated, thereby obscuring possible underlying mechanisms.

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 (ρ ≥ |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 (ρ = -0.75, P = 6.9 × 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.

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ß (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.

RNA Sequencing Reveals Interacting Key Determinants of Osteoarthritis Acting in Subchondral Bone and Articular Cartilage: Identification of IL11 and CHADL as Attractive Treatment Targets.

OBJECTIVE: To identify key determinants of the interactive pathophysiologic processes in subchondral bone and cartilage in osteoarthritis (OA). METHODS: We performed RNA sequencing on macroscopically preserved and lesional OA subchondral bone from patients in the Research Arthritis and Articular Cartilage study who underwent joint replacement surgery due to OA (n = 24 sample pairs: 6 hips and 18 knees). Unsupervised hierarchical clustering and differential expression analyses were conducted. Results were combined with data on previously identified differentially expressed genes in cartilage (partly overlapping samples) as well as data on recently identified OA risk genes. RESULTS: We identified 1,569 genes that were significantly differentially expressed between lesional and preserved subchondral bone, including CNTNAP2 (fold change [FC] 2.4, false discovery rate [FDR] 3.36 × 10-5 ) and STMN2 (FC 9.6, FDR 2.36 × 10-3 ). Among these 1,569 genes, 305 were also differentially expressed, and with the same direction of effect, in cartilage, including the recently recognized OA susceptibility genes IL11 and CHADL. Upon differential expression analysis with stratification for joint site, we identified 509 genes that were exclusively differentially expressed in subchondral bone of the knee, including KLF11 and WNT4. These genes that were differentially expressed exclusively in the knee were enriched for involvement in epigenetic processes, characterized by, e.g., HIST1H3J and HIST1H3H. CONCLUSION: IL11 and CHADL were among the most consistently differentially expressed genes OA pathophysiology-related genes in both bone and cartilage. As these genes were recently also identified as robust OA risk genes, they classify as attractive therapeutic targets acting on 2 OA-relevant tissues.

Identification of DIO2 as a new susceptibility locus for symptomatic osteoarthritis.

Osteoarthritis [MIM 165720] is a common late-onset articular joint disease for which no pharmaceutical intervention is available to attenuate the cartilage degeneration. To identify a new osteoarthritis susceptibility locus, a genome-wide linkage scan and combined linkage association analysis were applied to 179 affected siblings and four trios with generalized osteoarthritis (The GARP study). We tested, for confirmation by association, 1478 subjects who required joint replacement and 734 controls in a UK population. Additional replication was tested in 1582 population-based females from the Rotterdam study that contained 94 cases with defined hip osteoarthritis and in 267 Japanese females with symptomatic hip osteoarthritis and 465 controls. Suggested evidence for linkage in the GARP study was observed on chromosome 14q32.11 (log of odds = 3.03, P = 1.9 x 10(-4)). Genotyping tagging single-nucleotide polymorphisms covering three important candidate genes revealed a common coding variant (rs225014; Thr92Ala) in the iodothyronine-deiodinase enzyme type 2 (D2) gene (DIO2 [MIM 601413]) which significantly explained the linkage signal (P = 0.006). Confirmation and replication by association in the additional osteoarthritis studies indicated a common DIO2 haplotype, exclusively containing the minor allele of rs225014 and common allele of rs12885300, with a combined recessive odds ratio of 1.79, 95% confidence interval (CI) 1.37-2.34 with P = 2.02 x 10(-5) in female cases with advanced/symptomatic hip osteoarthritis. The gene product of this DIO2 converts intracellular pro-hormone-3,3',5,5'-tetraiodothyronine (T4) into the active thyroid hormone 3,3',5-triiodothyronine (T3) thereby regulating intracellular levels of active T3 in target tissues such as the growth plate. Our results indicate a new susceptibility gene (DIO2) conferring risk to osteoarthritis.

Elucidating Epigenetic Regulation by Identifying Functional cis-Acting Long Noncoding RNAs and Their Targets in Osteoarthritic Articular Cartilage.

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 indicates that intergenic and antisense lncRNAs play an important role in regulating the pathophysiology of OA.

Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage.

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.