DNA methylation changes of whole blood cells in response to active smoking exposure in adults: a systematic review of DNA methylation studies

Normal bronchial field basal cells show persistent methylome-wide impact of tobacco smoking, including in known cancer genes

Lung carcinogenesis is causally linked to cigarette smoking, in part by epigenetic changes. We tested whether accumulated epigenetic change in smokers is apparent in bronchial basal cells as cells of origin of squamous cell carcinoma. Using an EM-seq platform covering 53.8 million CpGs (96% of the entire genome) at an average of 7.5 sequencing reads per CpG site at a single base resolution, we evaluated cytology-normal basal cells bronchoscopically brushed from the in situ tobacco smoke-exposed 'bronchial epithelial field' and isolated by short-term primary culture from 54 human subjects. We found that mean methylation was globally lower in ever (former and current) smokers versus never smokers (p = 0.0013) across promoters, CpG shores, exons, introns, 3'-UTRs, and intergenic regions, but not in CpG islands. Among 6mers with dinucleotides flanking CpG, those containing CGCG showed no effect from smoking, while those flanked with TT and AA displayed the strongest effects. At the gene level, smoking-related differences in methylation level were observed in CDKL1, ARTN, EDC3, CYP1B1, FAM131A, and MAGI2. Among candidate cancer genes, smoking reduced the methylation level in KRAS, ROS1, CDKN1A, CHRNB4, and CADM1. We conclude that smoking reduces long-term epigenome-wide methylation in bronchial stem cells, is impacted by the flanking sequence, and persists indefinitely beyond smoking cessation.

DNA methylation-based semen age prediction using the markers identified in Koreans and Europeans

In the forensic field, sexual assaults have consistently been the important issue, with semen frequently serving as the primary evidence. When the suspect is unidentified, estimating the perpetrator's age using investigating semen can provide important information. The VISAGE consortium conducted research on the semen age prediction focused on European semen samples, but the age prediction model has remained undisclosed. Additionally, several studies have reported methylation differences across populations, indicating that the European semen age prediction model might not be broadly applicable to other groups. A study did explore semen age prediction in Koreans using Illumina's Infinium Methylation450K BeadChip array, however recent developments in technology could enhance this approach. To address this, we conducted a study on Korean males aged 18-70 years. We initially analyzed 49 samples utilizing Illumina's Infinium MethylationEPIC BeadChip array to identify age-related CpG sites. From this analysis, we identified 9 age-related CpG markers, excluding one due to difficulties in locus-specific analysis. As a result, we used 11 markers including 8 newly identified CpGs from the EPIC array and 3 CpG markers from previous research utilizing the SNaPshot assay. Furthermore, we incorporated 13 CpG markers from the European study to analyze a total of 159 semen samples using the Illumina Nextera MPS system. This approach enabled us to test age-related markers identified in Europeans within the Korean population and to construct a more accurate age prediction model using markers from both Korean and European sources.

Smoking-independent DNA methylation markers for lung cancer risk: External validation in a large population-based cohort study

Smoking-associated epigenetic changes have been linked to lung cancer (LC) risk; however, the role of epigenetic alterations independent of smoking is yet to be fully understood. This study aimed to validate 16 previously reported CpG sites that are independent of smoking yet associated with LC risk within a population-based prospective cohort. Using the Infinium Methylation EPIC BeadChip kit or the Infinium HumanMethylation450K BeadChip Assay, DNA methylation (DNAm) in whole blood was assessed in four subsets (n = 736, 1027, 997, and 312) of a population-based cohort from Germany. The DNAm levels of the 16 smoking-independent CpG sites were analyzed. Hazard ratios (HRs) and their 95% confidence intervals (95% CIs) were calculated to assess associations of DNAm at the 16 CpG sites with LC risk, adjusting for multiple covariates, including smoking habits and a smoking-associated DNAm score. Over 17 years of follow-up, a total of 199 LCs were observed. Among the 16 CpGs, cg02211449 showed a negative association with LC risk (HR [95% CI] per SD increase, = 0.70 [0.63-0.78]), while cg11385536 (1.04 [1.01-1.07]), cg09736286 (1.64 [1.10-2.44]), cg19907023 (1.64 [1.01-2.66]), and cg22032485 (1.52 [1.04-2.21]) displayed positive associations with LC risk. Five of the 16 suggested smoking-independent CpGs could be externally validated as predictors of LC risk. Further research should address their potential contribution to enhanced LC risk stratification.

Cigarette smoke and biological age induce degenerative heterogeneity in retinal pigment epithelium

Environmental exposure such as cigarette smoke induces epigenetic changes that can induce degenerative heterogeneity and accelerate aging. In early age-related macular degeneration (AMD), the leading worldwide cause of blindness among the elderly, retinal pigment epithelial (RPE) cell heterogeneity is a key change. Since smoking is the strongest environmental risk factor for AMD, we hypothesized that cigarette smoke induces degenerative RPE heterogeneity through epigenetic changes that are distinct from aging, and that with aging, the RPE becomes vulnerable to cigarette smoke insult. We administered cigarette smoke condensate (CSC) intravitreally to young and aged mice and performed snRNA-seq and snATAC-seq on the RPE/choroid. This analysis identified separate cell clusters corresponding to healthy and abnormal, dedifferentiated RPE in both aged vehicle-treated and young CSC-treated mice. The dedifferentiated RPE were characterized by a global decrease in chromatin accessibility and decreased expression of genes in functional categories that were linked to hallmarks of aging. Notably, young, dedifferentiated RPE also exhibited a compensatory upregulation of hallmarks of aging-related genes, specifically those related to mitochondrial function and proteostasis. In contrast, aged dedifferentiated RPE did not express these compensatory changes, and did not survive CSC treatment, as experimentally verified with TUNEL labeling. These changes are relevant to early AMD because we identified through scRNA-seq, similar dedifferentiated and healthy macular RPE clusters in a donor who smoked and another with early AMD, but not from a nonsmoker. Degenerative cellular heterogeneity can include an abnormal cluster that jeopardizes cell survival and may represent an additional hallmark of ocular aging.

Neighborhood environment associations with cognitive function and structural brain measures in older African Americans

BACKGROUND

Since older adults spend significant time in their neighborhood environment, environmental factors such as neighborhood socioeconomic disadvantage, high racial segregation, low healthy food availability, low access to recreation, and minimal social engagement may have adverse effects on cognitive function and increase susceptibility to dementia. DNA methylation, which is associated with neighborhood characteristics as well as cognitive function and white matter hyperintensity (WMH), may act as a mediator between neighborhood characteristics and neurocognitive outcomes.

METHODS

In this study, we examined whether DNA methylation in peripheral blood leukocytes mediates the relationship between neighborhood characteristics and cognitive function (N = 542) or WMH (N = 466) in older African American (AA) participants without preliminary evidence of dementia from the Genetic Epidemiology Network of Arteriopathy (GENOA).

RESULTS

For a 1-mile buffer around a participant's residence, each additional fast food destination or unfavorable food store with alcohol per square mile was nominally associated with a 0.05 (95%CI: 0.01, 0.09) and a 0.04 (0.00, 0.08) second improvement in visual conceptual tracking score, respectively. Also, each additional alcohol drinking place per square mile was nominally associated with a 0.62 (0.05, 1.19) word increase in delayed recall score, indicating better memory function (all p < 0.05). Neighborhood characteristics were not associated with WMH. We did not find evidence that DNA methylation mediates the observed associations between neighborhood characteristics and cognitive function.

CONCLUSIONS

The presence of fast food destinations and unfavorable food stores with alcohol was associated cognitive measures, possibly due to greater social interaction provided in these venues. However, replication of these findings is necessary. Further examination of the potential pathways between the neighborhood environment and cognitive function/WMH may allow the development of potential behavioral, infrastructural, and pharmaceutical interventions to facilitate aging in place and healthy brain aging in older adults, especially in marginal populations that are most at risk.

Unraveling the causal impact of smoking and its DNA methylation signatures on cardiovascular disease: Mendelian randomization and colocalization analysis

BACKGROUND

To explore the mechanisms linking smoking to cardiovascular diseases (CVDs) from an epigenetic perspective.

METHODS

