Search

Lower muscle strength in midlife predicts disability and mortality in later life. Blood-borne factors, including growth differentiation factor 11 (GDF11), have been linked to muscle regeneration in animal models. We aimed to identify gene transcripts associated with muscle strength in adults. Meta-analysis of whole blood gene expression (overall 17,534 unique genes measured by microarray) and hand-grip strength in four independent cohorts (n= 7,781, ages: 20–104 yr, weighted mean = 56), adjusted for age, sex, height, weight, and leukocyte subtypes. Separate analyses were performed in subsets (older/younger than 60, men/women). Expression levels of 221 genes were associated with strength after adjustment for cofactors and for multiple statistical testing, including ALAS2(rate-limiting enzyme in heme synthesis), PRF1(perforin, a cytotoxic protein associated with inflammation), IGF1R, and IGF2BP2(both insulin like growth factor related). We identified statistical enrichment for hemoglobin biosynthesis, innate immune activation, and the stress response. Ten genes were associated only in younger individuals, four in men only and one in women only. For example, PIK3R2(a negative regulator of PI3K/AKT growth pathway) was negatively associated with muscle strength in younger (<60 yr) individuals but not older (≥60 yr). We also show that 115 genes (52%) have not previously been linked to muscle in NCBI PubMed abstracts. This first large-scale transcriptome study of muscle strength in human adults confirmed associations with known pathways and provides new evidence for over half of the genes identified. There may be age- and sex-specific gene expression signatures in blood for muscle strength. [ABSTRACT FROM AUTHOR]

Background: Impaired pulmonary function carries significant risks for lung, cardiovascular, and metabolic disorders., Research Question: Can circulating protein biomarkers of pulmonary function provide insights into the pathophysiology of lung function impairment and links to comorbidities?, Study Design and Methods: We analyzed plasma levels of 2922 proteins in 32,493 UK Biobank (UKB) participants (53% women, age=57±8 years) to investigate their associations with spirometry measures of lung-function (1-second forced-expiratory volume [FEV1], forced-vital capacity [FVC], FEV1/FVC ratio), and with obstructive (N=4713) and restrictive (N=3886) spirometry patterns. Significant protein signatures were functionally annotated and externally validated in 740 Framingham Heart Study participants (FHS). We inferred causality using Mendelian randomization and examined colocalization of genetic signals of protein biomarkers with corresponding lung traits., Results: In UKB, we identified 1240 proteins significantly associated (P UKB <0.000017) with FEV1, 1310 with FVC, and 513 with FEV1/FVC. Of these, 44, 99, and 13 proteins, respectively, were nominally significant (P FHS <0.01) in FHS. Plasma levels of 737 proteins (7 with P FHS <0.01) differed in individuals with an obstructive spirometry pattern (OSP), and 811 proteins (38 with P FHS <0.01) differed in restrictive spirometry pattern (RSP) compared to normal spirometry in the UKB. Putatively causal relations to FEV1, FVC, FEV1/FVC, and OSP were observed for 55, 63, 28, and 14 proteins, respectively. Of note, several circulating decoy receptors, including interleukin-1 receptor-like 1, tumor necrosis factor receptor superfamily member-6B, and macrophage scavenger receptor-1 emerged as causal and protective biomarkers of lung function. Enrichment analysis suggested a connection between reduced lung function and systemic inflammation driven by adipose tissue dysfunction and gut-dysbiosis. Protein biomarkers associated with lung function also were enriched for susceptibility to cardiovascular conditions and cancers., Interpretation: This study identifies proteomic signatures of reduced lung function linked to comorbidities, paving the way for improved diagnostics and treatment for lung disease., (Copyright © 2024. Published by Elsevier Inc.)

