Search

Background/objectives: Human identical twins are born at a rate of 3-4 per 1000 live births. Many other mammals also occasionally produce monozygotic twins, referred to as sporadic polyembryony. The underlying mechanisms are unknown. Through epigenome-wide association studies (EWAS), we identified a robust DNA methylation signature in somatic tissues from human monozygotic (MZ) twins, comprising 834 differentially methylated positions (MZ-DMPs). The results point to a connection between monozygotic twinning and early genome programming and enable new angles to study monozygotic twinning., Methods: The mammalian methylation array (MMA) measures 38,608 CpGs focusing on regions that are well-conserved across many mammalian species, allowing for pan-mammalian comparative epigenomic studies. Here, we successfully map human MZ-DMPs to probes of the mammalian methylation array across 157 mammalian genomes., Results: As expected, based on the modest probe overlap between Illumina 450k/EPIC and mammalian methylation array probes, only a subset of MZ-DMPs reside in conserved regions covered by the mammalian methylation array. These include probes mapping to NPAS3 , KLHL35 , CASZ1 , and ATP2B2 . Re-analysis restricting the original EWAS in humans to conserved MMA regions yielded additional MZ-DMPs, suggesting that more loci may be detected by application of the mammalian array to monozygotic twins., Conclusions: In conclusion, the mammalian methylation array may prove to be a promising platform to study whether a shared DNA methylation signature of sporadic polyembryony exists across diverse mammalian species. This may potentially point to shared underlying mechanisms.

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., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. No commercial re-use. See rights and permissions. Published by BMJ.)

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling ( GNG7 , RGS3 ) and innate immune response ( SLC15A4 ), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm., Competing Interests: Conflicts of Interest Nothing to declare.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers ( 131 I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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., (© 2024. The Author(s).)

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., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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., Competing Interests: Declaration of interests DLM is a part-time employee of Optima Partners Ltd. ETGK received a grant from the Netherlands Lung Foundation. GHK received grants or contracts from ZON-MW, Vertex, Netherlands Lung Foundation, GSK, TEVA the Netherlands, and European Union; consulting fees from Astra Zeneca (money to institution); honoraria from Sanofi, Boehringer Ingelheim; and chairs the exquAlro Foundation. MCM received grants from the Research Council of Norway and European Research Council. REM is a scientific advisor to the Epigenetic Clock Development Foundation and Optima Partners. All other authors have nothing to disclose., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of interests P.N. reports investigator-initiated research grants from Allelica, Apple, Amgen, Boston Scientific, Genentech/Roche, and Novartis, personal fees from Allelica, Apple, AstraZeneca, Blackstone Life Sciences, Eli Lilly & Co, Foresite Labs, Genentech/Roche, GV, HeartFlow, Magnet Biomedicine, and Novartis, scientific advisory board membership of Esperion Therapeutics, Preciseli, and TenSixteen Bio, scientific co-founder of TenSixteen Bio, equity in MyOme, Preciseli, and TenSixteen Bio, and spousal employment at Vertex Pharmaceuticals, all unrelated to the present work. M.B.L. has received speaker and advisory fees from Otsuka, Reata, Bayer, and Sanofi Genzyme. G.P. has received speaker and advisory fees from Amgen, Bayer, Novartis, and Sanofi, and has received research grants from Bayer and Sanofi. L.M.R. is a consultant for the TOPMed Administrative Coordinating Center (through Westat)., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals., (© 2023. The Author(s).)

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.

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., (Copyright © 2022. Published by Elsevier España, S.L.U.)

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

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

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., (© 2022. The Author(s).)

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., (© 2022. The Author(s).)

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., (© 2021. The Author(s).)

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., (© 2022. Springer Science+Business Media, LLC, part of Springer Nature.)

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., (© 2021. The Author(s).)

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., (© 2021. The Author(s).)

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.

Exceptionally long-lived species, including many bats, rarely show overt signs of aging, making it difficult to determine why species differ in lifespan. Here, we use DNA methylation (DNAm) profiles from 712 known-age bats, representing 26 species, to identify epigenetic changes associated with age and longevity. We demonstrate that DNAm accurately predicts chronological age. Across species, longevity is negatively associated with the rate of DNAm change at age-associated sites. Furthermore, analysis of several bat genomes reveals that hypermethylated age- and longevity-associated sites are disproportionately located in promoter regions of key transcription factors (TF) and enriched for histone and chromatin features associated with transcriptional regulation. Predicted TF binding site motifs and enrichment analyses indicate that age-related methylation change is influenced by developmental processes, while longevity-related DNAm change is associated with innate immunity or tumorigenesis genes, suggesting that bat longevity results from augmented immune response and cancer suppression.

Purpose: To independently assess upper urinary tract Memokath (MMK-051) stent outcomes in a national tertiary referral center. Materials and Methods: Two researchers, completely independent to the treating team, reviewed electronic MMK-051 stent(s) patient management records. Outcomes included time to first complication, complication(s)-severity, MMK-051 stent lifespan and change incidence, salvage therapy, further surgical intervention, and mortality. Results and Limitations: One hundred patients received 162 MMK-051 stent(s) (59% with malignant and 63% with distal ureteral obstruction [UO]) with only three lost to follow-up (FU). At 5-year mean FU, only 25 patients had complication-free original MMK-051 stents (14 alive, 11 dead). Of the remaining 75 patients, 22 had other stents, 12 had major surgery (e.g., nephrectomy), 3 became dialysis dependent, and 14 stabilized without ureteral stenting after original MMK-051 removal. Malignant obstruction patients had greater original MMK-051 stent longevity ( p  < 0.02), but also 20 of the 21 deaths (95%). The 72% mean 5-year stent complication rate included migration (46%), blockage (34%), nonfunctioning kidney (8%), urosepsis needing intravenous antibiotics (8%), and others (6%), including one postoperative death, one ureteral injury, and two with intractable pain. Median time to first complication was 12.5 months. Conclusions: MMK-051 stents had optimal utility in managing malignant UO and in those unfit for corrective surgery. Longer independently assessed mean 5-year outcomes review revealed much higher complication rates (72%) than previously reported. Future international metallic ureteral stent guidelines should encourage clinicians to adopt patient-centered multidisciplinary assessment and selection, with counseling plus goal-setting, and harmonized long-term protocol-based reporting, for optimized future patient safety and outcomes.

Background: Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants., Methods: We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity., Findings: When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10 -11 ; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10 -9 ; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10 -9 ). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10 -9 , nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10 -9 , CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants., Interpretation: This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry., Competing Interests: Declaration of Competing Interest GRA is employed by Regeneron Pharmaceuticals and he owns stock and stock options for Regeneron Pharmaceuticals. BMP serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. BMP reports serving on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. Y-DIC, LRY, JCM, BDM, JIR, KDT, JPL, EB, JAS, GRA report grants from NIH during the conduct of the study. Remaining authors have nothing to disclose., (Copyright © 2020 The Author(s). Published by Elsevier B.V. All rights reserved.)