Mendelian Randomization (MR) analysis was performed to assess the causal effects of smoking behavior and DNA methylation levels at smoking-related CpG sites on nine CVDs, including aortic aneurysm, atrial fibrillation, coronary atherosclerosis, coronary heart disease, heart failure, intracerebral hemorrhage, ischemic stroke, myocardial infarction, subarachnoid hemorrhage. Colocalization analysis was used to further identify key smoking-related CpG sites from the MR causal estimates. Reactome enrichment analysis was used to elucidate the potential mechanisms.

RESULTS

MR analysis indicates that smoking behaviors are significantly associated with an increased risk of nine CVDs (OR > 1, P < 0.05). Through MR and colocalization analysis, five key smoking-related CpG sites were ultimately determined. DNA methylation alteration at cg25313468 (located in the TSS1500 region of REST) is simultaneously associated with the risk of atrial fibrillation, coronary atherosclerosis, coronary heart disease, and myocardial infarction. Additionally, cg21647257 (located in the TSS200 region of CLIP3) is associated with the risk of atrial fibrillation; cg06197751 (located in SGEF gene body) and cg07520810 (located in ARID5B gene body) are associated with the risk of coronary atherosclerosis; cg16822035 (located in MCF2L gene body) is associated with the risk of myocardial infarction. Enrichment analysis suggests that phosphatase and tensin homologue (PTEN) may be involved in the downstream mechanisms of cg25313468 (REST).

CONCLUSION

This study uncovers the relationship between smoking, DNA methylation, and CVDs, providing new insights into the pathogenic effect of smoking on CVDs from an epigenetic perspective.

Role of epigenetic mechanisms in the pathogenesis of chronic respiratory diseases and response to inhaled exposures: From basic concepts to clinical applications

Epigenetic modifications are chemical groups in our DNA (and chromatin) that determine which genes are active and which are shut off. Importantly, they integrate environmental signals to direct cellular function. Upon chronic environmental exposures, the epigenetic signature of lung cells gets altered, triggering aberrant gene expression programs that can lead to the development of chronic lung diseases. In addition to driving disease, epigenetic marks can serve as attractive lung disease biomarkers, due to early onset, disease specificity, and stability, warranting the need for more epigenetic research in the lung field. Despite substantial progress in mapping epigenetic alterations (mostly DNA methylation) in chronic lung diseases, the molecular mechanisms leading to their establishment are largely unknown. This review is meant as a guide for clinicians and lung researchers interested in epigenetic regulation with a focus on DNA methylation. It provides a short introduction to the main epigenetic mechanisms (DNA methylation, histone modifications and non-coding RNA) and the machinery responsible for their establishment and removal. It presents examples of epigenetic dysregulation across a spectrum of chronic lung diseases and discusses the current state of epigenetic therapies. Finally, it introduces the concept of epigenetic editing, an exciting novel approach to dissecting the functional role of epigenetic modifications. The promise of this emerging technology for the functional study of epigenetic mechanisms in cells and its potential future use in the clinic is further discussed.

Epigenome-Wide DNA Methylation in Unipolar Depression: Predictive Biomarker of Antidepressant Treatment Response?

BACKGROUND

Despite the well-documented efficacy of antidepressant agents for the treatment of major depressive disorder (MDD), initial treatment nonresponse rates are high. Recent years have seen an increase in research into predictive biomarkers toward improving diagnosis and individualized treatment. Among those, epigenetic mechanisms such as DNA methylation constitute promising candidate markers in predicting antidepressant treatment response in MDD. The present study sought to address epigenome-wide DNA methylation as a predictor of antidepressant treatment response in the largest sample to date of patients with MDD.

METHODS

Epigenome-wide DNA methylation was analyzed using the Infinium MethylationEPIC BeadChip in peripheral blood of n = 230 Caucasian patients with MDD receiving 6-week antidepressant treatment in a naturalistic in-patient setting as well as in a subsample of n = 107 patients primarily receiving continuous treatment with serotonin reuptake inhibitors or serotonin and norepinephrine reuptake inhibitors. Treatment response was assessed by means of the Hamilton Depression Scale.

RESULTS

No genome-wide significant hits were observed. Suggestive (P < 1E-5) epigenome-wide evidence was discerned for altered DNA methylation at 6 CpG sites (LOC102724467, LOC100506023, RSPO2, SAG, IL16, PRKCI) to predict response to naturalistic antidepressant treatment. In patients treated with serotonin reuptake inhibitors or serotonin and norepinephrine reuptake inhibitors, differential DNA methylation at 11 CpGs, for example, mapping to the TIMP2, VDAC1, or SORL1 genes, was suggestively associated with treatment response.

CONCLUSIONS

The present results provide preliminary evidence for altered DNA methylation patterns to be associated with antidepressant treatment response in MDD. Provided significant replication in independent and larger samples, the present findings might in the future aid in clinical decision-making toward more individualized and thus more efficacious treatments of MDD.

DNA methylation at AHRR as a master predictor of smoke exposure and a biomarker for sleep and exercise

BACKGROUND

DNA methylation profiling may provide a more accurate measure of the smoking status than self-report and may be useful in guiding clinical interventions and forensic investigations. In the current study, blood DNA methylation profiles of nearly 800 Polish individuals were assayed using Illuminia EPIC and the inference of smoking from epigenetic data was explored. In addition, we focused on the role of the AHRR gene as a top marker for smoking and investigated its responsiveness to other lifestyle behaviors.

RESULTS

We found > 450 significant CpGs associated with cigarette consumption, and overrepresented in various biological functions including cell communication, response to stress, blood vessel development, cell death, and atherosclerosis. The model consisting of cg05575921 in AHRR (p = 4.5 × 10-32) and three additional CpGs (cg09594361, cg21322436 in CNTNAP2 and cg09842685) was able to predict smoking status with a high accuracy of AUC = 0.8 in the test set. Importantly, a gradual increase in the probability of smoking was observed, starting from occasional smokers to regular heavy smokers. Furthermore, former smokers displayed the intermediate DNA methylation profiles compared to current and never smokers, and thus our results indicate the potential reversibility of DNA methylation after smoking cessation. The AHRR played a key role in a predictive analysis, explaining 21.5% of the variation in smoking. In addition, the AHRR methylation was analyzed for association with other modifiable lifestyle factors, and showed significance for sleep and physical activity. We also showed that the epigenetic score for smoking was significantly correlated with most of the epigenetic clocks tested, except for two first-generation clocks.

CONCLUSIONS

Our study suggests that a more rapid return to never-smoker methylation levels after smoking cessation may be achievable in people who change their lifestyle in terms of physical activity and sleep duration. As cigarette smoking has been implicated in the literature as a leading cause of epigenetic aging and AHRR appears to be modifiable by multiple exogenous factors, it emerges as a promising target for intervention and investment.

A systematic review of phenotypic and epigenetic clocks used for aging and mortality quantification in humans

Aging is the leading driver of disease in humans and has profound impacts on mortality. Biological clocks are used to measure the aging process in the hopes of identifying possible interventions. Biological clocks may be categorized as phenotypic or epigenetic, where phenotypic clocks use easily measurable clinical biomarkers and epigenetic clocks use cellular methylation data. In recent years, methylation clocks have attained phenomenal performance when predicting chronological age and have been linked to various age-related diseases. Additionally, phenotypic clocks have been proven to be able to predict mortality better than chronological age, providing intracellular insights into the aging process. This review aimed to systematically survey all proposed epigenetic and phenotypic clocks to date, excluding mitotic clocks (i.e., cancer risk clocks) and those that were modeled using non-human samples. We reported the predictive performance of 33 clocks and outlined the statistical or machine learning techniques used. We also reported the most influential clinical measurements used in the included phenotypic clocks. Our findings provide a systematic reporting of the last decade of biological clock research and indicate possible avenues for future research.

