Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

A functional genomics approach reveals suggestive quantitative trait loci associated with combined TLR4 and BCP crystal-induced inflammation and osteoarthritis.

OBJECTIVE: Basic calcium phosphate (BCP) crystals can activate the NLRP3 inflammasome and are potentially involved in the pathogenesis of osteoarthritis (OA). In order to elucidate relevant inflammatory mechanisms in OA, we used a functional genomics approach to assess genetic variation influencing BCP crystal-induced cytokine production. METHOD: Peripheral blood mononuclear cells (PBMCs) were isolated from healthy volunteers who were previously genotyped and stimulated with BCP crystals and/or lipopolysaccharide (LPS) after which cytokines release was assessed. Cytokine quantitative trait locus (cQTL) mapping was performed. For in vitro validation of the cQTL located in anoctamin 3 (ANO3), PBMCs were incubated with Tamoxifen and Benzbromarone prior to stimulation. Additionally, we performed co-localisation analysis of our top cQTLs with the most recent OA meta-analysis of genome-wide association studies (GWAS). RESULTS: We observed that BCP crystals and LPS synergistically induce IL-1ß in human PBMCs. cQTL analysis revealed several suggestive loci influencing cytokine release upon stimulation, among which are quantitative trait locus annotated to ANO3 and GLIS3. As functional validation, anoctamin inhibitors reduced IL-1ß release in PBMCs after stimulation. Co-localisation analysis showed that the GLIS3 locus was shared between LPS/BCP crystal-induced IL-1ß and genetic association with Knee OA. CONCLUSIONS: We identified and functionally validated a new locus, ANO3, associated with LPS/BCP crystal-induced inflammation in PBMCs. Moreover, the cQTL in the GLIS3 locus co-localises with the previously found locus associated with Knee OA, suggesting that this Knee OA locus might be explained through an inflammatory mechanism. These results form a basis for further exploration of inflammatory mechanisms in OA.

A systematic review of the sex differences in risk factors for knee osteoarthritis.

OBJECTIVES: Previous systematic reviews focused on the evidence of common risk factors for knee OA (KOA); however, the effect and strength of association between risk factors and KOA might be different between the two sexes. The aim of the present systematic review was to determine the current evidence on sex differences in the association between risk factors and KOA and their prevalence. METHODS: We searched the following electronic bibliographic databases: MEDLINE (PubMed), EMBASE and Web of Science. A methodological quality assessment was conducted independently by two researchers according to an adapted version of the standardized set of criteria known as the Newcastle-Ottawa Quality Assessment Scale (NOS). The NOS, a star system, was converted to three categories of quality. RESULTS: In total, 27 studies reported sex-specific risk estimates on several risk factors for KOA. Out of the 22 longitudinal cohort studies (except one nested case-control), 12 were of good quality and 10 were of fair quality. The five cross-sectional studies consisted of one of good, three of fair and one of poor quality. There was an indication of sex differences in risk factors leading to higher risk of KOA: high BMI, alcohol consumption, atherosclerosis and high vitamin E levels in women, and high physical activity, soft drink consumption and abdominal obesity in men. Knee injury, high blood pressure and low step rate seem to affect both women and men. CONCLUSION: More good quality studies are needed to assess sex differences in risk factors for KOA, especially for symptomatic/clinical OA.

Plasma proteomics identifies CRTAC1 as a biomarker for osteoarthritis severity and progression.

OBJECTIVES: The aim of this study was to identify biomarkers for radiographic OA severity and progression acting within the inflammation and metabolic pathways. METHODS: For 3517 Rotterdam Study participants, 184 plasma protein levels were measured using Olink inflammation and cardiometabolic panels. We studied associations with severity and progression of knee, hip and hand OA and a composite overall OA burden score by multivariable regression models, adjusting for age, sex, cell counts and BMI. RESULTS: We found 18 significantly associated proteins for overall OA burden, of which 5 stayed significant after multiple testing correction: circulating cartilage acidic protein 1 (CRTAC1), cartilage oligomeric matrix protein (COMP), thrombospondin 4, IL-18 receptor 1 (IL-18R1) and TNF ligand superfamily member 14. These proteins were also associated with progression of knee OA, with the exception of IL-18R1. The strongest association was found for the level of CRTAC1, with 1 s.d. increase in protein level resulting in an increase of 0.09 (95% CI 0.06, 0.12) in the overall OA Kellgren-Lawrence sum score (P = 2.9 × 10-8) in the model adjusted for age, sex, BMI and cell counts. This association was also present with the severity of OA in all three joints and progression of knee OA and was independent of BMI. We observed a stronger association for CRTAC1 with OA than for the well-known OA biomarker COMP. CONCLUSION: We identified several compelling biomarkers reflecting the overall OA burden and the increased risk for OA progression. CRTAC1 was the most compelling and robust biomarker for OA severity and progression. Such a biomarker may be used for disease monitoring.

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly.

