Epigenetic Signatures of Cigarette Smoking

DNA Methylation Patterns Associated with Tinnitus in Young Adults-A Pilot Study

PURPOSE

Tinnitus, the perception of sound without any external sound source, is a prevalent hearing health concern. Mounting evidence suggests that a confluence of genetic, environmental, and lifestyle factors can influence the pathogenesis of tinnitus. We hypothesized that alteration in DNA methylation, an epigenetic modification that occurs at cytosines of cytosine-phosphate-guanine (CpG) dinucleotide sites, where a methyl group from S-adenyl methionine gets transferred to the fifth carbon of the cytosine, could contribute to tinnitus. DNA methylation patterns are tissue-specific, but the tissues involved in tinnitus are not easily accessible in humans. This pilot study used saliva as a surrogate tissue to identify differentially methylated CpG regions (DMRs) associated with tinnitus. The study was conducted on healthy young adults reporting bilateral continuous chronic tinnitus to limit the influence of age-related confounding factors and health-related comorbidities.

METHODS

The present study evaluated the genome-wide methylation levels from saliva-derived DNA samples from 24 healthy young adults with bilateral continuous chronic tinnitus (> 1 year) and 24 age, sex, and ethnicity-matched controls with no tinnitus. Genome-wide DNA methylation was evaluated for > 850,000 CpG sites using the Infinium Human Methylation EPIC BeadChip. The association analysis used the Bumphunter algorithm on 23 cases and 20 controls meeting the quality control standards. The methylation level was expressed as the area under the curve of CpG sites within DMRs.The FDR-adjusted p-value threshold of 0.05 was used to identify statistically significant DMRs associated with tinnitus.

RESULTS

We obtained 25 differentially methylated regions (DMRs) associated with tinnitus. Genes within or in the proximity of the hypermethylated DMRs related to tinnitus included LCLAT1, RUNX1, RUFY1, NUDT12, TTC23, SLC43A2, C4orf27 (STPG2), and EFCAB4B. Genes within or in the proximity of hypomethylated DMRs associated with tinnitus included HLA-DPB2, PM20D1, TMEM18, SNTG2, MUC4, MIR886, MIR596, TXNRD1, EID3, SDHAP3, HLA-DPB2, LASS3 (CERS3), C10orf11 (LRMDA), HLA-DQB1, NADK, SZRD1, MFAP2, NUP210L, TPM3, INTS9, and SLC2A14. The burden of genetic variation could explain the differences in the methylation levels for DMRs involving HLA-DPB2, HLA-DQB1, and MUC4, indicating the need for replication in large independent cohorts.

CONCLUSION

Consistent with the literature on comorbidities associated with tinnitus, we identified genes within or close to DMRs involved in auditory functions, chemical dependency, cardiovascular diseases, psychiatric conditions, immune disorders, and metabolic syndromes. These results indicate that epigenetic mechanisms could influence tinnitus, and saliva can be a good surrogate for identifying the epigenetic underpinnings of tinnitus in humans. Further research with a larger sample size is needed to identify epigenetic biomarkers and investigate their influence on the phenotypic expression of tinnitus.

Maternal smoking DNA methylation risk score associated with health outcomes in offspring of European and South Asian ancestry

Background

Maternal smoking has been linked to adverse health outcomes in newborns but the extent to which it impacts newborn health has not been quantified through an aggregated cord blood DNA methylation (DNAm) score. Here, we examine the feasibility of using cord blood DNAm scores leveraging large external studies as discovery samples to capture the epigenetic signature of maternal smoking and its influence on newborns in White European and South Asian populations.

Methods

We first examined the association between individual CpGs and cigarette smoking during pregnancy, and smoking exposure in two White European birth cohorts (n=744). Leveraging established CpGs for maternal smoking, we constructed a cord blood epigenetic score of maternal smoking that was validated in one of the European-origin cohorts (n=347). This score was then tested for association with smoking status, secondary smoking exposure during pregnancy, and health outcomes in offspring measured after birth in an independent White European (n=397) and a South Asian birth cohort (n=504).

Results

Several previously reported genes for maternal smoking were supported, with the strongest and most consistent association signal from the GFI1 gene (6 CpGs with p<5 × 10-5). The epigenetic maternal smoking score was strongly associated with smoking status during pregnancy (OR = 1.09 [1.07, 1.10], p=5.5 × 10-33) and more hours of self-reported smoking exposure per week (1.93 [1.27, 2.58], p=7.8 × 10-9) in White Europeans. However, it was not associated with self-reported exposure (p>0.05) among South Asians, likely due to a lack of smoking in this group. The same score was consistently associated with a smaller birth size (-0.37±0.12 cm, p=0.0023) in the South Asian cohort and a lower birth weight (-0.043±0.013 kg, p=0.0011) in the combined cohorts.

Conclusions

This cord blood epigenetic score can help identify babies exposed to maternal smoking and assess its long-term impact on growth. Notably, these results indicate a consistent association between the DNAm signature of maternal smoking and a small body size and low birth weight in newborns, in both White European mothers who exhibited some amount of smoking and in South Asian mothers who themselves were not active smokers.

Funding

This study was funded by the Canadian Institutes of Health Research Metabolomics Team Grant: MWG-146332.

Proteomic aging clock predicts mortality and risk of common age-related diseases in diverse populations

Circulating plasma proteins play key roles in human health and can potentially be used to measure biological age, allowing risk prediction for age-related diseases, multimorbidity and mortality. Here we developed a proteomic age clock in the UK Biobank (n = 45,441) using a proteomic platform comprising 2,897 plasma proteins and explored its utility to predict major disease morbidity and mortality in diverse populations. We identified 204 proteins that accurately predict chronological age (Pearson r = 0.94) and found that proteomic aging was associated with the incidence of 18 major chronic diseases (including diseases of the heart, liver, kidney and lung, diabetes, neurodegeneration and cancer), as well as with multimorbidity and all-cause mortality risk. Proteomic aging was also associated with age-related measures of biological, physical and cognitive function, including telomere length, frailty index and reaction time. Proteins contributing most substantially to the proteomic age clock are involved in numerous biological functions, including extracellular matrix interactions, immune response and inflammation, hormone regulation and reproduction, neuronal structure and function and development and differentiation. In a validation study involving biobanks in China (n = 3,977) and Finland (n = 1,990), the proteomic age clock showed similar age prediction accuracy (Pearson r = 0.92 and r = 0.94, respectively) compared to its performance in the UK Biobank. Our results demonstrate that proteomic aging involves proteins spanning multiple functional categories and can be used to predict age-related functional status, multimorbidity and mortality risk across geographically and genetically diverse populations.

A Bivariate Twin Study of Lifetime cannabis Initiation and Lifetime Regular Tobacco Smoking Across Three Different Countries

Regular cigarette smoking and cannabis consumption are strongly positively related to each other, yet few studies explore their underlying variation and covariation. We evaluated the genetic and environmental decomposition of variance and covariance of these two traits in twin data from three countries with different social norms and legislation. Data from the Netherlands Twin Register, FinnTwin12/16, and the Minnesota Center for Twin Family Research (total N = 21,617) were analyzed in bivariate threshold models of lifetime regular smoking initiation (RSI) and lifetime cannabis initiation (CI). We ran unstratified models and models stratified by sex and country. Prevalence of RSI was lowest in the Netherlands and prevalence of CI was highest in Minnesota. In the unstratified model, genetic (A) and common environmental factors (C) contributed substantially to the liabilities of RSI (A = 0.47, C = 0.34) and CI (A = 0.28, C = 0.51). The two liabilities were significantly phenotypically (rP = 0.56), genetically (rA = 0.74), and environmentally correlated in the unstratified model (rC = 0.47and rE = 0.48, representing correlations between common and unique environmental factors). The magnitude of phenotypic correlation between liabilities varied by country but not sex (Minnesota rP ~ 0.70, Netherlands rP ~ 0.59, Finland rP ~ 0.45). Comparisons of decomposed correlations could not be reliably tested in the stratified models. The prevalence and association of RSI and CI vary by sex and country. These two behaviors are correlated because there is genetic and environmental overlap between their underlying latent liabilities. There is heterogeneity in the genetic architecture of these traits across country.

