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.

DNA Methylation in Newborns and Maternal Smoking in Pregnancy: Genome-wide Consortium Meta-analysis.

Epigenetic modifications, including DNA methylation, represent a potential mechanism for environmental impacts on human disease. Maternal smoking in pregnancy remains an important public health problem that impacts child health in a myriad of ways and has potential lifelong consequences. The mechanisms are largely unknown, but epigenetics most likely plays a role. We formed the Pregnancy And Childhood Epigenetics (PACE) consortium and meta-analyzed, across 13 cohorts (n = 6,685), the association between maternal smoking in pregnancy and newborn blood DNA methylation at over 450,000 CpG sites (CpGs) by using the Illumina 450K BeadChip. Over 6,000 CpGs were differentially methylated in relation to maternal smoking at genome-wide statistical significance (false discovery rate, 5%), including 2,965 CpGs corresponding to 2,017 genes not previously related to smoking and methylation in either newborns or adults. Several genes are relevant to diseases that can be caused by maternal smoking (e.g., orofacial clefts and asthma) or adult smoking (e.g., certain cancers). A number of differentially methylated CpGs were associated with gene expression. We observed enrichment in pathways and processes critical to development. In older children (5 cohorts, n = 3,187), 100% of CpGs gave at least nominal levels of significance, far more than expected by chance (p value < 2.2 × 10(-16)). Results were robust to different normalization methods used across studies and cell type adjustment. In this large scale meta-analysis of methylation data, we identified numerous loci involved in response to maternal smoking in pregnancy with persistence into later childhood and provide insights into mechanisms underlying effects of this important exposure.

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.

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.

Cell Specific eQTL Analysis without Sorting Cells.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn's disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.

A meta-analysis of gene expression signatures of blood pressure and hypertension.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

Tobacco smoking is associated with DNA methylation of diabetes susceptibility genes.

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.

Determinants for Quantitative Sensory Testing and the Association with Chronic Musculoskeletal Pain in the General Elderly Population.

OBJECTIVE: Chronic musculoskeletal pain is accompanied by central sensitization, which can be determined with quantitative sensory testing (QST). In this study, we aim to investigate whether central sensitization, as measured by thermal QST, is detectable in community-dwelling elderly individuals suffering from self-reported chronic pain and identify determinants influencing thermal QST measurement analyses and interpretation. METHODS: In 3,936 participants of the Rotterdam Study, cold and warmth sensitivity and heat pain thresholds were determined using the thermo-sensory analyzer TSA II (Medoc Advanced Medical Systems, Durham, NC, U.S.A.). Using Cox regression, associations were studied with chronic pain and potential determinants (body mass index [BMI], reaction speed, systolic and diastolic blood pressure, skin color, skin temperature, seasonal influence, depression, anxiety, atopic eczema, age at menarche, years since menopause, hormone replacement therapy (HRT) use during menopause, and reproductive lifespan). RESULTS: In addition to the effect of age and gender on thermal sensitivity, darker skin color and the presence of atopic eczema were associated with higher sensitivity for heat pain. Cold sensitivity and warmth sensitivity thresholds were both influenced by BMI, reaction speed, skin temperature, season, depression, dark skin color, years since menopause, and reproductive lifespan. The presence of chronic pain was associated with 0.2 degrees lower heat pain threshold in all participants, and 0.3 degrees lower in individuals with chronic pain in more than 2 sites. CONCLUSION: Higher sensitivity for heat pain, one feature of central sensitization, is present in community-dwelling elderly with chronic pain. Additional determinants should be considered when analyzing and interpreting QST measurements.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Genome-wide association and functional studies identify the DOT1L gene to be involved in cartilage thickness and hip osteoarthritis.

Hip osteoarthritis (HOA) is one of the most disabling and common joint disorders with a large genetic component that is, however, still ill-defined. To date, genome-wide association studies (GWAS) in osteoarthritis (OA) and specifically in HOA have yielded only few loci, which is partly explained by heterogeneity in the OA definition. Therefore, we here focused on radiographically measured joint-space width (JSW), a proxy for cartilage thickness and an important underlying intermediate trait for HOA. In a GWAS of 6,523 individuals on hip-JSW, we identified the G allele of rs12982744 on chromosome 19p13.3 to be associated with a 5% larger JSW (P = 4.8 × 10(-10)). The association was replicated in 4,442 individuals from three United Kingdom cohorts with an overall meta-analysis P value of 1.1 × 10(-11). The SNP was also strongly associated with a 12% reduced risk for HOA (P = 1 × 10(-4)). The SNP is located in the DOT1L gene, which is an evolutionarily conserved histone methyltransferase, recently identified as a potentially dedicated enzyme for Wnt target-gene activation in leukemia. Immunohistochemical staining of the DOT1L protein in mouse limbs supports a role for DOT1L in chondrogenic differentiation and adult articular cartilage. DOT1L is also expressed in OA articular chondrocytes. Silencing of Dot1l inhibited chondrogenesis in vitro. Dot1l knockdown reduces proteoglycan and collagen content, and mineralization during chondrogenesis. In the ATDC5 chondrogenesis model system, DOT1L interacts with TCF and Wnt signaling. These data are a further step to better understand the role of Wnt-signaling during chondrogenesis and cartilage homeostasis. DOT1L may represent a therapeutic target for OA.

