Search

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.

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.

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.

Objectives 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 3 D 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 eight changed significantly, with BMP6 showing the 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.

Description, Fibronectin binding to collagen type II is crucial for articular chondrocytes to stay in a healthy steady state., 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.

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.

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 &gt, 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.

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.

ObjectiveTo uncover the microRNA (miRNA) interactome of the osteoarthritis (OA) pathophysiological process in the cartilage.MethodsWe 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.ResultsWe 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.ConclusionsBy 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.

Background Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. Methods Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. Results Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99–1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96–1.04; p value 0.59). Conclusions Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

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.

The blood metabolome incorporates cues from the environment and the host's genetic background, potentially offering a holistic view of an individual's health status.We have compiled a vast resource of proton nuclear magnetic resonance metabolomics and phenotypic data encompassing over 25 000 samples derived from 26 community and hospital-based cohorts.Using this resource, we constructed a metabolomics-based age predictor (metaboAge) to calculate an individual's biological age. Exploration in independent cohorts demonstrates that being judged older by one's metabolome, as compared with one's chronological age, confers an increased risk on future cardiovascular disease, mortality, and functionality in older individuals. A web-based tool for calculating metaboAge (metaboage.researchlumc.nl) allows easy incorporation in other epidemiological studies. Access to data can be requested at bbmri.nl/samples-images-data.In summary, we present a vast resource of metabolomics data and illustrate its merit by constructing a metabolomics-based score for biological age that captures aspects of current and future cardiometabolic health.

Background: The blood metabolome incorporates cues from the environment and the host’s genetic background, potentially offering a holistic view of an individual’s health status. Methods: We have compiled a vast resource of proton nuclear magnetic resonance metabolomics and phenotypic data encompassing over 25 000 samples derived from 26 community and hospital-based cohorts. Results: Using this resource, we constructed a metabolomics-based age predictor (metaboAge) to calculate an individual’s biological age. Exploration in independent cohorts demonstrates that being judged older by one’s metabolome, as compared with one’s chronological age, confers an increased risk on future cardiovascular disease, mortality, and functionality in older individuals. A web-based tool for calculating metaboAge ( metaboage.researchlumc.nl ) allows easy incorporation in other epidemiological studies. Access to data can be requested at bbmri.nl/samples-images-data . Conclusions: In summary, we present a vast resource of metabolomics data and illustrate its merit by constructing a metabolomics-based score for biological age that captures aspects of current and future cardiometabolic health.

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.

Introduction: High mortality rates of approximately 20% within 1 year after treatment are observed for patients with proximal femoral fractures. This preliminary study explores the prognostic value of a previously constructed mortality risk score based on a set of 14 metabolites for the survival and functional recovery in patients with proximal femoral fractures. Materials and Methods: A prospective observational cohort study was conducted including patients admitted with a proximal femoral fracture. The primary outcome was patient survival, and the recovery of independence in activities of daily living was included as a secondary outcome. The mortality risk score was constructed for each patient and its prognostic value was tested for the whole population. Results: Data was available form 136 patients. The mean age of all patients was 82.1 years, with a median follow-up of 6 months. Within this period, 19.0% of all patients died and 51.1% recovered to their prefracture level of independence. The mortality score was significantly associated with mortality (HR, 2.74; 95% CI, 1.61-4.66; P < 0.001), but showed only a fair prediction accuracy (AUC = 0.68) and a borderline significant comparison of the mortality score tertile groups in survival analyses (P = 0.049). No decisive associations were found in any of the analyses for the functional recovery of patients. Discussion: These findings support the previously determined prognostic value of the mortality risk score. However, the independent prognostic value when adjusted for potential confounding factors is yet to be assessed. Also, a risk score constructed for this specific patient population might achieve higher accuracies for the prediction of survival and functional recovery. Conclusions: A modest prediction accuracy was observed for the mortality risk score in this population. More elaborate studies are needed to validate these findings and develop a tailored model for clinical purposes in this patient population.