BACKGROUND: Looking older for one's chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. OBJECTIVES: To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. METHODS: We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5-87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). RESULTS: We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one's age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62-0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77-0.95), less age-related hearing loss (model 2; B = -0.76, 95% CI -1.35 to -0.17) and fewer cataracts (OR 0.84, 95% CI 0.73-0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04-0.10). CONCLUSIONS: PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

Towards sex-specific osteoarthritis risk models: evaluation of risk factors for knee osteoarthritis in males and females.

OBJECTIVES: The aim of this study was to identify sex-specific prevalence and strength of risk factors for the incidence of radiographic knee OA (incRKOA). METHODS: Our study population consisted of 10 958 Rotterdam Study participants free of knee OA in one or both knees at baseline. One thousand and sixty-four participants developed RKOA after a median follow-up time of 9.6 years. We estimated the association between each available risk factor and incRKOA using sex stratified multivariate regression models with generalized estimating equations. Subsequently, we statistically tested sex differences between risk estimates and calculated the population attributable fractions (PAFs) for modifiable risk factors. RESULTS: The prevalence of the investigated risk factors was, in general, higher in women compared with men, except that alcohol intake and smoking were higher in men and high BMI showed equal prevalence. We found significantly different risk estimates between men and women: high level of physical activity [relative risk (RR) 1.76 (95% CI: 1.29-2.40)] or a Kellgren and Lawrence score 1 at baseline [RR 5.48 (95% CI: 4.51-6.65)] was higher in men. Among borderline significantly different risk estimates was BMI ≥27, associated with higher risk for incRKOA in women [RR 2.00 (95% CI: 1.74-2.31)]. The PAF for higher BMI was 25.6% in women and 19.3% in men. CONCLUSION: We found sex-specific differences in both presence and relative risk of several risk factors for incRKOA. Especially BMI, a modifiable risk factor, impacts women more strongly than men. These risk factors can be used in the development of personalized prevention strategies and in building sex-specific prediction tools to identify high risk profile patients.

Shorter periconception maternal telomere length and the risk of congenital cardiac outflow defects in the offspring.

BACKGROUND: Congenital cardiac outflow defects (COD) are the largest group of congenital heart defects, with ventricular septal defect (VSD) as the most prevalent phenotype. Increased maternal age, excessive oxidative stress and inflammation are involved in the pathophysiology of COD and enhance telomere length (TL) shortening. We investigated the association between periconception maternal TL and the risk of having a child with COD. METHODS: From a multicentre case-control trial, 306 case mothers of a child with COD and 424 control mothers of a child without a congenital malformation were selected. Relative TL was measured by qPCR. Multivariable logistic regression was used to compute crude and adjusted odds ratios, per standard deviation decrease, between maternal T/S ratio and COD and VSD risk. Adjustments were made for maternal age. Additional adjustments were made in a second model. RESULTS: Shorter maternal relative TL was significantly associated with an OR of 1.29 (95% CI 1.04-1.61), p = .02, for the risk of VSD in offspring, which remained significant after an adjustment for maternal age (adjOR 1.25(95% CI 1.01-1.55), p = .04). No association between maternal TL and the risk of overall COD in offspring was observed. CONCLUSION: Shorter maternal relative TL is associated with an approximately 1.3-OR for the risk, per SD in relative TL shortening, of VSD in the offspring. These findings need further confirmation in other studies on the predictive value of maternal TL.

