Association analysis between an epigenetic risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Association analysis between an epigenetic alcohol risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Expression quantitative trait methylation analysis elucidates gene regulatory effects of DNA methylation: the Framingham Heart Study.

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.

Epigenome-wide DNA methylation association study of circulating IgE levels identifies novel targets for asthma.

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.

Whole genome DNA and RNA sequencing of whole blood elucidates the genetic architecture of gene expression underlying a wide range of diseases.

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.

Elucidating the genetic architecture of DNA methylation to identify promising molecular mechanisms of disease.

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.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases.

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 < 5x10 - 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 7,233 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.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases.

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 7,233 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.

A genomic approach identifies sRAGE as a putatively causal protein for asthma.

BACKGROUND: Asthma is a complex respiratory condition caused by environmental and genetic factors. Although lower concentrations of the anti-inflammatory protein soluble receptor for advanced glycation end products (sRAGE) have been associated with asthma in humans and mouse models, it is uncertain whether sRAGE plays a causal role in asthma. OBJECTIVE: We designed a 2-stage study of sRAGE in relation to asthma with association analysis in FHS participants as well as causal inference testing using Mendelian randomization (MR). METHODS: We measured plasma levels of sRAGE and performed cross-sectional analysis to examine the association between plasma sRAGE concentration and asthma status in 6546 FHS participants. We then used sRAGE protein advanced glycation end products (pQTLs) derived from a genome-wide association study of plasma sRAGE levels in ∼7000 FHS participants with UK Biobank asthma genome-wide association study in MR to consider sRAGE as a putatively causal protein for asthma. We also performed replication MR using an externally derived sRAGE pQTL from the INTERVAL study. Last, we conducted colocalization using cis-pQTL variants at the advanced glycosylation end-product specific receptor (AGER) locus with variants from the UK Biobank asthma genome-wide association study. RESULTS: Association analysis revealed that each 1 SD increment in sRAGE concentration was associated with a 14% lower odds of asthma in FHS participants (95% CI 0.76-0.96). MR identified sRAGE as putatively causal for and protective against asthma on the basis of self-reported (odds ratio [per 1 SE increment in inverse-rank-normalized sRAGE] = 0.97, 95% CI 0.95-0.99; P = .005) and doctor-diagnosed asthma (odds ratio = 0.97, 95% CI 0.95-0.99; P = .011). CONCLUSION: Through this genomic approach, we identified sRAGE as a putatively causal, biologically important, and protective protein in relation to asthma. Functional studies in cell/animal models are needed to confirm our findings.

Redox regulation of protein kinase C by selenometabolites and selenoprotein thioredoxin reductase limits cancer prevention by selenium.

The cancer-preventive mechanism of selenium should address the way low concentrations of selenometabolites react with cellular targets without being diffused from the sites of generation, the way selenium selectively kills tumor cells, and the intriguing U-shaped curve that is seen with dietary supplementation of selenium and cancer prevention. Protein kinase C (PKC), a receptor for tumor promoters, is well suited for this mechanism. Due to the catalytic redox cycle, low concentrations of methylselenol, a postulated active metabolite of selenium, react with the tumor-promoting lipid hydroperoxide bound to PKC to form methylseleninic acid (MSA), which selectively reacts with thiol residues present within the vicinity of the PKC catalytic domain to inactivate it. Given that lipid hydroperoxide levels are high in promoting cells, PKC inactivation selectively leads to death in these cells. A biphasic effect of MSA in inducing cell death was observed in certain prostate cancer cell lines; lower concentrations of MSA induced cell death, while higher concentrations failed to do so. Lower concentrations of selenium inactivate more sensitive antiapoptotic isoenzymes of PKC (ε and α), sparing less sensitive proapoptotic isoenzymes (PKCδ and PKCζ). Higher concentrations of selenium also inactivate proapoptotic isoenzymes and consequently make tumor cells resistant to apoptosis. Due to a high-affinity binding of thioredoxin to the PKC catalytic domain, this thiol oxidation is explicitly reversed by thioredoxin reductase (TXNRD), a selenoprotein. Therefore, overexpression of TXNRD in advanced tumor cells could make them resistant to selenium-induced death. Conceivably, this mechanism, at least in part, explains why selenium prevents cancer only in certain cases.

Laminin-1 induces endocytosis of 67KDa laminin receptor and protects Neuroscreen-1 cells against death induced by serum withdrawal.

Although the function of laminin in the basement membrane is known, the function of soluble "neuronal" laminin is unknown. Since laminin is neuroprotective, we determined whether the soluble laminin-1 induces signaling for neuroprotection via its 67KDa laminin-1 receptor (67LR). Treatment of Neuroscreen-1 (NS-1) cells with laminin-1 or YIGSR peptide, which corresponds to a sequence in laminin-1 ß1 chain that binds to 67LR, induced a decrease in the cell-surface expression of 67LR and caused its internalization. Furthermore, intracellular cAMP-elevating agents, dibutyryl-cAMP, forskolin, and rolipram, also induced this internalization. Both soluble laminin-1 and YIGSR induced a sustained elevation of intracellular cAMP under defined conditions, suggesting a causal role of cAMP in the endocytosis of 67LR. This endocytosis was not observed in cells deficient in protein kinase A (PKA) nor in cells treated with either SQ 22536, an inhibitor for adenylyl cyclase, or ESI-09, an inhibitor for the exchange protein directly activated by cAMP (Epac). In addition, when internalization occurred in NS-1 cells, 67LR and adenylyl cyclase were localized in early endosomes. Under conditions in which endocytosis had occurred, both laminin-1 and YIGSR protected NS-1 cells from cell death induced by serum withdrawal. However, under conditions in which endocytosis did not occur, neither laminin-1 nor YIGSR protected these cells. Conceivably, the binding of laminin-1 to 67LR causes initial signaling through PKA and Epac, which causes the internalization of 67LR, along with signaling enzymes, such as adenylyl cyclase, into early endosomes. This causes sustained signaling for protection against cell death induced by serum withdrawal.