Aim: UVA radiation drives skin photoaging in the dermis, plausibly via persistent changes to DNA methylation in dermal fibroblasts. Methods: Genome-wide DNA methylation changes after five repeated daily UVA doses were determined at 48 h (transitionary) and 1 week (recovery) post final irradiation. Results: Differential methylation was found at the transitionary time point in active chromatin states near genes that are highly expressed in fibroblasts and are involved in cellular defensive mechanisms; the majority of these methylation differences were restored to control levels after 7 day recovery. At the recovery time point, new differential methylation occurred at repressed regions near developmental genes, normally weakly expressed in fibroblasts. Conclusion: UVA irradiation induces transitionary and recovery-associated DNA methylation responses in fibroblasts with contrasting functional characteristics.

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.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Progress in high-throughput metabolic profiling provides unprecedented opportunities to obtain insights into the effects of drugs on human metabolism. The Biobanking BioMolecular Research Infrastructure of the Netherlands has constructed an atlas of drug-metabolite associations for 87 commonly prescribed drugs and 150 clinically relevant plasma-based metabolites assessed by proton nuclear magnetic resonance. The atlas includes a meta-analysis of ten cohorts (18,873 persons) and uncovers 1,071 drug-metabolite associations after evaluation of confounders including co-treatment. We show that the effect estimates of statins on metabolites from the cross-sectional study are comparable to those from intervention and genetic observational studies. Further data integration links proton pump inhibitors to circulating metabolites, liver function, hepatic steatosis and the gut microbiome. Our atlas provides a tool for targeted experimental pharmaceutical research and clinical trials to improve drug efficacy, safety and repurposing. We provide a web-based resource for visualization of the atlas (http://bbmri.researchlumc.nl/atlas/).

Telomere shortening has been associated with multiple age-related diseases such as cardiovascular disease, diabetes, and dementia. However, the biological mechanisms responsible for these associations remain largely unknown. In order to gain insight into the metabolic processes driving the association of leukocyte telomere length (LTL) with age-related diseases, we investigated the association between LTL and serum metabolite levels in 7,853 individuals from seven independent cohorts. LTL was determined by quantitative polymerase chain reaction and the levels of 131 serum metabolites were measured with mass spectrometry in biological samples from the same blood draw. With partial correlation analysis, we identified six metabolites that were significantly associated with LTL after adjustment for multiple testing: lysophosphatidylcholine acyl C17:0 (lysoPC a C17:0, p-value = 7.1 × 10−6), methionine (p-value = 9.2 × 10−5), tyrosine (p-value = 2.1 × 10−4), phosphatidylcholine diacyl C32:1 (PC aa C32:1, p-value = 2.4 × 10−4), hydroxypropionylcarnitine (C3-OH, p-value = 2.6 × 10−4), and phosphatidylcholine acyl-alkyl C38:4 (PC ae C38:4, p-value = 9.0 × 10−4). Pathway analysis showed that the three phosphatidylcholines and methionine are involved in homocysteine metabolism and we found supporting evidence for an association of lipid metabolism with LTL. In conclusion, we found longer LTL associated with higher levels of lysoPC a C17:0 and PC ae C38:4, and with lower levels of methionine, tyrosine, PC aa C32:1, and C3-OH. These metabolites have been implicated in inflammation, oxidative stress, homocysteine metabolism, and in cardiovascular disease and diabetes, two major drivers of morbidity and mortality.

Metabolomics examines the small molecules involved in cellular metabolism. Approximately 50% of total phenotypic differences in metabolite levels is due to genetic variance, but heritability estimates differ across metabolite classes and lipid species. We performed a review of all genetic association studies, and identified > 800 class-specific metabolite loci that influence metabolite levels. In a twin-family cohort (N= 5,117), these metabolite loci were leveraged to simultaneously estimate total heritability (h2total), and the proportion of heritability captured by known metabolite loci (h2Metabolite-hits) for 309 lipids and 52 organic acids. Our study revealed significant differences inh2Metabolite-hitsamong different classes of lipids and organic acids. Furthermore, phosphatidylcholines with a high degree of unsaturation had higherh2Metabolite-hitsestimates than phosphatidylcholines with a low degree of unsaturation. This study highlights the importance of common genetic variants for metabolite levels, and elucidates the genetic architecture of metabolite classes and lipid species.

Skeletal muscles control posture, mobility and strength, affecting whole-body metabolism. The muscle is composed of different types of myofibers, which are defined by their contractile properties, and are marked by the expression of sarcomeric proteins. Three myosin heavy-chain (MyHC) isoforms, MyHC-1, MyHC-2A and MyHC-2X characterize myofiber properties in humans, yet characterization of myofibers is often performed considering only a single MyHC isoform expression. We aimed to study heterogeneity in muscle tissue composition, by detailed analysis of the expression of the three MyHC isoforms in consortium. We created a dataset with high-throughput MyHC immunofluorescence imaging in muscle tissues of 56 older adults. Our analyses revealed six distinct myofiber clusters, some of which could not be detected or are misclassified by the traditional visual myofiber assessment. Transcriptome analyses revealed that the abundance of these myofiber clusters was correlated with the expression of distinct modules of sarcomeric genes. Interestingly, one of the novel clusters, expressing all three MyHC isoforms, correlated to multiple histological measures of muscle health. Collectively, we present a data-driven procedure for a deeper characterization of muscle tissue composition, that opens new options to unravel the relations between alterations in muscle composition and changes in a variety of health conditions and aging.

Contains fulltext : 209008.pdf (Publisher’s version ) (Open Access) BACKGROUND: Previous studies have reported associations between attention-deficit/hyperactivity disorder symptoms and DNA methylation in children. We report the first epigenome-wide association study meta-analysis of adult attention-deficit/hyperactivity disorder symptoms, based on peripheral blood DNA methylation (Infinium HumanMethylation450K array) in three population-based adult cohorts. METHODS: An epigenome-wide association study was performed in the Netherlands Twin Register (N = 2258, mean age 37 years), Dunedin Multidisciplinary Health and Development Study (N = 800, age 38 years), and Environmental Risk Longitudinal Twin Study (N = 1631, age 18 years), and results were combined through meta-analysis (total sample size N = 4689). Region-based analyses accounting for the correlation between nearby methylation sites were also performed. RESULTS: One epigenome-wide significant differentially methylated position was detected in the Dunedin study, but meta-analysis did not detect differentially methylated positions that were robustly associated across cohorts. In region-based analyses, six significant differentially methylation regions (DMRs) were identified in the Netherlands Twin Register, 19 in the Dunedin study, and none in the Environmental Risk Longitudinal Twin Study. Of these DMRs, 92% were associated with methylation quantitative trait loci, and 68% showed moderate to large blood-brain correlations for DNA methylation levels. DMRs included six nonoverlapping DMRs (three in the Netherlands Twin Register, three in the Dunedin study) in the major histocompatibility complex, which were associated with expression of genes in the major histocompatibility complex, including C4A and C4B, previously implicated in schizophrenia. CONCLUSIONS: Our findings point at new candidate loci involved in immune and neuronal functions that await further replication. Our work also illustrates the need for further research to examine to what extent epigenetic associations with psychiatric traits depend on characteristics such as age, comorbidities, exposures, and genetic background.

Contains fulltext : 171367.pdf (Publisher’s version ) (Open Access) Genome-wide association and fine-mapping studies in 14 autoimmune diseases (AID) have implicated more than 250 loci in one or more of these diseases. As more than 90% of AID-associated SNPs are intergenic or intronic, pinpointing the causal genes is challenging. We performed a systematic analysis to link 460 SNPs that are associated with 14 AID to causal genes using transcriptomic data from 629 blood samples. We were able to link 71 (39%) of the AID-SNPs to two or more nearby genes, providing evidence that for part of the AID loci multiple causal genes exist. While 54 of the AID loci are shared by one or more AID, 17% of them do not share candidate causal genes. In addition to finding novel genes such as ULK3, we also implicate novel disease mechanisms and pathways like autophagy in celiac disease pathogenesis. Furthermore, 42 of the AID SNPs specifically affected the expression of 53 non-coding RNA genes. To further understand how the non-coding genome contributes to AID, the SNPs were linked to functional regulatory elements, which suggest a model where AID genes are regulated by network of chromatin looping/non-coding RNAs interactions. The looping model also explains how a causal candidate gene is not necessarily the gene closest to the AID SNP, which was the case in nearly 50% of cases.

Genome-wide association studies have identified hundreds of genetic variants associated with blood pressure (BP), but sequence variation accounts for a small fraction of the phenotypic variance. Epigenetic changes may alter the expression of genes involved in BP regulation and explain part of the missing heritability. We therefore conducted a two-stage meta-analysis of the cross-sectional associations of systolic and diastolic BP with blood-derived genome-wide DNA methylation measured on the Infinium HumanMethylation450 BeadChip in 17,010 individuals of European, African American, and Hispanic ancestry. Of 31 discovery-stage cytosine-phosphate-guanine (CpG) dinucleotides, 13 replicated after Bonferroni correction (discovery: N = 9,828, p < 1.0× 10-7; replication: N = 7,182, p 30%) and independent of known BP genetic variants, explaining an additional 1.4% and 2.0% of the interindividual variation in systolic and diastolic BP, respectively. Bidirectional Mendelian randomization among up to 4,513 individuals of European ancestry from 4 cohorts suggested that methylation at cg08035323 (TAF1B-YWHAQ) influences BP, while BP influences methylation at cg00533891 (ZMIZ1), cg00574958 (CPT1A), and cg02711608 (SLC1A5). Gene expression analyses further identified six genes (TSPAN2, SLC7A11, UNC93B1, CPT1A, PTMS, and LPCAT3) with evidence of triangular associations between methylation, gene expression, and BP. Additional integrative Mendelian randomization analyses of gene expression and DNA methylation suggested that the expression of TSPAN2 is a putative mediator of association between DNA methylation at cg23999170 and BP. These findings suggest that heritable DNA methylation plays a role in regulating BP independently of previously known genetic variants.

Determining cell identity and maturation status of differentiated pluripotent stem cells (PSCs) requires knowledge of the transcriptional and epigenetic trajectory of organs during development. Here, we generate a transcriptional and DNA methylation atlas covering 21 organs during human fetal development. Analysis of multiple isogenic organ sets shows that organ-specific DNA methylation patterns are highly dynamic between week 9 (W9) and W22 of gestation. We investigate the impact of reprogramming on organ-specific DNA methylation by generating human induced pluripotent stem cell (hiPSC) lines from six isogenic organs. All isogenic hiPSCs acquire DNA methylation patterns comparable to existing hPSCs. However, hiPSCs derived from fetal brain retain brain-specific DNA methylation marks that seem sufficient to confer higher propensity to differentiate to neural derivatives. This systematic analysis of human fetal organs during development and associated isogenic hiPSC lines provides insights in the role of DNA methylation in lineage commitment and epigenetic reprogramming in humans., While DNA methylation and gene expression data are widely available for animal models, comprehensive data from human development is rarer. Here, the authors generated transcriptional and DNA methylation data from 21 organs during human development and 6 isogenic induced pluripotent stem cell lines.

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) 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.

