Background: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential ...methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS). Methods: We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C). Findings: Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4–71 CpGs may be modified by sex or dietary intake. Nearly half (35–50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs. Interpretation: Many smoking-related methylation sites were identified with Illumina’s EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases. Funding: Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.
Research for the use of biomarkers in osteoarthritis (OA) is promising, however, adequate discrimination between patients and controls may be hampered due to innate differences. We set out to ...identify loci influencing levels of serum cartilage oligomeric protein (sCOMP) and urinary C-telopeptide of type II collagen (uCTX-II).
Meta-analysis of genome-wide association studies was applied to standardised residuals of sCOMP (N=3316) and uCTX-II (N=4654) levels available in 6 and 7 studies, respectively, from TreatOA. Effects were estimated using a fixed-effects model. Six promising signals were followed up by de novo genotyping in the Cohort Hip and Cohort Knee study (N = 964). Subsequently, their role in OA susceptibility was investigated in large-scale genome-wide association studies meta-analyses for OA. Differential expression of annotated genes was assessed in cartilage.
Genome-wide significant association with sCOMP levels was found for a SNP within MRC1 (rs691461, p = 1.7 × 10(-12)) and a SNP within CSMD1 associated with variation in uCTX-II levels with borderline genome-wide significance (rs1983474, p = 8.5 × 10(-8)). Indication for association with sCOMP levels was also found for a locus close to the COMP gene itself (rs10038, p = 7.1 × 10(-6)). The latter SNP was subsequently found to be associated with hip OA whereas COMP expression appeared responsive to the OA pathophysiology in cartilage.
We have identified genetic loci affecting either uCTX-II or sCOMP levels. The genome wide significant association of MRC1 with sCOMP levels was found likely to act independent of OA subtypes. Increased sensitivity of biomarkers with OA may be accomplished by taking genetic variation into account.
The MTHFR C677T polymorphism is associated with mildly elevated homocysteine levels when folate and/or riboflavin status is low. Furthermore, a mildly elevated homocysteine level is a risk factor for ...osteoporotic fractures. We studied whether dietary intake of riboflavin and folate modifies the effects of the MTHFR C677T variant on fracture risk in 5035 men and women from the Rotterdam Study. We found that the MTHFR C677T variant interacts with dietary riboflavin intake to influence fracture risk in women.
Introduction: The MTHFR C677T polymorphism is associated with mildly elevated homocysteine (Hcy) levels in the presence of low folate and/or riboflavin status. A mildly elevated Hcy level was recently identified as a modifiable risk factor for osteoporotic fracture. We studied whether dietary intake of riboflavin and folate modifies the effects of the MTHFR C677T polymorphism on BMD and fracture risk.
Materials and Methods: We studied 5035 individuals from the Rotterdam Study, ≥55 yr of age, who had data available on MTHFR, nutrient intake, and fracture risk. We performed analysis on Hcy levels in a total of 666 individuals, whereas BMD data were present for 4646 individuals (2692women).
Results: In the total population, neither the MTHFR C677T polymorphism nor low riboflavin intake was associated with fracture risk and BMD. However, in the lowest quartile of riboflavin intake, female 677‐T homozygotes had a 1.8 (95% CI: 1.1‐2.9, p = 0.01) times higher risk for incident osteoporotic fractures and a 2.6 (95% CI: 1.3‐5.1, p = 0.01) times higher risk for fragility fractures compared with the 677‐CC genotype (interaction, p = 0.0002). This effect was not seen for baseline BMD in both men and women. No significant influence was found for dietary folate intake on the association between the MTHFR C677T genotype and fracture risk or BMD. In the lowest quartile of dietary riboflavin intake, T‐homozygous individuals (men and women combined) had higher (22.5%) Hcy levels compared with C‐homozygotes (mean difference = 3.44 νM, p = 0. 01; trend, p = 0.02).
Conclusions: In this cohort of elderly whites, the MTHFR C677T variant interacts with dietary riboflavin intake to influence fracture risk in women.
We tested whether DNA-methylation profiles account for inter-individual variation in body mass index (BMI) and height and whether they predict these phenotypes over and above genetic factors. Genetic ...predictors were derived from published summary results from the largest genome-wide association studies on BMI (n ∼ 350,000) and height (n ∼ 250,000) to date. We derived methylation predictors by estimating probe-trait effects in discovery samples and tested them in external samples. Methylation profiles associated with BMI in older individuals from the Lothian Birth Cohorts (LBCs, n = 1,366) explained 4.9% of the variation in BMI in Dutch adults from the LifeLines DEEP study (n = 750) but did not account for any BMI variation in adolescents from the Brisbane Systems Genetic Study (BSGS, n = 403). Methylation profiles based on the Dutch sample explained 4.9% and 3.6% of the variation in BMI in the LBCs and BSGS, respectively. Methylation profiles predicted BMI independently of genetic profiles in an additive manner: 7%, 8%, and 14% of variance of BMI in the LBCs were explained by the methylation predictor, the genetic predictor, and a model containing both, respectively. The corresponding percentages for LifeLines DEEP were 5%, 9%, and 13%, respectively, suggesting that the methylation profiles represent environmental effects. The differential effects of the BMI methylation profiles by age support previous observations of age modulation of genetic contributions. In contrast, methylation profiles accounted for almost no variation in height, consistent with a mainly genetic contribution to inter-individual variation. The BMI results suggest that combining genetic and epigenetic information might have greater utility for complex-trait prediction.
Although it is assumed that epigenetic mechanisms, such as changes in DNA methylation (DNAm), underlie the relationship between adverse intrauterine conditions and adult metabolic health, evidence ...from human studies remains scarce. Therefore, we evaluated whether DNAm in whole blood mediated the association between prenatal famine exposure and metabolic health in 422 individuals exposed to famine in utero and 463 (sibling) controls. We implemented a two-step analysis, namely, a genome-wide exploration across 342,596 cytosine-phosphate-guanine dinucleotides (CpGs) for potential mediators of the association between prenatal famine exposure and adult body mass index (BMI), serum triglycerides (TG), or glucose concentrations, which was followed by formal mediation analysis. DNAm mediated the association of prenatal famine exposure with adult BMI and TG but not with glucose. DNAm at
(cg09349128), a gene involved in energy metabolism, mediated 13.4% 95% confidence interval (CI), 5 to 28% of the association between famine exposure and BMI. DNAm at six CpGs, including
(cg19693031), influencing β cell function, and
(cg07397296), affecting lipid metabolism, together mediated 80% (95% CI, 38.5 to 100%) of the association between famine exposure and TG. Analyses restricted to those exposed to famine during early gestation identified additional CpGs mediating the relationship with TG near
(glycolysis) and
(adipogenesis). DNAm at the CpGs involved was associated with gene expression in an external data set and correlated with DNAm levels in fat depots in additional postmortem data. Our data are consistent with the hypothesis that epigenetic mechanisms mediate the influence of transient adverse environmental factors in early life on long-term metabolic health. The specific mechanism awaits elucidation.
Recent genome-wide association studies (GWASs) identified the first genetic loci associated with attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD). The next step is ...to use these results to increase our understanding of the biological mechanisms involved. Most of the identified variants likely influence gene regulation. The aim of the current study is to shed light on the mechanisms underlying the genetic signals and prioritize genes by integrating GWAS results with gene expression and DNA methylation (DNAm) levels.
We applied summary-data–based Mendelian randomization to integrate ADHD and ASD GWAS data with fetal brain expression and methylation quantitative trait loci, given the early onset of these disorders. We also analyzed expression and methylation quantitative trait loci datasets of adult brain and blood, as these provide increased statistical power. We subsequently used summary-data–based Mendelian randomization to investigate if the same variant influences both DNAm and gene expression levels.
We identified multiple gene expression and DNAm levels in fetal brain at chromosomes 1 and 17 that were associated with ADHD and ASD, respectively, through pleiotropy at shared genetic variants. The analyses in brain and blood showed additional associated gene expression and DNAm levels at the same and additional loci, likely because of increased statistical power. Several of the associated genes have not been identified in ADHD and ASD GWASs before.
Our findings identified the genetic variants associated with ADHD and ASD that likely act through gene regulation. This facilitates prioritization of candidate genes for functional follow-up studies.
Human iris color was one of the first traits for which Mendelian segregation was established. To date, the genetics of iris color is still not fully understood and is of interest, particularly in ...view of forensic applications. In three independent genome-wide association (GWA) studies of a total of 1406 persons and a genome-wide linkage study of 1292 relatives, all from the Netherlands, we found that the 15q13.1 region is the predominant region involved in human iris color. There were no other regions showing consistent genome-wide evidence for association and linkage to iris color. Single nucleotide polymorphisms (SNPs) in the
HERC2 gene and, to a lesser extent, in the neighboring
OCA2 gene were independently associated to iris color variation.
OCA2 has been implicated in iris color previously. A replication study within two populations confirmed that the
HERC2 gene is a new and significant determinant of human iris color variation, in addition to
OCA2. Furthermore,
HERC2 rs916977 showed a clinal allele distribution across 23 European populations, which was significantly correlated to iris color variation. We suggest that genetic variants regulating expression of the
OCA2 gene exist in the
HERC2 gene or, alternatively, within the 11.7 kb of sequence between
OCA2 and
HERC2, and that most iris color variation in Europeans is explained by those two genes. Testing markers in the
HERC2-OCA2 region may be useful in forensic applications to predict eye color phenotypes of unknown persons of European genetic origin.
DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA ...methylation patterns in blood using large-scale population genomics data.
By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g., NFKBIE, CDCA7(L), and NLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation.
We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation.
Abstract 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.
Many workers are daily exposed to occupational agents like gases/fumes, mineral dust or biological dust, which could induce adverse health effects. Epigenetic mechanisms, such as DNA methylation, ...have been suggested to play a role. We therefore aimed to identify differentially methylated regions (DMRs) upon occupational exposures in never-smokers and investigated if these DMRs associated with gene expression levels. To determine the effects of occupational exposures independent of smoking, 903 never-smokers of the LifeLines cohort study were included. We performed three genome-wide methylation analyses (Illumina 450 K), one per occupational exposure being gases/fumes, mineral dust and biological dust, using robust linear regression adjusted for appropriate confounders. DMRs were identified using comb-p in Python. Results were validated in the Rotterdam Study (233 never-smokers) and methylation-expression associations were assessed using Biobank-based Integrative Omics Study data (n = 2802). Of the total 21 significant DMRs, 14 DMRs were associated with gases/fumes and 7 with mineral dust. Three of these DMRs were associated with both exposures (RPLP1 and LINC02169 (2×)) and 11 DMRs were located within transcript start sites of gene expression regulating genes. We replicated two DMRs with gases/fumes (VTRNA2-1 and GNAS) and one with mineral dust (CCDC144NL). In addition, nine gases/fumes DMRs and six mineral dust DMRs significantly associated with gene expression levels. Our data suggest that occupational exposures may induce differential methylation of gene expression regulating genes and thereby may induce adverse health effects. Given the millions of workers that are exposed daily to occupational exposures, further studies on this epigenetic mechanism and health outcomes are warranted.