The interrelatedness of chronic cough and chronic pain.

Since chronic cough has common neurobiological mechanisms and pathophysiology with chronic pain, both clinical disorders might be interrelated. Hence, we examined the association between chronic cough and chronic pain in adult subjects in the Rotterdam Study, a large prospective population-based cohort study.Using a standardised questionnaire, chronic pain was defined as pain lasting up to 6 months and grouped into a frequency of weekly/monthly or daily pain. Chronic cough was described as daily coughing for at least 3 months duration. The longitudinal and cross-sectional associations were investigated bi-directionally.Of 7141 subjects in the study, 54% (n=3888) reported chronic pain at baseline. The co-prevalence of daily chronic pain and chronic cough was 4.4%. Chronic cough was more prevalent in subjects with daily and weekly/monthly chronic pain compared with those without chronic pain (13.8% and 10.3% versus 8.2%; p<0.001). After adjustment for potential confounders, prevalent chronic pain was significantly associated with incident chronic cough (OR 1.47, 95% CI 1.08-1.99). The association remained significant in subjects with daily chronic pain (OR 1.49, 95% CI 1.06-2.11) with a similar effect estimate, albeit non-significant in those with weekly/monthly chronic pain (OR 1.43, 95% CI 0.98-2.10). After adjustment for covariables, subjects with chronic cough had a significant risk of developing chronic pain (OR 1.63, 95% CI 1.02-2.62) compared with those without chronic cough.Chronic cough and chronic pain confer risk on each other among adult subjects, indicating that both conditions might share common risk factors and/or pathophysiologic mechanisms.

Mendelian randomization study on vitamin D levels and osteoarthritis risk: a concise report.

OBJECTIVE: The role of vitamin D in OA is unclear and previous epidemiological studies have provided inconsistent results. We conducted a two-sample Mendelian randomization (MR) study to investigate the causal relationship between genetically determined serum vitamin D levels and hip/knee OA. METHODS: Six single-nucleotide polymorphisms (SNPs) associated with vitamin D levels in the Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits Consortium were selected as instrumental variables. Summary statistics of the SNPs effects on OA were derived from the Iceland and UK Biobank, comprising 23 877 knee OA cases, 17 151 hip OA cases and >562 000 controls. The control samples match the OA cases in age, sex and county of origin. RESULTS: The MR analyses showed no causal association between genetically determined vitamin D levels and knee OA [odds ratio (OR) 1.03 (95% CI 0.84, 1.26)] or hip OA [OR 1.06 (95% CI 0.83, 1.35)]. CONCLUSION: Genetic variations associated with low vitamin D serum levels are not associated with increased risk of hip or knee OA in community-dwelling older adults, suggesting that vitamin D levels are not causally linked to OA. It is therefore unlikely that vitamin D supplementation protects against hip or knee OA.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression.

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

DNA methylation age is associated with an altered hemostatic profile in a multiethnic meta-analysis.

Many hemostatic factors are associated with age and age-related diseases; however, much remains unknown about the biological mechanisms linking aging and hemostatic factors. DNA methylation is a novel means by which to assess epigenetic aging, which is a measure of age and the aging processes as determined by altered epigenetic states. We used a meta-analysis approach to examine the association between measures of epigenetic aging and hemostatic factors, as well as a clotting time measure. For fibrinogen, we performed European and African ancestry-specific meta-analyses which were then combined via a random effects meta-analysis. For all other measures we could not estimate ancestry-specific effects and used a single fixed effects meta-analysis. We found that 1-year higher extrinsic epigenetic age as compared with chronological age was associated with higher fibrinogen (0.004 g/L/y; 95% confidence interval, 0.001-0.007; P = .01) and plasminogen activator inhibitor 1 (PAI-1; 0.13 U/mL/y; 95% confidence interval, 0.07-0.20; P = 6.6 × 10-5) concentrations, as well as lower activated partial thromboplastin time, a measure of clotting time. We replicated PAI-1 associations using an independent cohort. To further elucidate potential functional mechanisms, we associated epigenetic aging with expression levels of the PAI-1 protein encoding gene (SERPINE1) and the 3 fibrinogen subunit-encoding genes (FGA, FGG, and FGB) in both peripheral blood and aorta intima-media samples. We observed associations between accelerated epigenetic aging and transcription of FGG in both tissues. Collectively, our results indicate that accelerated epigenetic aging is associated with a procoagulation hemostatic profile, and that epigenetic aging may regulate hemostasis in part via gene transcription.

