Genome-wide identification of DNA methylation QTLs in whole blood highlights pathways for cardiovascular disease

Detecting methylation quantitative trait loci using a methylation random field method

DNA methylation may be regulated by genetic variants within a genomic region, referred to as methylation quantitative trait loci (mQTLs). The changes of methylation levels can further lead to alterations of gene expression, and influence the risk of various complex human diseases. Detecting mQTLs may provide insights into the underlying mechanism of how genotypic variations may influence the disease risk. In this article, we propose a methylation random field (MRF) method to detect mQTLs by testing the association between the methylation level of a CpG site and a set of genetic variants within a genomic region. The proposed MRF has two major advantages over existing approaches. First, it uses a beta distribution to characterize the bimodal and interval properties of the methylation trait at a CpG site. Second, it considers multiple common and rare genetic variants within a genomic region to identify mQTLs. Through simulations, we demonstrated that the MRF had improved power over other existing methods in detecting rare variants of relatively large effect, especially when the sample size is small. We further applied our method to a study of congenital heart defects with 83 cardiac tissue samples and identified two mQTL regions, MRPS10 and PSORS1C1, which were colocalized with expression QTL in cardiac tissue. In conclusion, the proposed MRF is a useful tool to identify novel mQTLs, especially for studies with limited sample sizes.

Integrating Genome and Methylome Data to Identify Candidate DNA Methylation Biomarkers for Pancreatic Cancer Risk

BACKGROUND

The role of methylation in pancreatic cancer risk remains unclear. We integrated genome and methylome data to identify CpG sites (CpG) with the genetically predicted methylation to be associated with pancreatic cancer risk. We also studied gene expression to understand the identified associations.

METHODS

Using genetic data and white blood cell methylation data from 1,595 subjects of European descent, we built genetic models to predict DNA methylation levels. After internal and external validation, we applied prediction models with satisfactory performance to the genetic data of 8,280 pancreatic cancer cases and 6,728 controls of European ancestry to investigate the associations of predicted methylation with pancreatic cancer risk. For associated CpGs, we compared their measured levels in pancreatic tumor versus benign tissue.

RESULTS

We identified 45 CpGs at nine loci showing an association with pancreatic cancer risk, including 15 CpGs showing an association independent from identified risk variants. We observed significant correlations between predicted methylation of 16 of the 45 CpGs and predicted expression of eight adjacent genes, of which six genes showed associations with pancreatic cancer risk. Of the 45 CpGs, we were able to compare measured methylation of 16 in pancreatic tumor versus benign pancreatic tissue. Of them, six showed differentiated methylation.

CONCLUSIONS

We identified methylation biomarker candidates associated with pancreatic cancer using genetic instruments and added additional insights into the role of methylation in regulating gene expression in pancreatic cancer development.

IMPACT

A comprehensive study using genetic instruments identifies 45 CpG sites at nine genomic loci for pancreatic cancer risk.

Meta-analysis of epigenome-wide association studies of carotid intima-media thickness

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.

DNA methylation signatures of incident coronary heart disease: findings from epigenome-wide association studies

Coronary heart disease (CHD) is a type of cardiovascular disease (CVD) that affects the coronary arteries, which provide oxygenated blood to the heart. It is a major cause of mortality worldwide. Various prediction methods have been developed to assess the likelihood of developing CHD, including those based on clinical features and genetic variation. Recent epigenome-wide studies have identified DNA methylation signatures associated with the development of CHD, indicating that DNA methylation may play a role in predicting future CHD. This narrative review summarises recent findings from DNA methylation studies of incident CHD (iCHD) events from epigenome-wide association studies (EWASs). The results suggest that DNA methylation signatures may identify new mechanisms involved in CHD progression and could prove a useful adjunct for the prediction of future CHD.

Genome-wide sequencing-based identification of methylation quantitative trait loci and their role in schizophrenia risk

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.

The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Connecting the epigenome, metabolome and proteome for a deeper understanding of disease

Epigenome-wide association studies (EWAS) identify genes that are dysregulated by the studied clinical endpoints, thereby indicating potential new diagnostic biomarkers, drug targets and therapy options. Combining EWAS with deep molecular phenotyping, such as approaches enabled by metabolomics and proteomics, allows further probing of the underlying disease-associated pathways. For instance, methylation of the TXNIP gene is associated robustly with prevalent type 2 diabetes and further with metabolites that are short-term markers of glycaemic control. These associations reflect TXNIP's function as a glucose uptake regulator by interaction with the major glucose transporter GLUT1 and suggest that TXNIP methylation can be used as a read-out for the organism's exposure to glucose stress. Another case is the association between DNA methylation of the AHRR and F2RL3 genes with smoking and a protein that is involved in the reprogramming of the bronchial epithelium. These examples show that associations between DNA methylation and intermediate molecular traits can open new windows into how the body copes with physiological challenges. This knowledge, if carefully interpreted, may indicate novel therapy options and, together with monitoring of the methylation state of specific methylation sites, may in the future allow the early diagnosis of impending disease. It is essential for medical practitioners to recognize the potential that this field holds in translating basic research findings to clinical practice. In this review, we present recent advances in the field of EWAS with metabolomics and proteomics and discuss both the potential and the challenges of translating epigenetic associations, with deep molecular phenotypes, to biomedical applications.

Novel DNA methylation loci and genes showing pleiotropic association with Alzheimer's dementia: a network Mendelian randomization analysis

Previous genome-wide association studies (GWAS) have identified potential genetic variants involved in the risk of Alzheimer's dementia, but their underlying biological interpretation remains largely unclear. In addition, the effects of DNA methylation and gene expression on Alzheimer's dementia are not well understood. A network summary data-based Mendelian randomization (SMR) analysis was performed integrating cis- DNA methylation quantitative trait loci (mQTL) /cis- gene expression QTL (eQTL) data in the brain and blood, as well as GWAS summarized data for Alzheimer's dementia to evaluate the pleiotropic associations of DNA methylation and gene expression with Alzheimer's dementia and to explore the complex mechanisms underpinning Alzheimer's dementia. After correction for multiple testing (false discovery rate [FDR] P < 0.05) and filtering using the heterogeneity in dependent instruments (HEIDI) test (P_HEIDI>0.01), we identified dozens of DNA methylation sites and genes showing pleiotropic associations with Alzheimer's dementia. We found 22 and 16 potentially causal pathways of Alzheimer's dementia (i.e., SNP→DNA methylation→Gene expression→Alzheimer's dementia) in the brain and blood, respectively. Approximately two-thirds of the identified DNA methylation sites had an influence on gene expression and the expression of almost all the identified genes was regulated by DNA methylation. Our network SMR analysis provided evidence supporting the pleiotropic association of some novel DNA methylation sites and genes with Alzheimer's dementia and revealed possible causal pathways underlying the pathogenesis of Alzheimer's dementia. Our findings shed light on the role of DNA methylation in gene expression and in the development of Alzheimer's dementia.