Background: The manipulation of pregnancy diets in animals can lead to changes in DNA methylation with phenotypic consequences in the offspring. Human studies have concentrated on the effects of nutrition during early gestation. Lacking in humans is an epigenome-wide association study of DNA methylation in relation to perturbations in nutrition across all gestation periods. Methods: We used the quasi-experimental setting of the Dutch famine of 1944-45 to evaluate the impact of famine exposure during specific 10-week gestation periods, or during any time in gestation, on genome-wide DNA methylation levels at age ~59 years. In addition, we evaluated the impact of exposure during a shorter pre-and post-conception period. DNA methylation was assessed using the Illumina 450k array in whole blood among 422 individuals with prenatal famine exposure and 463 time-or sibling-controls without prenatal famine exposure. Results: Famine exposure during gestation weeks 1-10, but not weeks 11-20, 21-30 or 31-delivery, was associated with an increase in DNA methylation of CpG dinucleotides cg20823026 (FAM150B), cg10354880 (SLC38A2) and cg27370573 (PPAP2C) and a decrease of cg11496778 (OSBPL5/MRGPRG) (P-7, PFDR-7, PFDR=0.034) and a decrease of cg23989336 (TMEM105) after exposure around conception. These changes represent a shift of 0.3-0.6 standard deviations and are linked to genes involved in growth, development and metabolism. Conclusion: Early gestation, and not mid or late gestation, is identified as a critical timeperiod for adult DNA methylation changes in whole blood after prenatal exposure to famine.

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.