Artificial intelligence-based epigenomic, transcriptomic and histologic signatures of tobacco use in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) biomarker studies rarely employ multi-omic biomarker strategies and pertinent clinicopathologic characteristics to predict mortality. In this study we determine for the first time a combined epigenetic, gene expression, and histology signature that differentiates between patients with different tobacco use history (heavy tobacco use with ≥10 pack years vs. no tobacco use). Using The Cancer Genome Atlas (TCGA) cohort (n = 257) and an internal cohort (n = 40), we identify 3 epigenetic markers (GPR15, GNG12, GDNF) and 13 expression markers (IGHA2, SCG5, RPL3L, NTRK1, CD96, BMP6, TFPI2, EFEMP2, RYR3, DMTN, GPD2, BAALC, and FMO3), which are dysregulated in OSCC patients who were never smokers vs. those who have a ≥ 10 pack year history. While mortality risk prediction based on smoking status and clinicopathologic covariates alone is inaccurate (c-statistic = 0.57), the combined epigenetic/expression and histologic signature has a c-statistic = 0.9409 in predicting 5-year mortality in OSCC patients.

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.

The effect of maternal prenatal tobacco smoking on offspring academic achievement: A systematic review and meta-analysis

BACKGROUND

Previous epidemiological studies examining the prospective association between maternal prenatal tobacco smoking and offspring academic achievement have reported conflicting results. Therefore, this systematic review and meta-analysis was conducted to examine the magnitude and consistency of association reported by those studies.

METHODS

This systematic review and meta-analysis was guided by the PRISMA protocol. Relevant epidemiological studies on the topic were extracted from four main databases (PubMed/Medline, Embase, PsycINFO, and Scopus). The Newcastle-Ottawa Scale (NOS) was used to appraise the methodological quality of the included studies. We conducted a narrative assessment of the studies that did not report effect estimates. Inverse variance-weighted random effect meta-analysis was used to combine studies reporting effect sizes to estimate pooled adjusted odds ratio with 95% confidence intervals (95% CI). The review was prospectively registered in PROSPERO (CRD42022350901).

RESULTS

Nineteen observational studies, published between 1973 and 2021 with a total of 1.25 million study participants were included in the final review. Of these, fifteen studies (79 %) reported reduced academic achievement in offspring exposed to maternal prenatal tobacco smoking. The eight primary studies (sample size = 723,877) included in the meta-analysis together suggested a 49 % higher risk of reduced academic achievement in offspring exposed to maternal prenatal tobacco smoking when compared to non-exposed offspring (Pooled odds ratio = 1.49, 95 % CI:1.17-1.91).

CONCLUSION

Our review found a positive association between maternal prenatal tobacco smoking and offspring reduced academic achievement. However, variation in the adjustment of potential confounders and significant heterogeneity across included studies limited more conclusive inference. Mechanistic studies to identify causal pathways and specific academic impacts are needed to inform targeted developmental programs to assist child learning and academic performance.

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.

No association between peripheral serotonin-gene-related DNA methylation and brain serotonin neurotransmission in the healthy and depressed state

BACKGROUND

Methylation of serotonin-related genes has been proposed as a plausible gene-by-environment link which may mediate environmental stress, depressive and anxiety symptoms. DNA methylation is often measured in blood cells, but little is known about the association between this peripheral epigenetic modification and brain serotonergic architecture. Here, we evaluated the association between whole-blood-derived methylation of four CpG sites in the serotonin transporter (SLC6A4) and six CpG sites of the tryptophan hydroxylase 2 (TPH2) gene and in-vivo brain levels of serotonin transporter (5-HTT) and serotonin 4 receptor (5-HT4) in a cohort of healthy individuals (N = 254) and, for 5-HT4, in a cohort of unmedicated patients with depression (N = 90). To do so, we quantified SLC6A4/TPH2 methylation using bisulfite pyrosequencing and estimated brain 5-HT4 and 5-HTT levels using positron emission tomography. In addition, we explored the association between SLC6A4 and TPH2 methylation and measures of early life and recent stress, depressive and anxiety symptoms on 297 healthy individuals.

RESULTS

We found no statistically significant association between peripheral DNA methylation and brain markers of serotonergic neurotransmission in patients with depression or in healthy individuals. In addition, although SLC6A4 CpG2 (chr17:30,236,083) methylation was marginally associated with the parental bonding inventory overprotection score in the healthy cohort, statistical significance did not remain after accounting for blood cell heterogeneity.

CONCLUSIONS

We suggest that findings on peripheral DNA methylation in the context of brain serotonin-related features should be interpreted with caution. More studies are needed to rule out a role of SLC6A4 and TPH2 methylation as biomarkers for environmental stress, depressive or anxiety symptoms.

Inference of tobacco and alcohol consumption habits from DNA methylation analysis of blood

DNA methylation has become a biomarker of great interest in the forensic and clinical fields. In criminal investigations, the study of this epigenetic marker has allowed the development of DNA intelligence tools providing information that can be useful for investigators, such as age prediction. Following a similar trend, when the origin of a sample in a criminal scenario is unknown, the inference of an individual's lifestyle such as tobacco use and alcohol consumption could provide relevant information to help in the identification of DNA donors at the crime scene. At the same time, in the clinical domain, prediction of these trends of consumption could allow the identification of people at risk or better identification of the causes of different pathologies. In the present study, DNA methylation data from the UK AIRWAVE study was used to build two binomial logistic models for the inference of smoking and drinking status. A total of 348 individuals (116 non-smokers, 116 former smokers and 116 smokers) plus a total of 237 individuals (79 non-drinkers, 79 moderate drinkers and 79 drinkers) were used for development of tobacco and alcohol consumption prediction models, respectively. The tobacco prediction model was composed of two CpGs (cg05575921 in AHRR and cg01940273) and the alcohol prediction model three CpGs (cg06690548 in SLC7A11, cg0886875 and cg21294714 in MIR4435-2HG), providing correct classifications of 86.49% and 74.26%, respectively. Validation of the models was performed using leave-one-out cross-validation. Additionally, two independent testing sets were also assessed for tobacco and alcohol consumption. Considering that the consumption of these substances could underlie accelerated epigenetic ageing patterns, the effect of these lifestyles on the prediction of age was evaluated. To do that, a quantile regression model based on previous studies was generated, and the potential effect of tobacco and alcohol consumption with the epigenetic age was assessed. The Wilcoxon test was used to evaluate the residuals generated by the model and no significant differences were observed between the categories analyzed.

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

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.

An Emerging Approach of Age-Related Hearing Loss Research: Application of Integrated Multi-Omics Analysis

As one of the most common otologic diseases in the elderly, age-related hearing loss (ARHL) usually characterized by hearing loss and cognitive disorders, which have a significant impact on the elderly's physical and mental health and quality of life. However, as a typical disease of aging, it is unclear why aging causes widespread hearing impairment in the elderly. As molecular biological experiments have been conducted for research recently, ARHL is gradually established at various levels with the application and development of integrated multi-omics analysis in the studies of ARHL. Here, the recent progress in the application of multi-omics analysis in the molecular mechanisms of ARHL development and therapeutic regimens, including the combined analysis of different omics, such as transcriptome, proteome, and metabolome, to screen for risk sites, risk genes, and differences in lipid metabolism, etc., is outlined and the integrated histological data further promote the profound understanding of the disease process as well as physiological mechanisms of ARHL. The advantages and disadvantages of multi-omics analysis in disease research are also discussed and the authors speculate on the future prospects and applications of this part-to-whole approach, which may provide more comprehensive guidance for ARHL and aging disease prevention and treatment.

A validated restriction enzyme ddPCR cg05575921 (AHRR) assay to accurately assess smoking exposure

BACKGROUND & METHODS

In this study, a novel restriction enzyme (RE) digestion-based droplet digital polymerase chain reaction (ddPCR) assay was designed for cg005575921 within the AHRR gene body and compared with matching results obtained by bisulfite conversion (BIS) ddPCR and Illumina DNA methylation array.

RESULTS

The RE ddPCR cg05575921 assay appeared concordant with BIS ddPCR (r2 = 0.94, P < 0.0001) and, when compared with the Illumina array, had significantly better smoking status classification performance for current versus never smoked (AUC 0.96 versus 0.93, P < 0.04) and current versus ex-smoker (AUC 0.88 versus 0.83, P < 0.04) comparisons.

CONCLUSIONS

The RE ddPCR cg05575921 assay accurately predicts smoking status and could be a useful component of 'precision-medicine' chronic disease risk screening tools.

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.

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.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Accelerated epigenetic age, inflammation, and gene expression in lung: comparisons of smokers and vapers with non-smokers