Meta-analysis of epigenome-wide association studies of cognitive abilities.

Cognitive functions are important correlates of health outcomes across the life-course. Individual differences in cognitive functions are partly heritable. Epigenetic modifications, such as DNA methylation, are susceptible to both genetic and environmental factors and may provide insights into individual differences in cognitive functions. Epigenome-wide meta-analyses for blood-based DNA methylation levels at ~420,000 CpG sites were performed for seven measures of cognitive functioning using data from 11 cohorts. CpGs that passed a Bonferroni correction, adjusting for the number of CpGs and cognitive tests, were assessed for: longitudinal change; being under genetic control (methylation QTLs); and associations with brain health (structural MRI), brain methylation and Alzheimer's disease pathology. Across the seven measures of cognitive functioning (meta-analysis n range: 2557-6809), there were epigenome-wide significant (P < 1.7 × 10-8) associations for global cognitive function (cg21450381, P = 1.6 × 10-8), and phonemic verbal fluency (cg12507869, P = 2.5 × 10-9). The CpGs are located in an intergenic region on chromosome 12 and the INPP5A gene on chromosome 10, respectively. Both probes have moderate correlations (~0.4) with brain methylation in Brodmann area 20 (ventral temporal cortex). Neither probe showed evidence of longitudinal change in late-life or associations with white matter brain MRI measures in one cohort with these data. A methylation QTL analysis suggested that rs113565688 was a cis methylation QTL for cg12507869 (P = 5 × 10-5 and 4 × 10-13 in two lookup cohorts). We demonstrate a link between blood-based DNA methylation and measures of phonemic verbal fluency and global cognitive ability. Further research is warranted to understand the mechanisms linking genomic regulatory changes with cognitive function to health and disease.

Validated inference of smoking habits from blood with a finite DNA methylation marker set.

Inferring a person's smoking habit and history from blood is relevant for complementing or replacing self-reports in epidemiological and public health research, and for forensic applications. However, a finite DNA methylation marker set and a validated statistical model based on a large dataset are not yet available. Employing 14 epigenome-wide association studies for marker discovery, and using data from six population-based cohorts (N = 3764) for model building, we identified 13 CpGs most suitable for inferring smoking versus non-smoking status from blood with a cumulative Area Under the Curve (AUC) of 0.901. Internal fivefold cross-validation yielded an average AUC of 0.897 ± 0.137, while external model validation in an independent population-based cohort (N = 1608) achieved an AUC of 0.911. These 13 CpGs also provided accurate inference of current (average AUCcrossvalidation 0.925 ± 0.021, AUCexternalvalidation0.914), former (0.766 ± 0.023, 0.699) and never smoking (0.830 ± 0.019, 0.781) status, allowed inferring pack-years in current smokers (10 pack-years 0.800 ± 0.068, 0.796; 15 pack-years 0.767 ± 0.102, 0.752) and inferring smoking cessation time in former smokers (5 years 0.774 ± 0.024, 0.760; 10 years 0.766 ± 0.033, 0.764; 15 years 0.767 ± 0.020, 0.754). Model application to children revealed highly accurate inference of the true non-smoking status (6 years of age: accuracy 0.994, N = 355; 10 years: 0.994, N = 309), suggesting prenatal and passive smoking exposure having no impact on model applications in adults. The finite set of DNA methylation markers allow accurate inference of smoking habit, with comparable accuracy as plasma cotinine use, and smoking history from blood, which we envision becoming useful in epidemiology and public health research, and in medical and forensic applications.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis.

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10-8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

An Epigenome-Wide Association Study of Obesity-Related Traits.