Epigenetic effects of cannabis: A systematic scoping review of behavioral and emotional symptoms associated with cannabis use and exocannabinoid exposure

BACKGROUND

Recent research suggests that epigenetic modifications may mediate the behavioral effects of cannabis, influencing exocannabinnoids' long term effects in cognitive function and its role in the emergence of psychotic symptoms.

BASIC PROCEDURES

In this systematic scoping review, we assessed the current evidence of epigenetic effects associated with the use of cannabis or exocannabinoid administration and their relationship with behavioral and emotional symptoms. We searched PubMed, Cochrane CENTRAL, and Web of Science, up to January 2022, using the terms "cannabis" and "epigenetics." The search yielded 178 articles, of which 43 underwent full article revision; 37 articles were included in the review.

MAIN FINDINGS

The gathered evidence included observational cross-sectional studies conducted on human subjects and experimental designs using animal models that conveyed disparity in administration dosage, methods of cannabis use assessment and targeted epigenetic mechanisms. Nine studies performed epigenome-wide analysis with identification of differentially methylated sites; most of these studies found a global hypomethylation, and enrichment in genes related to cellular survival and neurodevelopment. Other studies assessed methylation at specific genes and found that cannabis exposure was associated with reduced methylation at Cg05575921, DNMT1, DRD2, COMT, DLGAP2, Arg1, STAT3, MGMT, and PENK, while hypermethylation was found at DNMT3a/b, NCAM1, and AKT1.

CONCLUSIONS

The review found evidence of an exocannabinoid-induced epigenetic changes that modulate depressive-anxious, psychotic, and addictive behavioural phenotypes. Further studies will require dosage exposure/administration uniformization and a customized pool of genes to assess their suitability as biomarkers for psychiatric diseases.

Non-cigarette tobacco products, aryl-hydrocarbon receptor repressor gene methylation and smoking-related health outcomes

INTRODUCTION

Cigarette smoking leads to altered DNA methylation at the aryl-hydrocarbon receptor repressor (AHRR) gene. However, it remains unknown whether pipe or cigar smoking is associated with AHRR methylation. We evaluated associations of non-cigarette tobacco use with AHRR methylation and determined if AHRR methylation was associated with smoking-related health outcomes.

METHODS

Data were pooled across four population-based cohorts that enrolled participants from 1985 to 2002. Tobacco exposures were evaluated using smoking questionnaires. AHRR cg05575921 methylation was measured in peripheral blood leucocyte DNA. Spirometry and respiratory symptoms were evaluated at the time of methylation measurements and in subsequent visits. Vital status was monitored using the National Death Index.

RESULTS

Among 8252 adults (mean age 56.7±10.3 years, 58.1% women, 40.6% black), 4857 (58.9%) participants used cigarettes and 634 (7.7%) used non-cigarette tobacco products. Exclusive use of non-cigarette tobacco products was independently associated with lower AHRR methylation (-2.44 units, 95% CI -4.42 to -0.45), though to a lesser extent than exclusive use of cigarettes (-6.01 units, 95% CI -6.01 to -4.10). Among participants who exclusively used non-cigarette tobacco products, reduced AHRR methylation was associated with increased respiratory symptom burden (OR 1.60, 95% CI 1.03 to 2.68) and higher all-cause mortality (log-rank p=0.02).

CONCLUSION

Pipe and cigar smoking were independently associated with lower AHRR methylation in a multiethnic cohort of US adults. Among users of non-cigarette tobacco products, lower AHRR methylation was associated with poor respiratory health outcomes and increased mortality. AHRR methylation may identify non-cigarette tobacco users with an increased risk of adverse smoking-related health outcomes.

A novel hypothesis-generating approach for detecting phenotypic associations using epigenetic data

Aim: Hypotheses about what phenotypes to include in causal analyses, that in turn can have clinical and policy implications, can be guided by hypothesis-free approaches leveraging the epigenome, for example. Materials & methods: Minimally adjusted epigenome-wide association studies (EWAS) using ALSPAC data were performed for example conditions, dysmenorrhea and heavy menstrual bleeding (HMB). Differentially methylated CpGs were searched in the EWAS Catalog and associated traits identified. Traits were compared between those with and without the example conditions in ALSPAC. Results: Seven CpG sites were associated with dysmenorrhea and two with HMB. Smoking and adverse childhood experience score were associated with both conditions in the hypothesis-testing phase. Conclusion: Hypothesis-generating EWAS can help identify associations for future analyses.

Epigenome-wide DNA Methylation Association Study of CHIP Provides Insight into Perturbed Gene Regulation

With age, hematopoietic stem cells can acquire somatic mutations in leukemogenic genes that confer a proliferative advantage in a phenomenon termed "clonal hematopoiesis of indeterminate potential" (CHIP). How these mutations confer a proliferative advantage and result in increased risk for numerous age-related diseases remains poorly understood. We conducted a multiracial meta-analysis of epigenome-wide association studies (EWAS) of CHIP and its subtypes in four cohorts (N=8196) to elucidate the molecular mechanisms underlying CHIP and illuminate how these changes influence cardiovascular disease risk. The EWAS findings were functionally validated using human hematopoietic stem cell (HSC) models of CHIP. A total of 9615 CpGs were associated with any CHIP, 5990 with DNMT3A CHIP, 5633 with TET2 CHIP, and 6078 with ASXL1 CHIP (P <1×10-7). CpGs associated with CHIP subtypes overlapped moderately, and the genome-wide DNA methylation directions of effect were opposite for TET2 and DNMT3A CHIP, consistent with their opposing effects on global DNA methylation. There was high directional concordance between the CpGs shared from the meta-EWAS and human edited CHIP HSCs. Expression quantitative trait methylation analysis further identified transcriptomic changes associated with CHIP-associated CpGs. Causal inference analyses revealed 261 CHIP-associated CpGs associated with cardiovascular traits and all-cause mortality (FDR adjusted p-value <0.05). Taken together, our study sheds light on the epigenetic changes impacted by CHIP and their associations with age-related disease outcomes. The novel genes and pathways linked to the epigenetic features of CHIP may serve as therapeutic targets for preventing or treating CHIP-mediated diseases.

Map of epigenetic age acceleration: A worldwide analysis

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Effect of smoking on thrombotic antiphospholipid syndrome: a 10-year prospective cohort study

OBJECTIVES

Cigarette smoking is an established risk factor for autoimmune diseases. However, whether smoking plays a clear role in thrombotic APS (TAPS) has not been determined. We aimed to investigate the effects of smoking on the clinical characteristics and prognosis of TAPS.

METHODS

This was a prospective cohort study from 2013 to 2022. During the study period, 297 patients were diagnosed with TAPS, including 82 smokers and 215 non-smokers. After propensity score matching, 57 smokers and 57 non-smokers matched by age and sex were analysed.

RESULTS

