African ancestry-derived APOL1 risk genotypes show proximal epigenetic associations.

Apolipoprotein L1 (APOL1) coding variants, termed G1 and G2, are established genetic risk factors for a growing spectrum of diseases, including kidney disease, in individuals of African ancestry. Evidence suggests that the risk variants, which show a recessive mode of inheritance, lead to toxic gain-of-function changes of the APOL1 protein. Disease occurrence and presentation vary, likely due to modifiers or second hits. To understand the role of the epigenetic landscape in relation to APOL1 risk variants, we performed methylation quantitative trait locus (meQTL) analysis to identify differentially methylated CpGs influenced by APOL1 risk variants in 611 African American individuals. We identified five CpGs that were significantly associated with APOL1 risk alleles in discovery and replication studies, and one CpG-APOL1 association was independent of other genomic variants. Our study highlights proximal DNA methylation alterations that may help explain the variable disease risk and clinical manifestation of APOL1 variants.

Epigenome-wide association study of lung cancer among never smokers in two prospective cohorts in Shanghai, China.

BACKGROUND: The aetiology of lung cancer among individuals who never smoked remains elusive, despite 15% of lung cancer cases in men and 53% in women worldwide being unrelated to smoking. Epigenetic alterations, particularly DNA methylation (DNAm) changes, have emerged as potential drivers. Yet, few prospective epigenome-wide association studies (EWAS), primarily focusing on peripheral blood DNAm with limited representation of never smokers, have been conducted. METHODS: We conducted a nested case-control study of 80 never-smoking incident lung cancer cases and 83 never-smoking controls within the Shanghai Women's Health Study and Shanghai Men's Health Study. DNAm was measured in prediagnostic oral rinse samples using Illumina MethylationEPIC array. Initially, we conducted an EWAS to identify differentially methylated positions (DMPs) associated with lung cancer in the discovery sample of 101 subjects. The top 50 DMPs were further evaluated in a replication sample of 62 subjects, and results were pooled using fixed-effect meta-analysis. RESULTS: Our study identified three DMPs significantly associated with lung cancer at the epigenome-wide significance level of p<8.22×10-8. These DMPs were identified as cg09198866 (MYH9; TXN2), cg01411366 (SLC9A10) and cg12787323. Furthermore, examination of the top 1000 DMPs indicated significant enrichment in epithelial regulatory regions and their involvement in small GTPase-mediated signal transduction pathways. Additionally, GrimAge acceleration was identified as a risk factor for lung cancer (OR=1.19 per year; 95% CI 1.06 to 1.34). CONCLUSIONS: While replication in a larger sample size is necessary, our findings suggest that DNAm patterns in prediagnostic oral rinse samples could provide novel insights into the underlying mechanisms of lung cancer in never smokers.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits.

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Genome-wide study investigating effector genes and polygenic prediction for kidney function in persons with ancestry from Africa and the Americas.

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p < 5 × 10-8), of which two have not been previously reported in any ancestry group. We integrated fine-mapped loci with epigenomic and transcriptomic resources to highlight potential effector genes relevant to kidney physiology and disease, and reveal key regulatory elements and pathways involved in renal function and development. We demonstrate the varying but increased predictive power offered by a multi-ancestry polygenic score for eGFR and highlight the importance of population diversity in GWASs and multi-omics resources to enhance opportunities for clinical translation for all.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases.

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation.

BACKGROUND: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS). METHODS: We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C). FINDINGS: Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs. INTERPRETATION: Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases. FUNDING: Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

DNA methylation networks underlying mammalian traits.