BACKGROUND

Cigarette smoking and aging are the main risk factors for pulmonary diseases, including cancer. Epigenetic aging may explain the relationship between smoking, electronic cigarette vaping, and pulmonary health. No study has examined smoking and vaping-related epigenetic aging in relation to lung biomarkers.

METHODS

Lung epigenetic aging measured by DNA methylation (mAge) and its acceleration (mAA) was assessed in young (age 21-30) electronic cigarette vapers (EC, n = 14, including 3 never-smoking EC), smokers (SM, n = 16), and non-EC/non-SM (NS, n = 39). We investigated relationships of mAge estimates with chronological age (Horvath-mAge), lifespan/mortality (Grim-mAge), telomere length (TL-mAge), smoking/EC history, urinary biomarkers, lung cytokines, and transcriptome.

RESULTS

Compared to NS, EC and SM had significantly older Grim-mAge, shorter TL-mAge, significantly accelerated Grim-mAge and decelerated TL-mAge. Among SM, Grim-mAA was associated with nicotine intake and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL). For EC, Horvath-mAA was significantly correlated with puffs per day. Overall, cytokines (IL-1β, IL-6, and IL-8) and 759 transcripts (651 unique genes) were significantly associated with Grim-mAA. Grim-mAA-associated genes were highly enriched in immune-related pathways and genes that play a role in the morphology and structures of cells/tissues.

CONCLUSIONS

Faster lung mAge for SM is consistent with prior studies of blood. Faster lung mAge for EC compared to NS indicates possible adverse pulmonary effects of EC on biological aging. Our findings support further research, particularly on epigenetic markers, on effects of smoking and vaping on pulmonary health. Given that most EC are former smokers, further study is needed to understand unique effects of electronic cigarettes on biological aging.

Dissecting the pathogenic effects of smoking and its hallmarks in blood DNA methylation on colorectal cancer risk

BACKGROUND

Tobacco smoking is suggested as a risk factor for colorectal cancer (CRC), but the complex relationship and the potential pathway are not fully understood.

METHODS

We performed two-sample Mendelian randomisation (MR) analyses with genetic instruments for smoking behaviours and related DNA methylation in blood and summary-level GWAS data of colorectal cancer to disentangle the relationship. Colocalization analyses and prospective gene-environment interaction analyses were also conducted as replication.

RESULTS

Convincing evidence was identified for the pathogenic effect of smoking initiation on CRC risk and suggestive evidence was observed for the protective effect of smoking cessation in the univariable MR analyses. Multivariable MR analysis revealed that these associations were independent of other smoking phenotypes and alcohol drinking. Genetically predicted methylation at CpG site cg17823346 [ZMIZ1] were identified to decrease CRC risk; while genetically predicted methylation at cg02149899 would increase CRC risk. Colocalization and gene-environment interaction analyses added further evidence to the relationship between epigenetic modification at cg17823346 [ZMIZ1] as well as cg02149899 and CRC risk.

DISCUSSION

Our study confirms the significant association between tobacco smoking, DNA methylation and CRC risk and yields a novel insight into the pathogenic effect of tobacco smoking on CRC risk.

Effects of smoking on genome-wide DNA methylation profiles: A study of discordant and concordant monozygotic twin pairs

Background

Smoking-associated DNA methylation levels identified through epigenome-wide association studies (EWASs) are generally ascribed to smoking-reactive mechanisms, but the contribution of a shared genetic predisposition to smoking and DNA methylation levels is typically not accounted for.

Methods

We exploited a strong within-family design, that is, the discordant monozygotic twin design, to study reactiveness of DNA methylation in blood cells to smoking and reversibility of methylation patterns upon quitting smoking. Illumina HumanMethylation450 BeadChip data were available for 769 monozygotic twin pairs (mean age = 36 years, range = 18-78, 70% female), including pairs discordant or concordant for current or former smoking.

Results

In pairs discordant for current smoking, 13 differentially methylated CpGs were found between current smoking twins and their genetically identical co-twin who never smoked. Top sites include multiple CpGs in CACNA1D and GNG12, which encode subunits of a calcium voltage-gated channel and G protein, respectively. These proteins interact with the nicotinic acetylcholine receptor, suggesting that methylation levels at these CpGs might be reactive to nicotine exposure. All 13 CpGs have been previously associated with smoking in unrelated individuals and data from monozygotic pairs discordant for former smoking indicated that methylation patterns are to a large extent reversible upon smoking cessation. We further showed that differences in smoking level exposure for monozygotic twins who are both current smokers but differ in the number of cigarettes they smoke are reflected in their DNA methylation profiles.

Conclusions

In conclusion, by analysing data from monozygotic twins, we robustly demonstrate that DNA methylation level in human blood cells is reactive to cigarette smoking.

Funding

We acknowledge funding from the National Institute on Drug Abuse grant DA049867, the Netherlands Organization for Scientific Research (NWO): Biobanking and Biomolecular Research Infrastructure (BBMRI-NL, NWO 184.033.111) and the BBRMI-NL-financed BIOS Consortium (NWO 184.021.007), NWO Large Scale infrastructures X-Omics (184.034.019), Genotype/phenotype database for behaviour genetic and genetic epidemiological studies (ZonMw Middelgroot 911-09-032); Netherlands Twin Registry Repository: researching the interplay between genome and environment (NWO-Groot 480-15-001/674); the Avera Institute, Sioux Falls (USA), and the National Institutes of Health (NIH R01 HD042157-01A1, MH081802, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995); epigenetic data were generated at the Human Genomics Facility (HuGe-F) at ErasmusMC Rotterdam. Cotinine assaying was sponsored by the Neuroscience Campus Amsterdam. DIB acknowledges the Royal Netherlands Academy of Science Professor Award (PAH/6635).

A novel quantitative targeted analysis of X-chromosome inactivation (XCI) using nanopore sequencing

X-chromosome inactivation (XCI) analyses often assist in diagnostics of X-linked traits, however accurate assessment remains challenging with current methods. We developed a novel strategy using amplification-free Cas9 enrichment and Oxford nanopore technologies sequencing called XCI-ONT, to investigate and rigorously quantify XCI in human androgen receptor gene (AR) and human X-linked retinitis pigmentosa 2 gene (RP2). XCI-ONT measures methylation over 116 CpGs in AR and 58 CpGs in RP2, and separate parental X-chromosomes without PCR bias. We show the usefulness of the XCI-ONT strategy over the PCR-based golden standard XCI technique that only investigates one or two CpGs per gene. The results highlight the limitations of using the golden standard technique when the XCI pattern is partially skewed and the advantages of XCI-ONT to rigorously quantify XCI. This study provides a universal XCI-method on DNA, which is highly valuable in clinical and research framework of X-linked traits.

Short-term smoking cessation leads to a universal decrease in whole blood genomic DNA methylation in patients with a smoking history

BACKGROUND

Epigenetics is involved in various human diseases. Smoking is one of the most common environmental factors causing epigenetic changes. The DNA methylation changes and mechanisms after quitting smoking have yet to be defined. The present study examined the changes in DNA methylation levels before and after short-term smoking cessation and explored the potential mechanism.

METHODS

Whole blood and clinical data were collected from 8 patients before and after short-term smoking cessation, DNA methylation was assessed, and differentially methylated sites were analyzed, followed by a comprehensive analysis of the differentially methylated sites with clinical data. GO/KEGG enrichment and protein-protein interaction (PPI) network analyses identified the hub genes. The differentially methylated sites between former and current smokers in GSE50660 from the GEO database were detected by GEO2R. Then, a Venn analysis was carried out using the differentially methylated sites. GO/KEGG enrichment analysis was performed on the genes corresponding to the common DNA methylation sites, the PPI network was constructed, and hub genes were predicted. The enriched genes associated with the cell cycle were selected, and the pan-cancer gene expression and clinical significance in lung cancer were analyzed based on the TCGA database.

RESULTS

Most genes showed decreased DNA methylation levels after short-term smoking cessation; 694 upregulated methylation CpG sites and 3184 downregulated methylation CpG sites were identified. The DNA methylation levels were altered according to the clinical data (body weight, expiratory, and tobacco dependence score). Enrichment analysis, construction of the PPI network, and pan-cancer analysis suggested that smoking cessation may affect various biological processes.