A genetic variant in the seed region of miR-4513 shows pleiotropic effects on lipid and glucose homeostasis, blood pressure, and coronary artery disease.

MicroRNAs (miRNA) play a crucial role in the regulation of diverse biological processes by post-transcriptional modulation of gene expression. Genetic polymorphisms in miRNA-related genes can potentially contribute to a wide range of phenotypes. The effect of such variants on cardiometabolic diseases has not yet been defined. We systematically investigated the association of genetic variants in the seed region of miRNAs with cardiometabolic phenotypes, using the thus far largest genome-wide association studies on 17 cardiometabolic traits/diseases. We found that rs2168518:G>A, a seed region variant of miR-4513, associates with fasting glucose, low-density lipoprotein-cholesterol, total cholesterol, systolic and diastolic blood pressure, and risk of coronary artery disease. We experimentally showed that miR-4513 expression is significantly reduced in the presence of the rs2168518 mutant allele. We sought to identify miR-4513 target genes that may mediate these associations and revealed five genes (PCSK1, BNC2, MTMR3, ANK3, and GOSR2) through which these effects might be taking place. Using luciferase reporter assays, we validated GOSR2 as a target of miR-4513 and further demonstrated that the miRNA-mediated regulation of this gene is changed by rs2168518. Our findings indicate a pleiotropic effect of miR-4513 on cardiometabolic phenotypes and may improve our understanding of the pathophysiology of cardiometabolic diseases.

Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region.

BACKGROUND AND OBJECTIVES: Chronic widespread pain (CWP) is a common disorder affecting ∼10% of the general population and has an estimated heritability of 48-52%. In the first large-scale genome-wide association study (GWAS) meta-analysis, we aimed to identify common genetic variants associated with CWP. METHODS: We conducted a GWAS meta-analysis in 1308 female CWP cases and 5791 controls of European descent, and replicated the effects of the genetic variants with suggestive evidence for association in 1480 CWP cases and 7989 controls. Subsequently, we studied gene expression levels of the nearest genes in two chronic inflammatory pain mouse models, and examined 92 genetic variants previously described associated with pain. RESULTS: The minor C-allele of rs13361160 on chromosome 5p15.2, located upstream of chaperonin-containing-TCP1-complex-5 gene (CCT5) and downstream of FAM173B, was found to be associated with a 30% higher risk of CWP (minor allele frequency=43%; OR=1.30, 95% CI 1.19 to 1.42, p=1.2×10(-8)). Combined with the replication, we observed a slightly attenuated OR of 1.17 (95% CI 1.10 to 1.24, p=4.7×10(-7)) with moderate heterogeneity (I2=28.4%). However, in a sensitivity analysis that only allowed studies with joint-specific pain, the combined association was genome-wide significant (OR=1.23, 95% CI 1.14 to 1.32, p=3.4×10(-8), I2=0%). Expression levels of Cct5 and Fam173b in mice with inflammatory pain were higher in the lumbar spinal cord, not in the lumbar dorsal root ganglions, compared to mice without pain. None of the 92 genetic variants previously described were significantly associated with pain (p>7.7×10(-4)). CONCLUSIONS: We identified a common genetic variant on chromosome 5p15.2 associated with joint-specific CWP in humans. This work suggests that CCT5 and FAM173B are promising targets in the regulation of pain.

Identification of genetic loci associated with Helicobacter pylori serologic status.