Using DNA methylation profiles (n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

Potentially causal associations between placental DNA methylation and schizophrenia and other neuropsychiatric disorders.

Increasing evidence supports the role of placenta in neurodevelopment and potentially, in the later onset of neuropsychiatric disorders. Recently, methylation quantitative trait loci (mQTL) and interaction QTL (iQTL) maps have proven useful to understand SNP-genome wide association study (GWAS) relationships, otherwise missed by conventional expression QTLs. In this context, we propose that part of the genetic predisposition to complex neuropsychiatric disorders acts through placental DNA methylation (DNAm). We constructed the first public placental cis-mQTL database including nearly eight million mQTLs calculated in 368 fetal placenta DNA samples from the INMA project, ran cell type- and gestational age-imQTL models and combined those data with the summary statistics of the largest GWAS on 10 neuropsychiatric disorders using Summary-based Mendelian Randomization (SMR) and colocalization. Finally, we evaluated the influence of the DNAm sites identified on placental gene expression in the RICHS cohort. We found that placental cis-mQTLs are highly enriched in placenta-specific active chromatin regions, and useful to map the etiology of neuropsychiatric disorders at prenatal stages. Specifically, part of the genetic burden for schizophrenia, bipolar disorder and major depressive disorder confers risk through placental DNAm. The potential causality of several of the observed associations is reinforced by secondary association signals identified in conditional analyses, regional pleiotropic methylation signals associated to the same disorder, and cell type-imQTLs, additionally associated to the expression levels of relevant immune genes in placenta. In conclusion, the genetic risk of several neuropsychiatric disorders could operate, at least in part, through DNAm and associated gene expression in placenta.

Predictors of postoperative complications and readmissions in laparoscopic pancreas resection: Results of a cohort 105 consecutive cases. A retrospective study.

INTRODUCTION: Laparoscopic resection of the pancreas (LRP) has been implemented to a varying degree because it is technically demanding and requires a long learning curve. In the present study we analyze the risk factors for complications and hospital readmissions in a single center study of 105 consecutive LRPs. METHODS: We conducted a retrospective study using a prospective database. Data were collected on age, gender, BMI, ASA score, type of surgery, histologic type, operative time, hospital stay, postoperative complications, degree of severity and hospital readmission. RESULTS: The cohort included 105 patients, 63 females and 42 males with a median age and BMI of 58 (53-70) and 25.5 (22,2-27.9) respectively. Eighteen (17%) central pancreatectomies, 5 (4.8%) enucleations, 81 (77.6%) distal pancreatectomies and one total pancreatectomy were performed. Fifty-six patients (53.3%) experienced some type of complication, of which 13 (12.3%) were severe (Clavien-Dindo > IIIb) and 11 (10.5%) patients were readmitted in the first 30 days after surgery. In the univariate analysis, age, male gender, ASA score, central pancreatectomy and operative time were significantly associated with the development of complications (P <0.05). In the multivariate analysis, male gender (OR 7.97; 95% CI 1.08-58.88)), severe complications (OR 59.40; 95% CI, 7.69-458.99), and the development of intrabdominal collections (OR 8.97; 95% CI, 1.28-63.02)) were associated with hospital readmission. CONCLUSIONS: Age, male gender, ASA score, operative time and central pancreatectomy are associated with a higher incidence of complications. Male gender, severe complications and intraabdominal collections are associated with more hospital readmissions.

Diversity in EWAS: current state, challenges, and solutions.

Here, we report a lack of diversity in epigenome-wide association studies (EWAS) and DNA methylation (DNAm) data, discuss current challenges, and propose solutions for EWAS and DNAm research in diverse populations. The strategies we propose include fostering community involvement, new data generation, and cost-effective approaches such as locus-specific analysis and ancestry variable region analysis.

Cell Type-Specific Signal Analysis in Epigenome-Wide Association Studies.

Hundreds of epigenome-wide association studies (EWAS) have been performed, successfully identifying replicated epigenomic signals in processes such as aging and smoking. Despite this progress, it remains a major challenge in EWAS to detect both cell type-specific and cell type confounding effects impacting study results. One way to identify these effects is through eFORGE (experimentally derived Functional element Overlap analysis of ReGions from EWAS), a published tool that uses 815 datasets from large-scale mapping studies to detect enriched tissues, cell types, and genomic regions. Here, I show that eFORGE analysis can be extended to EWAS differentially variable positions (DVPs), identifying target cell types and tissues. In addition, I also show that eFORGE tissue-specific enrichment can be detected for sites below EWAS significance threshold. I develop on these and other analysis examples, extending our knowledge of eFORGE cell type- and tissue-specific enrichment results for different EWAS.

A pan-tissue DNA methylation atlas enables in silico decomposition of human tissue methylomes at cell-type resolution.

Bulk-tissue DNA methylomes represent an average over many different cell types, hampering our understanding of cell-type-specific contributions to disease development. As single-cell methylomics is not scalable to large cohorts of individuals, cost-effective computational solutions are needed, yet current methods are limited to tissues such as blood. Here we leverage the high-resolution nature of tissue-specific single-cell RNA-sequencing datasets to construct a DNA methylation atlas defined for 13 solid tissue types and 40 cell types. We comprehensively validate this atlas in independent bulk and single-nucleus DNA methylation datasets. We demonstrate that it correctly predicts the cell of origin of diverse cancer types and discovers new prognostic associations in olfactory neuroblastoma and stage 2 melanoma. In brain, the atlas predicts a neuronal origin for schizophrenia, with neuron-specific differential DNA methylation enriched for corresponding genome-wide association study risk loci. In summary, the DNA methylation atlas enables the decomposition of 13 different human tissue types at a high cellular resolution, paving the way for an improved interpretation of epigenetic data.

A mammalian methylation array for profiling methylation levels at conserved sequences.

Infinium methylation arrays are not available for the vast majority of non-human mammals. Moreover, even if species-specific arrays were available, probe differences between them would confound cross-species comparisons. To address these challenges, we developed the mammalian methylation array, a single custom array that measures up to 36k CpGs per species that are well conserved across many mammalian species. We designed a set of probes that can tolerate specific cross-species mutations. We annotate the array in over 200 species and report CpG island status and chromatin states in select species. Calibration experiments demonstrate the high fidelity in humans, rats, and mice. The mammalian methylation array has several strengths: it applies to all mammalian species even those that have not yet been sequenced, it provides deep coverage of conserved cytosines facilitating the development of epigenetic biomarkers, and it increases the probability that biological insights gained in one species will translate to others.