Cell type heterogeneity presents a challenge to the interpretation of epigenome data, compounded by the difficulty in generating reliable single-cell DNA methylomes for large numbers of cells and samples. We present EPISCORE, a computational algorithm that performs virtual microdissection of bulk tissue DNA methylation data at single cell-type resolution for any solid tissue. EPISCORE applies a probabilistic epigenetic model of gene regulation to a single-cell RNA-seq tissue atlas to generate a tissue-specific DNA methylation reference matrix, allowing quantification of cell-type proportions and cell-type-specific differential methylation signals in bulk tissue data. We validate EPISCORE in multiple epigenome studies and tissue types.

Summary: It is well recognized that cell-type heterogeneity hampers the interpretation of Epigenome-Wide Association Studies (EWAS). Many tools have emerged to address this issue, including several R/Bioconductor packages that infer cell-type composition. Here we present a web application for cell-type deconvolution, which offers the functionality of our EpiDISH Bioconductor/R package in a user-friendly GUI environment. Users can upload their data to infer cell-type composition and differentially methylated cytosines in individual cell-types (DMCTs) for a range of different tissues., Availability and Implementation: EpiDISH web server is implemented with Shiny in R, and is freely available at https://www.biosino.org/EpiDISH/., (© The Author(s) 2019. Published by Oxford University Press.)

Summary: The Illumina Infinium EPIC BeadChip is a new high-throughput array for DNA methylation analysis, extending the earlier 450k array by over 400 000 new sites. Previously, a method named eFORGE was developed to provide insights into cell type-specific and cell-composition effects for 450k data. Here, we present a significantly updated and improved version of eFORGE that can analyze both EPIC and 450k array data. New features include analysis of chromatin states, transcription factor motifs and DNase I footprints, providing tools for epigenome-wide association study interpretation and epigenome editing., Availability and Implementation: eFORGE v2.0 is implemented as a web tool available from https://eforge.altiusinstitute.org and https://eforge-tf.altiusinstitute.org/., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2019. Published by Oxford University Press.)