Thousands of CpGs Show DNA Methylation Differences in ACPA-Positive Individuals

High levels of anti-citrullinated protein antibodies (ACPA) are often observed prior to a diagnosis of rheumatoid arthritis (RA). We undertook a replication study to confirm CpG sites showing evidence of differential methylation in subjects positive vs. negative for ACPA, in a new subset of 112 individuals sampled from the population cohort and biobank CARTaGENE in Quebec, Canada. Targeted custom capture bisulfite sequencing was conducted at approximately 5.3 million CpGs located in regulatory or hypomethylated regions from whole blood; library and protocol improvements had been instituted between the original and this replication study, enabling better coverage and additional identification of differentially methylated regions (DMRs). Using binomial regression models, we identified 19,472 ACPA-associated differentially methylated cytosines (DMCs), of which 430 overlapped with the 1909 DMCs reported by the original study; 814 DMRs of relevance were clustered by grouping adjacent DMCs into regions. Furthermore, we performed an additional integrative analysis by looking at the DMRs that overlap with RA related loci published in the GWAS Catalog, and protein-coding genes associated with these DMRs were enriched in the biological process of cell adhesion and involved in immune-related pathways.

Challenges and Opportunities in Understanding Genetics of Fungal Diseases: Towards a Functional Genomics Approach

Infectious diseases are a leading cause of morbidity and mortality worldwide, and human pathogens have long been recognized as one of the main sources of evolutionary pressure, resulting in a high variable genetic background in immune-related genes. The study of the genetic contribution to infectious diseases has undergone tremendous advances over the last decades. Here, focusing on genetic predisposition to fungal diseases, we provide an overview of the available approaches for studying human genetic susceptibility to infections, reviewing current methodological and practical limitations. We describe how the classical methods available, such as family-based studies and candidate gene studies, have contributed to the discovery of crucial susceptibility factors for fungal infections. We will also discuss the contribution of novel unbiased approaches to the field, highlighting their success but also their limitations for the fungal immunology field. Finally, we show how a systems genomics approach can overcome those limitations and can lead to efficient prioritization and identification of genes and pathways with a critical role in susceptibility to fungal diseases. This knowledge will help to stratify at-risk patient groups and, subsequently, develop early appropriate prophylactic and treatment strategies.

Multi-omics analysis identifies CpGs near G6PC2 mediating the effects of genetic variants on fasting glucose

AIMS/HYPOTHESIS

An elevated fasting glucose level in non-diabetic individuals is a key predictor of type 2 diabetes. Genome-wide association studies (GWAS) have identified hundreds of SNPs for fasting glucose but most of their functional roles in influencing the trait are unclear. This study aimed to identify the mediation effects of DNA methylation between SNPs identified as significant from GWAS and fasting glucose using Mendelian randomisation (MR) analyses.

METHODS

We first performed GWAS analyses for three cohorts (Taiwan Biobank with 18,122 individuals, the Healthy Aging Longitudinal Study in Taiwan with 1989 individuals and the Stanford Asia-Pacific Program for Hypertension and Insulin Resistance with 416 individuals) with individuals of Han Chinese ancestry in Taiwan, followed by a meta-analysis for combining the three GWAS analysis results to identify significant and independent SNPs for fasting glucose. We determined whether these SNPs were methylation quantitative trait loci (meQTLs) by testing their associations with DNA methylation levels at nearby CpG sites using a subsample of 1775 individuals from the Taiwan Biobank. The MR analysis was performed to identify DNA methylation with causal effects on fasting glucose using meQTLs as instrumental variables based on the 1775 individuals. We also used a two-sample MR strategy to perform replication analysis for CpG sites with significant MR effects based on literature data.

RESULTS

Our meta-analysis identified 18 significant (p < 5 × 10-8) and independent SNPs for fasting glucose. Interestingly, all 18 SNPs were meQTLs. The MR analysis identified seven CpGs near the G6PC2 gene that mediated the effects of a significant SNP (rs2232326) in the gene on fasting glucose. The MR effects for two CpGs were replicated using summary data based on the European population, using an exonic SNP rs2232328 in G6PC2 as the instrument.

CONCLUSIONS/INTERPRETATION

Our analysis results suggest that rs2232326 and rs2232328 in G6PC2 may affect DNA methylation at CpGs near the gene and that the methylation may have downstream effects on fasting glucose. Therefore, SNPs in G6PC2 and CpGs near G6PC2 may reside along the pathway that influences fasting glucose levels. This is the first study to report CpGs near G6PC2, an important gene for regulating insulin secretion, mediating the effects of GWAS-significant SNPs on fasting glucose.

Genetic drivers of m6A methylation in human brain, lung, heart and muscle

The most prevalent post-transcriptional mRNA modification, N⁶-methyladenosine (m⁶A), plays diverse RNA-regulatory roles, but its genetic control in human tissues remains uncharted. Here we report 129 transcriptome-wide m⁶A profiles, covering 91 individuals and 4 tissues (brain, lung, muscle and heart) from GTEx/eGTEx. We integrate these with interindividual genetic and expression variation, revealing 8,843 tissue-specific and 469 tissue-shared m⁶A quantitative trait loci (QTLs), which are modestly enriched in, but mostly orthogonal to, expression QTLs. We integrate m⁶A QTLs with disease genetics, identifying 184 GWAS-colocalized m⁶A QTL, including brain m⁶A QTLs underlying neuroticism, depression, schizophrenia and anxiety; lung m⁶A QTLs underlying expiratory flow and asthma; and muscle/heart m⁶A QTLs underlying coronary artery disease. Last, we predict novel m⁶A regulators that show preferential binding in m⁶A QTLs, protein interactions with known m⁶A regulators and expression correlation with the m⁶A levels of their targets. Our results provide important insights and resources for understanding both cis and trans regulation of epitranscriptomic modifications, their interindividual variation and their roles in human disease.

The epigenetic etiology of cardiovascular disease in a longitudinal Swedish twin study

BACKGROUND

Studies on DNA methylation have the potential to discover mechanisms of cardiovascular disease (CVD) risk. However, the role of DNA methylation in CVD etiology remains unclear.

RESULTS