We conducted an epigenome-wide association study on obesity-related traits. We used data from 2 prospective, population-based cohort studies: the Rotterdam Study (RS) (2006-2013) and the Atherosclerosis Risk in Communities (ARIC) Study (1990-1992). We used the RS (n = 1,450) as the discovery panel and the ARIC Study (n = 2,097) as the replication panel. Linear mixed-effect models were used to assess the cross-sectional associations between genome-wide DNA methylation in leukocytes and body mass index (BMI) and waist circumference (WC), adjusting for sex, age, smoking, leukocyte proportions, array number, and position on array. The latter 2 variables were modeled as random effects. Fourteen 5'-C-phosphate-G-3' (CpG) sites were associated with BMI and 26 CpG sites with WC in the RS after Bonferroni correction (P < 1.07 × 10-7), of which 12 and 13 CpGs were replicated in the ARIC Study, respectively. The most significant novel CpGs were located on the Musashi RNA binding protein 2 gene (MSI2; cg21139312) and the leucyl-tRNA synthetase 2, mitochondrial gene (LARS2; cg18030453) and were associated with both BMI and WC. CpGs at BRDT, PSMD1, IFI44L, MAP1A, and MAP3K5 were associated with BMI. CpGs at LGALS3BP, MAP2K3, DHCR24, CPSF4L, and TMEM49 were associated with WC. We report novel associations between methylation at MSI2 and LARS2 and obesity-related traits. These results provide further insight into mechanisms underlying obesity-related traits, which can enable identification of new biomarkers in obesity-related chronic diseases.

A whole-blood transcriptome meta-analysis identifies gene expression signatures of cigarette smoking.

Cigarette smoking is a leading modifiable cause of death worldwide. We hypothesized that cigarette smoking induces extensive transcriptomic changes that lead to target-organ damage and smoking-related diseases. We performed a meta-analysis of transcriptome-wide gene expression using whole blood-derived RNA from 10,233 participants of European ancestry in six cohorts (including 1421 current and 3955 former smokers) to identify associations between smoking and altered gene expression levels. At a false discovery rate (FDR) <0.1, we identified 1270 differentially expressed genes in current vs. never smokers, and 39 genes in former vs. never smokers. Expression levels of 12 genes remained elevated up to 30 years after smoking cessation, suggesting that the molecular consequence of smoking may persist for decades. Gene ontology analysis revealed enrichment of smoking-related genes for activation of platelets and lymphocytes, immune response, and apoptosis. Many of the top smoking-related differentially expressed genes, including LRRN3 and GPR15, have DNA methylation loci in promoter regions that were recently reported to be hypomethylated among smokers. By linking differential gene expression with smoking-related disease phenotypes, we demonstrated that stroke and pulmonary function show enrichment for smoking-related gene expression signatures. Mediation analysis revealed the expression of several genes (e.g. ALAS2) to be putative mediators of the associations between smoking and inflammatory biomarkers (IL6 and C-reactive protein levels). Our transcriptomic study provides potential insights into the effects of cigarette smoking on gene expression in whole blood and their relations to smoking-related diseases. The results of such analyses may highlight attractive targets for treating or preventing smoking-related health effects.

Epigenome-Wide Association Study Identifies Methylation Sites Associated With Liver Enzymes and Hepatic Steatosis.

BACKGROUND & AIMS: Epigenetic mechanisms might be involved in the regulation of liver enzyme level. We aimed to identify CpG sites at which DNA methylation levels are associated with blood levels of liver enzymes and hepatic steatosis. METHODS: We conducted an epigenome-wide association study in whole blood for liver enzyme levels, including gamma-glutamyl transferase (GGT), alanine aminotransferase (ALT), and aspartate aminotransferase (AST), among a discovery set of 731 participants of the Rotterdam Study and sought replication in a non-overlapping sample of 719 individuals. Significant DNA methylation changes were further analyzed to evaluate their relation with hepatic steatosis. Expression levels of the top identified gene were measured in 9 human liver cell lines and compared with expression profiles of its potential targets associated with lipid traits. The candidate gene was subsequently knocked down in human hepatoma cells using lentiviral vectors expressing small hairpin RNAs. RESULTS: Eight probes annotated to SLC7A11, SLC1A5, SLC43A1, PHGDH, PSORS1C1, SREBF1, ANKS3 were associated with GGT and 1 probe annotated to SLC7A11 was associated with ALT after Bonferroni correction (1.0 × 10-7). No probe was identified for AST levels. Four probes for GGT levels including cg06690548 (SLC7A11), cg11376147 (SLC43A1), cg22304262 (SLC1A5), and cg14476101 (PHGDH), and 1 for ALT cg06690548 (SLC7A11) were replicated. DNA methylation at SLC7A11 was associated with reduced risk of hepatic steatosis in participants (odds ratio, 0.69; 95% CI= 0.55-0.93; P value: 2.7 × 10-3). In functional experiments, SLC7A11 was highly expressed in human liver cells; its expression is positively correlated with expression of a panel of lipid-associated genes, indicating a role of SLC7A11 in lipid metabolism. CONCLUSIONS: Our results provide new insights into epigenetic mechanisms associated with markers of liver function and hepatic steatosis, laying the groundwork for future diagnostic and therapeutic applications.