Background: While smoking is known to associate with development of multiple diseases, the underlying mechanisms are still poorly understood. Tobacco smoking can modify the chemical integrity of DNA leading to changes in transcriptional activity, partly through an altered epigenetic state. We aimed to investigate the impact of smoking on lung cells collected from bronchoalveolar lavage (BAL)., Methods: We profiled changes in DNA methylation (5mC) and its oxidised form hydroxymethylation (5hmC) using conventional bisulphite (BS) treatment and oxidative bisulphite treatment with Illumina Infinium MethylationEPIC BeadChip, and examined gene expression by RNA-seq in healthy smokers., Findings: We identified 1667 total 5mC + 5hmC, 1756 5mC and 67 5hmC differentially methylated positions (DMPs) between smokers and non-smokers (FDR-adjusted P <.05, absolute Δβ >0.15). Both 5mC DMPs and to a lesser extent 5mC + 5hmC were predominantly hypomethylated. In contrast, almost all 5hmC DMPs were hypermethylated, supporting the hypothesis that smoking-associated oxidative stress can lead to DNA demethylation, via the established sequential oxidation of which 5hmC is the first step. While we confirmed differential methylation of previously reported smoking-associated 5mC + 5hmC CpGs using former generations of BeadChips in alveolar macrophages, the large majority of identified DMPs, 5mC + 5hmC (1639/1667), 5mC (1738/1756), and 5hmC (67/67), have not been previously reported. Most of these novel smoking-associating sites are specific to the EPIC BeadChip and, interestingly, many of them are associated to FANTOM5 enhancers. Transcriptional changes affecting 633 transcripts were consistent with DNA methylation profiles and converged to alteration of genes involved in migration, signalling and inflammatory response of immune cells., Interpretation: Collectively, these findings suggest that tobacco smoke exposure epigenetically modifies BAL cells, possibly involving a continuous active demethylation and subsequent increased activity of inflammatory processes in the lungs. FUND: The study was supported by the Swedish Research Council, the Swedish Heart-Lung Foundation, the Stockholm County Council (ALF), the King Gustav's and Queen Victoria's Freemasons' Foundation, Knut and Alice Wallenberg Foundation, Neuro Sweden, and the Swedish MS foundation., (Copyright © 2019. Published by Elsevier B.V.)

The molecular alterations that occur in cells before cancer is manifest are largely uncharted. Lung carcinoma in situ (CIS) lesions are the pre-invasive precursor to squamous cell carcinoma. Although microscopically identical, their future is in equipoise, with half progressing to invasive cancer and half regressing or remaining static. The cellular basis of this clinical observation is unknown. Here, we profile the genomic, transcriptomic, and epigenomic landscape of CIS in a unique patient cohort with longitudinally monitored pre-invasive disease. Predictive modeling identifies which lesions will progress with remarkable accuracy. We identify progression-specific methylation changes on a background of widespread heterogeneity, alongside a strong chromosomal instability signature. We observed mutations and copy number changes characteristic of cancer and chart their emergence, offering a window into early carcinogenesis. We anticipate that this new understanding of cancer precursor biology will improve early detection, reduce overtreatment, and foster preventative therapies targeting early clonal events in lung cancer.

Introduction: Even though ocular refractive state is highly heritable and under strong genetic control, the identification of susceptibility genes remains a challenge. Several HGF (hepatocyte growth factor) gene variants have been associated with ocular refractive errors and corneal pathology., Purpose: Here, we assess the association of an HGF gene variant, previously reported as associated with hyperopia, and ocular biometric parameters in a multicenter Spanish cohort., Methods: An observational prospective multicenter cross-sectional study was designed, including a total of 403 unrelated subjects comprising 188 hyperopic children (5 to 17 years) and 2 control groups: 52 emmetropic adolescents (13 to 17 years) and 163 emmetropic young adults (18 to 28 years). Each individual underwent a comprehensive eye examination including cycloplegic refraction, and topographic and ocular biometric analysis. Genomic DNA was extracted from oral swabs. HGF single nucleotide polymorphism (SNP) rs12536657 was genotyped. Genotypic, allelic, and logistic regression analyses were performed comparing the different groups. A quantitative trait association test analyzing several biometric parameters was also performed using generalized estimating equations (GEEs) adjusting for age and gender., Results: No association between rs12536657 and hyperopia was found through gender-adjusted logistic regression comparing the hyperopic children with either of the two control groups. Significant associations between mean topographic corneal curvature and rs12536657 for G/A (slope = +0.32; CI 95%: 0.04-0.60; p =0.023) and A/A (slope = +0.76; CI 95%: 0.12-1.40; p =0.020) genotypes were observed with the age- and gender-adjusted univariate GEE model. Both flat and steep corneal topographic meridians were also significantly associated with rs12536657 for the G/A and A/A genotypes. No association was found between rs12536657 and any other topographic or biometric measurements., Conclusions: Our results support a possible role for HGF gene variant rs12536657 in corneal curvature in our population. To our knowledge, this is the first multicenter quantitative trait association study of HGF genotypes and ocular biometric parameters comprising a pediatric cohort.