We performed an epigenome-wide association study (EWAS) on CVD in a longitudinal sample of Swedish twins (535 individuals). We selected CpGs reaching the Bonferroni-corrected significance level (2 [Formula: see text] 10-7) or the top-ranked 20 CpGs with the lowest P values if they did not reach this significance level in EWAS analysis associated with non-stroke CVD, overall stroke, and ischemic stroke, respectively. We further applied a bivariate autoregressive latent trajectory model with structured residuals (ALT-SR) to evaluate the cross-lagged effect between DNA methylation of these CpGs and cardiometabolic traits (blood lipids, blood pressure, and body mass index). Furthermore, mediation analysis was performed to evaluate whether the cross-lagged effects had causal impacts on CVD. In the EWAS models, none of the CpGs we selected reached the Bonferroni-corrected significance level. The ALT-SR model showed that DNA methylation levels were more likely to predict the subsequent level of cardiometabolic traits rather than the other way around (numbers of significant cross-lagged paths of methylation → trait/trait → methylation were 84/4, 45/6, 66/1 for the identified three CpG sets, respectively). Finally, we demonstrated significant indirect effects from DNA methylation on CVD mediated by cardiometabolic traits.

CONCLUSIONS

We present evidence for a directional association from DNA methylation on cardiometabolic traits and CVD, rather than the opposite, highlighting the role of epigenetics in CVD development.

Higher diet quality relates to decelerated epigenetic aging

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.

Multi-Omics Approaches in Immunological Research

The immune system plays a vital role in health and disease, and is regulated through a complex interactive network of many different immune cells and mediators. To understand the complexity of the immune system, we propose to apply a multi-omics approach in immunological research. This review provides a complete overview of available methodological approaches for the different omics data layers relevant for immunological research, including genetics, epigenetics, transcriptomics, proteomics, metabolomics, and cellomics. Thereafter, we describe the various methods for data analysis as well as how to integrate different layers of omics data. Finally, we discuss the possible applications of multi-omics studies and opportunities they provide for understanding the complex regulatory networks as well as immune variation in various immune-related diseases.

Framingham Heart Study: JACC Focus Seminar, 1/8

The Framingham Heart Study is the longest-running cardiovascular epidemiological study, starting in 1948. This paper gives an overview of the various cohorts, collected data, and most important research findings to date. In brief, the Framingham Heart Study, funded by the National Institutes of Health and managed by Boston University, spans 3 generations of well phenotyped White persons and 2 cohorts comprised of racial and ethnic minority groups. These cohorts are densely phenotyped, with extensive longitudinal follow-up, and they continue to provide us with important information on human cardiovascular and noncardiovascular physiology over the lifespan, as well as to identify major risk factors for cardiovascular disease. This paper also summarizes some of the more recent progress in molecular epidemiology and discusses the future of the study.

Statistical methods for mediation analysis in the era of high-throughput genomics: Current successes and future challenges

Mediation analysis investigates the intermediate mechanism through which an exposure exerts its influence on the outcome of interest. Mediation analysis is becoming increasingly popular in high-throughput genomics studies where a common goal is to identify molecular-level traits, such as gene expression or methylation, which actively mediate the genetic or environmental effects on the outcome. Mediation analysis in genomics studies is particularly challenging, however, thanks to the large number of potential mediators measured in these studies as well as the composite null nature of the mediation effect hypothesis. Indeed, while the standard univariate and multivariate mediation methods have been well-established for analyzing one or multiple mediators, they are not well-suited for genomics studies with a large number of mediators and often yield conservative p-values and limited power. Consequently, over the past few years many new high-dimensional mediation methods have been developed for analyzing the large number of potential mediators collected in high-throughput genomics studies. In this work, we present a thorough review of these important recent methodological advances in high-dimensional mediation analysis. Specifically, we describe in detail more than ten high-dimensional mediation methods, focusing on their motivations, basic modeling ideas, specific modeling assumptions, practical successes, methodological limitations, as well as future directions. We hope our review will serve as a useful guidance for statisticians and computational biologists who develop methods of high-dimensional mediation analysis as well as for analysts who apply mediation methods to high-throughput genomics studies.

Epigenome-wide scan identifies differentially methylated regions for lung cancer using pre-diagnostic peripheral blood

BACKGROUND

DNA methylation markers have been associated with lung cancer risk and may identify aetiologically relevant genomic regions, or alternatively, be markers of disease risk factors or biological processes associated with disease development.

METHODS

In a nested case-control study, we measured blood leukocyte DNA methylation levels in pre-diagnostic samples collected from 430 participants (208 cases; 222 controls) in the 1989 CLUE II cohort. We compared DNA methylation levels with case/control status to identify novel genomic regions, both single CpG sites and differentially methylated regions (DMRs), while controlling for known DNA methylation changes associated with smoking using a previously described pack-years-based smoking methylation score. Stratification analyses were conducted over time from blood draw to diagnosis, histology, and smoking status.

RESULTS

We identified 16 single CpG sites and 40 DMRs significantly associated with lung cancer risk (q < 0.05). The identified genomic regions were associated with genes including H19, HOXA3/HOXA4, RUNX3, BRICD5, PLXNB2, and RP13. For the single CpG sites, the strongest association was noted for cg09736286 in the DIABLO gene (OR [for 1 SD] = 2.99, 95% CI: 1.95-4.59, P-value = 4.81 × 10-7). We found that CpG sites in the HOXA3/HOXA4 region were hypermethylated in cases compared to controls.

CONCLUSION

The single CpG sites and DMRs that we identified represented significant measurable differences in lung cancer risk, providing potential biomarkers for lung cancer risk stratification. Future studies will need to examine whether these regions are causally related to lung cancer.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

Distinct Associations of BMI and Fatty Acids With DNA Methylation in Fasting and Postprandial States in Men