Overall, smokers with TAPS had more cardiovascular risk factors than non-smokers, including hypertension (36.59% vs 14.42%, P < 0.001), obesity (15.85% vs 7.44%, P = 0.029), dyslipidaemia (64.63% vs 48.37%, P = 0.012) and hyperhomocysteinaemia (62.20% vs 36.28%, P < 0.001). Arterial thrombotic events were more common in smokers at diagnosis (62.20% vs 46.05%, P = 0.013), especially myocardial infarction, visceral thrombosis and peripheral vascular thrombosis. After matching, smokers showed balanced cardiovascular risk factors with non-smokers at baseline, but retained a higher prevalence of arterial thrombosis (59.65% vs 33.33%, P = 0.005), mainly distributed in cerebral vascular, cardiovascula and retinal vascular territories. During follow-up, smokers presented a tendency for more recurrent arterial thrombosis and less recurrent venous thrombosis. Smokers had significantly poorer outcomes for organ damage with higher Damage Index for APS score (median 2.00 vs 1.00, P = 0.008), especially in the cardiovascular (26.32% vs 3.51%, P = 0.001), gastrointestinal (15.79% vs 1.75%, P = 0.016) and ophthalmologic (10.53% vs 00.00%, P = 0.027) systems.

CONCLUSION

Smoking is related to increased arterial events and poor prognosis in TAPS patients. Patients with TAPS should be fully encouraged to avoid smoking.

penalizedclr: an R package for penalized conditional logistic regression for integration of multiple omics layers

BACKGROUND

The matched case-control design, up until recently mostly pertinent to epidemiological studies, is becoming customary in biomedical applications as well. For instance, in omics studies, it is quite common to compare cancer and healthy tissue from the same patient. Furthermore, researchers today routinely collect data from various and variable sources that they wish to relate to the case-control status. This highlights the need to develop and implement statistical methods that can take these tendencies into account.

RESULTS

We present an R package penalizedclr, that provides an implementation of the penalized conditional logistic regression model for analyzing matched case-control studies. It allows for different penalties for different blocks of covariates, and it is therefore particularly useful in the presence of multi-source omics data. Both L1 and L2 penalties are implemented. Additionally, the package implements stability selection for variable selection in the considered regression model.

CONCLUSIONS

The proposed method fills a gap in the available software for fitting high-dimensional conditional logistic regression models accounting for the matched design and block structure of predictors/features. The output consists of a set of selected variables that are significantly associated with case-control status. These variables can then be investigated in terms of functional interpretation or validation in further, more targeted studies.

Twin-based Mendelian Randomization Analyses Highlight Smoking's Effects on Blood DNA Methylation, with Putative Reverse Causation

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.

The influence of epigenetic biological age on key complications and outcomes in aneurysmal subarachnoid haemorrhage

BACKGROUND

We aimed to investigate the association between DNA-methylation biological age (B-age) calculated as age acceleration (ageAcc) and key aneurysmal subarachnoid haemorrhage (aSAH) complications such as vasospasm, delayed cerebral ischaemia (DCI), poor outcome, and mortality.

METHODS

We conducted a prospective study involving 277 patients with aSAH. B-age was determined in whole blood samples using five epigenetic clocks: Hannum's, Horvath's, Levine's and both versions of Zhang's clocks. Age acceleration was calculated as the residual obtained from regressing out the effect of C-age on the mismatch between C-age and B-age. We then tested the association between ageAcc and vasospasm, DCI and 12-month poor outcome (mRS 3-5) and mortality using linear regression models adjusted for confounders.

RESULTS

Average C-age was 55.0 years, with 66.8% being female. Vasospasm occurred in 143 cases (51.6%), DCI in 70 (25.3%) and poor outcomes in 99 (35.7%), with a mortality rate of 20.6%. Lower ageAcc was linked to vasospasm in Horvath's and Levine's clocks, whereas increased ageAcc was associated with 12-month mortality in Hannum's clock. No significant differences in ageAcc were found for DCI or poor outcome at 12 months with other clocks.

CONCLUSIONS

Our study indicates that B-age is independently associated with vasospasm and 12-month mortality in patients with aSAH. These findings underscore the potential role of epigenetics in understanding the pathophysiology of aSAH-related complications and outcomes.

Smoking-Induced DNA Hydroxymethylation Signature Is Less Pronounced than True DNA Methylation: The Population-Based KORA Fit Cohort

Despite extensive research on 5-methylcytosine (5mC) in relation to smoking, there has been limited exploration into the interaction between smoking and 5-hydroxymethylcytosine (5hmC). In this study, total DNA methylation (5mC+5hmC), true DNA methylation (5mC) and hydroxymethylation (5hmC) levels were profiled utilizing conventional bisulphite (BS) and oxidative bisulphite (oxBS) treatment, measured with the Illumina Infinium Methylation EPIC BeadChip. An epigenome-wide association study (EWAS) of 5mC+5hmC methylation revealed a total of 38,575 differentially methylated positions (DMPs) and 2023 differentially methylated regions (DMRs) associated with current smoking, along with 82 DMPs and 76 DMRs associated with former smoking (FDR-adjusted p < 0.05). Additionally, a focused examination of 5mC identified 33 DMPs linked to current smoking and 1 DMP associated with former smoking (FDR-adjusted p < 0.05). In the 5hmC category, eight DMPs related to current smoking and two DMPs tied to former smoking were identified, each meeting a suggestive threshold (p < 1 × 10-5). The substantial number of recognized DMPs, including 5mC+5hmC (7069/38,575, 2/82), 5mC (0/33, 1/1), and 5hmC (2/8, 0/2), have not been previously reported. Our findings corroborated previously established methylation positions and revealed novel candidates linked to tobacco smoking. Moreover, the identification of hydroxymethylated CpG sites with suggestive links provides avenues for future research.

Epigenetic Age Mediates the Association of Life's Essential 8 With Cardiovascular Disease and Mortality

BACKGROUND

Life's Essential 8 (LE8) is an enhanced metric for cardiovascular health. The interrelations among LE8, biomarkers of aging, and disease risks are unclear.

METHODS AND RESULTS

LE8 score was calculated for 5682 Framingham Heart Study participants. We implemented 4 DNA methylation-based epigenetic age biomarkers, with older epigenetic age hypothesized to represent faster biological aging, and examined whether these biomarkers mediated the associations between the LE8 score and cardiovascular disease (CVD), CVD-specific mortality, and all-cause mortality. We found that a 1 SD increase in the LE8 score was associated with a 35% (95% CI, 27-41; P=1.8E-15) lower risk of incident CVD, a 36% (95% CI, 24-47; P=7E-7) lower risk of CVD-specific mortality, and a 29% (95% CI, 22-35; P=7E-15) lower risk of all-cause mortality. These associations were partly mediated by epigenetic age biomarkers, particularly the GrimAge and the DunedinPACE scores. The potential mediation effects by epigenetic age biomarkers tended to be more profound in participants with higher genetic risk for older epigenetic age, compared with those with lower genetic risk. For example, in participants with higher GrimAge polygenic scores (greater than median), the mean proportion of mediation was 39%, 39%, and 78% for the association of the LE8 score with incident CVD, CVD-specific mortality, and all-cause mortality, respectively. No significant mediation was observed in participants with lower GrimAge polygenic score.

CONCLUSIONS

DNA methylation-based epigenetic age scores mediate the associations between the LE8 score and incident CVD, CVD-specific mortality, and all-cause mortality, particularly in individuals with higher genetic predisposition for older epigenetic age.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Cigarette Smoking and E-cigarette Use Induce Shared DNA Methylation Changes Linked to Carcinogenesis