MicroRNAs, crucial in regulating protein-coding gene expression, are implicated in various diseases. We performed a genome-wide association study of plasma miRNAs (ex-miRNAs) in 3,743 Framingham Heart Study (FHS) participants and identified 1,027 cis- ex-miRNA-eQTLs ( cis- exQTLs) for 37 ex-miRNAs, with 55% replication in an independent study. Colocalization analyses suggested potential genetic coregulation of ex-miRNAs with whole blood mRNAs. Mendelian randomization indicated 29 ex-miRNAs potentially influencing 35 traits. Notably, the chromosome 14q23 and 14q32 miRNA clusters emerged as the top signal, contributing over 50% of the significant cis- exQTL results, and were associated with a diverse range of traits including platelet count. Correlations of 10 ex-miRNAs (such as miR-376c-3p) in 14q32 with platelet count and volume were confirmed in FHS participants. These findings shed light on the genetic basis of ex-miRNA expression and their involvement in complex traits., Competing Interests: The authors declare no conflict of interest.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n  = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel ( CADM3 , NALCN, NLRC5, ZNF792 , and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level ( p  < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5 , NOD2 , and AIM2 . Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

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., Competing Interests: Competing interests F.P. is an employee of Biomodal. L.M.R serves as a consultant for the NHLBI TOPMed Administrative Coordinating Center (through Westat). P.N. reports research grants from Allelica, Amgen, Apple, Boston Scientific, Genentech / Roche, and Novartis, personal fees from Allelica, Apple, AstraZeneca, Blackstone Life Sciences, Creative Education Concepts, CRISPR Therapeutics, Eli Lilly & Co, Foresite Labs, Genentech / Roche, GV, HeartFlow, Magnet Biomedicine, Merck, and Novartis, scientific advisory board membership of Esperion Therapeutics, Preciseli, TenSixteen Bio, and Tourmaline Bio, scientific co-founder of TenSixteen Bio, equity in MyOme, Preciseli, and TenSixteen Bio, and spousal employment at Vertex Pharmaceuticals, all unrelated to the present work. Psaty serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. No other authors have competing interests. The Jimma University Internal Review Board approved the study ensuring all procedures involving human participation adhere to ethical guidelines established by both the Institution and National Research Committee before any research began.

Machine learning methods have been used in identifying omics markers for a variety of phenotypes. We aimed to examine whether a supervised machine learning algorithm can improve identification of alcohol-associated transcriptomic markers. In this study, we analysed array-based, whole-blood derived expression data for 17 873 gene transcripts in 5508 Framingham Heart Study participants. By using the Boruta algorithm, a supervised random forest (RF)-based feature selection method, we selected twenty-five alcohol-associated transcripts. In a testing set (30 % of entire study participants), AUC (area under the receiver operating characteristics curve) of these twenty-five transcripts were 0·73, 0·69 and 0·66 for non-drinkers v . moderate drinkers, non-drinkers v . heavy drinkers and moderate drinkers v . heavy drinkers, respectively. The AUC of the selected transcripts by the Boruta method were comparable to those identified using conventional linear regression models, for example, AUC of 1958 transcripts identified by conventional linear regression models (false discovery rate < 0·2) were 0·74, 0·66 and 0·65, respectively. With Bonferroni correction for the twenty-five Boruta method-selected transcripts and three CVD risk factors (i.e. at P < 6·7e-4), we observed thirteen transcripts were associated with obesity, three transcripts with type 2 diabetes and one transcript with hypertension. For example, we observed that alcohol consumption was inversely associated with the expression of DOCK4 , IL4R , and SORT1 , and DOCK4 and SORT1 were positively associated with obesity, and IL4R was inversely associated with hypertension. In conclusion, using a supervised machine learning method, the RF-based Boruta algorithm, we identified novel alcohol-associated gene transcripts.

Background: Extracellular microRNAs (miRNAs) are a class of noncoding RNAs that remain stable in the extracellular milieu, where they contribute to various physiological and pathological processes by facilitating intercellular signaling. Previous studies have reported associations between miRNAs and cardiovascular diseases (CVDs); however, the plasma miRNA signatures of CVD and its risk factors have not been fully elucidated at the population level., Methods and Results: Plasma miRNA levels were measured in 4440 FHS (Framingham Heart Study) participants. Linear regression analyses were conducted to test the cross-sectional associations of each miRNA with 8 CVD risk factors. Prospective analyses of the associations of miRNAs with new-onset obesity, hypertension, type 2 diabetes, CVD, and all-cause mortality were conducted using proportional hazards regression. Replication was carried out in 1999 RS (Rotterdam Study) participants. Pathway enrichment analyses were conducted and target genes were predicted for miRNAs associated with ≥5 risk factors in the FHS. In the FHS, 6 miRNAs (miR-193b-3p, miR-122-5p, miR-365a-3p, miR-194-5p, miR-192-5p, and miR-193a-5p) were associated with ≥5 risk factors. This miRNA signature was enriched for pathways associated with CVD and several genes annotated to these pathways were predicted targets of the identified miRNAs. Furthermore, miR-193b-3p, miR-194-5p, and miR-193a-5p were each associated with ≥2 risk factors in the RS. Prospective analysis revealed 8 miRNAs associated with all-cause mortality in the FHS., Conclusions: These findings highlight associations between miRNAs and CVD risk factors that may provide valuable insights into the underlying pathogenesis of CVD.

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.

Background: Dietary intake of antioxidants such as vitamins C and E protect against oxidative stress, and may also be associated with altered DNA methylation patterns., Methods: We meta-analysed epigenome-wide association study (EWAS) results from 11,866 participants across eight population-based cohorts to evaluate the association between self-reported dietary and supplemental intake of vitamins C and E with DNA methylation. EWAS were adjusted for age, sex, BMI, caloric intake, blood cell type proportion, smoking status, alcohol consumption, and technical covariates. Significant results of the meta-analysis were subsequently evaluated in gene set enrichment analysis (GSEA) and expression quantitative trait methylation (eQTM) analysis., Results: In meta-analysis, methylation at 4,656 CpG sites was significantly associated with vitamin C intake at FDR ≤ 0.05. The most significant CpG sites associated with vitamin C (at FDR ≤ 0.01) were enriched for pathways associated with systems development and cell signalling in GSEA, and were associated with downstream expression of genes enriched in the immune response in eQTM analysis. Furthermore, methylation at 160 CpG sites was significantly associated with vitamin E intake at FDR ≤ 0.05, but GSEA and eQTM analysis of the top most significant CpG sites associated with vitamin E did not identify significant enrichment of any biological pathways investigated., Conclusions: We identified significant associations of many CpG sites with vitamin C and E intake, and our results suggest that vitamin C intake may be associated with systems development and the immune response.

Background: Metabolite signatures of long-term alcohol consumption are lacking. To better understand the molecular basis linking alcohol drinking and cardiovascular disease (CVD), we investigated circulating metabolites associated with long-term alcohol consumption and examined whether these metabolites were associated with incident CVD., Methods: Cumulative average alcohol consumption (g/day) was derived from the total consumption of beer, wine, and liquor on average of 19 years in 2428 Framingham Heart Study Offspring participants (mean age 56 years, 52% women). We used linear mixed models to investigate the associations of alcohol consumption with 211 log-transformed plasma metabolites, adjusting for age, sex, batch, smoking, diet, physical activity, BMI, and familial relationship. Cox models were used to test the association of alcohol-related metabolite scores with fatal and nonfatal incident CVD (myocardial infarction, coronary heart disease, stroke, and heart failure)., Results: We identified 60 metabolites associated with cumulative average alcohol consumption (p < 0.05/211 ≈ 0.00024). For example, 1 g/day increase of alcohol consumption was associated with higher levels of cholesteryl esters (e.g., CE 16:1, beta = 0.023 ± 0.002, p = 6.3e - 45) and phosphatidylcholine (e.g., PC 32:1, beta = 0.021 ± 0.002, p = 3.1e - 38). Survival analysis identified that 10 alcohol-associated metabolites were also associated with a differential CVD risk after adjusting for age, sex, and batch. Further, we built two alcohol consumption weighted metabolite scores using these 10 metabolites and showed that, with adjustment age, sex, batch, and common CVD risk factors, the two scores had comparable but opposite associations with incident CVD, hazard ratio 1.11 (95% CI = [1.02, 1.21], p = 0.02) vs 0.88 (95% CI = [0.78, 0.98], p = 0.02)., Conclusions: We identified 60 long-term alcohol consumption-associated metabolites. The association analysis with incident CVD suggests a complex metabolic basis between alcohol consumption and CVD., (© 2023. The Author(s).)

Background: Insulin resistance (IR) is a major risk factor for Alzheimer's disease (AD) dementia. The mechanisms by which IR predisposes to AD are not well-understood. Epigenetic studies may help identify molecular signatures of IR associated with AD, thus improving our understanding of the biological and regulatory mechanisms linking IR and AD., Methods: We conducted an epigenome-wide association study of IR, quantified using the homeostatic model assessment of IR (HOMA-IR) and adjusted for body mass index, in 3,167 participants from the Framingham Heart Study (FHS) without type 2 diabetes at the time of blood draw used for methylation measurement. We identified DNA methylation markers associated with IR at the genome-wide level accounting for multiple testing (P < 1.1 × 10 -7 ) and evaluated their association with neurological traits in participants from the FHS (N = 3040) and the Religious Orders Study/Memory and Aging Project (ROSMAP, N = 707). DNA methylation profiles were measured in blood (FHS) or dorsolateral prefrontal cortex (ROSMAP) using the Illumina HumanMethylation450 BeadChip. Linear regressions (ROSMAP) or mixed-effects models accounting for familial relatedness (FHS) adjusted for age, sex, cohort, self-reported race, batch, and cell type proportions were used to assess associations between DNA methylation and neurological traits accounting for multiple testing., Results: We confirmed the strong association of blood DNA methylation with IR at three loci (cg17901584-DHCR24, cg17058475-CPT1A, cg00574958-CPT1A, and cg06500161-ABCG1). In FHS, higher levels of blood DNA methylation at cg00574958 and cg17058475 were both associated with lower IR (P = 2.4 × 10 -11 and P = 9.0 × 10 -8 ), larger total brain volumes (P = 0.03 and P = 9.7 × 10 -4 ), and smaller log lateral ventricular volumes (P = 0.07 and P = 0.03). In ROSMAP, higher levels of brain DNA methylation at the same two CPT1A markers were associated with greater risk of cognitive impairment (P = 0.005 and P = 0.02) and higher AD-related indices (CERAD score: P = 5 × 10 -4 and 0.001; Braak stage: P = 0.004 and P = 0.01)., Conclusions: Our results suggest potentially distinct epigenetic regulatory mechanisms between peripheral blood and dorsolateral prefrontal cortex tissues underlying IR and AD at CPT1A locus., (© 2023. The Author(s).)

The alcohol-associated biological aging remains to be studied across adulthood. We conducted linear regression analyses to investigate the associations between alcohol consumption and two DNA methylation-based biological age acceleration metrics in 3823 Framingham Heart Study participants (24-92 years and 53.8% women) adjusting for covariates. We also investigated whether the two epigenetic aging metrics mediated the association of alcohol consumption with hypertension. We found that higher long-term average alcohol consumption was significantly associated with biological age acceleration assessed by GrimAge acceleration (GAA) and PhenoAge acceleration (PAA) in middle-aged (45-64 years, n = 1866) and older (65-92 years, n = 1267) participants while not in young participants (24-44 years, n = 690). For example, one additional standard drink of alcohol (~14 grams of ethanol per day) was associated with a 0.71 ± 0.15-year ( p = 2.1e-6) and 0.60 ± 0.18-year ( p = 7.5e-4) increase in PAA in middle-aged and older participants, respectively, but the association was not significant in young participants ( p = 0.23). One additional standard serving of liquor (~14 grams of ethanol) was associated with a greater increase in GAA (0.82-year, p = 4.8e-4) and PAA (1.45-year, p = 7.4e-5) than beer (GAA: 0.45-year, p = 5.2e-4; PAA: 0.48-year, p = 0.02) and wine (GAA: 0.51-year, p = 0.02; PAA: 0.91-year, p = 0.008) in middle-aged participant group. We observed that up to 28% of the association between alcohol consumption and hypertension was mediated by GAA or PAA in the pooled sample. Our findings suggest that alcohol consumption is associated with greater biological aging quantified by epigenetic aging metrics, which may mediate the association of alcohol consumption with quantitative traits, such as hypertension.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs., Competing Interests: Declaration of interests P.N. reports research grants from Allelica, Apple, Amgen, Boston Scientific, Genentech/Roche, and Novartis; personal fees from Allelica, Apple, AstraZeneca, Blackstone Life Sciences, Eli Lilly & Co, Foresite Labs, Genentech/Roche, GV, HeartFlow, Magnet Biomedicine, and Novartis; scientific advisory board membership of Esperion Therapeutics, Preciseli, and TenSixteen Bio; scientific co-founder of TenSixteen Bio; equity in MyOme, Preciseli, and TenSixteen Bio; and spousal employment at Vertex Pharmaceuticals, all unrelated to the present work. B.M.P. serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. L.M.R., S.S.R., and R.M. are consultants for the TOPMed Administrative Coordinating Center (through Westat). M.E.M. receives funding from Regeneron Pharmaceutical Inc. unrelated to this work. X. Lin is a consultant of AbbVie Pharmaceuticals and Verily Life Sciences. P.T.E. receives sponsored research support from Bayer AG, IBM Research, Bristol Myers Squibb, Pfizer, and Novo Nordisk; he has also served on advisory boards or consulted for Bayer AG, MyoKardia, and Novartis. A.P.C. previously received investigator-initiated grant support from Amgen, Inc. unrelated to the present work., (Copyright © 2023 American Society of Human Genetics. All rights reserved.)

Altered DNA methylation (DNAm) might be a biological intermediary in the pathway from smoking to lung cancer. In this study, we investigated the contribution of differential blood DNAm to explain the association between smoking and lung cancer incidence. Blood DNAm was measured in 2321 Strong Heart Study (SHS) participants. Incident lung cancer was assessed as time to event diagnoses. We conducted mediation analysis, including validation with DNAm and paired gene expression data from the Framingham Heart Study (FHS). In the SHS, current versus never smoking and pack-years single-mediator models showed, respectively, 29 and 21 differentially methylated positions (DMPs) for lung cancer with statistically significant mediated effects (14 of 20 available, and five of 14 available, positions, replicated, respectively, in FHS). In FHS, replicated DMPs showed gene expression downregulation largely in trans, and were related to biological pathways in cancer. The multimediator model identified that DMPs annotated to the genes AHRR and IER3 jointly explained a substantial proportion of lung cancer. Thus, the association of smoking with lung cancer was partly explained by differences in baseline blood DNAm at few relevant sites. Experimental studies are needed to confirm the biological role of identified eQTMs and to evaluate potential implications for early detection and control of lung cancer., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Ana Navas-Acien reports financial support was provided by National Heart Lung and Blood Institute. Ana Navas-Acien reports financial support was provided by National Institute of Environmental Health Sciences. Daniel Levy reports financial support was provided by National Institutes of Health. Arce Domingo-Relloso reports financial support was provided by La Caixa Foundation. Maria Tellez-Plaza reports financial support provided by AstraZeneca Spain, the Instituto de Salud Carlos III and the Spanish Ministry of Science and Innovation and co-funded with European Funds for Regional Development., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Background: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases., Methods: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables., Findings: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P < 7.94E-04) that can be explored as potential therapeutic targets., Interpretation: By integrating eQTM and MR analyses in general and clinical asthma populations, our findings provide a deeper understanding of the multidimensional inter-relations of DNA methylation, gene expression, and IgE levels., Funding: US NIH/NHLBI grants: P01HL132825, K99HL159234. N01-HC-25195 and HHSN268201500001I., Competing Interests: Declaration of interests DLD has received grants from Bayer and honoraria from Novartis. STW has received royalties from UpToDate and served on the Board of Histolix, a digital pathology company. Rest authors declared no conflicts of interest., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p < 1E-7 where the top most significant eGenes (i.e., gene transcripts associated with a CpG) were enriched in biological pathways related to cell signaling, and for 1208 clinical traits (enrichment false discovery rate [FDR] ≤ 0.05). We also identified 246,667 significant trans CpG-transcript pairs at p < 1E-14 where the top most significant eGenes were enriched in biological pathways related to activation of the immune response, and for 1191 clinical traits (enrichment FDR ≤ 0.05). Independent and external replication of the top 1000 significant cis and trans CpG-transcript pairs was completed in the Women's Health Initiative and Jackson Heart Study cohorts. Using significant cis CpG-transcript pairs, we identified significant mediation of the association between CpG sites and cardiometabolic traits through gene expression and identified shared genetic regulation between CpGs and transcripts associated with cardiometabolic traits. In conclusion, we developed a robust and powerful resource of whole blood eQTM CpG-transcript pairs that can help inform future functional studies that seek to understand the molecular basis of disease., (© 2023. Springer Nature Limited.)

Introduction: We investigated associations of obesity with the expression of Alzheimer's disease (AD)-related genes in a large community-based cohort., Methods: The sample consisted of 5619 participants from the Framingham Heart Study. Obesity metrics included body mass index (BMI) and waist-to-hip ratio (WHR). Gene expression was measured for a set of 74 AD-related genes, derived by integrating genome-wide association study results with functional genomics data., Results: Obesity metrics were associated with the expression of 21 AD-related genes. The strongest associations were observed with CLU, CD2AP, KLC3, and FCER1G. Unique associations were noted with TSPAN14, SLC24A4 for BMI, and ZSCAN21, BCKDK for WHR. After adjustment for cardiovascular risk factors, 13 associations remained significant for BMI and 8 for WHR. Dichotomous obesity metrics exhibited unique associations with EPHX2 for BMI, and with TSPAN14 for WHR., Discussion: Obesity was associated with AD-related gene expression; these findings shed light on the molecular pathways linking obesity to AD., (© 2023 the Alzheimer's Association.)