IMPORTANCE: Helicobacter pylori is a major cause of gastritis and gastroduodenal ulcer disease and can cause cancer. H. pylori prevalence is as high as 90% in some developing countries but 10% of a given population is never colonized, regardless of exposure. Genetic factors are hypothesized to confer H. pylori susceptibility. OBJECTIVE: To identify genetic loci associated with H. pylori seroprevalence in 2 independent population-based cohorts and to determine their putative pathophysiological role by whole-blood RNA gene expression profiling. DESIGN, SETTING, AND PARTICIPANTS: Two independent genome-wide association studies (GWASs) and a subsequent meta-analysis were conducted for anti-H. pylori IgG serology in the Study of Health in Pomerania (SHIP) (recruitment, 1997-2001 [n = 3830]) as well as the Rotterdam Study (RS-I) (recruitment, 1990-1993) and RS-II (recruitment, 2000-2001 [n = 7108]) populations. Whole-blood RNA gene expression profiles were analyzed in RS-III (recruitment, 2006-2008 [n = 762]) and SHIP-TREND (recruitment, 2008-2012 [n = 991]), and fecal H. pylori antigen in SHIP-TREND (n = 961). MAIN OUTCOMES AND MEASURES: H. pylori seroprevalence. RESULTS: Of 10,938 participants, 6160 (56.3%) were seropositive for H. pylori. GWASs identified the toll-like receptor (TLR) locus (4p14; top-ranked single-nucleotide polymorphism (SNP), rs10004195; P = 1.4 × 10(-18); odds ratio, 0.70 [95% CI, 0.65 to 0.76]) and the FCGR2A locus (1q23.3; top-ranked SNP, rs368433; P = 2.1 × 10(-8); odds ratio, 0.73 [95% CI, 0.65 to 0.81]) as associated with H. pylori seroprevalence. Among the 3 TLR genes at 4p14, only TLR1 was differentially expressed per copy number of the minor rs10004195-A allele (ß = -0.23 [95% CI, -0.34 to -0.11]; P = 2.1 × 10(-4)). Individuals with high fecal H. pylori antigen titers (optical density >1) also exhibited the highest 25% of TLR1 expression levels (P = .01 by χ2 test). Furthermore, TLR1 exhibited an Asn248Ser substitution in the extracellular domain strongly linked to the rs10004195 SNP. CONCLUSIONS AND RELEVANCE: GWAS meta-analysis identified an association between TLR1 and H. pylori seroprevalence, a finding that requires replication in nonwhite populations. If confirmed, genetic variations in TLR1 may help explain some of the observed variation in individual risk for H. pylori infection.

An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis.

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

Caenorhabditis elegans orthologs of human genes differentially expressed with age are enriched for determinants of longevity.

We report a systematic RNAi longevity screen of 82 Caenorhabditis elegans genes selected based on orthology to human genes differentially expressed with age. We find substantial enrichment in genes for which knockdown increased lifespan. This enrichment is markedly higher than published genomewide longevity screens in C. elegans and similar to screens that preselected candidates based on longevity-correlated metrics (e.g., stress resistance). Of the 50 genes that affected lifespan, 46 were previously unreported. The five genes with the greatest impact on lifespan (>20% extension) encode the enzyme kynureninase (kynu-1), a neuronal leucine-rich repeat protein (iglr-1), a tetraspanin (tsp-3), a regulator of calcineurin (rcan-1), and a voltage-gated calcium channel subunit (unc-36). Knockdown of each gene extended healthspan without impairing reproduction. kynu-1(RNAi) alone delayed pathology in C. elegans models of Alzheimer's disease and Huntington's disease. Each gene displayed a distinct pattern of interaction with known aging pathways. In the context of published work, kynu-1, tsp-3, and rcan-1 are of particular interest for immediate follow-up. kynu-1 is an understudied member of the kynurenine metabolic pathway with a mechanistically distinct impact on lifespan. Our data suggest that tsp-3 is a novel modulator of hypoxic signaling and rcan-1 is a context-specific calcineurin regulator. Our results validate C. elegans as a comparative tool for prioritizing human candidate aging genes, confirm age-associated gene expression data as valuable source of novel longevity determinants, and prioritize select genes for mechanistic follow-up.

Blood RNA expression profiles undergo major changes during the seventh decade.

Genome-wide alterations in RNA expression profiles are age-associated. Yet the rate and temporal pattern of those alterations are poorly understood. We investigated temporal changes in RNA expression profiles in blood from population-based studies using a quadratic regression model. Comparative analysis between two independent studies was carried out after sample-weighting that downsized differences in sample distribution over age between the datasets. We show that age-associated expression profiles are clustered into two major inclinations and transcriptional alternations occur predominantly from the seventh decade onwards. The age-associated genes in blood are enriched in functional groups of the translational machinery and the immune system. The results are highly consistent between the two population-based studies indicating that our analysis overcomes potential confounders in population-based studies. We suggest that the critical age when major transcriptional alterations occur could help understanding aging and disease risk during adulthood.

Associations between joint effusion in the knee and gene expression levels in the circulation: a meta-analysis.