Tobacco use is a major modifiable risk factor for adverse health outcomes, including cancer, and elicits profound epigenetic changes thought to be associated with long-term cancer risk. While electronic cigarettes (e-cigarettes) have been advocated as harm reduction alternatives to tobacco products, recent studies have revealed potential detrimental effects, highlighting the urgent need for further research into the molecular and health impacts of e-cigarettes. Here, we applied computational deconvolution methods to dissect the cell- and tissue-specific epigenetic effects of tobacco or e-cigarette use on DNA methylation (DNAme) in over 3,500 buccal/saliva, cervical, or blood samples, spanning epithelial and immune cells at directly and indirectly exposed sites. The 535 identified smoking-related DNAme loci [cytosine-phosphate-guanine sites (CpG)] clustered into four functional groups, including detoxification or growth signaling, based on cell type and anatomic site. Loci hypermethylated in buccal epithelial cells of smokers associated with NOTCH1/RUNX3/growth factor receptor signaling also exhibited elevated methylation in cancer tissue and progressing lung carcinoma in situ lesions, and hypermethylation of these sites predicted lung cancer development in buccal samples collected from smokers up to 22 years prior to diagnosis, suggesting a potential role in driving carcinogenesis. Alarmingly, these CpGs were also hypermethylated in e-cigarette users with a limited smoking history. This study sheds light on the cell type-specific changes to the epigenetic landscape induced by smoking-related products.

SIGNIFICANCE

The use of both cigarettes and e-cigarettes elicits cell- and exposure-specific epigenetic effects that are predictive of carcinogenesis, suggesting caution when broadly recommending e-cigarettes as aids for smoking cessation.

High-dimensional mediation analysis for continuous outcome with confounders using overlap weighting method in observational epigenetic study

BACKGROUND

Mediation analysis is a powerful tool to identify factors mediating the causal pathway of exposure to health outcomes. Mediation analysis has been extended to study a large number of potential mediators in high-dimensional data settings. The presence of confounding in observational studies is inevitable. Hence, it's an essential part of high-dimensional mediation analysis (HDMA) to adjust for the potential confounders. Although the propensity score (PS) related method such as propensity score regression adjustment (PSR) and inverse probability weighting (IPW) has been proposed to tackle this problem, the characteristics with extreme propensity score distribution of the PS-based method would result in the biased estimation.

METHODS

In this article, we integrated the overlapping weighting (OW) technique into HDMA workflow and proposed a concise and powerful high-dimensional mediation analysis procedure consisting of OW confounding adjustment, sure independence screening (SIS), de-biased Lasso penalization, and joint-significance testing underlying the mixture null distribution. We compared the proposed method with the existing method consisting of PS-based confounding adjustment, SIS, minimax concave penalty (MCP) variable selection, and classical joint-significance testing.

RESULTS

Simulation studies demonstrate the proposed procedure has the best performance in mediator selection and estimation. The proposed procedure yielded the highest true positive rate, acceptable false discovery proportion level, and lower mean square error. In the empirical study based on the GSE117859 dataset in the Gene Expression Omnibus database using the proposed method, we found that smoking history may lead to the estimated natural killer (NK) cell level reduction through the mediation effect of some methylation markers, mainly including methylation sites cg13917614 in CNP gene and cg16893868 in LILRA2 gene.

CONCLUSIONS

The proposed method has higher power, sufficient false discovery rate control, and precise mediation effect estimation. Meanwhile, it is feasible to be implemented with the presence of confounders. Hence, our method is worth considering in HDMA studies.

Differential Impact of Smoking on Methylome and Transcriptome in Crohn's Disease and Ulcerative Colitis

BACKGROUND

Smoking is an environmental factor that differentially impacts Crohn's disease (CD) and ulcerative colitis (UC). The mechanism of impact of smoking on disease risk and clinical outcomes remains to be established.

METHODS

This study used a prospective cohort of patients with CD or UC. Self-reported smoking status was validated using serum cotinine measurement. We profiled methylation changes in peripheral blood using the Illumina Methylation BeadChip. Transcriptomic profiling was performed on ileal and colonic tissue using an Illumina TruSeq platform. We compared the methylation and transcriptional changes in current, former, and never smokers stratified by disease type.

RESULTS

Our cohort included 200 patients with CD or UC with methylation profiles and 160 with transcriptomic data. The mean serum cotinine level was higher in current compared with former or never smokers. Epigenetic changes common to both CD and UC included hypomethylation at AHRR. Smoking-associated MGAT3 hypomethylation was associated with severe disease course only in UC, while IER3 hypomethylation was associated with worse course only in CD. Smoking downregulated several inflammatory pathways in UC. Current smoking in CD but not in UC was associated with upregulation of several genes mediating Paneth cell function. Genes with opposite direction of effects in CD and UC include HSD3B2 and GSTA1.

CONCLUSIONS

Our findings suggest both common and differential effects of cigarette smoking on CD and UC. Paneth cell dysfunction may mediate adverse impact of smoking on CD. Bile acid and oxidative stress pathways may be relevant for the differential effect of smoking on CD and UC.

Methods for the Analysis of Multiple Epigenomic Mediators in Environmental Epidemiology

PURPOSE OF REVIEW

Epigenetic changes can be highly influenced by environmental factors and have in turn been proposed to influence chronic disease. Being able to quantify to which extent epigenomic processes are mediators of the association between environmental exposures and diseases is of interest for epidemiologic research. In this review, we summarize the proposed mediation analysis methods with applications to epigenomic data.

RECENT FINDINGS

The ultra-high dimensionality and high correlations that characterize omics data have hindered the precise quantification of mediated effects. Several methods have been proposed to deal with mediation in high-dimensional settings, including methods that incorporate dimensionality reduction techniques to the mediation algorithm. Although important methodological advances have been conducted in the previous years, key challenges such as the development of sensitivity analyses, dealing with mediator-mediator interactions, including environmental mixtures as exposures, or the integration of different omic data should be the focus of future methodological developments for epigenomic mediation analysis.

Quantifying the stochastic component of epigenetic aging

DNA methylation clocks can accurately estimate chronological age and, to some extent, also biological age, yet the process by which age-associated DNA methylation (DNAm) changes are acquired appears to be quasi-stochastic, raising a fundamental question: how much of an epigenetic clock's predictive accuracy could be explained by a stochastic process of DNAm change? Here, using DNAm data from sorted immune cells, we build realistic simulation models, subsequently demonstrating in over 22,770 sorted and whole-blood samples from 25 independent cohorts that approximately 66-75% of the accuracy underpinning Horvath's clock could be driven by a stochastic process. This fraction increases to 90% for the more accurate Zhang's clock, but is lower (63%) for the PhenoAge clock, suggesting that biological aging is reflected by nonstochastic processes. Confirming this, we demonstrate that Horvath's age acceleration in males and PhenoAge's age acceleration in severe coronavirus disease 2019 cases and smokers are not driven by an increased rate of stochastic change but by nonstochastic processes. These results significantly deepen our understanding and interpretation of epigenetic clocks.

Demographic characteristics and epigenetic biological aging among post-9/11 veterans: Associations of DunedinPACE with sex, race, and age

Measures of epigenetic aging derived from DNA methylation (DNAm) have enabled the assessment of biological aging in new populations and cohorts. In the present study, we used an epigenetic measure of aging, DunedinPACE, to examine rates of aging across demographic groups in a sample of 2,309 United States military veterans from the VISN 6 MIRECC's Post-Deployment Mental Health Study. As assessed by DunedinPACE, female veterans were aging faster than male veterans (β = 0.39, 95 % CI [0.29, 0.48], p < .001), non-Hispanic Black veterans were aging faster than non-Hispanic White veterans (β = 0.58, 95 % CI [0.50, 0.66], p < .001), and older veterans were biologically aging faster than younger veterans (β = 0.21, 95 % CI [0.18, 0.25], p < .001). In secondary analyses, these differences in rates of aging were not explained by a variety of biopsychosocial covariates. In addition, the percentage of European genetic admixture in non-Hispanic Black veterans was not associated with DunedinPACE. Our findings suggest that female and non-Hispanic Black veterans are at greater risk of accelerated aging among post-9/11 veterans. Interventions that slow aging might provide relatively greater benefit among veterans comprising these at-risk groups.