CONCLUSIONS

Smoking cessation leads to epigenetic changes, mainly decreased in the decline of most DNA methylation levels. Bioinformatics further identified the biologically relevant changes after short-term smoking cessation.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Genome-Wide Methylation Profiling in 229 Patients With Crohn's Disease Requiring Intestinal Resection: Epigenetic Analysis of the Trial of Prevention of Post-operative Crohn's Disease (TOPPIC)

BACKGROUND & AIMS

DNA methylation alterations may provide important insights into gene-environment interaction in cancer, aging, and complex diseases, such as inflammatory bowel disease (IBD). We aim first to determine whether the circulating DNA methylome in patients requiring surgery may predict Crohn's disease (CD) recurrence following intestinal resection; and second to compare the circulating methylome seen in patients with established CD with that we had reported in a series of inception cohorts.

METHODS

TOPPIC was a placebo-controlled, randomized controlled trial of 6-mercaptopurine at 29 UK centers in patients with CD undergoing ileocolic resection between 2008 and 2012. Genomic DNA was extracted from whole blood samples from 229 of the 240 patients taken before intestinal surgery and analyzed using 450KHumanMethylation and Infinium Omni Express Exome arrays (Illumina, San Diego, CA). Coprimary objectives were to determine whether methylation alterations may predict clinical disease recurrence; and to assess whether the epigenetic alterations previously reported in newly diagnosed IBD were present in the patients with CD recruited into the TOPPIC study. Differential methylation and variance analysis was performed comparing patients with and without clinical evidence of recurrence. Secondary analyses included investigation of methylation associations with smoking, genotype (MeQTLs), and chronologic age. Validation of our previously published case-control observation of the methylome was performed using historical control data (CD, n = 123; Control, n = 198).

RESULTS

CD recurrence in patients following surgery is associated with 5 differentially methylated positions (Holm P < .05), including probes mapping to WHSC1 (P = 4.1 × 10-9, Holm P = .002) and EFNA3 (P = 4.9 × 10-8, Holm P = .02). Five differentially variable positions are demonstrated in the group of patients with evidence of disease recurrence including a probe mapping to MAD1L1 (P = 6.4 × 10-5). DNA methylation clock analyses demonstrated significant age acceleration in CD compared with control subjects (GrimAge + 2 years; 95% confidence interval, 1.2-2.7 years), with some evidence for accelerated aging in patients with CD with disease recurrence following surgery (GrimAge +1.04 years; 95% confidence interval, -0.04 to 2.22). Significant methylation differences between CD cases and control subjects were seen by comparing this cohort in conjunction with previously published control data, including validation of our previously described differentially methylated positions (RPS6KA2 P = 1.2 × 10-19, SBNO2 = 1.2 × 10-11) and regions (TXK [false discovery rate, P = 3.6 × 10-14], WRAP73 [false discovery rate, P = 1.9 × 10-9], VMP1 [false discovery rate, P = 1.7 × 10-7], and ITGB2 [false discovery rate, P = 1.4 × 10-7]).

CONCLUSIONS

We demonstrate differential methylation and differentially variable methylation in patients developing clinical recurrence within 3 years of surgery. Moreover, we report replication of the CD-associated methylome, previously characterized only in adult and pediatric inception cohorts, in patients with medically refractory disease needing surgery.

Patients with schizophrenia and bipolar disorder display a similar global gene expression signature in whole blood that reflects elevated proportion of immature neutrophil cells with association to lipid changes

Schizophrenia (SCZ) and bipolar disorder (BD) share clinical characteristics, genetic susceptibility, and immune alterations. We aimed to identify differential transcriptional patterns in peripheral blood cells of patients with SCZ or BD versus healthy controls (HC). We analyzed microarray-based global gene expression data in whole blood from a cohort of SCZ (N = 329), BD (N = 203) and HC (N = 189). In total, 65 genes were significantly differentially expressed in SCZ and 125 in BD, as compared to HC, with similar ratio of up- and downregulated genes in both disorders. Among the top differentially expressed genes, we found an innate immunity signature that was shared between SCZ and BD, consisting of a cluster of upregulated genes (e.g., OLFM4, ELANE, BPI and MPO) that indicate an increased fraction of immature neutrophils. Several of these genes displayed sex differences in the expression pattern, and post-hoc analysis demonstrated a positive correlation with triglyceride and a negative correlation with HDL cholesterol. We found that many of the downregulated genes in SCZ and BD were associated with smoking. These findings of neutrophil granulocyte-associated transcriptome signatures in both SCZ and BD point at altered innate immunity pathways with association to lipid changes and potential for clinical translation.

ADHD genetic burden associates with older epigenetic age: mediating roles of education, behavioral and sociodemographic factors among older adults

BACKGROUND

Shortened lifespans are associated with having Attention Deficit Hyperactivity Disorder (ADHD), which is likely mediated by related behavioral and sociodemographic factors that are also associated with accelerated physiological aging. Such factors include exhibiting more depressive symptoms, more cigarette smoking, higher body mass index, lower educational attainment, lower income in adulthood, and more challenges with cognitive processes compared to the general population. A higher polygenic score for ADHD (ADHD-PGS) is associated with having more characteristic features of ADHD. The degree to which (1) the ADHD-PGS associates with an epigenetic biomarker developed to predict accelerated aging and earlier mortality is unknown, as are whether (2) an association would be mediated by behavioral and sociodemographic correlates of ADHD, or (3) an association would be mediated first by educational attainment, then by behavioral and sociodemographic correlates. We evaluated these relationships in a population-based sample from the US Health and Retirement Study, among N = 2311 adults age 50 and older, of European-ancestry, with blood-based epigenetic and genetic data. The ADHD-PGS was calculated from a prior genomewide meta-analysis. Epigenome-wide DNA methylation levels that index biological aging and earlier age of mortality were quantified by a blood-based biomarker called GrimAge. We used a structural equation modeling approach to test associations with single and multi-mediation effects of behavioral and contextual indicators on GrimAge, adjusted for covariates.

RESULTS

The ADHD-PGS was significantly and directly associated with GrimAge when adjusting for covariates. In single mediation models, the effect of the ADHD-PGS on GrimAge was partially mediated via smoking, depressive symptoms, and education. In multi-mediation models, the effect of the ADHD-PGS on GrimAge was mediated first through education, then smoking, depressive symptoms, BMI, and income.

CONCLUSIONS

Findings have implications for geroscience research in elucidating lifecourse pathways through which ADHD genetic burden and symptoms can alter risks for accelerated aging and shortened lifespans, when indexed by an epigenetic biomarker. More education appears to play a central role in attenuating negative effects on epigenetic aging from behavioral and sociodemographic risk factors related to ADHD. We discuss implications for the potential behavioral and sociodemographic mediators that may attenuate negative biological system effects.

Targeted DNA methylation analysis and prediction of smoking habits in blood based on massively parallel sequencing

Tobacco smoking is a frequent habit sustained by > 1.3 billion people in 2020 and the leading preventable factor for health risk and premature mortality worldwide. In the forensic context, predicting smoking habits from biological samples may allow broadening DNA phenotyping. In this study, we aimed to implement previously published smoking habit classification models based on blood DNA methylation at 13 CpGs. First, we developed a matching lab tool based on bisulfite conversion and multiplex PCR followed by amplification-free library preparation and targeted paired-end massively parallel sequencing (MPS). Analysis of six technical duplicates revealed high reproducibility of methylation measurements (Pearson correlation of 0.983). Artificially methylated standards uncovered marker-specific amplification bias, which we corrected via bi-exponential models. We then applied our MPS tool to 232 blood samples from Europeans of a wide age range, of which 90 were current, 71 former and 71 never smokers. On average, we obtained 189,000 reads/sample and 15,000 reads/CpG, without marker drop-out. Methylation distributions per smoking category roughly corresponded to previous microarray analysis, showcasing large inter-individual variation but with technology-driven bias. Methylation at 11 out of 13 smoking-CpGs correlated with daily cigarettes in current smokers, while solely one was weakly correlated with time since cessation in former smokers. Interestingly, eight smoking-CpGs correlated with age, and one displayed weak but significant sex-associated methylation differences. Using bias-uncorrected MPS data, smoking habits were relatively accurately predicted using both two- (current/non-current) and three- (never/former/current) category model, but bias correction resulted in worse prediction performance for both models. Finally, to account for technology-driven variation, we built new, joint models with inter-technology corrections, which resulted in improved prediction results for both models, with or without PCR bias correction (e.g. MPS cross-validation F₁-score > 0.8; 2-categories). Overall, our novel assay takes us one step closer towards the forensic application of viable smoking habit prediction from blood traces. However, future research is needed towards forensically validating the assay, especially in terms of sensitivity. We also need to further shed light on the employed biomarkers, particularly on the mechanistics, tissue specificity and putative confounders of smoking epigenetic signatures.