Study Objectives: (1) To determine the prevalence of snoring, breathing pauses during sleep and obstructive sleep apnoea syndrome in the United Kingdom (UK) and determine the relation between these events and obesity and other sociodemographic variables. (2) To compare and integrate this data with published UK population data., Methods: A total of 664 women and 575 men aged 18-100 years who formed a representative sample of the non-institutionalised UK population participated in an online interview survey directed by a previously validated computerised system., Results: Overall, 38% of men and 30.4% of women report that they snore at night. Furthermore, 8.7% of men and 5.6% of women state that they stop breathing at night. Comparing our data to published data from the 1990s, this study observes a highly significant increase in the rates of reported breathing pauses during sleep (sleep apnoea) in the UK over the last 20 years (p < 0.0001). In addition, we observe a highly significant increase in the prevalence of obesity (BMI>30) in the UK population between 1994 and 2015 (p < 0.0001). Integration of our data with NHS and public health England data on obesity confirms this increase., Conclusions: Our data demonstrate a significant increase in the rates of reported breathing pauses during sleep (sleep apnoea) and obesity in the UK over the last 20 years. Sociodemographic and behavioural changes have likely contributed to this. Moreover, our data also suggests that sleep disordered breathing (SDB) is widely underdiagnosed in the UK., (Crown Copyright © 2018. Published by Elsevier B.V. All rights reserved.)

An outstanding challenge of epigenome-wide association studies (EWASs) performed in complex tissues is the identification of the specific cell type(s) responsible for the observed differential DNA methylation. Here we present a statistical algorithm called CellDMC ( https://github.com/sjczheng/EpiDISH ), which can identify differentially methylated positions and the specific cell type(s) driving the differential methylation. We validated CellDMC on in silico mixtures of DNA methylation data generated with different technologies, as well as on real mixtures from epigenome-wide association and cancer epigenome studies. CellDMC achieved over 90% sensitivity and specificity in scenarios where current state-of-the-art methods did not identify differential methylation. By applying CellDMC to an EWAS performed in buccal swabs, we identified smoking-associated differentially methylated positions occurring in the epithelial compartment, which we validated in smoking-related lung cancer. CellDMC may be useful in the identification of causal DNA-methylation alterations in disease.

Background: DNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA., Results: Good-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1-10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10 -16 ) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10 -16 ). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3., Conclusions: We conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.

Background: Intra-sample cellular heterogeneity presents numerous challenges to the identification of biomarkers in large Epigenome-Wide Association Studies (EWAS). While a number of reference-based deconvolution algorithms have emerged, their potential remains underexplored and a comparative evaluation of these algorithms beyond tissues such as blood is still lacking., Results: Here we present a novel framework for reference-based inference, which leverages cell-type specific DNAse Hypersensitive Site (DHS) information from the NIH Epigenomics Roadmap to construct an improved reference DNA methylation database. We show that this leads to a marginal but statistically significant improvement of cell-count estimates in whole blood as well as in mixtures involving epithelial cell-types. Using this framework we compare a widely used state-of-the-art reference-based algorithm (called constrained projection) to two non-constrained approaches including CIBERSORT and a method based on robust partial correlations. We conclude that the widely-used constrained projection technique may not always be optimal. Instead, we find that the method based on robust partial correlations is generally more robust across a range of different tissue types and for realistic noise levels. We call the combined algorithm which uses DHS data and robust partial correlations for inference, EpiDISH (Epigenetic Dissection of Intra-Sample Heterogeneity). Finally, we demonstrate the added value of EpiDISH in an EWAS of smoking., Conclusions: Estimating cell-type fractions and subsequent inference in EWAS may benefit from the use of non-constrained reference-based cell-type deconvolution methods.