microRNAs, oxidative stress, and genotoxicity as the main inducers in the pathobiology of cancer development

Cancer is one of the most serious leading causes of death in the world. Many eclectic factors are involved in cancer progression including genetic and epigenetic alongside environmental ones. In this account, the performance and fluctuations of microRNAs are significant in cancer diagnosis and treatment, particularly as diagnostic biomarkers in oncology. So, microRNAs manage and control the gene expression after transcription by mRNA degradation, or also they can inhibit their translation. Conspicuously, these molecular structures take part in controlling the cellular, physiological and pathological functions, which many of them can accomplish as tumor inhibitors or oncogenes. Relatively, Oxidative stress is defined as the inequality between the creation of reactive oxygen species (ROS) and the body's ability to detoxify the reactive mediators or repair the resulting injury. ROS and microRNAs have been recognized as main cancer promoters and possible treatment targets. Importantly, genotoxicity has been established as the primary reason for many diseases as well as several malignancies. The procedures have no obvious link with mutagenicity and influence the organization, accuracy of the information, or fragmentation of DNA. Conclusively, mutations in these patterns can lead to carcinogenesis. In this review article, we report the impressive and practical roles of microRNAs, oxidative stress, and genotoxicity in the pathobiology of cancer development in conjunction with their importance as reliable cancer biomarkers and their association with circulating miRNA, exosomes and exosomal miRNAs, RNA remodeling, DNA methylation, and other molecular elements in oncology.

Development and Validation of a Novel Placental DNA Methylation Biomarker of Maternal Smoking during Pregnancy in the ECHO Program

BACKGROUND

Maternal cigarette smoking during pregnancy (MSDP) is associated with numerous adverse health outcomes in infants and children with potential lifelong consequences. Negative effects of MSDP on placental DNA methylation (DNAm), placental structure, and function are well established.

OBJECTIVE

Our aim was to develop biomarkers of MSDP using DNAm measured in placentas (N = 96 ), collected as part of the Vitamin C to Decrease the Effects of Smoking in Pregnancy on Infant Lung Function double-blind, placebo-controlled randomized clinical trial conducted between 2012 and 2016. We also aimed to develop a digital polymerase chain reaction (PCR) assay for the top ranking cytosine-guanine dinucleotide (CpG) so that large numbers of samples can be screened for exposure at low cost.

METHODS

We compared the ability of four machine learning methods [logistic least absolute shrinkage and selection operator (LASSO) regression, logistic elastic net regression, random forest, and gradient boosting machine] to classify MSDP based on placental DNAm signatures. We developed separate models using the complete EPIC array dataset and on the subset of probes also found on the 450K array so that models exist for both platforms. For comparison, we developed a model using CpGs previously associated with MSDP in placenta. For each final model, we used model coefficients and normalized beta values to calculate placental smoking index (PSI) scores for each sample. Final models were validated in two external datasets: the Extremely Low Gestational Age Newborn observational study, N = 426 ; and the Rhode Island Children's Health Study, N = 237 .

RESULTS

Logistic LASSO regression demonstrated the highest performance in cross-validation testing with the lowest number of input CpGs. Accuracy was greatest in external datasets when using models developed for the same platform. PSI scores in smokers only (n = 72 ) were moderately correlated with maternal plasma cotinine levels. One CpG (cg27402634), with the largest coefficient in two models, was measured accurately by digital PCR compared with measurement by EPIC array (R 2 = 0.98 ).

DISCUSSION

To our knowledge, we have developed the first placental DNAm-based biomarkers of MSDP with broad utility to studies of prenatal disease origins. https://doi.org/10.1289/EHP13838.

Integrative analysis of genetics, epigenetics and RNA expression data reveal three susceptibility loci for smoking behavior in Chinese Han population

Despite numerous studies demonstrate that genetics and epigenetics factors play important roles on smoking behavior, our understanding of their functional relevance and coordinated regulation remains largely unknown. Here we present a multiomics study on smoking behavior for Chinese smoker population with the goal of not only identifying smoking-associated functional variants but also deciphering the pathogenesis and mechanism underlying smoking behavior in this under-studied ethnic population. After whole-genome sequencing analysis of 1329 Chinese Han male samples in discovery phase and OpenArray analysis of 3744 samples in replication phase, we discovered that three novel variants located near FOXP1 (rs7635815), and between DGCR6 and PRODH (rs796774020), and in ARVCF (rs148582811) were significantly associated with smoking behavior. Subsequently cis-mQTL and cis-eQTL analysis indicated that these variants correlated significantly with the differential methylation regions (DMRs) or differential expressed genes (DEGs) located in the regions where these variants present. Finally, our in silico multiomics analysis revealed several hub genes, like DRD2, PTPRD, FOXP1, COMT, CTNNAP2, to be synergistic regulated each other in the etiology of smoking.

When smoke meets gut: deciphering the interactions between tobacco smoking and gut microbiota in disease development

Tobacco smoking is a prevalent and detrimental habit practiced worldwide, increasing the risk of various diseases, including chronic obstructive pulmonary disease (COPD), cardiovascular disease, liver disease, and cancer. Although previous research has explored the detrimental health effects of tobacco smoking, recent studies suggest that gut microbiota dysbiosis may play a critical role in these outcomes. Numerous tobacco smoke components, such as nicotine, are found in the gastrointestinal tract and interact with gut microbiota, leading to lasting impacts on host health and diseases. This review delves into the ways tobacco smoking and its various constituents influence gut microbiota composition and functionality. We also summarize recent advancements in understanding how tobacco smoking-induced gut microbiota dysbiosis affects host health. Furthermore, this review introduces a novel perspective on how changes in gut microbiota following smoking cessation may contribute to withdrawal syndrome and the degree of health improvements in smokers.

Associations of Epigenetic Age Acceleration With CVD Risks Across the Lifespan: The Bogalusa Heart Study

Although epigenetic age acceleration (EAA) might serve as a molecular signature of childhood cardiovascular disease (CVD) risk factors and further promote midlife subclinical CVD, few studies have comprehensively examined these life course associations. This study sought to test whether childhood CVD risk factors predict EAA in adulthood and whether EAA mediates the association between childhood CVD risks and midlife subclinical disease. Among 1,580 Bogalusa Heart Study participants, we estimated extrinsic EAA, intrinsic EAA, PhenoAge acceleration (PhenoAgeAccel), and GrimAge acceleration (GrimAgeAccel) during adulthood. We tested prospective associations of longitudinal childhood body mass index (BMI), blood pressure, lipids, and glucose with EAAs using linear mixed effects models. After confirming EAAs with midlife carotid intima-media thickness and carotid plaque, structural equation models examined mediating effects of EAAs on associations of childhood CVD risk factors with subclinical CVD measures. After stringent multiple testing corrections, each SD increase in childhood BMI was significantly associated with 0.6-, 0.9-, and 0.5-year increases in extrinsic EAA, PhenoAgeAccel, and GrimAgeAccel, respectively (P < 0.001 for all 3 associations). Likewise, each SD increase in childhood log-triglycerides was associated with 0.5- and 0.4-year increases in PhenoAgeAccel and GrimAgeAccel (P < 0.001 for both), respectively, whereas each SD increase in childhood high-density lipoprotein cholesterol was associated with a 0.3-year decrease in GrimAgeAccel (P = 0.002). Our findings indicate that PhenoAgeAccel mediates an estimated 27.4% of the association between childhood log-triglycerides and midlife carotid intima-media thickness (P = 0.022). Our data demonstrate that early life CVD risk factors may accelerate biological aging and promote subclinical atherosclerosis.

Blood DNA Methylation Analysis Reveals a Distinctive Epigenetic Signature of Vasospasm in Aneurysmal Subarachnoid Hemorrhage