We have previously shown that blood global DNA methylation (DNAm) differs between postprandial state (PS) and fasting state (FS) and is associated with BMI and polyunsaturated fatty acid (PUFA) (negatively and positively, respectively) in 12 metabolically healthy adult Mexican men (AMM cohort) equally distributed among conventional BMI classes. Here, we detailed those associations at CpG dinucleotide level by exploiting the Infinium methylation EPIC array (Illumina). We sought differentially methylated CpG (dmCpG) that were (1) associated with BMI (BMI-dmCpG) and/or fatty acids (FA) (FA-dmCpG) in FS or PS and (2) different across FS and PS within a BMI class. BMI-dmCpG and FA-dmCpG were more numerous in FS compared to PS and largely prandial state-specific. For saturated and monounsaturated FA, dmCpG overlap was higher across than within the respective saturation group. Several BMI- and FA-dmCpG mapped to genes involved in metabolic disease and in some cases matched published experimental data sets. Notably, SETDB1 and MTHFS promoter dmCpG could explain the previously observed associations between global DNAm, PUFA content, and BMI in FS. Surprisingly, overlap between BMI-dmCpG and FA-dmCpG was limited and the respective dmCpG were differentially distributed across functional genomic elements. BMI-dmCpG showed the highest overlap with dmCpG of the saturated FA palmitate, monounsaturated C20:1 and PUFA C20:2. Of these, selected promoter BMI-dmCpG showed opposite associations with palmitate compared to C20:1 and C20:2. As for the comparison between FS and PS within BMI classes, dmCpG were strikingly more abundant and variably methylated in overweight relative to normoweight or obese subjects (∼70–139-fold, respectively). Overweight-associated dmCpG-hosting genes were significantly enriched in targets for E47, SREBP1, and RREB1 transcription factors, which are known players in obesity and lipid homeostasis, but none overlapped with BMI-dmCpG. We show for the first time that the association of BMI and FA with methylation of disease-related genes is distinct in FS and PS and that limited overlap exists between BMI- and FA-dmCpG within and across prandial states. Our study also identifies a transcriptional regulation circuitry in overweight that might contribute to adaptation to that condition or to transition to obesity. Further work is necessary to define the pathophysiological implications of these findings.

The Interaction Analysis of SNP Variants and DNA Methylation Identifies Novel Methylated Pathogenesis Genes in Congenital Heart Diseases

Congenital heart defect (CHD) is a rare and complicated disease with a high mortality rate. Its etiology remains unclear and includes many aspects. DNA methylation has been indicated to be involved in heart development in the early stage of life, and aberrant methylation level was related to CHDs. This study provides the first evidence of the cross talk of SNP variants and DNA methylation in clarifying CHD underlying genomic cause. We gathered whole exome sequencing (WES) data for Group 1 consisting of patients with PA (n = 78), TOF (n = 20), TAPVC (n = 78), and PDA (n = 40), and 100 healthy children as control group. Rare non-synonymous mutations and novel genes were found and highlighted. Meanwhile, we carried out the second analysis of DNA methylation data from patients with PA (n = 3), TAPVC (n = 3), TOF (n = 3), and PDA (n = 2), and five healthy controls using 850 K array in Group 2. DNA methylation was linked to WES data, and we explored an obvious overlap of hyper/hypomethylated genes. Next, we identified some candidate genes by Fisher’s exact test and Burden analysis; then, those methylated genes were figured out by the criteria of the mutation located in the CpG islands of the genome, differential methylation sites (DMS), and DNA methylation quantitative trait loci (meQTLs) in the database, respectively. Also, the interaction of differentially methylated candidate genes with known CHD pathogenetic genes was depicted in a molecular network. Taken together, our findings show that nine novel genes (ANGPTL4, VEGFA, PAX3, MUC4, HLA-DRB1, TJP2, BCR, PKD1, and HK2) in methylation level are critical to CHD and reveal a new insight into the molecular pathogenesis of CHD.

Causal Inference Methods to Integrate Omics and Complex Traits

Major biotechnological advances have facilitated a tremendous boost to the collection of (gen-/transcript-/prote-/methyl-/metabol-)omics data in very large sample sizes worldwide. Coordinated efforts have yielded a deluge of studies associating diseases with genetic markers (genome-wide association studies) or with molecular phenotypes. Whereas omics-disease associations have led to biologically meaningful and coherent mechanisms, the identified (non-germline) disease biomarkers may simply be correlates or consequences of the explored diseases. To move beyond this realm, Mendelian randomization provides a principled framework to integrate information on omics- and disease-associated genetic variants to pinpoint molecular traits causally driving disease development. In this review, we show the latest advances in this field, flag up key challenges for the future, and propose potential solutions.

African genetic diversity and adaptation inform a precision medicine agenda

The deep evolutionary history of African populations, since the emergence of modern humans more than 300,000 years ago, has resulted in high genetic diversity and considerable population structure. Selected genetic variants have increased in frequency due to environmental adaptation, but recent exposures to novel pathogens and changes in lifestyle render some of them with properties leading to present health liabilities. The unique discoverability potential from African genomic studies promises invaluable contributions to understanding the genomic and molecular basis of health and disease. Globally, African populations are understudied, and precision medicine approaches are largely based on data from European and Asian-ancestry populations, which limits the transferability of findings to the continent of Africa. Africa needs innovative precision medicine solutions based on African data that use knowledge and implementation strategies aligned to its climatic, cultural, economic and genomic diversity.

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci

BACKGROUND

DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach.

METHODS

The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses.

RESULTS

We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development.

CONCLUSIONS

We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Multi-Omics Approaches to Define Calcific Aortic Valve Disease Pathogenesis

Calcific aortic valve disease sits at the confluence of multiple world-wide epidemics of aging, obesity, diabetes, and renal dysfunction, and its prevalence is expected to nearly triple over the next 3 decades. This is of particularly dire clinical relevance, as calcific aortic valve disease can progress rapidly to aortic stenosis, heart failure, and eventually premature death. Unlike in atherosclerosis, and despite the heavy clinical toll, to date, no pharmacotherapy has proven effective to halt calcific aortic valve disease progression, with invasive and costly aortic valve replacement representing the only treatment option currently available. This substantial gap in care is largely because of our still-limited understanding of both normal aortic valve biology and the key regulatory mechanisms that drive disease initiation and progression. Drug discovery is further hampered by the inherent intricacy of the valvular microenvironment: a unique anatomic structure, a complex mixture of dynamic biomechanical forces, and diverse and multipotent cell populations collectively contributing to this currently intractable problem. One promising and rapidly evolving tactic is the application of multiomics approaches to fully define disease pathogenesis. Herein, we summarize the application of (epi)genomics, transcriptomics, proteomics, and metabolomics to the study of valvular heart disease. We also discuss recent forays toward the omics-based characterization of valvular (patho)biology at single-cell resolution; these efforts promise to shed new light on cellular heterogeneity in healthy and diseased valvular tissues and represent the potential to efficaciously target and treat key cell subpopulations. Last, we discuss systems biology- and network medicine-based strategies to extract meaning, mechanisms, and prioritized drug targets from multiomics datasets.

Genetic impacts on DNA methylation: research findings and future perspectives

Multiple recent studies highlight that genetic variants can have strong impacts on a significant proportion of the human DNA methylome. Methylation quantitative trait loci, or meQTLs, allow for the exploration of biological mechanisms that underlie complex human phenotypes, with potential insights for human disease onset and progression. In this review, we summarize recent milestones in characterizing the human genetic basis of DNA methylation variation over the last decade, including heritability findings and genome-wide identification of meQTLs. We also discuss challenges in this field and future areas of research geared to generate insights into molecular processes underlying human complex traits.