OBJECTIVE: To identify molecular biomarkers for early knee osteoarthritis (OA), we examined whether joint effusion in the knee associated with different gene expression levels in the circulation. MATERIALS AND METHODS: Joint effusion grades measured with magnetic resonance (MR) imaging and gene expression levels in blood were determined in women of the Rotterdam Study (N=135) and GARP (N=98). Associations were examined using linear regression analyses, adjusted for age, fasting status, RNA quality, technical batch effects, blood cell counts, and BMI. To investigate enriched pathways and protein-protein interactions, we used the DAVID and STRING webtools. RESULTS: In a meta-analysis, we identified 257 probes mapping to 189 unique genes in blood that were nominally significantly associated with joint effusion grades in the knee. Several compelling genes were identified such as C1orf38 and NFATC1. Significantly enriched biological pathways were: response to stress, gene expression, negative regulation of intracellular signal transduction, and antigen processing and presentation of exogenous pathways. CONCLUSION: Meta-analyses and subsequent enriched biological pathways resulted in interesting candidate genes associated with joint effusion that require further characterization. Associations were not transcriptome-wide significant most likely due to limited power. Additional studies are required to replicate our findings in more samples, which will greatly help in understanding the pathophysiology of OA and its relation to inflammation, and may result in biomarkers urgently needed to diagnose OA at an early stage.

Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns.

Folate is vital for fetal development. Periconceptional folic acid supplementation and food fortification are recommended to prevent neural tube defects. Mechanisms whereby periconceptional folate influences normal development and disease are poorly understood: epigenetics may be involved. We examine the association between maternal plasma folate during pregnancy and epigenome-wide DNA methylation using Illumina's HumanMethyl450 Beadchip in 1,988 newborns from two European cohorts. Here we report the combined covariate-adjusted results using meta-analysis and employ pathway and gene expression analyses. Four-hundred forty-three CpGs (320 genes) are significantly associated with maternal plasma folate levels during pregnancy (false discovery rate 5%); 48 are significant after Bonferroni correction. Most genes are not known for folate biology, including APC2, GRM8, SLC16A12, OPCML, PRPH, LHX1, KLK4 and PRSS21. Some relate to birth defects other than neural tube defects, neurological functions or varied aspects of embryonic development. These findings may inform how maternal folate impacts the developing epigenome and health outcomes in offspring.

Peripheral Blood Transcriptomic Signatures of Fasting Glucose and Insulin Concentrations.

Genome-wide association studies (GWAS) have successfully identified genetic loci associated with glycemic traits. However, characterizing the functional significance of these loci has proven challenging. We sought to gain insights into the regulation of fasting insulin and fasting glucose through the use of gene expression microarray data from peripheral blood samples of participants without diabetes in the Framingham Heart Study (FHS) (n = 5,056), the Rotterdam Study (RS) (n = 723), and the InCHIANTI Study (Invecchiare in Chianti) (n = 595). Using a false discovery rate q <0.05, we identified three transcripts associated with fasting glucose and 433 transcripts associated with fasting insulin levels after adjusting for age, sex, technical covariates, and complete blood cell counts. Among the findings, circulating IGF2BP2 transcript levels were positively associated with fasting insulin in both the FHS and RS. Using 1000 Genomes-imputed genotype data, we identified 47,587 cis-expression quantitative trait loci (eQTL) and 6,695 trans-eQTL associated with the 433 significant insulin-associated transcripts. Of note, we identified a trans-eQTL (rs592423), where the A allele was associated with higher IGF2BP2 levels and with fasting insulin in an independent genetic meta-analysis comprised of 50,823 individuals. We conclude that integration of genomic and transcriptomic data implicate circulating IGF2BP2 mRNA levels associated with glucose and insulin homeostasis.

DNA methylation-based measures of biological age: meta-analysis predicting time to death.

Estimates of biological age based on DNA methylation patterns, often referred to as "epigenetic age", "DNAm age", have been shown to be robust biomarkers of age in humans. We previously demonstrated that independent of chronological age, epigenetic age assessed in blood predicted all-cause mortality in four human cohorts. Here, we expanded our original observation to 13 different cohorts for a total sample size of 13,089 individuals, including three racial/ethnic groups. In addition, we examined whether incorporating information on blood cell composition into the epigenetic age metrics improves their predictive power for mortality. All considered measures of epigenetic age acceleration were predictive of mortality (p≤8.2x10-9), independent of chronological age, even after adjusting for additional risk factors (p<5.4x10-4), and within the racial/ethnic groups that we examined (non-Hispanic whites, Hispanics, African Americans). Epigenetic age estimates that incorporated information on blood cell composition led to the smallest p-values for time to death (p=7.5x10-43). Overall, this study a) strengthens the evidence that epigenetic age predicts all-cause mortality above and beyond chronological age and traditional risk factors, and b) demonstrates that epigenetic age estimates that incorporate information on blood cell counts lead to highly significant associations with all-cause mortality.