Vasospasm is a potentially preventable cause of poor prognosis in patients with aneurysmal subarachnoid hemorrhage (aSAH). Epigenetics might provide insight on its molecular mechanisms. We aimed to analyze the association between differential DNA methylation (DNAm) and development of vasospasm. We conducted an epigenome-wide association study in 282 patients with aSAH admitted to our hospital. DNAm was assessed with the EPIC Illumina chip (> 850 K CpG sites) in whole-blood samples collected at hospital admission. We identified differentially methylated positions (DMPs) at the CpG level using Cox regression models adjusted for potential confounders, and then we used the DMP results to find differentially methylated regions (DMRs) and enriched biological pathways. A total of 145 patients (51%) experienced vasospasm. In the DMP analysis, we identified 31 CpGs associated with vasospasm at p-value < 10-5. One of them (cg26189827) was significant at the genome-wide level (p-value < 10-8), being hypermethylated in patients with vasospasm and annotated to SUGCT gene, mainly expressed in arteries. Region analysis revealed 13 DMRs, some of them annotated to interesting genes such as POU5F1, HLA-DPA1, RUFY1, and CYP1A1. Functional enrichment analysis showed the involvement of biological processes related to immunity, inflammatory response, oxidative stress, endothelial nitric oxide, and apoptosis. Our findings show, for the first time, a distinctive epigenetic signature of vasospasm in aSAH, establishing novel links with essential biological pathways, including inflammation, immune responses, and oxidative stress. Although further validation is required, our results provide a foundation for future research into the complex pathophysiology of vasospasm.

On epigenetic stochasticity, entropy and cancer risk

Epigenetic changes are known to accrue in normal cells as a result of ageing and cumulative exposure to cancer risk factors. Increasing evidence points towards age-related epigenetic changes being acquired in a quasi-stochastic manner, and that they may play a causal role in cancer development. Here, I describe the quasi-stochastic nature of DNA methylation (DNAm) changes in ageing cells as well as in normal cells at risk of neoplastic transformation, discussing the implications of this stochasticity for developing cancer risk prediction strategies, and in particular, how it may require a conceptual paradigm shift in how we select cancer risk markers. I also describe the mounting evidence that a significant proportion of DNAm changes in ageing and cancer development are related to cell proliferation, reflecting tissue-turnover and the opportunity this offers for predicting cancer risk via the development of epigenetic mitotic-like clocks. Finally, I describe how age-associated DNAm changes may be causally implicated in cancer development via an irreversible suppression of tissue-specific transcription factors that increases epigenetic and transcriptomic entropy, promoting a more plastic yet aberrant cancer stem-cell state. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Association analysis between an epigenetic risk score and blood pressure

Background

Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases.

Results

We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN.

Conclusions

Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Association analysis between an epigenetic alcohol risk score and blood pressure

Background

Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases.

Results

We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN.

Conclusions

Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

APOBEC shapes tumor evolution and age at onset of lung cancer in smokers

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

A prospective study of smoking-related white blood cell DNA methylation markers and risk of bladder cancer

Bladder cancer, a common neoplasm, is primarily caused by tobacco smoking. Epigenetic alterations including DNA methylation have the potential to be used as prospective markers of increased risk, particularly in at-risk populations such as smokers. We aimed to investigate the potential of smoking-related white blood cell (WBC) methylation markers to contribute to an increase in bladder cancer risk prediction over classical questionnaire-based smoking metrics (i.e., duration, intensity, packyears) in a nested case-control study within the prospective prostate, lung, colorectal, and ovarian (PLCO) Cancer Screening Trial and the alpha-tocopherol, beta-carotene cancer (ATBC) Prevention Study (789 cases; 849 controls). We identified 200 differentially methylated sites associated with smoking status and 28 significantly associated (after correction for multiple testing) with bladder cancer risk among 2670 previously reported smoking-related cytosine-phosphate-guanines sites (CpGs). Similar patterns were observed across cohorts. Receiver operating characteristic (ROC) analyses indicated that cg05575921 (AHHR), the strongest smoking-related association we identified for bladder cancer risk, alone yielded similar predictive performance (AUC: 0.60) than classical smoking metrics (AUC: 0.59-0.62). Best prediction was achieved by including the first principal component (PC1) from the 200 smoking-related CpGs alongside smoking metrics (AUC: 0.63-0.65). Further, PC1 remained significantly associated with elevated bladder cancer risk after adjusting for smoking metrics. These findings suggest DNA methylation profiles reflect aspects of tobacco smoke exposure in addition to those captured by smoking duration, intensity and packyears, and/or individual susceptibility relevant to bladder cancer etiology, warranting further investigation.

Maternal smoking during pregnancy and cortical structure in children with attention-deficit/hyperactivity disorder

Maternal smoking during pregnancy (MSDP) is considered a risk factor for ADHD. While the mechanisms underlying this association are not well understood, MSDP may impact the developing brain in ways that lead to ADHD. Here, we investigated the effect of prenatal smoking exposure on cortical brain structures in children with ADHD using two methods of assessing prenatal exposure: maternal recall and epigenetic typing. Exposure groups were defined according to: (1) maternal recall (+MSDP: n = 24; -MSDP: n = 85) and (2) epigenetic markers (EM) (+EM: n = 14 -EM: n = 21). CIVET-1.1.12 and RMINC were used to acquire cortical brain measurements and perform statistical analyses, respectively. The vertex with highest significance was tested for association with Continuous Performance Test (CPT) dimensions. While no differences of brain structures were identified between +MSDP and -MSDP, +EM children (n = 10) had significantly smaller surface area in the right orbitofrontal cortex (ROFc), middle temporal cortex (RTc) and parahippocampal gyrus (RPHg) (15% FDR) compared to -EM children (n = 20). Cortical surface area in the RPHg significantly correlated with CPT commission errors T-scores. This study suggests that molecular markers may better define exposure to environmental risks, as compared to human recall.

Epigenome-wide association study of total nicotine equivalents in multiethnic current smokers from three prospective cohorts

The impact of tobacco exposure on health varies by race and ethnicity and is closely tied to internal nicotine dose, a marker of carcinogen uptake. DNA methylation is strongly responsive to smoking status and may mediate health effects, but study of associations with internal dose is limited. We performed a blood leukocyte epigenome-wide association study (EWAS) of urinary total nicotine equivalents (TNEs; a measure of nicotine uptake) and DNA methylation measured using the MethylationEPIC v1.0 BeadChip (EPIC) in six racial and ethnic groups across three cohort studies. In the Multiethnic Cohort Study (discovery, n = 1994), TNEs were associated with differential methylation at 408 CpG sites across >250 genomic regions (p < 9 × 10-8). The top significant sites were annotated to AHRR, F2RL3, RARA, GPR15, PRSS23, and 2q37.1, all of which had decreasing methylation with increasing TNEs. We identified 45 novel CpG sites, of which 42 were unique to the EPIC array and eight annotated to genes not previously linked with smoking-related DNA methylation. The most significant signal in a novel gene was cg03748458 in MIR383;SGCZ. Fifty-one of the 408 discovery sites were validated in the Singapore Chinese Health Study (n = 340) and the Southern Community Cohort Study (n = 394) (Bonferroni corrected p < 1.23 × 10-4). Significant heterogeneity by race and ethnicity was detected for CpG sites in MYO1G and CYTH1. Furthermore, TNEs significantly mediated the association between cigarettes per day and DNA methylation at 15 sites (average 22.5%-44.3% proportion mediated). Our multiethnic study highlights the transethnic and ethnic-specific methylation associations with internal nicotine dose, a strong predictor of smoking-related morbidities.