Animal lifespan is regulated by conserved metabolic signalling pathways and specific transcription factors, but whether these pathways affect common downstream mechanisms remains largely elusive. Here we show that NCL-1/TRIM2/Brat tumour suppressor extends lifespan and limits nucleolar size in the major C. elegans longevity pathways, as part of a convergent mechanism focused on the nucleolus. Long-lived animals representing distinct longevity pathways exhibit small nucleoli, and decreased expression of rRNA, ribosomal proteins, and the nucleolar protein fibrillarin, dependent on NCL-1. Knockdown of fibrillarin also reduces nucleolar size and extends lifespan. Among wildtype C. elegans, individual nucleolar size varies, but is highly predictive for longevity. Long-lived dietary restricted fruit flies and insulin-like-peptide mutants exhibit small nucleoli and fibrillarin expression, as do long-lived dietary restricted and IRS1 knockout mice. Furthermore, human muscle biopsies from individuals who underwent modest dietary restriction coupled with exercise also display small nucleoli. We suggest that small nucleoli are a cellular hallmark of longevity and metabolic health conserved across taxa., Animal lifespan is plastic and is regulated by conserved signalling pathways. Here, Tiku et al. show that longevity-enhancing mutations or interventions are associated with reduced nucleolar size in worms, flies, mice and humans, and that nucleolar size can predict life-expectancy in individual worms.

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-1131-9) contains supplementary material, which is available to authorized users.