Methylation quantitative trait locus analysis of chronic postsurgical pain uncovers epigenetic mediators of genetic risk

Background: Overlap of pathways enriched by single nucleotide polymorphisms and DNA-methylation underlying chronic postsurgical pain (CPSP), prompted pilot study of CPSP-associated methylation quantitative trait loci (meQTL). Materials & methods: Children undergoing spine-fusion were recruited prospectively. Logistic-regression for genome- and epigenome-wide CPSP association and DNA-methylation-single nucleotide polymorphism association/mediation analyses to identify meQTLs were followed by functional genomics analyses. Results: CPSP (n = 20/58) and non-CPSP groups differed in pain-measures. Of 2753 meQTLs, DNA-methylation at 127 cytosine-guanine dinucleotides mediated association of 470 meQTLs with CPSP (p < 0.05). At PARK16 locus, CPSP risk meQTLs were associated with decreased DNA-methylation at RAB7L1 and increased DNA-methylation at PM20D1. Corresponding RAB7L1/PM20D1 blood eQTLs (GTEx) and cytosine-guanine dinucleotide-loci enrichment for histone marks, transcription factor binding sites and ATAC-seq peaks suggest altered transcription factor-binding. Conclusion: CPSP-associated meQTLs indicate epigenetic mechanisms mediate genetic risk. Clinical trial registration: NCT01839461, NCT01731873 (ClinicalTrials.gov).

Evidence for CAT gene being functionally involved in the susceptibility of COVID-19

Novel coronary pneumonia (COVID-19) is a respiratory distress syndrome caused by a new type of coronavirus. Understanding the genetic basis of susceptibility and prognosis to COVID-19 is of great significance to disease prevention, molecular typing, prognosis, and treatment. However, so far, there have been only two genome-wide association studies (GWASs) on the susceptibility of COVID-19. Starting with these reported DNA variants, we found the genes regulated by these variants through cis-eQTL and cis-meQTL acting. We further did a series of bioinformatics analysis on these potential risk genes. The analysis shows that the genetic variants on EHF regulate the expression of its neighbor CAT gene via cis-eQTL. There was significant evidence that CAT and the SARS-CoV-2-related S protein binding protein ACE2 interact with each other. Intracellular localization results showed that CAT and ACE2 proteins both exists in the cell membrane and extracellular area and their interaction could have an impact on the cell invasion ability of S protein. In addition, the expression of these three genes showed a significant positive correlation in the lungs. Based on these results, we propose that CAT plays a crucial intermediary role in binding effectiveness of ACE2, thereby affecting the susceptibility to COVID-19.

Genome-wide association analysis reveals regulation of at-risk loci by DNA methylation in prostate cancer

Epigenetic changes are potentially important for the ontogeny and progression of tumors but are not usually studied because of the complexity of analyzing transcript regulation resulting from epigenetic alterations. Prostate cancer (PCa) is characterized by variable clinical manifestations and frequently unpredictable outcomes. We performed an expression quantitative trait loci (eQTL) analysis to identify the genomic regions that regulate gene expression in PCa and identified a relationship between DNA methylation and clinical information. Using multi-level information published in The Cancer Genome Atlas, we performed eQTL-based analyses on DNA methylation and gene expression. To better interpret these data, we correlated loci and clinical indexes to identify the important loci for both PCa development and progression. Our data demonstrated that although only a small proportion of genes are regulated via DNA methylation in PCa, these genes are enriched in important cancer-related groups. In addition, single nucleotide polymorphism analysis identified the locations of CpG sites and genes within at-risk loci, including the 19q13.2–q13.43 and 16q22.2–q23.1 loci. Further, an epigenetic association study of clinical indexes detected risk loci and pyrosequencing for site validation. Although DNA methylation-regulated genes across PCa samples are a small proportion, the associated genes play important roles in PCa carcinogenesis.

DNA methylation of blood cells is associated with prevalent type 2 diabetes in a meta-analysis of four European cohorts

BACKGROUND

Type 2 diabetes (T2D) is a heterogeneous disease with well-known genetic and environmental risk factors contributing to its prevalence. Epigenetic mechanisms related to changes in DNA methylation (DNAm), may also contribute to T2D risk, but larger studies are required to discover novel markers, and to confirm existing ones.

RESULTS

We performed a large meta-analysis of individual epigenome-wide association studies (EWAS) of prevalent T2D conducted in four European studies using peripheral blood DNAm. Analysis of differentially methylated regions (DMR) was also undertaken, based on the meta-analysis results. We found three novel CpGs associated with prevalent T2D in Europeans at cg00144180 (HDAC4), cg16765088 (near SYNM) and cg24704287 (near MIR23A) and confirmed three CpGs previously identified (mapping to TXNIP, ABCG1 and CPT1A). We also identified 77 T2D associated DMRs, most of them hypomethylated in T2D cases versus controls. In adjusted regressions among diabetic-free participants in ALSPAC, we found that all six CpGs identified in the meta-EWAS were associated with white cell-types. We estimated that these six CpGs captured 11% of the variation in T2D, which was similar to the variation explained by the model including only the common risk factors of BMI, sex, age and smoking (R2 = 10.6%).

CONCLUSIONS

This study identifies novel loci associated with T2D in Europeans. We also demonstrate associations of the same loci with other traits. Future studies should investigate if our findings are generalizable in non-European populations, and potential roles of these epigenetic markers in T2D etiology or in determining long term consequences of T2D.

The Molecular Basis of Predicting Atherosclerotic Cardiovascular Disease Risk

Atherosclerotic cardiovascular disease (ASCVD) proceeds through a series of stages: initiation, progression (or regression), and complications. By integrating known biology regarding molecular signatures of each stage with recent advances in high-dimensional molecular data acquisition platforms (to assay the genome, epigenome, transcriptome, proteome, metabolome, and gut microbiome), snapshots of each phase of atherosclerotic cardiovascular disease development can be captured. In this review, we will summarize emerging approaches for assessment of atherosclerotic cardiovascular disease risk in humans using peripheral blood molecular signatures and molecular imaging approaches. We will then discuss the potential (and challenges) for these snapshots to be integrated into a personalized movie providing dynamic readouts of an individual's atherosclerotic cardiovascular disease risk status throughout the life course.

Flexible and efficient perovskite quantum dot solar cells via hybrid interfacial architecture