Rheumatoid arthritis (RA) is a risk factor for atherosclerotic cardiovascular diseases (CVD). Given the critical roles of the immune system and inflammatory signals in the pathogenesis of CVD, we hypothesized that interrogation of CVD-related proteins using integrative genomics might provide new insights into the pathophysiology of RA. We utilized two-sample Mendelian randomization (MR) for causal inference between circulating protein levels and RA by incorporating genetic variants, followed by colocalization to characterize the causal associations. Genetic variants from three sources were obtained: those associated with 71 CVD-related proteins measured in nearly 7000 Framingham Heart Study participants, a published genome-wide association study (GWAS) of RA (19 234 cases, 61 565 controls), and GWAS of rheumatoid factor (RF) levels from the UK Biobank (n = 30 565). We identified the soluble receptor for advanced glycation end products (sRAGE), a critical inflammatory pathway protein, as putatively causal and protective for both RA (odds ratio per 1-standard deviation increment in inverse-rank normalized sRAGE level = 0.364; 95% confidence interval 0.342-0.385; P = 6.40 × 10 -241 ) and RF levels (β [change in RF level per sRAGE increment] = - 1.318; SE = 0.434; P = 0.002). Using an integrative genomic approach, we highlight the AGER/RAGE axis as a putatively causal and promising therapeutic target for RA., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Measurement of circulating immunoglobulin E (IgE) concentration is helpful for diagnosing and treating asthma and allergic diseases. Identifying gene expression signatures associated with IgE might elucidate novel pathways for IgE regulation. To this end, we performed a discovery transcriptome-wide association study to identify differentially expressed genes associated with circulating IgE levels in whole-blood derived RNA from 5,345 participants in the Framingham Heart Study across 17,873 mRNA gene-level transcripts. We identified 216 significant transcripts at a false discovery rate <0.05. We conducted replication using the meta-analysis of two independent external studies: the Childhood Asthma Management Program (n=610) and the Genetic Epidemiology of Asthma in Costa Rica Study (n=326); we then reversed the discovery and replication cohorts, which revealed 59 significant genes that replicated in both directions. Gene ontology analysis revealed that many of these genes were implicated in immune function pathways, including defense response, inflammatory response, and cytokine production. Mendelian randomization (MR) analysis revealed four genes ( CLC , CCDC21 , S100A13 , and GCNT1 ) as putatively causal ( p <0.05) regulators of IgE levels. GCNT1 (beta=1.5, p =0.01)-which is a top result in the MR analysis of expression in relation to asthma and allergic diseases-plays a role in regulating T helper type 1 cell homing, lymphocyte trafficking, and B cell differentiation. Our findings build upon prior knowledge of IgE regulation and provide a deeper understanding of underlying molecular mechanisms. The IgE-associated genes that we identified-particularly those implicated in MR analysis-can be explored as promising therapeutic targets for asthma and IgE-related diseases., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Recto, Huan, Lee, Lee, Gereige, Yao, Hwang, Joehanes, Kelly, Lasky-Su, O’Connor and Levy.)