Most disease associated genetic risk factors are non-coding, making it challenging to design experiments to understand their functional consequences. Identification of expression quantitative trait loci (eQTLs) has been a powerful approach to infer downstream effects of disease variants but the large majority remains unexplained.. The analysis of DNA methylation, a key component of the epigenome, offers highly complementary data on the regulatory potential of genomic regions. However, a large-scale, combined analysis of methylome and transcriptome data to infer downstream effects of disease variants is lacking. Here, we show that disease variants have wide-spread effects on DNA methylation in trans that likely reflect the downstream effects on binding sites of cis-regulated transcription factors. Using data on 3,841 Dutch samples, we detected 272,037 independent cis-meQTLs (FDR < 0.05) and identified 1,907 trait-associated SNPs that affect methylation levels of 10,141 different CpG sites in trans (FDR < 0.05), an eight-fold increase in the number of downstream effects that was known from trans-eQTL studies. Trans-meQTL CpG sites are enriched for active regulatory regions, being correlated with gene expression and overlap with Hi-C determined interchromosomal contacts. We detected many trans-meQTL SNPs that affect expression levels of nearby transcription factors (including NFKB1, CTCF and NKX2-3), while the corresponding trans-meQTL CpG sites frequently coincide with its respective binding site. Trans-meQTL mapping therefore provides a strategy for identifying and better understanding downstream functional effects of many disease-associated variants.

The stress response is regulated by the mineralocorticoid receptor (MR) and the glucocorticoid receptor (GR). When the balance between GR and MR signalling is disturbed, one's capacity to cope with a stressful event is diminished. In this pilot study, we tested the hypothesis that an interaction between common variants in the MR (rs5522) or GR gene (rs41423247) and stressful life events influences perfectionism levels in a group of patients with an eating disorder (ED; n = 113). Patients carrying the minor G allele of rs5522 had a higher perfectionism score if more stressful life events were experienced [β = 0.95, t(109) = 3.75, p