All-inorganic CsPbI₃ perovskite quantum dots have received substantial research interest for photovoltaic applications because of higher efficiency compared to solar cells using other quantum dots materials and the various exciting properties that perovskites have to offer. These quantum dot devices also exhibit good mechanical stability amongst various thin-film photovoltaic technologies. We demonstrate higher mechanical endurance of quantum dot films compared to bulk thin film and highlight the importance of further research on high-performance and flexible optoelectronic devices using nanoscale grains as an advantage. Specifically, we develop a hybrid interfacial architecture consisting of CsPbI₃ quantum dot/PCBM heterojunction, enabling an energy cascade for efficient charge transfer and mechanical adhesion. The champion CsPbI₃ quantum dot solar cell has an efficiency of 15.1% (stabilized power output of 14.61%), which is among the highest report to date. Building on this strategy, we further demonstrate a highest efficiency of 12.3% in flexible quantum dot photovoltaics.

Mendelian randomization analysis identified genes pleiotropically associated with the risk and prognosis of COVID-19

OBJECTIVES

COVID-19 has caused a large global pandemic. Patients with COVID-19 exhibited considerable variation in disease behavior. Pervious genome-wide association studies have identified potential genetic variants involved in the risk and prognosis of COVID-19, but the underlying biological interpretation remains largely unclear.

METHODS

We applied the summary data-based Mendelian randomization (SMR) method to identify genes that were pleiotropically associated with the risk and various outcomes of COVID-19, including severe respiratory confirmed COVID-19 and hospitalized COVID-19.

RESULTS

In blood, we identified 2 probes, ILMN_1765146 and ILMN_1791057 tagging IFNAR2, that showed pleiotropic association with hospitalized COVID-19 (β [SE]=0.42 [0.09], P = 4.75 × 10-06 and β [SE]=-0.48 [0.11], P = 6.76 × 10-06, respectively). Although no other probes were significant after correction for multiple testing in both blood and lung, multiple genes as tagged by the top 5 probes were involved in inflammation or antiviral immunity, and several other tagged genes, such as PON2 and HPS5, were involved in blood coagulation.

CONCLUSIONS

We identified IFNAR2 and other potential genes that could be involved in the susceptibility or prognosis of COVID-19. These findings provide important leads to a better understanding of the mechanisms of cytokine storm and venous thromboembolism in COVID-19 and potential therapeutic targets for the effective treatment of COVID-19.

DNA Methylation in Ovarian Cancer Susceptibility

Epigenetic alterations are somatically acquired over the lifetime and during neoplastic transformation but may also be inherited as widespread 'constitutional' alterations in normal tissues that can cause cancer predisposition. Epithelial ovarian cancer (EOC) has an established genetic susceptibility and mounting epidemiological evidence demonstrates that DNA methylation (DNAm) intermediates as well as independently contributes to risk. Targeted studies of known EOC susceptibility genes (CSGs) indicate rare, constitutional BRCA1 promoter methylation increases familial and sporadic EOC risk. Blood-based epigenome-wide association studies (EWAS) for EOC have detected a total of 2846 differentially methylated probes (DMPs) with 71 genes replicated across studies despite significant heterogeneity. While EWAS detect both symptomatic and etiologic DMPs, adjustments and analytic techniques may enrich risk associations, as evidenced by the detection of dysregulated methylation of BNC2-a known CSG identified by genome-wide associations studies (GWAS). Integrative genetic-epigenetic approaches have mapped methylation quantitative trait loci (meQTL) to EOC risk, revealing DNAm variations that are associated with nine GWAS loci and, further, one novel risk locus. Increasing efforts to mapping epigenome variation across populations and cell types will be key to decoding both the genomic and epigenomic causal pathways to EOC.

A simple, rapid, interpretable, actionable and implementable digital PCR based mortality index

Mortality assessments are conducted for both civil and commercial purposes. Recent advances in epigenetics have resulted in DNA methylation tools to assess risk and aid in this task. However, widely available array-based algorithms are not readily translatable into clinical tools and do not provide a good foundation for clinical recommendations. Further, recent work shows evidence of heritability and possible racial bias in these indices. Using a publicly available array data set, the Framingham Heart Study (FHS), we develop and test a five-locus mortality-risk algorithm using only previously validated methylation biomarkers that have been shown to be free of racial bias, and that provide specific assessments of smoking, alcohol consumption, diabetes and heart disease. We show that a model using age, sex and methylation measurements at these five loci outperforms the 513 probe Levine index and approximates the predictive power of the 1030 probe GrimAge index. We then show each of the five loci in our algorithm can be assessed using a more powerful, reference-free digital PCR approach, further demonstrating that it is readily clinically translatable. Finally, we show the loci do not reflect ethnically specific variation. We conclude that this algorithm is a simple, yet powerful tool for assessing mortality risk. We further suggest that the output from this or similarly derived algorithms using either array or digital PCR can be used to provide powerful feedback to patients, guide recommendations for additional medical assessments, and help monitor the effect of public health prevention interventions.

A Transcriptome-Wide Association Study Identifies Candidate Susceptibility Genes for Pancreatic Cancer Risk