A systematic analysis highlights multiple long non-coding RNAs associated with cardiometabolic disorders.

Genome-wide association studies (GWAS) have identified many susceptibility loci for cardiometabolic disorders. Most of the associated variants reside in non-coding regions of the genome including long non-coding RNAs (lncRNAs), which are thought to play critical roles in diverse biological processes. Here, we leveraged data from the available GWAS meta-analyses on lipid and obesity-related traits, blood pressure, type 2 diabetes, and coronary artery disease and identified 179 associated single-nucleotide polymorphisms (SNPs) in 102 lncRNAs (p-value < 2.3 × 10-7). Of these, 55 SNPs, either the lead SNP or in strong linkage disequilibrium with the lead SNP in the related loci, were selected for further investigations. Our in silico predictions and functional annotations of the SNPs as well as expression and DNA methylation analysis of their lncRNAs demonstrated several lncRNAs that fulfilled predefined criteria for being potential functional targets. In particular, we found evidence suggesting that LOC157273 (at 8p23.1) is involved in regulating serum lipid-cholesterol. Our results showed that rs4841132 in the second exon and cg17371580 in the promoter region of LOC157273 are associated with lipids; the lncRNA is expressed in liver and associates with the expression of its nearby coding gene, PPP1R3B. Collectively, we highlight a number of loci associated with cardiometabolic disorders for which the association may act through lncRNAs.

Epigenetic changes in umbilical cord mesenchymal stromal cells upon stimulation and culture expansion.

BACKGROUND: Mesenchymal stromal cells (MSCs) are studied for their immunotherapeutic potential. Prior to therapeutic use, MSCs are culture expanded to obtain the required cell numbers and, to improve their efficacy, MSCs may be primed in vitro. Culture expansion and priming induce phenotypical and functional changes in MSCs and thus standardisation and quality control measurements come in need. We investigated the impact of priming and culturing on MSC DNA methylation and examined the use of epigenetic profiling as a quality control tool. METHODS: Human umbilical cord-derived MSCs (ucMSCs) were cultured for 3 days with interferon (IFN)γ, transforming growth factor (TGF)ß or a multi-factor combination (MC; IFNγ, TGFß and retinoic acid). In addition, ucMSCs were culture expanded for 14 days. Phenotypical changes and T-cell proliferation inhibition capacity were examined. Genome-wide DNA methylation was measured with Infinium MethylationEPIC Beadchip. RESULTS: Upon priming, ucMSCs exhibited a different immunophenotype and ucMSC(IFNγ) and ucMSC(MC) had an increased capacity to inhibit T-cell proliferation. DNA methylation patterns were minimally affected by priming, with only one significantly differentially methylated site (DMS) in IFNγ- and MC-primed ucMSCs associated with autophagy activity. In contrast, 14 days after culture expansion, ucMSCs displayed minor phenotypical and functional changes but showed >4000 significantly DMSs, mostly concerning genes involved in membrane composition, cell adhesion and transmembrane signalling. DISCUSSION: These data show that DNA methylation of MSCs is only marginally affected by priming, whereas culture expansion and subsequent increased cellular interactions have a large impact on methylation. On account of this study, we suggest that DNA methylation analysis is a useful quality control tool for culture expanded therapeutic MSCs.