Background: The relations of alcohol consumption and gene expression remain to be elucidated., Materials and Methods: We examined cross-sectional associations between alcohol consumption and whole blood derived gene expression levels and between alcohol-associated genes and obesity, hypertension, and diabetes in 5531 Framingham Heart Study (FHS) participants., Results: We identified 25 alcohol-associated genes. We further showed cross-sectional associations of 16 alcohol-associated genes with obesity, nine genes with hypertension, and eight genes with diabetes at P < 0.002. For example, we observed decreased expression of PROK2 (β = -0.0018; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) and PAX5 (β = -0.0014; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) per 1 g/day increase in alcohol consumption. Consistent with our previous observation on the inverse association of alcohol consumption with obesity and positive association of alcohol consumption with hypertension, we found that PROK2 was positively associated with obesity (OR = 1.42; 95%CI: 1.17, 1.72; P = 4.5e - 4) and PAX5 was negatively associated with hypertension (OR = 0.73; 95%CI: 0.59, 0.89; P = 1.6e - 3). We also observed that alcohol consumption was positively associated with expression of ABCA13 (β = 0.0012; 95%CI: 0.0007, 0.0017; P = 1.3e - 6) and ABCA13 was positively associated with diabetes (OR = 2.57; 95%CI: 1.73, 3.84; P = 3.5e - 06); this finding, however, was inconsistent with our observation of an inverse association between alcohol consumption and diabetes., Conclusions: We showed strong cross-sectional associations between alcohol consumption and expression levels of 25 genes in FHS participants. Nonetheless, complex relationships exist between alcohol-associated genes and CVD risk factors., (© The Author(s) 2022. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Purpose: Assessing an individual's systemic inflammatory state is vital to understand inflammation's role in cardiometabolic diseases and identify those at the greatest risk of disease. We generated global inflammation scores and investigated their associations with cardiometabolic risk factors and adverse outcomes., Patients and Methods: Aggregate Inflammation Scores (AIS) and Principal Component Analysis (PCA) scores were generated for 7287 Framingham Heart Study participants using up to 26 inflammation-related proteins, with higher scores reflecting a pro-inflammatory milieu. Multivariable regression and proportional hazards analyses were conducted to investigate the associations of inflammation with cardiometabolic risk factors and outcomes. The primary outcomes for cross-sectional analyses included age, cigarette smoking, fasting lipid and glucose levels, blood pressure, body mass index (BMI), and hypertension, diabetes, and obesity. For prospective analyses, new-onset hypertension, diabetes, obesity, cardiovascular disease and all-cause mortality were investigated., Results: Higher inflammation scores were associated with smoking and older age, higher BMI, systolic blood pressure, lipids, and glucose levels, and with greater odds of hypertension and diabetes after adjusting for age, sex, cohort, and BMI (all p < 0.001). Higher baseline scores were associated with greater odds of new-onset hypertension after adjusting for traditional risk factors (OR [95% CI] per one standard deviation [1-SD] increase, AIS: 1.33 [1.21-1.47], PCA score: 1.26 [1.12-1.42], p < 0.001). The AIS also was associated with new-onset diabetes (1.32 [1.14-1.52], p < 0.001). Proportional hazards analyses revealed greater risk of new-onset cardiovascular disease events and all-cause mortality (HR [95% CI] per 1-SD, AIS: 1.25 [1.14-1.37] and 1.32 [1.23-1.42], PCA score: 1.22 [1.13-1.33] and 1.40 [1.31-1.49], p < 0.001)., Conclusion: Global inflammation scores encompassing an array of pro- and anti-inflammatory proteins and pathways may enhance risk assessment for cardiometabolic diseases. The AIS and PCA scores provide further opportunities to investigate the mechanisms of inflammation-related risk of disease., Competing Interests: Potential conflict of interest: The co-author C. Yao aided in the completion of this project, while an employee of the National Heart, Lung, and Blood Institute; Dr Yao is currently an employee of Bristol Myers Squibb. Jennifer E Ho reports grants from NIH, grants from Bayer, AG, outside the submitted work. The authors declare that they have no other conflicts of interest., (© 2022 Tejada et al.)

Background: As omics measurements profiled on different molecular layers are interconnected, integrative approaches that incorporate the regulatory effect from multi-level omics data are needed. When the multi-level omics data are from the same individuals, gene expression (GE) clusters can be identified using information from regulators like genetic variants and DNA methylation. When the multi-level omics data are from different individuals, the choice of integration approaches is limited., Methods: We developed an approach to improve GE clustering from microarray data by integrating regulatory data from different but partially overlapping sets of individuals. We achieve this through (1) decomposing gene expression into the regulated component and the other component that is not regulated by measured factors, (2) optimizing the clustering goodness-of-fit objective function. We do not require the availability of different omics measurements on all individuals. A certain amount of individual overlap between GE data and the regulatory data is adequate for modeling the regulation, thus improving GE clustering., Results: A simulation study shows that the performance of the proposed approach depends on the strength of the GE-regulator relationship, degree of missingness, data dimensionality, sample size, and the number of clusters. Across the various simulation settings, the proposed method shows competitive performance in terms of accuracy compared to the alternative K-means clustering method, especially when the clustering structure is due mostly to the regulated component, rather than the unregulated component. We further validate the approach with an application to 8,902 Framingham Heart Study participants with data on up to 17,873 genes and regulation information of DNA methylation and genotype from different but partially overlapping sets of participants. We identify clustering structures of genes associated with pulmonary function while incorporating the predicted regulation effect from the measured regulators. We further investigate the over-representation of these GE clusters in pathways of other diseases that may be related to lung function and respiratory health., Conclusion: We propose a novel approach for clustering GE with the assistance of regulatory data that allowed for different but partially overlapping sets of individuals to be included in different omics data., (© 2022. The Author(s).)

Aim: To identify differential methylation related to prescribed opioid use. Methods: This study examined whether blood DNA methylation, measured using Illumina arrays, differs by recent opioid medication use in four population-based cohorts. We meta-analyzed results (282 users; 10,560 nonusers) using inverse-variance weighting. Results: Differential methylation (false discovery rate <0.05) was observed at six CpGs annotated to the following genes: KIAA0226 , CPLX2 , TDRP , RNF38 , TTC23 and GPR179 . Integrative epigenomic analyses linked implicated loci to regulatory elements in blood and/or brain. Additionally, 74 CpGs were differentially methylated in males or females. Methylation at significant CpGs correlated with gene expression in blood and/or brain. Conclusion: This study identified DNA methylation related to opioid medication use in general populations. The results could inform the development of blood methylation biomarkers of opioid use.

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis-eQTL variant-gene transcript (eGene) pairs at p < 5 × 10 -8 (2,855,111 unique cis-eQTL variants and 15,982 unique eGenes) and 1,469,754 trans-eQTL variant-eGene pairs at p < 1e-12 (526,056 unique trans-eQTL variants and 7233 unique eGenes). In addition, 442,379 cis-eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis-eGenes are enriched for immune functions (FDR < 0.05). The cis-eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases., (© 2022. The Author(s).)

DNA methylation commonly occurs at cytosine-phosphate-guanine sites (CpGs) that can serve as biomarkers for many diseases. We analyzed whole genome sequencing data to identify DNA methylation quantitative trait loci (mQTLs) in 4126 Framingham Heart Study participants. Our mQTL mapping identified 94,362,817 cis-mQTLvariant-CpG pairs (for 210,156 unique autosomal CpGs) at P < 1e-7 and 33,572,145 trans-mQTL variant-CpG pairs (for 213,606 unique autosomal CpGs) at P < 1e-14. Using cis-mQTL variants for 1258 CpGs associated with seven cardiovascular disease (CVD) risk factors, we found 104 unique CpGs that colocalized with at least one CVD trait. For example, cg11554650 (PPP1R18) colocalized with type 2 diabetes, and was driven by a single nucleotide polymorphism (rs2516396). We performed Mendelian randomization (MR) analysis and demonstrated 58 putatively causal relations of CVD risk factor-associated CpGs to one or more risk factors (e.g., cg05337441 [APOB] with LDL; MR P = 1.2e-99, and 17 causal associations with coronary artery disease (e.g. cg08129017 [SREBF1] with coronary artery disease; MR P = 5e-13). We also showed that three CpGs, e.g., cg14893161 (PM20D1), are putatively causally associated with COVID-19 severity. To assist in future analyses of the role of DNA methylation in disease pathogenesis, we have posted a comprehensive summary data set in the National Heart, Lung, and Blood Institute's BioData Catalyst., (© 2022. The Author(s).)

Cigarette smoking is associated with epigenetic changes that may be reversible following smoking cessation. Whole blood DNA methylation was evaluated in Framingham Heart Study Offspring (n = 169) and Third Generation (n = 30) cohort participants at two study visits 6 years apart and in Atherosclerosis Risk in Communities (ARIC) study (n = 222) participants at two study visits 20 years apart. Changes in DNA methylation (delta β values) at 483,565 cytosine-phosphate-guanine (CpG) sites and differentially methylated regions (DMRs) were compared between participants who were current, former, or never smokers at both visits (current-current, former-former, never-never, respectively), versus those who quit in the interim (current-former). Interim quitters had more hypermethylation at four CpGs annotated to AHRR , one CpG annotated to F2RL3 , and one intergenic CpG (cg21566642) compared with current-current smokers (FDR < 0.02 for all), and two significant DMRs were identified. While there were no significant differentially methylated CpGs in the comparison of interim quitters and former-former smokers, 106 DMRs overlapping with small nucleolar RNA were identified. As compared with all non-smokers, current-current smokers additionally had more hypermethylation at two CpG sites annotated to HIVEP3 and TMEM126A , respectively, and another intergenic CpG (cg14339116). Gene transcripts associated with smoking cessation were implicated in immune responses, cell homoeostasis, and apoptosis. Smoking cessation is associated with early reversion of blood DNA methylation changes at CpG sites annotated to AHRR and F2RL3 towards those of never smokers. Associated gene expression suggests a role of longitudinal smoking-related DNA methylation changes in immune response processes.

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV 1 , FVC, and FEV 1 /FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV 1 , FVC, or FEV 1 /FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Background: Epigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD., Methods: Blood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE -/- ) mouse model of atherosclerosis., Results: A total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic., Conclusions: Differential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

DNA methylation (DNAm) has been reported to be associated with many diseases and with mortality. We hypothesized that the integration of DNAm with clinical risk factors would improve mortality prediction. We performed an epigenome-wide association study of whole blood DNAm in relation to mortality in 15 cohorts (n = 15,013). During a mean follow-up of 10 years, there were 4314 deaths from all causes including 1235 cardiovascular disease (CVD) deaths and 868 cancer deaths. Ancestry-stratified meta-analysis of all-cause mortality identified 163 CpGs in European ancestry (EA) and 17 in African ancestry (AA) participants at p < 1 × 10 -7 , of which 41 (EA) and 16 (AA) were also associated with CVD death, and 15 (EA) and 9 (AA) with cancer death. We built DNAm-based prediction models for all-cause mortality that predicted mortality risk after adjusting for clinical risk factors. The mortality prediction model trained by integrating DNAm with clinical risk factors showed an improvement in prediction of cancer death with 5% increase in the C-index in a replication cohort, compared with the model including clinical risk factors alone. Mendelian randomization identified 15 putatively causal CpGs in relation to longevity, CVD, or cancer risk. For example, cg06885782 (in KCNQ4) was positively associated with risk for prostate cancer (Beta = 1.2, P MR  = 4.1 × 10 -4 ) and negatively associated with longevity (Beta = -1.9, P MR  = 0.02). Pathway analysis revealed that genes associated with mortality-related CpGs are enriched for immune- and cancer-related pathways. We identified replicable DNAm signatures of mortality and demonstrated the potential utility of CpGs as informative biomarkers for prediction of mortality risk., (© 2022 The Authors. Aging Cell published by Anatomical Society and John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD., (© 2022. The Author(s).)

Rationale: It has been speculated that shared mechanisms underlie respiratory and cardiovascular diseases (CVD) including systemic inflammation or mutual risk factors. In this context, we sought to examine the associations of CVD-related plasma proteins with lung function as measured by spirometry in a large community-based cohort of adults., Methods: The study included 5777 Framingham Heart Study participants who had spirometry and measurement of 71 CVD-related plasma proteins. The association of plasma proteins with lung function was assessed cross-sectionally and longitudinally using models accounting for familial correlations. Linear mixed models were used for the following measurements: FEV1%predicted, FVC%predicted, and FEV1/FVC ratio with secondary analyses examining obstructive and restrictive physiology at baseline and their new onset during follow up., Measurements and Main Results: Among the 71 CVD-related plasma proteins, 13 proteins were associated in cross-sectional analyses with FEV1%predicted, 17 proteins were associated with FVC%predicted, and 1 protein was associated with FEV1/FVC. The proteins with the greatest inverse relations to FEV1%predicted and FVC%predicted included leptin, adrenomedullin, and plasminogen activator inhibitor-1; in contrast there were three proteins with positive relations to FEV1%predicted and FVC%predicted including insulin growth factor binding protein 2, tetranectin, and soluble receptor for advanced glycation end products. In longitudinal analyses, three proteins were associated with longitudinal change in FEV1 (ΔFEV1) and four with ΔFVC; no proteins were associated with ΔFEV1/FVC., Conclusion: Our findings highlight CVD-related plasma proteins that are associated with lung function including markers of inflammation, adiposity, and fibrosis, representing proteins that may contribute both to respiratory and CVD risk., Competing Interests: The authors have read the journal’s policy and have the following competing interests: Jennifer Ho received research support from Gilead Sciences and Bayer AG, and research supplies from EcoNugenics. The National Heart, Lung, and Blood Institute and Ionis Pharmaceuticals entered into a Cooperative Research and Development Agreement (CRADA) to conduct research targeting the AGER gene in relation to lung disease. Daniel Levy is the NHLBI principal investigator on the CRADA. Neither Daniel Levy nor the NHLBI is receiving any funding from Ionis Pharmaceuticals in conjunction with the CRADA. There are no patents, products in development or marketed products associated with this research to declare. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Background: DNA methylation-based epigenetic age measures have been used as biological aging markers and are associated with a healthy lifespan. Few population-based studies have examined the relation between diet and epigenetic age acceleration., Objectives: We aimed to investigate the relation between diet quality and epigenetic age acceleration., Methods: We analyzed data from 1995 participants (mean age, 67 years; 55% women) of the Framingham Heart Study Offspring Cohort. Cross-sectional associations between the Dietary Approaches to Stop Hypertension (DASH) score and 3 whole-blood DNA methylation-derived epigenetic age acceleration measures-Dunedin Pace of Aging Methylation (DunedinPoAm), GrimAge acceleration (GrimAA), and PhenoAge acceleration (PhenoAA)-were examined. A mediation analysis was conducted to assess the mediating role of epigenetic age acceleration in relation to DASH and all-cause mortality., Results: A higher DASH score was associated with lower levels of DunedinPoAm (β = -0.05; SE = 0.02; P = 0.007), GrimAA (β = -0.09; SE = 0.02; P < 0.001), and PhenoAA (β = -0.07; SE = 0.02; P = 0.001). All 3 epigenetic measures mediated the association between the DASH score and all-cause mortality, with mean proportions of 22.1% for DunedinPoAm (Pmediation = 0.04), 45.1% for GrimAA (Pmediation = 0.001), and 22.9% for PhenoAA (Pmediation = 0.03). An interaction was observed between the DASH score and smoking status in relation to the epigenetic aging markers. The association between the DASH score and epigenetic aging markers tended to be stronger in "ever-smokers" (former and current smokers) compared to "never-smokers." The proportions of mediation were 31.3% for DunedinPoAm, 46.8% for GrimAA, and 10.3% for PhenoAA in ever-smokers, whereas no significant mediation was observed in never-smokers., Conclusions: Higher diet quality is associated with slower epigenetic age acceleration, which partially explains the beneficial effect of diet quality on the lifespan. Our findings emphasize that adopting a healthy diet is crucial for maintaining healthy aging., (© The Author(s) 2021. Published by Oxford University Press on behalf of the American Society for Nutrition.)

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Chronic kidney disease is a major public health burden. Elevated urinary albumin-to-creatinine ratio is a measure of kidney damage, and used to diagnose and stage chronic kidney disease. To extend the knowledge on regulatory mechanisms related to kidney function and disease, we conducted a blood-based epigenome-wide association study for estimated glomerular filtration rate (n = 33,605) and urinary albumin-to-creatinine ratio (n = 15,068) and detected 69 and seven CpG sites where DNA methylation was associated with the respective trait. The majority of these findings showed directionally consistent associations with the respective clinical outcomes chronic kidney disease and moderately increased albuminuria. Associations of DNA methylation with kidney function, such as CpGs at JAZF1, PELI1 and CHD2 were validated in kidney tissue. Methylation at PHRF1, LDB2, CSRNP1 and IRF5 indicated causal effects on kidney function. Enrichment analyses revealed pathways related to hemostasis and blood cell migration for estimated glomerular filtration rate, and immune cell activation and response for urinary albumin-to-creatinineratio-associated CpGs., (© 2021. The Author(s).)

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits., (© 2021. The Author(s).)

Common carotid intima-media thickness (cIMT) is an index of subclinical atherosclerosis that is associated with ischemic stroke and coronary artery disease (CAD). We undertook a cross-sectional epigenome-wide association study (EWAS) of measures of cIMT in 6400 individuals. Mendelian randomization analysis was applied to investigate the potential causal role of DNA methylation in the link between atherosclerotic cardiovascular risk factors and cIMT or clinical cardiovascular disease. The CpG site cg05575921 was associated with cIMT (beta = -0.0264, p value = 3.5 × 10 -8 ) in the discovery panel and was replicated in replication panel (beta = -0.07, p value = 0.005). This CpG is located at chr5:81649347 in the intron 3 of the aryl hydrocarbon receptor repressor gene (AHRR). Our results indicate that DNA methylation at cg05575921 might be in the pathway between smoking, cIMT and stroke. Moreover, in a region-based analysis, 34 differentially methylated regions (DMRs) were identified of which a DMR upstream of ALOX12 showed the strongest association with cIMT (p value = 1.4 × 10 -13 ). In conclusion, our study suggests that DNA methylation may play a role in the link between cardiovascular risk factors, cIMT and clinical cardiovascular disease., (© 2021. The Author(s).)

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Multiple methods have been proposed to aggregate genetic variants in a gene or a region and jointly test their association with a trait of interest. However, these joint tests do not provide estimates of the individual effect of each variant. Moreover, few methods have evaluated the joint association of multiple variants with DNA methylation. We propose a method based on linear mixed models to estimate the joint and individual effect of multiple genetic variants on DNA methylation leveraging genomic annotations. Our approach is flexible, can incorporate covariates and annotation features, and takes into account relatedness and linkage disequilibrium (LD). Our method had correct Type-I error and overall high power for different simulated scenarios where we varied the number and specificity of functional annotations, number of causal and total genetic variants, frequency of genetic variants, LD, and genetic variant effect. Our method outperformed the family Sequence Kernel Association Test and had more stable estimations of effects than a classical single-variant linear mixed-effect model. Applied genome-wide to the Framingham Heart Study data, our method identified 921 DNA methylation sites influenced by at least one rare or low-frequency genetic variant located within 50 kilobases (kb) of the DNA methylation site., (© 2020 Wiley Periodicals LLC.)

Background: DNA methylation is a key epigenetic modification that can directly affect gene regulation. DNA methylation is highly influenced by environmental factors such as cigarette smoking, which is causally related to chronic obstructive pulmonary disease (COPD) and lung cancer. To date, there have been few large-scale, combined analyses of DNA methylation and gene expression and their interrelations with lung diseases., Results: We performed an epigenome-wide association study of whole blood gene expression in ~ 6000 individuals from four cohorts. We discovered and replicated numerous CpGs associated with the expression of cis genes within 500 kb of each CpG, with 148 to 1,741 cis CpG-transcript pairs identified across cohorts. We found that the closer a CpG resided to a transcription start site, the larger its effect size, and that 36% of cis CpG-transcript pairs share the same causal genetic variant. Mendelian randomization analyses revealed that hypomethylation and lower expression of CHRNA5, which encodes a smoking-related nicotinic receptor, are causally linked to increased risk of COPD and lung cancer. This putatively causal relationship was further validated in lung tissue data., Conclusions: Our results provide a large and comprehensive association study of whole blood DNA methylation with gene expression. Expression platform differences rather than population differences are critical to the replication of cis CpG-transcript pairs. The low reproducibility of trans CpG-transcript pairs suggests that DNA methylation regulates nearby rather than remote gene expression. The putatively causal roles of methylation and expression of CHRNA5 in relation to COPD and lung cancer provide evidence for a mechanistic link between patterns of smoking-related epigenetic variation and lung diseases, and highlight potential therapeutic targets for lung diseases and smoking cessation.

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-β-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

Background: DNA methylation patterns associated with habitual diet have not been well studied., Methods: Diet quality was characterized using a Mediterranean-style diet score and the Alternative Healthy Eating Index score. We conducted ethnicity-specific and trans-ethnic epigenome-wide association analyses for diet quality and leukocyte-derived DNA methylation at over 400 000 CpGs (cytosine-guanine dinucleotides) in 5 population-based cohorts including 6662 European ancestry, 2702 African ancestry, and 360 Hispanic ancestry participants. For diet-associated CpGs identified in epigenome-wide analyses, we conducted Mendelian randomization (MR) analysis to examine their relations to cardiovascular disease risk factors and examined their longitudinal associations with all-cause mortality., Results: We identified 30 CpGs associated with either Mediterranean-style diet score or Alternative Healthy Eating Index, or both, in European ancestry participants. Among these CpGs, 12 CpGs were significantly associated with all-cause mortality (Bonferroni corrected P <1.6×10 -3 ). Hypermethylation of cg18181703 ( SOCS3 ) was associated with higher scores of both Mediterranean-style diet score and Alternative Healthy Eating Index and lower risk for all-cause mortality ( P =5.7×10 -15 ). Ten additional diet-associated CpGs were nominally associated with all-cause mortality ( P <0.05). MR analysis revealed 8 putatively causal associations for 6 CpGs with 4 cardiovascular disease risk factors (body mass index, triglycerides, high-density lipoprotein cholesterol concentrations, and type 2 diabetes mellitus; Bonferroni corrected MR P <4.5×10 -4 ). For example, hypermethylation of cg11250194 ( FADS2 ) was associated with lower triglyceride concentrations (MR, P =1.5×10 -14 ).and hypermethylation of cg02079413 ( SNORA54 ; NAP1L4 ) was associated with body mass index (corrected MR, P =1×10 -6 )., Conclusions: Habitual diet quality was associated with differential peripheral leukocyte DNA methylation levels of 30 CpGs, most of which were also associated with multiple health outcomes, in European ancestry individuals. These findings demonstrate that integrative genomic analysis of dietary information may reveal molecular targets for disease prevention and treatment.

Background: Sleep Disordered Breathing (SDB) is associated with a wide range of pathophysiological changes due, in part, to hypoxemia during sleep. We sought to identify gene transcription associations with measures of SDB and hypoxemia during sleep, and study their response to treatment., Methods: In two discovery cohorts, Framingham Offspring Study (FOS; N = 571) and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 580), we studied gene expression in peripheral blood mononuclear cells in association with three measures of SDB: Apnea Hypopnea Index (AHI); average oxyhemoglobin saturation (avgO2) during sleep; and minimum oxyhemoglobin saturation (minO2) during sleep. Associated genes were used for analysis of gene expression in the blood of 15 participants with moderate or severe obstructive sleep apnea (OSA) from the Heart Biomarkers In Apnea Treatment (HeartBEAT) trial. These genes were studied pre- and post-treatment (three months) with continuous positive airway pressure (CPAP). We also performed Gene Set Enrichment Analysis (GSEA) on all traits and cohort analyses., Findings: Twenty-two genes were associated with SDB traits in both MESA and FOS. Of these, lower expression of CD1D and RAB20 was associated with lower avgO2 in MESA and FOS. CPAP treatment increased the expression of these genes in HeartBEAT participants. Immunity and inflammation pathways were up-regulated in subjects with lower avgO2; i.e., in those with a more severe SDB phenotype (MESA), whereas immuno-inflammatory processes were down-regulated following CPAP treatment (HeartBEAT)., Interpretation: Low oxygen saturation during sleep is associated with alterations in gene expression and transcriptional programs that are partially reversed by CPAP treatment., Competing Interests: Declarations of Competing Interest Dr. Rotter reports grants from NIH, during the conduct of this study. Dr. Redline reports grants from NIH,  during the conduct of the study. Dr. Patel reports grants from Bayer Pharmaceuticals, Philips Respironics, and Respicardia outside the submitted work. Dr. Mehra reports receiving National Institutes of Health funding support from the National Heart, Lung, and Blood Institute [U01HL125177, UG3HL140144] and the American Heart Association. Her institution has received positive airway pressure devices and equipment from Philips Respironics, ResMed, GE Healthcare, and Natus for research. R. M serves as a consultant for Respicardia, Enhale and Merck; received funds for service on the American Board of Medicine Sleep Medicine Exam test writing committee and received royalties from UpToDate. Dr. Gottlieb reports grants from NIH/NHLBI,  during the conduct of the study. Dr. Quan reports personal fees from Jazz Pharmaceuticals, personal fees from Best Doctors, other from Sleep Research Society, personal fees from American Academy of Sleep Medicine,  outside the submitted work., (Copyright © 2020 The Author(s). Published by Elsevier B.V. All rights reserved.)

DNA methylation (DNAm) and microRNAs (miRNAs) have been implicated in a wide-range of human diseases. While often studied in isolation, DNAm and miRNAs are not independent. We analyzed associations of expression of 283 miRNAs with DNAm at >400K CpG sites in whole blood obtained from 3565 individuals and identified 227 CpGs at which differential methylation was associated with the expression of 40 nearby miRNAs ( cis -miR-eQTMs) at FDR<0.01, including 91 independent CpG sites at r 2 < 0.2. cis- miR-eQTMs were enriched for CpGs in promoter and polycomb-repressed state regions, and 60% were inversely associated with miRNA expression. Bidirectional Mendelian randomization (MR) analysis further identified 58 cis -miR-eQTMCpG-miRNA pairs where DNAm changes appeared to drive miRNA expression changes and opposite directional effects were unlikely. Integration of genetic variants in joint analyses revealed an average partial between cis- miR-eQTM CpGs and miRNAs of 2% after conditioning on site-specific genetic variation, suggesting that DNAm is an important epigenetic regulator of miRNA expression. Finally, two-step MR analysis was performed to identify putatively causal CpGs driving miRNA expression in relation to human complex traits. We found that an imprinted region on 14q32 that was previously identified in relation to age at menarche is enriched with cis- miR-eQTMs. Nine CpGs and three miRNAs at this locus tested causal for age at menarche, reflecting novel epigenetic-driven molecular pathways underlying this complex trait. Our study sheds light on the joint genetic and epigenetic regulation of miRNA expression and provides insights into the relations of miRNAs to their targets and to complex phenotypes.