Pancreatic cancer is among the most well-characterized cancer types, yet a large proportion of the heritability of pancreatic cancer risk remains unclear. Here, we performed a large transcriptome-wide association study to systematically investigate associations between genetically predicted gene expression in normal pancreas tissue and pancreatic cancer risk. Using data from 305 subjects of mostly European descent in the Genotype-Tissue Expression Project, we built comprehensive genetic models to predict normal pancreas tissue gene expression, modifying the UTMOST (unified test for molecular signatures). These prediction models were applied to the genetic data of 8,275 pancreatic cancer cases and 6,723 controls of European ancestry. Thirteen genes showed an association of genetically predicted expression with pancreatic cancer risk at an FDR ≤ 0.05, including seven previously reported genes (INHBA, SMC2, ABO, PDX1, RCCD1, CFDP1, and PGAP3) and six novel genes not yet reported for pancreatic cancer risk [6q27: SFT2D1 OR (95% confidence interval (CI), 1.54 (1.25-1.89); 13q12.13: MTMR6 OR (95% CI), 0.78 (0.70-0.88); 14q24.3: ACOT2 OR (95% CI), 1.35 (1.17-1.56); 17q12: STARD3 OR (95% CI), 6.49 (2.96-14.27); 17q21.1: GSDMB OR (95% CI), 1.94 (1.45-2.58); and 20p13: ADAM33 OR (95% CI): 1.41 (1.20-1.66)]. The associations for 10 of these genes (SFT2D1, MTMR6, ACOT2, STARD3, GSDMB, ADAM33, SMC2, RCCD1, CFDP1, and PGAP3) remained statistically significant even after adjusting for risk SNPs identified in previous genome-wide association study. Collectively, this analysis identified novel candidate susceptibility genes for pancreatic cancer that warrant further investigation. SIGNIFICANCE: A transcriptome-wide association analysis identified seven previously reported and six novel candidate susceptibility genes for pancreatic cancer risk.

Mendelian randomization integrating GWAS and mQTL data identified novel pleiotropic DNA methylation loci for neuropathology of Alzheimer's disease

The pathogenesis of Alzheimer's disease (AD) remains largely unclear. Exploring the genetic/epigenetic loci showing pleiotropic association with the neuropathologies of AD may greatly enhance understanding of the mechanisms underlying the development of AD. In this study, using data from the Religious Orders Study and the Rush Memory and Aging Project, we undertook a Mendelian randomization approach integrating genome-wide association studies (GWASs) and DNA methylation quantitative trait locus data to explore pleiotropic epigenetic loci for AD neuropathologies, including amyloid-β (Aβ) load and tau-containing neurofibrillary tangle density. We performed GWASs of DNA methylation in brain tissues from 592 participants and mapped 60,595 cis-SNP-CpG pairs after correction for multiple testing. By linking cis-DNA methylation quantitative trait locus with GWAS results for Aβ load and tau tangles, we identified 47 CpGs showing pleiotropic association with Aβ load by the Mendelian randomization analysis. We then used gene expression data from 537 individuals and performed quantitative trait methylation analysis. We found that 18 of the 47 CpGs were in cis associated with 25 mRNAs/genes, comprising 41 unique CpG-mRNA/gene pairs. Our findings shed light on the role of DNA methylation in the pathogenesis of Aβ.

Genome-Wide Characterization of Host Transcriptional and Epigenetic Alterations During HIV Infection of T Lymphocytes

Background and methods: Host genomic alterations are closely related to dysfunction of CD4+ T lymphocytes in the HIV-host interplay. However, the roles of aberrant DNA methylation and gene expression in the response to HIV infection are not fully understood. We investigated the genome-wide DNA methylation and transcriptomic profiles in two HIV-infected T lymphocyte cell lines using high-throughput sequencing. Results: Based on DNA methylation data, we identified 3,060 hypomethylated differentially methylated regions (DMRs) and 2,659 hypermethylated DMRs in HIV-infected cells. Transcription-factor-binding motifs were significantly associated with methylation alterations, suggesting that DNA methylation modulates gene expression by affecting the binding to transcription factors during HIV infection. In support of this hypothesis, genes with promoters overlapping with DMRs were enriched in the biological function related to transcription factor activities. Furthermore, the analysis of gene expression data identified 1,633 upregulated genes and 2,142 downregulated genes on average in HIV-infected cells. These differentially expressed genes (DEGs) were significantly enriched in apoptosis-related pathways. Our results suggest alternative splicing as an additional mechanism that may contribute to T-cell apoptosis during HIV infection. We also demonstrated a genome-scale correlation between DNA methylation and gene expression in HIV-infected cells. We identified 831 genes with alterations in both DNA methylation and gene expression, which were enriched in apoptosis. Our results were validated using various experimental methods. In addition, consistent with our in silico results, a luciferase assay showed that the activity of the PDX1 and SMAD3 promoters was significantly decreased in the presence of HIV proteins, indicating the potential of these genes as genetic markers of HIV infection. Conclusions: Our results suggest important roles for DNA methylation and gene expression regulation in T-cell apoptosis during HIV infection. We propose a list of novel genes related to these processes for further investigation. This study also provides a comprehensive characterization of changes occurring at the transcriptional and epigenetic levels in T cells in response to HIV infection.

Mendelian randomization analysis identified genes pleiotropically associated with the risk and prognosis of COVID-19

OBJECTIVES

COVID-19 has caused a large global pandemic. Patients with COVID-19 exhibited considerable variation in disease behavior. Pervious genome-wide association studies have identified potential genetic variants involved in the risk and prognosis of COVID-19, but the underlying biological interpretation remains largely unclear.

METHODS

We applied the summary data-based Mendelian randomization (SMR) method to identify genes that were pleiotropically associated with the risk and various outcomes of COVID-19, including severe respiratory confirmed COVID-19 and hospitalized COVID-19.

RESULTS

In blood, we identified 2 probes, ILMN_1765146 and ILMN_1791057 tagging IFNAR2, that showed pleiotropic association with hospitalized COVID-19 (Beta; [SE]=0.42 [0.09], P=4.75E-06 and Beta; [SE]=-0.48 [0.11], P=6.76E-06, respectively). Although no other probes were significant after correction for multiple testing in both blood and lung, multiple genes as tagged by the top 5 probes were involved in inflammation or antiviral immunity, and several other tagged genes, such as PON2 and HPS5, were involved in blood coagulation.

CONCLUSIONS

We identified IFNAR2 and other potential genes that could be involved in the susceptibility or prognosis of COVID-19. These findings provide important leads to a better understanding of the mechanisms of cytokine storm and venous thromboembolism in COVID-19 and potential therapeutic targets for the effective treatment of COVID-19.

Genome-wide identification of genes regulating DNA methylation using genetic anchors for causal inference

BACKGROUND

DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA methylation patterns in blood using large-scale population genomics data.

RESULTS

By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g., NFKBIE, CDCA7(L), and NLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation.

CONCLUSION

We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation.

Genetics meets proteomics: perspectives for large population-based studies

Proteomic analysis of cells, tissues and body fluids has generated valuable insights into the complex processes influencing human biology. Proteins represent intermediate phenotypes for disease and provide insight into how genetic and non-genetic risk factors are mechanistically linked to clinical outcomes. Associations between protein levels and DNA sequence variants that colocalize with risk alleles for common diseases can expose disease-associated pathways, revealing novel drug targets and translational biomarkers. However, genome-wide, population-scale analyses of proteomic data are only now emerging. Here, we review current findings from studies of the plasma proteome and discuss their potential for advancing biomedical translation through the interpretation of genome-wide association analyses. We highlight the challenges faced by currently available technologies and provide perspectives relevant to their future application in large-scale biobank studies.

An integrative multi-omics analysis to identify candidate DNA methylation biomarkers related to prostate cancer risk

It remains elusive whether some of the associations identified in genome-wide association studies of prostate cancer (PrCa) may be due to regulatory effects of genetic variants on CpG sites, which may further influence expression of PrCa target genes. To search for CpG sites associated with PrCa risk, here we establish genetic models to predict methylation (N = 1,595) and conduct association analyses with PrCa risk (79,194 cases and 61,112 controls). We identify 759 CpG sites showing an association, including 15 located at novel loci. Among those 759 CpG sites, methylation of 42 is associated with expression of 28 adjacent genes. Among 22 genes, 18 show an association with PrCa risk. Overall, 25 CpG sites show consistent association directions for the methylation-gene expression-PrCa pathway. We identify DNA methylation biomarkers associated with PrCa, and our findings suggest that specific CpG sites may influence PrCa via regulating expression of candidate PrCa target genes.

Investigating asthma heterogeneity through shared and distinct genetics: Insights from genome-wide cross-trait analysis

Asthma is a heterogeneous respiratory disease reflecting distinct pathobiologic mechanisms. These mechanisms are based, at least partly, on different genetic factors shared by many other conditions, such as allergic diseases and obesity. Investigating the shared genetic effects enables better understanding of the mechanisms of phenotypic correlations and is less subject to confounding by environmental factors. The increasing availability of large-scale genome-wide association study (GWAS) for asthma has enabled researchers to examine the genetic contributions to the epidemiologic associations between asthma subtypes and those between coexisting diseases and/or traits and asthma. Studies have found not only shared but also distinct genetic components between asthma subtypes, indicating that the heterogeneity is related to distinct genetics. This review summarizes a recently compiled analytic approach-genome-wide cross-trait analysis-to determine shared and distinct genetic architecture. The genome-wide cross-trait analysis features in several analytic aspects: genetic correlation, cross-trait meta-analysis, Mendelian randomization, polygenic risk score, and functional analysis. In this article, we discuss in detail the scientific goals that can be achieved by these analyses, their advantages, and their limitations. We also make recommendations for future directions: (1) ethnicity-specific asthma GWASs and (2) application of cross-trait methods to multiomics data to dissect the heritability found in GWASs. Finally, these analytic approaches are also applicable to complex and heterogeneous traits beyond asthma.

DNA methylation and incident cardiovascular disease

PURPOSE OF REVIEW

DNA methylation is an epigenetic mark that reflects both genetic and environmental influences over the life course and has the potential to be a robust biomarker for cardiovascular disease (CVD) risk. However, standard association studies linking DNA methylation and CVD are susceptible to reverse causation and may not directly translate into useful biomarkers of future disease. Studies of incident CVD represent a crucial tool for improving this evidence base.

RECENT FINDINGS

Recent investigations have started to provide links between DNA methylation and incident CVD. Epigenome-wide association studies have suggested individual genetic loci in which differential methylation exists prior to disease onset, while multivariate predictive modeling approaches have made progress towards realizing the potential for DNA methylation as a predictive biomarker of CVD risk. Meanwhile, complementary analysis strategies such as Mendelian randomization have provided clues as to the causality of these epigenomic associations.

SUMMARY

Taken together, this wave of studies provides the basis for a better understanding of CVD pathophysiology and the development of more confident biomarkers of CVD risk.

Array-Based Epigenetic Aging Indices May Be Racially Biased

Epigenetic aging (EA) indices are frequently used as predictors of mortality and other important health outcomes. However, each of the commonly used array-based indices has significant heritable components which could tag ethnicity and potentially confound comparisons across racial and ethnic groups. To determine if this was possible, we examined the relationship of DNA methylation in cord blood from 203 newborns (112 African American (AA) and 91 White) at the 513 probes from the Levine PhenoAge Epigenetic Aging index to ethnicity. Then, we examined all sites significantly associated with race in the newborn sample to determine if they were also associated with an index of ethnic genetic heritage in a cohort of 505 AA adults. After Bonferroni correction, methylation at 50 CpG sites was significantly associated with ethnicity in the newborn cohort. The five most significant sites predicted ancestry with a receiver operator characteristic area under the curve of 0.97. Examination of the top 50 sites in the AA adult cohort showed that methylation status at 11 of those sites was also associated with percentage European ancestry. We conclude that the Levine PhenoAge Index is influenced by cryptic ethnic-specific genetic influences. This influence may extend to similarly constructed EA indices and bias cross-race comparisons.

A genome-wide multiphenotypic association analysis identified candidate genes and gene ontology shared by four common risky behaviors

BACKGROUND

Risky behaviors can lead to huge economic and health losses. However, limited efforts are paid to explore the genetic mechanisms of risky behaviors.

RESULT

MASH analysis identified a group of target genes for risky behaviors, such as APBB2, MAPT and DCC. For GO enrichment analysis, FUMA detected multiple risky behaviors related GO terms and brain related diseases, such as regulation of neuron differentiation (adjusted P value = 2.84×10^-5), autism spectrum disorder (adjusted P value =1.81×10^-27) and intelligence (adjusted P value =5.89×10^-15).

CONCLUSION

We reported multiple candidate genes and GO terms shared by the four risky behaviors, providing novel clues for understanding the genetic mechanism of risky behaviors.

METHODS

Multivariate Adaptive Shrinkage (MASH) analysis was first applied to the GWAS data of four specific risky behaviors (automobile speeding, drinks per week, ever-smoker, number of sexual partners) to detect the common genetic variants shared by the four risky behaviors. Utilizing genomic functional annotation data of SNPs, the SNPs detected by MASH were then mapped to target genes. Finally, gene set enrichment analysis of the identified candidate genes were conducted by the FUMA platform to obtain risky behaviors related gene ontology (GO) terms as well as diseases and traits, respectively.

Thioredoxin-interacting protein promotes activation and inflammation of monocytes with DNA demethylation in coronary artery disease

Numerous studies have demonstrated that thioredoxin‐interacting protein (TXNIP) expression of peripheral blood leucocytes is increased in coronary artery disease (CAD). However, the molecular mechanism of this phenomenon remained unclear. DNA methylation plays important roles in the regulation of gene expression. Therefore, we speculated there might be a close association between the expression of TXNIP and methylation. In this study, we found that compared with controls, DNA methylation at cg19693031 was decreased in CAD, while mRNA expressions of TXNIP and inflammatory factors, NLRP3, IL‐1β, IL‐18, were increased. Methylation at cg19693031 was negatively associated with TXNIP expression in the cohort, THP‐1 and macrophages/foam cells. Furthermore, Transwell assay and co‐cultured adhesion assay were performed to investigate functions of TXNIP on the migration of THP‐1 or the adhesion of THP‐1 on the surface of endothelial cells, respectively. Notably, overexpressed TXNIP promoted the migration and adhesion of THP‐1 cells and expressions of NLRP3, IL‐18 and IL‐1β. Oppositely, knock‐down TXNIP inhibited the migration and adhesion of THP‐1 and expressions of NLRP3, IL‐18. In conclusion, increased TXNIP expression, related to cg19693031 demethylation orientates monocytes towards an inflammatory status through the NLRP3 inflammasome pathway involved in the development of CAD.