Epigenome-wide association studies (EWAS) provide an alternative approach for studying human disease through consideration of non-genetic variants such as altered DNA methylation. To advance the complex interpretation of EWAS, we developed eFORGE (http://eforge.cs.ucl.ac.uk/), a new standalone and web-based tool for the analysis and interpretation of EWAS data. eFORGE determines the cell type-specific regulatory component of a set of EWAS-identified differentially methylated positions. This is achieved by detecting enrichment of overlap with DNase I hypersensitive sites across 454 samples (tissues, primary cell types, and cell lines) from the ENCODE, Roadmap Epigenomics, and BLUEPRINT projects. Application of eFORGE to 20 publicly available EWAS datasets identified disease-relevant cell types for several common diseases, a stem cell-like signature in cancer, and demonstrated the ability to detect cell-composition effects for EWAS performed on heterogeneous tissues. Our approach bridges the gap between large-scale epigenomics data and EWAS-derived target selection to yield insight into disease etiology., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: The concept of reproducibility is a foundation of the scientific method. With the arrival of fast and powerful computers over the last few decades, there has been an explosion of results based on complex computational analyses and simulations. The reproducibility of these results has been addressed mainly in terms of exact replicability or numerical equivalence, ignoring the wider issue of the reproducibility of conclusions through equivalent, extended or alternative methods., Results: We use case studies from our own research experience to illustrate how concepts of reproducibility might be applied in computational biology. Several fields have developed 'minimum information' checklists to support the full reporting of computational simulations, analyses and results, and standardised data formats and model description languages can facilitate the use of multiple systems to address the same research question. We note the importance of defining the key features of a result to be reproduced, and the expected agreement between original and subsequent results. Dynamic, updatable tools for publishing methods and results are becoming increasingly common, but sometimes come at the cost of clear communication. In general, the reproducibility of computational research is improving but would benefit from additional resources and incentives., Conclusions: We conclude with a series of linked recommendations for improving reproducibility in computational biology through communication, policy, education and research practice. More reproducible research will lead to higher quality conclusions, deeper understanding and more valuable knowledge.

The exact timing and contribution of epigenetic reprogramming to carcinogenesis are unclear. Women harbouring BRCA1/2 mutations demonstrate a 30-40-fold increased risk of high-grade serous extra-uterine Müllerian cancers (HGSEMC), otherwise referred to as 'ovarian carcinomas', which frequently develop from fimbrial cells but not from the proximal portion of the fallopian tube. Here we compare the DNA methylome of the fimbrial and proximal ends of the fallopian tube in BRCA1/2 mutation carriers and non-carriers. We show that the number of CpGs displaying significant differences in methylation levels between fimbrial and proximal fallopian tube segments are threefold higher in BRCA mutation carriers than in controls, correlating with overexpression of activation-induced deaminase in their fimbrial epithelium. The differentially methylated CpGs accurately discriminate HGSEMCs from non-serous subtypes. Epigenetic reprogramming is an early pre-malignant event integral to BRCA1/2 mutation-driven carcinogenesis. Our findings may provide a basis for cancer-preventative strategies.

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation.

Background: There is growing evidence that DNA methylation alterations contribute to carcinogenesis. While cancer tissue exhibits widespread DNA methylation changes, the proportion of tissue-specific versus tissue-independent DNA methylation alterations in cancer is unclear. In addition, it is unknown which factors determine the patterns of aberrant DNA methylation in cancer., Results: Using HumanMethylation450 BeadChips (450k), we here analyze genome-wide DNA methylation patterns of ten types of fetal tissue, in addition to matched normal-cancer data for corresponding tissue types, encompassing over 3000 samples. We demonstrate that the level of aberrant cancer DNA methylation in gene promoters and gene bodies is highly correlated between cancer types. We estimate that up to 60 % of the DNA methylation variation in a cancer genome of a given tissue type is explained by the corresponding variation in a cancer genome of another type, implying that much of the cancer DNA methylation landscape is tissue independent. We further show that histone marks in normal cells are better predictors of aberrant cancer DNA methylation than the corresponding signals in human embryonic stem cells. We build predictors of cancer DNA methylation patterns and show that although inclusion of three histone marks (H3K4me3, H3K27me3 and H3K36me3) improves model accuracy, the bivalent marks are the most predictive. Finally, we show that chromatin accessibility of gene promoters in normal tissue dictates the promoter's propensity to acquire aberrant DNA methylation in cancer in so far as it determines its level of DNA methylation in normal tissue., Conclusions: Our data show that a considerable fraction of the aberrant cancer DNA methylation landscape results from a mechanism that is largely tissue specific. Histone marks as specified in the normal cell of origin provide highly predictive models of aberrant cancer DNA methylation and outperform those derived from the same marks in hESCs.