Analysis of Systemic Epigenetic Alterations in Inflammatory Bowel Disease: Defining Geographical, Genetic and Immune-Inflammatory influences on the Circulating Methylome

BACKGROUND

Epigenetic alterations may provide valuable insights into gene-environment interactions in the pathogenesis of inflammatory bowel disease [IBD].

METHODS

Genome-wide methylation was measured from peripheral blood using the Illumina 450k platform in a case-control study in an inception cohort (295 controls, 154 Crohn's disease [CD], 161 ulcerative colitis [UC], 28 IBD unclassified [IBD-U)] with covariates of age, sex and cell counts, deconvoluted by the Houseman method. Genotyping was performed using Illumina HumanOmniExpressExome-8 BeadChips and gene expression using the Ion AmpliSeq Human Gene Expression Core Panel. Treatment escalation was characterized by the need for biological agents or surgery after initial disease remission.

RESULTS

A total of 137 differentially methylated positions [DMPs] were identified in IBD, including VMP1/MIR21 [p = 9.11 × 10-15] and RPS6KA2 [6.43 × 10-13], with consistency seen across Scandinavia and the UK. Dysregulated loci demonstrate strong genetic influence, notably VMP1 [p = 1.53 × 10-15]. Age acceleration is seen in IBD [coefficient 0.94, p < 2.2 × 10-16]. Several immuno-active genes demonstrated highly significant correlations between methylation and gene expression in IBD, in particular OSM: IBD r = -0.32, p = 3.64 × 10-7 vs non-IBD r = -0.14, p = 0.77]. Multi-omic integration of the methylome, genome and transcriptome also implicated specific pathways that associate with immune activation, response and regulation at disease inception. At follow-up, a signature of three DMPs [TAP1, TESPA1, RPTOR] were associated with treatment escalation to biological agents or surgery (hazard ratio of 5.19 [CI: 2.14-12.56], logrank p = 9.70 × 10-4).

CONCLUSION

These data demonstrate consistent epigenetic alterations at diagnosis in European patients with IBD, providing insights into the pathogenetic importance and translational potential of epigenetic mapping in complex disease.

A variability in response of osteoclasts to zoledronic acid is mediated by smoking-associated modification in the DNA methylome

BACKGROUND

Clinical trials have shown zoledronic acid as a potent bisphosphonate in preventing bone loss, but with varying potency between patients. Human osteoclasts ex vivo reportedly displayed a variable sensitivity to zoledronic acid > 200-fold, determined by the half-maximal inhibitory concentration (IC50), with cigarette smoking as one of the reported contributors to this variation. To reveal the molecular basis of the smoking-mediated variation on treatment sensitivity, we performed a DNA methylome profiling on whole blood cells from 34 healthy female blood donors. Multiple regression models were fitted to associate DNA methylation with ex vivo determined IC50 values, smoking, and their interaction adjusting for age and cell compositions.

RESULTS

We identified 59 CpGs displaying genome-wide significance (p < 1e-08) with a false discovery rate (FDR) < 0.05 for the smoking-dependent association with IC50. Among them, 3 CpGs have p < 1e-08 and FDR < 2e-03. By comparing with genome-wide association studies, 15 significant CpGs were locally enriched (within < 50,000 bp) by SNPs associated with bone and body size measures. Furthermore, through a replication analysis using data from a published multi-omics association study on bone mineral density (BMD), we could validate that 29 out of the 59 CpGs were in close vicinity of genomic sites significantly associated with BMD. Gene Ontology (GO) analysis on genes linked to the 59 CpGs displaying smoking-dependent association with IC50, detected 18 significant GO terms including cation:cation antiporter activity, extracellular matrix conferring tensile strength, ligand-gated ion channel activity, etc. CONCLUSIONS: Our results suggest that smoking mediates individual sensitivity to zoledronic acid treatment through epigenetic regulation. Our novel findings could have important clinical implications since DNA methylation analysis may enable personalized zoledronic acid treatment.

Refining epigenetic prediction of chronological and biological age

BACKGROUND

Epigenetic clocks can track both chronological age (cAge) and biological age (bAge). The latter is typically defined by physiological biomarkers and risk of adverse health outcomes, including all-cause mortality. As cohort sample sizes increase, estimates of cAge and bAge become more precise. Here, we aim to develop accurate epigenetic predictors of cAge and bAge, whilst improving our understanding of their epigenomic architecture.

METHODS