Epigenetic Reprogramming Drives Epithelial Disruption in Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) remains a major public health challenge that contributes greatly to mortality and morbidity worldwide. Although it has long been recognized that the epithelium is altered in COPD, there has been little focus on targeting it to modify the disease course. Therefore, mechanisms that disrupt epithelial cell function in patients with COPD are poorly understood. In this study, we sought to determine whether epigenetic reprogramming of the cell-cell adhesion molecule E-cadherin, encoded by the CDH1 gene, disrupts epithelial integrity. By reducing these epigenetic marks, we can restore epithelial integrity and rescue alveolar airspace destruction. We used differentiated normal and COPD-derived primary human airway epithelial cells, genetically manipulated mouse tracheal epithelial cells, and mouse and human precision-cut lung slices to assess the effects of epigenetic reprogramming. We show that the loss of CDH1 in COPD is due to increased DNA methylation site at the CDH1 enhancer D through the downregulation of the ten-eleven translocase methylcytosine dioxygenase (TET) enzyme TET1. Increased DNA methylation at the enhancer D region decreases the enrichment of RNA polymerase II binding. Remarkably, treatment of human precision-cut slices derived from patients with COPD with the DNA demethylation agent 5-aza-2'-deoxycytidine decreased cell damage and reduced air space enlargement in the diseased tissue. Here, we present a novel mechanism that targets epigenetic modifications to reverse the tissue remodeling in human COPD lungs and serves as a proof of concept for developing a disease-modifying target.

Association of methylation risk score with incident type 2 diabetes mellitus: A nested case-control study

AIMS

To investigate the association of methylation risk score (MRS) and its interactions with environmental factors with type 2 diabetes mellitus (T2DM) risk.

METHODS

We conducted a nested case-control study with 241 onset cases and 241 matched controls. Conditional logistic regression models were employed to identify risk CpG sites. Simple and weighted MRSs were constructed based on the methylation levels of ATP-binding cassette G1 gene, fat mass and obesity associated gene, potassium voltage-gated channel member 1 gene, and thioredoxin-interacting protein gene previously associated with T2DM to estimate the association of MRS with T2DM risk. Stratified analyses were used to investigate interactions between MRS and environmental factors.

RESULTS

A total of 10 CpG loci were identified from the aforementioned genes to calculate MRS. After controlling for potential confounding factors, taking tertile 1 as reference, the odds ratios (ORs) and 95% confidence intervals (CIs) for T2DM of tertile 3 was 2.39 (1.36-4.20) for simple MRS and 2.59 (1.45-4.63) for weighted MRS. With per SD score increment in MRS, the OR (95% CI) was 1.66 (1.29-2.14) and 1.60 (1.24-2.08) for simple and weighted MRSs, respectively. J-curved associations were observed between both simple and weighted MRSs and T2DM risks. Additionally, multiplication interactions for smoking and hypertension with simple MRS on the risk of T2DM were found, similarly for smoking and obesity with weighted MRS on the risk of T2DM (all Pinteraction  < .05).

CONCLUSION

Elevated simple and weighted MRSs were associated with increased risk of T2DM. Environmental risk factors may influence the association between MRS and T2DM.

Genome-wide DNA methylation analysis of cannabis use disorder in a veteran cohort enriched for posttraumatic stress disorder

Cannabis use has been increasing over the past decade, not only in the general US population, but particularly among military veterans. With this rise in use has come a concomitant increase in cannabis use disorder (CUD) among veterans. Here, we performed an epigenome-wide association study for lifetime CUD in an Iraq/Afghanistan era veteran cohort enriched for posttraumatic stress disorder (PTSD) comprising 2,310 total subjects (1,109 non-Hispanic black and 1,201 non-Hispanic white). We also investigated CUD interactions with current PTSD status and examined potential indirect effects of DNA methylation (DNAm) on the relationship between CUD and psychiatric diagnoses. Four CpGs were associated with lifetime CUD, even after controlling for the effects of current smoking (AHRR cg05575921, LINC00299 cg23079012, VWA7 cg22112841, and FAM70A cg08760398). Importantly, cg05575921, a CpG strongly linked to smoking, remained associated with lifetime CUD even when restricting the analysis to veterans who reported never smoking cigarettes. Moreover, CUD interacted with current PTSD to affect cg05575921 and cg23079012 such that those with both CUD and PTSD displayed significantly lower DNAm compared to the other groups. Finally, we provide preliminary evidence that AHRR cg05575921 helps explain the association between CUD and any psychiatric diagnoses, specifically mood disorders.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

Pan-cancer proteogenomics characterization of tumor immunity

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Former smoking associated with epigenetic modifications in human granulosa cells among women undergoing assisted reproduction

Smoking exposure during adulthood can disrupt oocyte development in women, contributing to infertility and possibly adverse birth outcomes. Some of these effects may be reflected in epigenome profiles in granulosa cells (GCs) in human follicular fluid. We compared the epigenetic modifications throughout the genome in GCs from women who were former (N = 15) versus never smokers (N = 44) undergoing assisted reproductive technologies (ART). This study included 59 women undergoing ART. Smoking history including time since quitting was determined by questionnaire. GCs were collected during oocyte retrieval and DNA methylation (DNAm) levels were profiled using the Infinium MethylationEPIC BeadChip. We performed an epigenome-wide association study with robust linear models, regressing DNAm level at individual loci on smoking status, adjusting for age, ovarian stimulation protocol, and three surrogate variables. We performed differentially methylated regions (DMRs) analysis and over-representation analysis of the identified CpGs and corresponding gene set. 81 CpGs were differentially methylated among former smokers compared to never smokers (FDR < 0.05). We identified 2 significant DMRs (KCNQ1 and RHBDD2). The former smoking-associated genes were enriched in oxytocin signaling, adrenergic signaling in cardiomyocytes, platelet activation, axon guidance, and chemokine signaling pathway. These epigenetic variations have been associated with inflammatory responses, reproductive outcomes, cancer development, neurodevelopmental disorder, and cardiometabolic health. Secondarily, we examined the relationships between time since quitting and DNAm at significant CpGs. We observed three CpGs in negative associations with the length of quitting smoking (p < 0.05), which were cg04254052 (KCNIP1), cg22875371 (OGDHL), and cg27289628 (LOC148145), while one in positive association, which was cg13487862 (PLXNB1). As a pilot study, we demonstrated epigenetic modifications associated with former smoking in GCs. The study is informative to potential biological pathways underlying the documented association between smoking and female infertility and biomarker discovery for smoking-associated reproductive outcomes.

Developmental origins of psycho-cardiometabolic multimorbidity in adolescence and their underlying pathways through methylation markers: a two-cohort study

Understanding the biological mechanisms behind multimorbidity patterns in adolescence is important as they may act as intermediary risk factor for long-term health. We aimed to explore relationship between prenatal exposures and adolescent's psycho-cardiometabolic intermediary traits mediated through epigenetic biomarkers, using structural equation modeling (SEM). We used data from mother-child dyads from pregnancy and adolescents at 16-17 years from two prospective cohorts: Northern Finland Birth Cohort 1986 (NFBC1986) and Raine Study from Australia. Factor analysis was applied to generate two different latent factor structures: (a) prenatal exposures and (b) adolescence psycho-cardiometabolic intermediary traits. Furthermore, three types of epigenetic biomarkers were included: (1) DNA methylation score for maternal smoking during pregnancy (DNAmMSS), (2) DNAm age estimate PhenoAge and (3) DNAm estimate for telomere length (DNAmTL). Similar factor structure was observed between both cohorts yielding three prenatal factors, namely BMI (Body Mass Index), SOP (Socio-Obstetric-Profile), and Lifestyle, and four adolescent factors: Anthropometric, Insulin-Triglycerides, Blood Pressure, and Mental health. In the SEM pathways, stronger direct effects of F1prenatal-BMI (NFBC1986 = β: 0.27; Raine = β: 0.39) and F2prenatal-SOP (β: -0.11) factors were observed on adolescent psycho-cardiometabolic multimorbidity. We observed an indirect effect of prenatal latent factors through epigenetic markers on a psycho-cardiometabolic multimorbidity factor in Raine study (P < 0.05). The present study exemplifies an evidence-based approach in two different birth cohorts to demonstrate similar composite structure of prenatal exposures and psycho-cardiometabolic traits (despite cultural, social, and genetic differences) and a common plausible pathway between them through underlying epigenetic markers.