Integrative analysis of 3604 GWAS reveals multiple novel cell type-specific regulatory associations.

BACKGROUND: Genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) are known to preferentially co-locate to active regulatory elements in tissues and cell types relevant to disease aetiology. Further characterisation of associated cell type-specific regulation can broaden our understanding of how GWAS signals may contribute to disease risk. RESULTS: To gain insight into potential functional mechanisms underlying GWAS associations, we developed FORGE2 ( https://forge2.altiusinstitute.org/ ), which is an updated version of the FORGE web tool. FORGE2 uses an expanded atlas of cell type-specific regulatory element annotations, including DNase I hotspots, five histone mark categories and 15 hidden Markov model (HMM) chromatin states, to identify tissue- and cell type-specific signals. An analysis of 3,604 GWAS from the NHGRI-EBI GWAS catalogue yielded at least one significant disease/trait-tissue association for 2,057 GWAS, including > 400 associations specific to epigenomic marks in immune tissues and cell types, > 30 associations specific to heart tissue, and > 60 associations specific to brain tissue, highlighting the key potential of tissue- and cell type-specific regulatory elements. Importantly, we demonstrate that FORGE2 analysis can separate previously observed accessible chromatin enrichments into different chromatin states, such as enhancers or active transcription start sites, providing a greater understanding of underlying regulatory mechanisms. Interestingly, tissue-specific enrichments for repressive chromatin states and histone marks were also detected, suggesting a role for tissue-specific repressed regions in GWAS-mediated disease aetiology. CONCLUSION: In summary, we demonstrate that FORGE2 has the potential to uncover previously unreported disease-tissue associations and identify new candidate mechanisms. FORGE2 is a transparent, user-friendly web tool for the integrative analysis of loci discovered from GWAS.

Comparison of phenotypes and outcomes following resection of incidental versus symptomatic pancreatic neuroendocrine tumors.

BACKGROUND AND AIM: fifty to 70 percent of pancreatic neuroendocrine tumors are diagnosed incidentally. The objective of this study was to compare the phenotype and oncological outcomes of incidental versus symptomatic pancreatic neuroendocrine tumors. METHODS: a retrospective study was performed, identifying all incidental and symptomatic tumors resected between 2000 and 2019. Baseline characteristics, symptoms, operative variables and pathological stage were all recorded. Patterns of recurrence and overall and disease-free survival were analyzed in both groups. RESULTS: fifty-one incidental and 45 symptomatic pancreatic tumor resections were performed. Symptomatic tumors were more frequent in females (29 vs 17; p = 0.005) and younger patients (median years; 50 vs 58; p = 0.012) and were detected at a more advanced stage (p = 0.027). There were no differences in location and most resections (n = 49; 51 %) were performed laparoscopically. There were no operative mortalities and 17 (17.7 %) severe complications (≥ IIIb on the Clavien-Dindo classification) were recorded with no differences between the two groups. With a median follow-up of 64.4 months (range 13.5-90), overall survival at five and ten years was 89.7 % and 72.8 % for the non-incidental tumors and 80.9 % and 54.6 % for the incidental tumors (p = ns), respectively. Disease-free survival in both groups (excluding M1a) was 71.2 % and 47.5 %, and 93.7 % and 78.1 %, respectively (p = ns). CONCLUSIONS: symptomatic tumors are more frequent in females and present at more advanced pathological stages. There were no significant differences in overall and disease-free survival between the two groups. Resection of incidental tumors ≥ 1.5-2 cm seems advisable, although each case should be assessed on an individual basis.

Epigenetic clock and methylation studies in the rhesus macaque.

Methylation levels at specific CpG positions in the genome have been used to develop accurate estimators of chronological age in humans, mice, and other species. Although epigenetic clocks are generally species-specific, the principles underpinning them appear to be conserved at least across the mammalian class. This is exemplified by the successful development of epigenetic clocks for mice and several other mammalian species. Here, we describe epigenetic clocks for the rhesus macaque (Macaca mulatta), the most widely used nonhuman primate in biological research. Using a custom methylation array (HorvathMammalMethylChip40), we profiled n = 281 tissue samples (blood, skin, adipose, kidney, liver, lung, muscle, and cerebral cortex). From these data, we generated five epigenetic clocks for macaques. These clocks differ with regard to applicability to different tissue types (pan-tissue, blood, skin), species (macaque only or both humans and macaques), and measure of age (chronological age versus relative age). Additionally, the age-based human-macaque clock exhibits a high age correlation (R = 0.89) with the vervet monkey (Chlorocebus sabaeus), another Old World species. Four CpGs within the KLF14 promoter were consistently altered with age in four tissues (adipose, blood, cerebral cortex, skin). Future studies will be needed to evaluate whether these epigenetic clocks predict age-related conditions in the rhesus macaque.

Identical twins carry a persistent epigenetic signature of early genome programming.

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.