Within the Netherlands a national network of biobanks has been established (Biobanking and Biomolecular Research Infrastructure-Netherlands (BBMRI-NL)) as a national node of the European BBMRI. One of the aims of BBMRI-NL is to enrich biobanks with different types of molecular and phenotype data. Here, we describe the Genome of the Netherlands (GoNL), one of the projects within BBMRI-NL. GoNL is a whole-genome-sequencing project in a representative sample consisting of 250 trio-families from all provinces in the Netherlands, which aims to characterize DNA sequence variation in the Dutch population. The parent-offspring trios include adult individuals ranging in age from 19 to 87 years (mean=53 years; SD=16 years) from birth cohorts 1910-1994. Sequencing was done on blood-derived DNA from uncultured cells and accomplished coverage was 14-15x. The family-based design represents a unique resource to assess the frequency of regional variants, accurately reconstruct haplotypes by family-based phasing, characterize short indels and complex structural variants, and establish the rate of de novo mutational events. GoNL will also serve as a reference panel for imputation in the available genome-wide association studies in Dutch and other cohorts to refine association signals and uncover population-specific variants. GoNL will create a catalog of human genetic variation in this sample that is uniquely characterized with respect to micro-geographic location and a wide range of phenotypes. The resource will be made available to the research and medical community to guide the interpretation of sequencing projects. The present paper summarizes the global characteristics of the project. © 2014 Macmillan Publishers Limited All rights reserved.

Over the past 25 years, the Adult Netherlands Twin Register (ANTR) has collected a wealth of information on physical and mental health, lifestyle, and personality in adolescents and adults. This article provides an overview of the sources of information available, the main research findings, and an outlook for the future. Between 1991 and 2012, longitudinal surveys were completed by twins, their parents, siblings, spouses, and offspring. Data are available for 33,957 participants, with most individuals having completed two or more surveys. Smaller projects provided in-depth phenotyping, including measurements of the autonomic nervous system, neurocognitive function, and brain imaging. For 46% of the ANTR participants, DNA samples are available and whole genome scans have been obtained in more than 11,000 individuals. These data have resulted in numerous studies on heritability, gene x environment interactions, and causality, as well as gene finding studies. In the future, these studies will continue with collection of additional phenotypes, such as metabolomic and telomere length data, and detailed genetic information provided by DNA and RNA sequencing. Record linkage to national registers will allow the study of morbidity and mortality, thus providing insight into the development of health, lifestyle, and behavior across the lifespan. Copyright © The Authors 2013.

Aims/hypothesis Tobacco smoking, a risk factor for diabetes, is an established modifier of DNA methylation. We hypothesised that tobacco smoking modifies DNA methylation of genes previously identified for diabetes. Methods We annotated CpG sites available on the Illumina Human Methylation 450K array to diabetes genes previously identified by genome-wide association studies (GWAS), and investigated them for an association with smoking by comparing current to never smokers. The discovery study consisted of 630 individuals (Bonferroni-corrected p = 1.4 × 10−5), and we sought replication in an independent sample of 674 individuals. The replicated sites were tested for association with nearby genetic variants and gene expression and fasting glucose and insulin levels. Results We annotated 3,620 CpG sites to the genes identified in the GWAS on type 2 diabetes. Comparing current smokers to never smokers, we found 12 differentially methylated CpG sites, of which five replicated: cg23161492 within ANPEP (p = 1.3 × 10−12); cg26963277 (p = 1.2 × 10−9), cg01744331 (p = 8.0 × 10−6) and cg16556677 (p = 1.2 × 10−5) within KCNQ1 and cg03450842 (p = 3.1 × 10−8) within ZMIZ1. The effect of smoking on DNA methylation at the replicated CpG sites attenuated after smoking cessation. Increased DNA methylation at cg23161492 was associated with decreased gene expression levels of ANPEP (p = 8.9 × 10−5). rs231356-T, which was associated with hypomethylation of cg26963277 (KCNQ1), was associated with a higher odds of diabetes (OR 1.06, p = 1.3 × 10−5). Additionally, hypomethylation of cg26963277 was associated with lower fasting insulin levels (p = 0.04). Conclusions/interpretation Tobacco smoking is associated with differential DNA methylation of the diabetes risk genes ANPEP, KCNQ1 and ZMIZ1. Our study highlights potential biological mechanisms connecting tobacco smoking to excess risk of type 2 diabetes. Electronic supplementary material The online version of this article (doi:10.1007/s00125-016-3872-0) contains peer-reviewed but unedited supplementary material, which is available to authorised users.

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation., Differential impact of genetic and environmental influences on DNA methylation may result in sex- and age-related physiological variation and disease susceptibility. By analysing DNA methylome of 2,603 individuals from twin families, here, the authors establish a catalogue of between-individual variation in DNA methylation.