Epigenome-wide association study on short-, intermediate- and long-term ozone exposure in Han Chinese, the NSPT study

Epidemiological and epigenetic studies have acknowledged ambient ozone exposure associated with inflammatory and cardiovascular disease. However, the molecular mechanisms still remained unclear, and epigenome-wide analysis in cohort were lacking, especially in Chinese. We included blood-derived DNA methylation for 3365 Chinese participants from the NSPT cohort and estimated individual ozone exposure level of short-, intermediate- and long-term, based on a validated prediction model. We performed epigenome-wide association studies which identified 59 CpGs and 30 DMRs at a strict genome-wide significance (P < 5 ×10-8). We also conducted comparison on the DNA methylation alteration corresponding to different time windows, and observed an enhanced differentiated methylation trend for intermediate- and long-term exposure, while the short-term exposure associated methylation changes did not retain. The targeted genes of methylation alteration were involved in mechanism related to aging, inflammation disease, metabolic syndrome, neurodevelopmental disorders, and oncogenesis. Underlying pathways were enriched in biological activities including telomere maintenance process, DNA damage response and megakaryocyte differentiation. In conclusion, our study is the first EWAS on ozone exposure conducted in large-scale Han Chinese cohort and identified associated DNA methylation change on CpGs and regions, as well as related gene functions and pathways.

Decoding the DNA methylome of central nervous system tumors: An emerging modality for integrated diagnosis

The definitive diagnosis and classification of individual cancers are crucial for patient care and cancer research. To achieve a robust diagnosis of central nervous system (CNS) tumors, a genotype-phenotype integrated diagnostic approach was introduced in recent versions of the World Health Organization classification, followed by the incorporation of a genome-wide DNA methylome-based classification. Microarray-based platforms are widely used to obtain DNA methylome data, and the German Cancer Research Center (Deutsches Krebsforschungszentrum [DKFZ]) has a webtool for a DNA methylation-based classifier (DKFZ classifier). Integration of DNA methylome will further enhance the precision of CNS tumor classification, especially in diagnostically challenging cases. However, in the clinical application of DNA methylome-based classification, challenges related to data interpretation persist, in addition to technical caveats, regulations, and limited accessibility. Dimensionality reduction (DMR) can complement integrated diagnosis by visualizing a profile and comparing it with other known samples. Therefore, DNA methylome-based classification is a highly useful research tool for auxiliary analysis in challenging diagnostic and rare disease cases, and for establishing novel tumor concepts. Decoding the DNA methylome, especially by DMR in addition to DKFZ classifier, emphasizes the capability of grasping the fundamental biological principles that provide new perspectives on CNS tumors.

Tobacco Cessation Following Laryngeal Cancer Diagnosis Predicts Response to Treatment and Laryngectomy-Free Survival

OBJECTIVE

To determine the effect of tobacco cessation following laryngeal cancer diagnosis on response to first-line therapy, laryngectomy-free survival, and overall survival in patients who were current smokers at the time of diagnosis.

STUDY DESIGN

Retrospective, case-control study.

SETTING

OU Stephenson Cancer Center, National Cancer Institute-Designated Cancer Center.

METHODS

We included 140 patients diagnosed with laryngeal squamous cell carcinoma, who were current smokers at the time of diagnosis, and were treated with first-line definitive radiation or chemo/radiation with the intent to cure. The association between patient characteristics and treatment response was assessed using the χ2 test and logistic regression analysis. Survival outcomes were analyzed using Kaplan-Meier methods and Cox proportional-hazards models.

RESULTS

Of the 140 current smokers, 61 patients (45%) quit smoking prior to treatment initiation. In adjusted logistic regression analysis, quitters had 3.7 times higher odds of achieving a complete response to first-line therapy than active smokers (odds ratio: 3.694 [1.575-8.661]; P = .003). In the adjusted Cox proportional-hazards model, quitters were 54% less likely to require salvage laryngectomy within 7 years of diagnosis than active smokers (hazard ratio: 0.456 [0.246-0.848]; P = .013). Quitters had a statistically significant increase in 7-year overall survival compared to active smokers (P = .02).

CONCLUSION

This is the first study to show that in newly diagnosed laryngeal cancer patients who are current smokers at the time of diagnosis, tobacco cessation significantly increases therapy response, laryngectomy-free survival, and overall survival. These data stress the importance of systematically incorporating tobacco cessation programs into laryngeal cancer treatment plans.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation

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.

Prenatal cocaine exposure and its influence on pediatric epigenetic clocks and epigenetic scores in humans

The investigation of the effects of prenatal cocaine exposure (PCE) on offspring has been inconsistent, with few studies investigating biological outcomes in humans. We profiled genome-wide DNA methylation (DNAm) of umbilical cord blood (UCB) from newborns with (n = 35) and without (n = 47) PCE. We used DNAm data to (1) assess pediatric epigenetic clocks at birth and (2) to estimate epigenetic scores (ES) for lifetime disorders. We generated gestational epigenetic age estimates (DNAmGA) based on Knight and Bohlin epigenetic clocks. We also investigated the association between DNAmGA and UCB serum brain-derived neurotrophic factor (BDNF) levels. Considering the large-scale DNAm data availability and existing evidence regarding PCE as a risk for health problems later in life, we generated ES for tobacco smoking, psychosis, autism, diabetes, and obesity. A gene ontology (GO) analysis on the CpGs included in the ES with group differences was performed. PCE was associated with lower DNAmGA in newborns, and this effect remained significant when controlling for potential confounders, such as blood cell type composition predicted by DNAm and obstetric data. DNAmGA was negatively correlated with BDNF levels in the serum of UCB. Higher tobacco smoking, psychosis, and diabetes ES were found in the PCE group. The GO analysis revealed GABAergic synapses as a potential pathway altered by PCE. Our findings of decelerated DNAmGA and ES for adverse phenotypes associated with PCE, suggest that the effects of gestational cocaine exposure on the epigenetic landscape of human newborns are detectable at birth.

Altered DNA methylation within DNMT3A, AHRR, LTA/TNF loci mediates the effect of smoking on inflammatory bowel disease

This work aims to investigate how smoking exerts effect on the development of inflammatory bowel disease (IBD). A prospective cohort study and a Mendelian randomization study are first conducted to evaluate the association between smoking behaviors, smoking-related DNA methylation and the risks of Crohn's disease (CD) and ulcerative colitis (UC). We then perform both genome-wide methylation analysis and co-localization analysis to validate the observed associations. Compared to never smoking, current and previous smoking habits are associated with increased CD (P = 7.09 × 10-10) and UC (P < 2 × 10-16) risk, respectively. DNA methylation alteration at cg17742416 [DNMT3A] is linked to both CD (P = 7.30 × 10-8) and UC (P = 1.04 × 10-4) risk, while cg03599224 [LTA/TNF] is associated with CD risk (P = 1.91 × 10-6), and cg14647125 [AHRR] and cg23916896 [AHRR] are linked to UC risk (P = 0.001 and 0.002, respectively). Our study identifies biological mechanisms and pathways involved in the effects of smoking on the pathogenesis of IBD.