First, we perform large-scale (N = 18,413) epigenome-wide association studies (EWAS) of chronological age and all-cause mortality. Next, to create a cAge predictor, we use methylation data from 24,674 participants from the Generation Scotland study, the Lothian Birth Cohorts (LBC) of 1921 and 1936, and 8 other cohorts with publicly available data. In addition, we train a predictor of time to all-cause mortality as a proxy for bAge using the Generation Scotland cohort (1214 observed deaths). For this purpose, we use epigenetic surrogates (EpiScores) for 109 plasma proteins and the 8 component parts of GrimAge, one of the current best epigenetic predictors of survival. We test this bAge predictor in four external cohorts (LBC1921, LBC1936, the Framingham Heart Study and the Women's Health Initiative study).

RESULTS

Through the inclusion of linear and non-linear age-CpG associations from the EWAS, feature pre-selection in advance of elastic net regression, and a leave-one-cohort-out (LOCO) cross-validation framework, we obtain cAge prediction with a median absolute error equal to 2.3 years. Our bAge predictor was found to slightly outperform GrimAge in terms of the strength of its association to survival (HRGrimAge = 1.47 [1.40, 1.54] with p = 1.08 × 10-52, and HRbAge = 1.52 [1.44, 1.59] with p = 2.20 × 10-60). Finally, we introduce MethylBrowsR, an online tool to visualise epigenome-wide CpG-age associations.

CONCLUSIONS

The integration of multiple large datasets, EpiScores, non-linear DNAm effects, and new approaches to feature selection has facilitated improvements to the blood-based epigenetic prediction of biological and chronological age.

Effect of hubble-bubble smoking on global DNA methylation and transcription levels of protamine and histone genes in human spermatozoa

This study was conducted to assess the impact of hubble-bubble smoking on global DNA methylation, DNA fragmentation; protamine deficiency of spermatozoa, and to determine whether the transcription levels of the protamine and histone genes are different in hubble-bubble smokers compared to nonsmokers. Five hundred semen samples were collected from males with an average age of 32.2 ± 6.1 years (300 hubble-bubble smokers "60%" and 200 nonsmokers "40%"). The nucleic acid was isolated from purified sperm, then ELISA and qPCR were used to evaluate the global DNA methylation and transcription level of protamine and histone, respectively. A significant elevation in global DNA methylation, protamine deficiency, and DNA fragmentation was found in hubble-bubble smokers compared to nonsmokers (P < 0.0001). A significant decline was shown in transcription levels of protamine and histone genes in hubble-bubble compared to nonsmokers (P < 0.0001). Additionally, a down-regulation in the transcription levels of protamine and histone was revealed in hubble-bubble compared to nonsmokers with fold change (0.0001 and 0.007, respectively). In conclusion, this study provided proof that hubble-bubble smoking has a negative impact on global DNA methylation, DNA fragmentation, protamine deficiency, and the transcription of protamine and histone genes in spermatozoa, and these findings influence negatively males' fecundity.

Mediation analysis method review of high throughput data

High-throughput technologies have made high-dimensional settings increasingly common, providing opportunities for the development of high-dimensional mediation methods. We aimed to provide useful guidance for researchers using high-dimensional mediation analysis and ideas for biostatisticians to develop it by summarizing and discussing recent advances in high-dimensional mediation analysis. The method still faces many challenges when extended single and multiple mediation analyses to high-dimensional settings. The development of high-dimensional mediation methods attempts to address these issues, such as screening true mediators, estimating mediation effects by variable selection, reducing the mediation dimension to resolve correlations between variables, and utilizing composite null hypothesis testing to test them. Although these problems regarding high-dimensional mediation have been solved to some extent, some challenges remain. First, the correlation between mediators are rarely considered when the variables are selected for mediation. Second, downscaling without incorporating prior biological knowledge makes the results difficult to interpret. In addition, a method of sensitivity analysis for the strict sequential ignorability assumption in high-dimensional mediation analysis is still lacking. An analyst needs to consider the applicability of each method when utilizing them, while a biostatistician could consider extensions and improvements in the methodology.

Epigenetic aging in Williams syndrome

BACKGROUND

Williams syndrome (WS) is a rare genetic disorder caused by a microdeletion at the 7q11.23 region and is characterized by diverse symptoms encompassing physical and cognitive features. WS was reported to be associated to altered DNA methylation (DNAm) patterns. However, due to the limited information from long-term studies, it remains unclear whether WS accelerates aging. Genome-wide DNAm profiles can serve as "epigenetic clocks" to help estimate biological aging along with age-related markers, such as plasma proteins and telomere length.

METHODS

We investigated GrimAge, DNAm-based telomere length (DNAmTL), and other epigenetic clocks in blood samples of 32 patients with WS and 32 healthy controls.

RESULTS

We observed a significant acceleration in GrimAge, DNAmTL, and other epigenetic clocks in patients with WS as compared with those of controls. In addition, several GrimAge components, such as adrenomedullin, growth differentiation factor-15, leptin and plasminogen activator inhibitor-1, were altered in patients with WS.

CONCLUSIONS

This study provides novel evidence supporting the hypothesis that WS may be associated to accelerated biological aging. A better understanding of the overall underlying biological effects of WS can provide new foundations for improved patient care; thus, long-term follow-up studies are still warranted.

A systematic assessment of cell type deconvolution algorithms for DNA methylation data

We performed systematic assessment of computational deconvolution methods that play an important role in the estimation of cell type proportions from bulk methylation data. The proposed framework methylDeConv (available as an R package) integrates several deconvolution methods for methylation profiles (Illumina HumanMethylation450 and MethylationEPIC arrays) and offers different cell-type-specific CpG selection to construct the extended reference library which incorporates the main immune cell subsets, epithelial cells and cell-free DNAs. We compared the performance of different deconvolution algorithms via simulations and benchmark datasets and further investigated the associations of the estimated cell type proportions to cancer therapy in breast cancer and subtypes in melanoma methylation case studies. Our results indicated that the deconvolution based on the extended reference library is critical to obtain accurate estimates of cell proportions in non-blood tissues.

SNP-by-CpG Site Interactions in ABCA7 Are Associated with Cognition in Older African Americans

SNPs in ABCA7 confer the largest genetic risk for Alzheimer's Disease (AD) in African Americans (AA) after APOE ε4. However, the relationship between ABCA7 and cognitive function has not been thoroughly examined. We investigated the effects of five known AD risk SNPs and 72 CpGs in ABCA7, as well as their interactions, on general cognitive function (cognition) in 634 older AA without dementia from Genetic Epidemiology Network of Arteriopathy (GENOA). Using linear mixed models, no SNP or CpG was associated with cognition after multiple testing correction, but five CpGs were nominally associated (p < 0.05). Four SNP-by-CpG interactions were associated with cognition (FDR q < 0.1). Contrast tests show that methylation is associated with cognition in some genotype groups (p < 0.05): a 1% increase at cg00135882 and cg22271697 is associated with a 0.68 SD decrease and 0.14 SD increase in cognition for those with the rs3764647 GG/AG (p = 0.004) and AA (p = 2 × 10-4) genotypes, respectively. In addition, a 1% increase at cg06169110 and cg17316918 is associated with a 0.37 SD decrease (p = 2 × 10-4) and 0.33 SD increase (p = 0.004), respectively, in cognition for those with the rs115550680 GG/AG genotype. While AD risk SNPs in ABCA7 were not associated with cognition in this sample, some have interactions with proximal methylation on cognition.

Derivation and validation of an epigenetic frailty risk score in population-based cohorts of older adults

DNA methylation (DNAm) patterns in peripheral blood have been shown to be associated with aging related health outcomes. We perform an epigenome-wide screening to identify CpGs related to frailty, defined by a frailty index (FI), in a large population-based cohort of older adults from Germany, the ESTHER study. Sixty-five CpGs are identified as frailty related methylation loci. Using LASSO regression, 20 CpGs are selected to derive a DNAm based algorithm for predicting frailty, the epigenetic frailty risk score (eFRS). The eFRS exhibits strong associations with frailty at baseline and after up to five-years of follow-up independently of established frailty risk factors. These associations are confirmed in another independent population-based cohort study, the KORA-Age study, conducted in older adults. In conclusion, we identify 65 CpGs as frailty-related loci, of which 20 CpGs are used to calculate the eFRS with predictive performance for frailty over long-term follow-up.

From the prodromal stage of multiple sclerosis to disease prevention

A prodrome is an early set of signs or symptoms that indicate the onset of a disease before more typical symptoms develop. Prodromal stages are well recognized in some neurological and immune-mediated diseases such as Parkinson disease, schizophrenia, type 1 diabetes mellitus and rheumatoid arthritis. Emerging evidence indicates that a prodromal stage exists in multiple sclerosis (MS), raising the possibility of intervention at this stage to delay or prevent the development of classical MS. However, much remains unclear about the prodromal stage of MS and considerable research is needed to fully characterize the prodrome and develop standardized criteria to reliably identify individuals with prodromal MS who are at high risk of progressing to a diagnosis of MS. In this Roadmap, we draw on work in other diseases to propose a disease framework for MS that incorporates the prodromal stage, and set out key steps and considerations needed in future research to fully characterize the MS prodrome, identify early disease markers and develop standardized criteria that will enable reliable identification of individuals with prodromal MS, thereby facilitating trials of interventions to slow or stop progression beyond the prodrome.

The association between DNA methylation of 6p21.33 and AHRR in blood and coronary heart disease in Chinese population

BACKGROUND

Early detection could significantly improve the prognosis of coronary heart disease (CHD). In-invitro diagnostic technique may provide a solution when sufficient biomarkers could be identified. Pertinent associations between blood-based aberrant DNA methylation and smoking, the pathogenesis of atherosclerosis, and CHD have been robustly demonstrated and replicated, but that studies in Chinese populations are rare. The blood-based methylation of aryl-hydrocarbon receptor repressor (AHRR) cg05575921 and 6p21.33 cg06126421 has been associated with cardiovascular mortality in Caucasians. Here, we aim to investigate whether the AHRR and 6p21.33 methylation in the blood is associated with CHD in the Chinese population.

METHODS

In this case-control study, 180 CHD patients recruited at their first registration in our study center, and 184 controls randomly selected from the people who participated in the annual health examination were enrolled. Methylation intensities of 19 CpG sites, including AHRR cg05575921, 6p21.33 cg06126421, and their flanking CpG sites, were quantified by mass spectrometry. The association between methylation intensities and CHD was estimated by logistic regression analyses adjusted for covariant.

RESULTS

Compared to the controls, lower methylation of 6p21.33_CpG_4.5/cg06126421 was independently associated with increased odds of being a CHD patient (OR per - 10% methylation = 1.42 after adjustment for age, gender, and batch effect; p = 0.032 by multiple testing corrections). No association between blood-based AHRR methylation and CHD was found.

CONCLUSIONS

6p21.33 methylation exhibits a significant association with CHD. The combination of 6p21.33 methylation and conventional risk factors might be an intermediate step towards the early detection of CHD.

Social Epigenomics: Conceptualizations and Considerations for Oral Health

Advances in high-throughput technologies and the generation of multiomics, such as genomic, epigenomic, transcriptomic, and metabolomic data, are paving the way for the biological risk stratification and prediction of oral diseases. When integrated with electronic health records, survey, census, and/or epidemiologic data, multiomics are anticipated to facilitate data-driven precision oral health, or the delivery of the right oral health intervention to the right individuals/populations at the right time. Meanwhile, multiomics may be modified by a multitude of social exposures, cumulatively along the life course and at various time points from conception onward, also referred to as the socio-exposome. For example, adverse exposures, such as precarious social and living conditions and related psychosocial stress among others, have been linked to specific genes being switched "on and off" through epigenetic mechanisms. These in turn are associated with various health conditions in different age groups and populations. This article argues that considering the impact of the socio-exposome in the biological profiling for precision oral health applications is necessary to ensure that definitions of biological risk do not override social ones. To facilitate the uptake of the socio-exposome in multiomics oral health studies and subsequent interventions, 3 pertinent facets are discussed. First, a summary of the epigenetic landscape of oral health is presented. Next, findings from the nondental literature are drawn on to elaborate the pathways and mechanisms that link the socio-exposome with gene expression-or the biological embedding of social experiences through epigenetics. Then, methodological considerations for implementing social epigenomics into oral health research are highlighted, with emphasis on the implications for study design and interpretation. The article concludes by shedding light on some of the current and prospective opportunities for social epigenomics research applied to the study of life course oral epidemiology.

Predicting risk of lung function impairment and all-cause mortality using a DNA methylation-based classifier of tobacco smoke exposure

BACKGROUND

The Epigenetic Smoking Status Estimator (EpiSmokEr) predicts smoking phenotypes based on DNA methylation at 121 CpG sites.

OBJECTIVE

Evaluate associations of EpiSmokEr-predicted versus self-reported smoking phenotypes with lung function and all-cause mortality in a cohort of older adults.

METHODS

The prospective Normative Aging Study collected DNA methylation measurements from 1999 to 2012 with follow-up through 2016. The R package EpiSmokEr derived predicted smoking phenotypes based on DNA methylation levels assayed by the Illumina HumanMethylation450 Beadchip. Spirometry was collected every 3-5 years. Airflow limitation was defined as forced expiratory volume in 1 s/forced vital capacity <0.7. Vital status was monitored through periodic mailings.

RESULTS

Among 784 participants contributing 5414 person-years of follow-up, the EpiSmokEr-predicted smoking phenotypes matched the self-reported phenotypes for 228 (97%) never smokers and 22 (71%) current smokers. In contrast, EpiSmokEr classified 407 (79%) self-reported former smokers as never smokers. Nonetheless, the EpiSmokEr-predicted former smoking phenotype was more strongly associated with incident airflow limitation (hazard ratio [HR] = 3.15, 95% confidence interval [CI] = 1.50-6.59) and mortality (HR = 2.11, 95% CI = 1.56-2.85) compared to the self-reported former smoking phenotype (airflow limitation: HR = 2.21, 95% CI = 1.13-4.33; mortality: HR = 1.08, 95% CI = 0.86-1.36). Risk of airflow limitation and death did not differ among self-reported never smokers and former smokers who were classified as never smokers. The discriminative accuracy of EpiSmokEr-predicted phenotypes for incident airflow limitation and mortality was improved compared to self-reported phenotypes.

CONCLUSIONS

The DNA methylation-based EpiSmokEr classifier may be a useful surrogate of smoking-induced lung damage and may identify former smokers most at risk of adverse smoking-related health effects.

Cognition, Aryl Hydrocarbon Receptor Repressor Methylation, and Abstinence Duration-Associated Multimodal Brain Networks in Smoking and Long-Term Smoking Cessation

Cigarette smoking and smoking cessation are associated with changes in cognition and DNA methylation; however, the neurobiological correlates of these effects have not been fully elucidated, especially in long-term cessation. Cognitive performance, percent methylation of the aryl hydrocarbon receptor repressor (AHRR) gene, and abstinence duration were used as references to supervise a multimodal fusion analysis of functional, structural, and diffusion magnetic resonance imaging (MRI) data, in order to identify associated brain networks in smokers and ex-smokers. Correlations among these networks and with smoking-related measures were performed. Cognition-, methylation-, and abstinence duration-associated networks discriminated between smokers and ex-smokers and correlated with differences in fractional amplitude of low frequency fluctuations (fALFF) values, gray matter volume (GMV), and fractional anisotropy (FA) values. Long-term smoking cessation was associated with more accurate cognitive performance, as well as lower fALFF and more GMV in the hippocampus complex. The methylation- and abstinence duration-associated networks positively correlated with smoking-related measures of abstinence duration and percent methylation, respectively, suggesting they are complementary measures. This analysis revealed structural and functional co-alterations linked to smoking abstinence and cognitive performance in brain regions including the insula, frontal gyri, and lingual gyri. Furthermore, AHRR methylation, a promising epigenetic biomarker of smoking recency, may provide an important complement to self-reported abstinence duration.

Smoke-induced SAV1 Gene Promoter Hypermethylation Disrupts YAP Negative Feedback and Promotes Malignant Progression of Non-small Cell Lung Cancer

YAP (gene symbol YAP1) as a potential oncoprotein, is positively correlated with the malignancy of various tumors. However, overexpression of YAP alone in multiple normal tissue cells has failed to induce tumor formation and the underlying mechanism is poorly understood. Herein, we show that YAP activation directly induces transcription of its negative regulator, SAV1, to constitute a negative feedback loop, which plays a vital role in maintaining lung epithelial cell homeostasis and was dysregulated in non-small cell lung cancer (NSCLC). Notably, smoking promotes the hypermethylation of the SAV1 promoter region, which disrupts YAP negative feedback by inactivating the Hippo pathway. Besides, exogenous overexpression of SAV1 can act as a traffic protein, activating the Hippo signaling and concurrently inhibiting the WNT pathway to decrease cancer cell growth. Furthermore, using the lung cancer organoids, we found that lentivirus-mediated SAV1 gene transfer combined with methylation inhibitor and YAP-TEAD inhibitor is a potential feasible clinical medication regimen for the lung cancer patient, especially among the smoking population. Thus, this SAV1 mediated feedback loop provides an efficient mechanism to establish the robustness and homeostasis of YAP regulation and as a potential target of gene therapy for the smoking NSCLC population.

Effect of menopausal hormone therapy on methylation levels in early and late postmenopausal women

BACKGROUND

Cardiovascular disease (CVD) remains the leading cause of death among postmenopausal women but standard primary prevention strategies in women are not as effective as in men. By comparison, the Early versus Late Intervention Trial with Estradiol (ELITE) study demonstrated that hormone therapy (HT) was associated with significant reduction in atherosclerosis progression in women who were within six years of menopause compared to those who were 10 or more years from menopause. These findings are consistent with other studies showing significant reductions in all-cause mortality and CVD with HT, particularly when initiated in women younger than 60 years of age or within 10 years since menopause. To explore the biological mechanisms underlying the age-related atheroprotective effects of HT, we investigated changes in methylation of blood cells of postmenopausal women who participated in ELITE.

RESULTS

We first validated the epigenetic data generated from blood leukocytes of ELITE participants by replicating previously known associations between smoking and methylation levels at previously identified CpG sites, such as cg05575921 at the AHRR locus. An epigenome-wide association study (EWAS) evaluating changes in methylation through interactions with time-since-menopause and HT revealed two significantly associated CpG sites on chromosomes 12 (cg19552895; p = 1.1 × 10^-9) and 19 (cg18515510; p = 2.4 × 10^-8). Specifically, HT resulted in modest, but significant, increases in methylation levels at both CpGs but only in women who were 10 or more years since menopause and randomized to HT. Changes in carotid artery intima-media thickness (CIMT) from baseline to 36 months after HT were not significantly correlated with changes in methylation levels at either cg19552895 or cg18515510. Evaluation of other previously identified CpG sites at which methylation levels in either blood or vascular tissue were associated with atherosclerosis also did not reveal any differences in methylation as a function of HT and time-since-menopause or with changes in CIMT.

CONCLUSIONS

We identified specific methylation differences in blood in response to HT among women who were 10 or more years since menopause. The functional consequence of these change with respect to atherosclerosis progression and protective effects of HT remains to be determined and will require additional studies.