To the editor: Recent large-scale sequencing studies report recurrent somatic mutations in the blood of elderly individuals in genes previously linked to clonal expansion of hematopoietic stem cells.[1][1][⇓][2][⇓][3]-[4][4] Particularly for DNMT3A and TET2 , a steep age-associated increase in

In this study, we investigated 102 single-nucleotide polymorphisms (SNPs) covering the common genetic variation in 16 genes recurrently regarded as candidates for human longevity: APOE; ACE; CETP; HFE; IL6; IL6R; MTHFR; TGFB1; APOA4; APOC3; SIRTs 1, 3, 6; and HSPAs 1A, 1L, 14. In a case-control study of 1,089 oldest-old (ages 92-93) and 736 middle-aged Danes, the minor allele frequency (MAF) of rs769449 (APOE) was significantly decreased in the oldest-old, while the MAF of rs9923854 (CETP) was significantly enriched. These effects were supported when investigating 1,613 oldest-old (ages 95-110) and 1,104 middle-aged Germans. rs769449 was in modest linkage equilibrium (R 2 = 0.55) with rs429358 of the APOE-ε4 haplotype and adjusting for rs429358 eliminated the association of rs769449, indicating that the association likely reflects the well-known effect of rs429358. Gene-based analysis confirmed the effects of variation in APOE and CETP and furthermore pointed to HSPA14 as a longevity gene. In a longitudinal study with 11 years of follow-up on survival in the oldest-old Danes, only one SNP, rs2069827 (IL6), was borderline significantly associated with survival from age 92 (P-corrected = 0.064). This advantageous effect of the minor allele was supported when investigating a Dutch longitudinal cohort (N = 563) of oldest-old (age 85+). Since rs2069827 was located in a putative transcription factor binding site, quantitative RNA expression studies were conducted. However, no difference in IL6 expression was observed between rs2069827 genotype groups. In conclusion, we here support and expand the evidence suggesting that genetic variation in APOE, CETP, and IL6, and possible HSPA14, is associated with human longevity.

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

Human height is a highly heritable trait, with genetic factors explaining up to 90% of phenotypic variation. Vitamin D levels are known to influence several physiological processes, including skeletal growth. The vitamin D receptor (VDR) gene has been reported as contributing to variation in height. A meta-analysis of 13607 adult individuals found a small but significant association with the rs1544410 (BsmI) polymorphism. In contrast, the meta-analysis found no effect in a sample of 550 children. Two recent studies reported variants with large effect on height elsewhere in VDR (rs10735810 [FokI] and rs7139166 [-1521] polymorphisms). We genotyped large Caucasian samples from Australia (N= 3906) and the Netherlands (N= 1689) for polymorphisms in VDR. The Australian samples were twin families with height measures from 3 time points throughout adolescence. The Dutch samples were adult twins. We use the available family data to perform both within and between family tests of association. We found no significant associations for any of the genotyped variants after multiple testing correction. The (non-significant) effect of rs1544410 in the Australian adolescent cohort was in the same direction and of similar magnitude (additive effect 0.3cm) to the effect observed in the published adult meta-analysis. An effect of this size explains ~0.1% of the phenotypic variance in height — this implies that many, probably hundreds, of such variants are responsible for the observed genetic variation. Our results did not support any role for two other regions (rs10735810, rs7139166) of VDR in explaining variation in height.

