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.
Epigenetic remodeling is one of the major features of the aging process. We recently demonstrated that DNA methylation of ELOVL2 and FHL2 CpG islands is highly correlated with age in whole blood. ...Here we investigated several aspects of age-associated hypermethylation of ELOVL2 and FHL2. We showed that ELOVL2 methylation is significantly different in primary dermal fibroblast cultures from donors of different ages. Using epigenomic data from public resources, we demonstrated that most of the tissues show ELOVL2 and FHL2 hypermethylation with age. Interestingly, ELOVL2 hypermethylation was not found in tissues with very low replication rate. We demonstrated that ELOVL2 hypermethylation is associated with in vitro cell replication rather than with senescence. We confirmed intra-individual hypermethylation of ELOVL2 and FHL2 in longitudinally assessed participants from the Doetinchem Cohort Study. Finally we showed that, although the methylation of the two loci is not associated with longevity/mortality in the Leiden Longevity Study, ELOVL2 methylation is associated with cytomegalovirus status in nonagenarians, which could be informative of a higher number of replication events in a fraction of whole-blood cells. Collectively, these results indicate that ELOVL2 methylation is a marker of cell divisions occurring during human aging.
X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained ...by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants.
Human leukocyte telomere length (LTL) decreases with age and shorter LTL has previously been associated with increased prospective mortality. However, it is not clear whether LTL merely marks the ...health status of an individual by its association with parameters of immune function, for example, or whether telomere shortening also contributes causally to lifespan variation in humans.
We measured LTL in 870 nonagenarian siblings (mean age 93 years), 1580 of their offspring and 725 spouses thereof (mean age 59 years) from the Leiden Longevity Study (LLS).
We found that shorter LTL is associated with increased prospective mortality in middle (30-80 years; hazard ratio (HR)=0.75, P=0.001) and highly advanced age (≥90 years; HR=0.92, P=0.028), and show that this association cannot be explained by the association of LTL with the immune-related markers insulin-like growth factor 1 to insulin-like growth factor binding protein 3 molar ratio, C-reactive protein, interleukin 6, cytomegalovirus serostatus or white blood cell counts. We found no difference in LTL between the middle-aged LLS offspring and their spouses (β=0.006, P=0.932). Neither did we observe an association of LTL-associated genetic variants with mortality in a prospective meta-analysis of multiple cohorts (n=8165).
We confirm LTL to be a marker of prospective mortality in middle and highly advanced age and additionally show that this association could not be explained by the association of LTL with various immune-related markers. Furthermore, the approaches performed here do not further support the hypothesis that LTL variation contributes to the genetic propensity for 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.
Objective
To identify osteoarthritis (OA) progression–modulating pathways in articular cartilage and their respective regulatory epigenetic and genetic determinants in end‐stage disease.
Methods
...Transcriptional activity of CpG was assessed using gene expression data and DNA methylation data for preserved and lesional articular cartilage samples. Disease‐responsive transcriptionally active CpG were identified by means of differential methylation between preserved and lesional cartilage. Transcriptionally relevant genetic determinants were addressed by means of single‐nucleotide polymorphisms (SNPs) proximal to the OA‐responsive transcriptionally active CpG. Statistical analyses were corrected for age, sex, joint, and technical covariates. A random effect was included to correct for possible correlations between paired samples.
Results
Of 9,838 transcribed genes in articular cartilage, 2,324 correlated with the methylation status of 3,748 transcriptionally active CpG; both negative (n = 1,741) and positive (n = 2,007) correlations were observed. Hypomethylation and hypermethylation (false discovery rate of <0.05, |Δβ| > 0.05) were observed for 62 and 25 transcriptionally active CpG, respectively, covering 70 unique genes. Enrichment for developmental and extracellular matrix maintenance pathways indicated possible reactivation of endochondral ossification. Finally, we observed 31 and 26 genes for which methylation and expression, respectively, were additionally affected by genetic variation.
Conclusion
We identified tissue‐specific genes involved in OA disease progression, reflected by genetic and pathologic epigenetic regulation of transcription, primarily at genes involved in development. Therefore, transcriptionally active SNPs near these genes may serve as putative susceptibility alleles. Our results constitute an important step toward understanding the reported widespread epigenetic changes occurring in OA articular cartilage and toward subsequent development of treatments targeting disease‐driving pathways.
To identify FN1 transcripts associated with OA pathophysiology and investigate the downstream effects of modulating FN1 expression and relative transcript ratio.
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.
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.
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.
DNA methylation has been recognized as a key mechanism in cell differentiation. Various studies have compared tissues to characterize epigenetically regulated genomic regions, but due to differences ...in study design and focus there still is no consensus as to the annotation of genomic regions predominantly involved in tissue-specific methylation. We used a new algorithm to identify and annotate tissue-specific differentially methylated regions (tDMRs) from Illumina 450k chip data for four peripheral tissues (blood, saliva, buccal swabs and hair follicles) and six internal tissues (liver, muscle, pancreas, subcutaneous fat, omentum and spleen with matched blood samples).
The majority of tDMRs, in both relative and absolute terms, occurred in CpG-poor regions. Further analysis revealed that these regions were associated with alternative transcription events (alternative first exons, mutually exclusive exons and cassette exons). Only a minority of tDMRs mapped to gene-body CpG islands (13%) or CpG islands shores (25%) suggesting a less prominent role for these regions than indicated previously. Implementation of ENCODE annotations showed enrichment of tDMRs in DNase hypersensitive sites and transcription factor binding sites. Despite the predominance of tissue differences, inter-individual differences in DNA methylation in internal tissues were correlated with those for blood for a subset of CpG sites in a locus- and tissue-specific manner.
We conclude that tDMRs preferentially occur in CpG-poor regions and are associated with alternative transcription. Furthermore, our data suggest the utility of creating an atlas cataloguing variably methylated regions in internal tissues that correlate to DNA methylation measured in easy accessible peripheral tissues.
Summary
By studying the loci that contribute to human longevity, we aim to identify mechanisms that contribute to healthy aging. To identify such loci, we performed a genome‐wide association study ...(GWAS) comparing 403 unrelated nonagenarians from long‐living families included in the Leiden Longevity Study (LLS) and 1670 younger population controls. The strongest candidate SNPs from this GWAS have been analyzed in a meta‐analysis of nonagenarian cases from the Rotterdam Study, Leiden 85‐plus study, and Danish 1905 cohort. Only one of the 62 prioritized SNPs from the GWAS analysis (P < 1 × 10−4) showed genome‐wide significance with survival into old age in the meta‐analysis of 4149 nonagenarian cases and 7582 younger controls OR = 0.71 (95% CI 0.65–0.77), P = 3.39 × 10−17. This SNP, rs2075650, is located in TOMM40 at chromosome 19q13.32 close to the apolipoprotein E (APOE) gene. Although there was only moderate linkage disequilibrium between rs2075650 and the ApoE ε4 defining SNP rs429358, we could not find an APOE‐independent effect of rs2075650 on longevity, either in cross‐sectional or in longitudinal analyses. As expected, rs429358 associated with metabolic phenotypes in the offspring of the nonagenarian cases from the LLS and their partners. In addition, we observed a novel association between this locus and serum levels of IGF‐1 in women (P = 0.005). In conclusion, the major locus determining familial longevity up to high age as detected by GWAS was marked by rs2075650, which tags the deleterious effects of the ApoE ε4 allele. No other major longevity locus was found.
EpiTYPER® is a mass spectrometry-based bisulfite sequencing method that enables region-specific DNA methylation analysis in a quantitative and high-throughput fashion. The technology targets genomic ...regions of 100-600 base pairs and results in the quantitative measurement of DNA methylation levels largely at single-nucleotide resolution. It is particularly suitable for larger scale efforts to study candidate regions or to validate regions from genome-wide DNA methylation studies. Here, we describe in detail how to design and perform EpiTYPER measurements and preprocess the data, providing details for high quality measurements not provided in the standard EpiTYPER protocol.