Remodelling the methylome is a hallmark of mammalian development and cell differentiation. However, current knowledge of DNA methylation dynamics in human tissue specification and organ development largely stems from the extrapolation of studies in vitro and animal models. Here, we report on the DNA methylation landscape using the 450k array of four human tissues (amnion, muscle, adrenal and pancreas) during the first and second trimester of gestation (9,18 and 22 weeks). We show that a tissue-specific signature, constituted by tissue-specific hypomethylated CpG sites, was already present at 9 weeks of gestation (W9). Furthermore, we report large-scale remodelling of DNA methylation from W9 to W22. Gain of DNA methylation preferentially occurred near genes involved in general developmental processes, whereas loss of DNA methylation mapped to genes with tissue-specific functions. Dynamic DNA methylation was associated with enhancers, but not promoters. Comparison of our data with external fetal adrenal, brain and liver revealed striking similarities in the trajectory of DNA methylation during fetal development. The analysis of gene expression data indicated that dynamic DNA methylation was associated with the progressive repression of developmental programs and the activation of genes involved in tissue-specific processes. The DNA methylation landscape of human fetal development provides insight into regulatory elements that guide tissue specification and lead to organ functionality., Author Summary Methylation of DNA is a key epigenetic mark. Adult tissues have highly distinct genome-wide DNA methylation signatures. How these signatures arise during human fetal development is largely unknown. Here, we studied DNA methylation profiles of four tissues (amnion, muscle, adrenal, pancreas) during first and second trimester of human fetal development. Already in the first trimester, a tissue-specific signature was found in each of the tissues. However, during the first and second trimester, a substantial number of genomic regions were found to gain and lose DNA methylation. Genomic regions that gained methylation were associated with the shut-down of developmental processes, while genomic regions that lose methylation were associated with the activation of tissue-specific functions. These findings on the DNA methylation landscape of human fetal development are important as they provide insight into regulatory elements that guide tissue specification and lead to organ functionality.

In order to meaningfully analyze common and rare genetic variants, results from genome-wide association studies (GWASASs) of multiple cohorts need to be combined in a meta-analysis in order to obtain enough power. This requires all cohorts to have the same single-nucleotide polymorphisms (SNPSNPSNPs) in their GWASASs. To this end, genotypes that have not been measured in a given cohort can be imputed on the basis of a set of reference haplotypes. This protocol provides guidelines for performing imputations with two widely used tools: minimac and IMPUTEPUTEPUTEPUTE2. These guidelines were developed and used by the Genome of the Netherlands (GoNL) consortium, which has created a population-specific reference panel for genetic imputations and used this reference to impute various Dutch biobanks. We also describe several factors that might influence the final imputation quality. This protocol, which has been used by the largest Dutch biobanks, should take approximately several days, depending on the sample size of the biobank and the computer resources available.

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.

To identify osteoarthritis (OA) progression-modulating pathways in articular cartilage and their respective regulatory epigenetic and genetic determinants in end-stage disease.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.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 of0.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.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.

The stress response is regulated by the mineralocorticoid receptor (MR) and the glucocorticoid receptor (GR). When the balance between GR and MR signalling is disturbed, one's capacity to cope with a stressful event is diminished. In this pilot study, we tested the hypothesis that an interaction between common variants in the MR (rs5522) or GR gene (rs41423247) and stressful life events influences perfectionism levels in a group of patients with an eating disorder (ED; n = 113). Patients carrying the minor G allele of rs5522 had a higher perfectionism score if more stressful life events were experienced [β = 0.95, t(109) = 3.75, p 0.01]. This effect was not found for patients carrying the AA genotype. These results suggest that rs5522 G allele carriers might be vulnerable to stressful life events. When patients with an ED are carriers and experience multiple life events, this might fuel their insecurity, which in turn may engender higher levels of perfectionism. Further studies are necessary to replicate and expand our findings.

Background: 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. Methods: 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). Results: 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 (a parts per thousand yen90 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 (beta = 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). Conclusions: 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.

Contains fulltext : 137213.pdf (Publisher’s version ) (Closed access) Whole-genome sequencing enables complete characterization of genetic variation, but geographic clustering of rare alleles demands many diverse populations be studied. Here we describe the Genome of the Netherlands (GoNL) Project, in which we sequenced the whole genomes of 250 Dutch parent-offspring families and constructed a haplotype map of 20.4 million single-nucleotide variants and 1.2 million insertions and deletions. The intermediate coverage ( approximately 13x) and trio design enabled extensive characterization of structural variation, including midsize events (30-500 bp) previously poorly catalogued and de novo mutations. We demonstrate that the quality of the haplotypes boosts imputation accuracy in independent samples, especially for lower frequency alleles. Population genetic analyses demonstrate fine-scale structure across the country and support multiple ancient migrations, consistent with historical changes in sea level and flooding. The GoNL Project illustrates how single-population whole-genome sequencing can provide detailed characterization of genetic variation and may guide the design of future population studies.